Question: Circular Genome Alignment
0
gravatar for isaac.rstefan
6 months ago by
isaac.rstefan0 wrote:

Hello all,

I'm aligning sequencing data to the human mitochondrial genome with Bowtie2. To handle the fact that the genome is circular, I'm doing two alignments (in which one has a shifted origin) as aligners treat the genome as linear. My goal was to replace the reads from one to the other where the origin is. However, when I do this, regions that should have the same number of reads in each alignment are significantly different. Further, the total number of reads that align is also VERY different. Does anyone have any suggestions on what might be going wrong or what to do?

Best, Stefan

alignment • 362 views
ADD COMMENTlink modified 5 months ago by Biostar ♦♦ 20 • written 6 months ago by isaac.rstefan0

A common technique is adding a bit of "over-hangs" from the other side of the break at the extremes of the linearized circular contig - this would allow the reads to map at the extremities, and is simpler than the approach you are taking.

ADD REPLYlink written 6 months ago by h.mon24k

bwa has the option to apply soft and hard clipping, it may be usefull, or hisat2 by default performs spliced alignments.

ADD REPLYlink written 5 months ago by Buffo1.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1388 users visited in the last hour