Question: Mapping to circular genome (mtDNA)
0
gravatar for Emilie
20 months ago by
Emilie0
Emilie0 wrote:

Hi,

I am working on detection of mitochondrial DNA mutations. I have performed sequencing of mtDNA from mice, and am using MapSplice2 for further analysis.

I was getting some unexpected results (or rather, I was lacking some expected large deletions) with MapSplice2. I looked at the coverage across the mt genome, and found that the sequence coverage is dropping considerably at the ends of the mtDNA. Which of course in retrospect makes sense, as my circular genome is being handled as a linear one, which is probably causing problems in the mapping (MapSplice2 uses Bowtie)....

But what are the options for mappers that can handle circularity? I am familiar with STAR, but did not find any information on it supporting circular genomes.

Thanks!

mapsplice alignment genome • 802 views
ADD COMMENTlink written 20 months ago by Emilie0

Hi, I'm also not aware of any mappers that can take circular references as input.

A simple workaround could be that you just cut your mt genome in two and reverse paste the ends back to each other , creating a 'new genome' , do the mapping with STAR and combine the results (if need-be)

ADD REPLYlink modified 20 months ago • written 20 months ago by lieven.sterck6.0k
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: 1025 users visited in the last hour