Question: Circular Genome Alignment
gravatar for isaac.rstefan
12 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 • 578 views
ADD COMMENTlink modified 11 months ago by Biostar ♦♦ 20 • written 12 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 12 months ago by h.mon27k

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

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


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