Question: Circular Genome Alignment
gravatar for isaac.rstefan
22 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 • 976 views
ADD COMMENTlink modified 8 months ago by Jonathanjacobs250 • written 22 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 22 months ago by h.mon30k

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

ADD REPLYlink written 21 months ago by Buffo1.8k
gravatar for Jonathanjacobs
8 months ago by
Rockville, MD
Jonathanjacobs250 wrote:

FWIW - CLC Genomics Workbench's read mapper supports circularized reference genomes and will map reads across origins without any issue.

ADD COMMENTlink written 8 months ago by Jonathanjacobs250
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: 920 users visited in the last hour