Question: Which aligner is best suited for circular human mitochondrial genome?
2
gravatar for nalandaatmi
3.8 years ago by
nalandaatmi70
United States
nalandaatmi70 wrote:

Dear All,

To perform a mitochondrial genome alignment. Which aligner is best suited for circular human mitochondrial genome?

I found out some of them bwa-mem, bowtie2.

Can you guys explain which is better and why?

assembly bwa alignment bowtie2 • 2.6k views
ADD COMMENTlink modified 3.8 years ago by Brian Bushnell16k • written 3.8 years ago by nalandaatmi70
1
gravatar for h.mon
3.8 years ago by
h.mon27k
Brazil
h.mon27k wrote:

Are you trying to align two mitochondrial genomes, next-gen DNAseq reads to a mt genome, or RNAseq to a mt genome? For DNAseq, BWA will be fine, for RNAseq, STAR or Subread should be fine. If you are building an index with just the mtDNA, you may have to adjust the index to account for the small genome. And if your reads are from whole cells / tissues, you will probably align NuMT reads to the mtDNA, this problem would probably be smaller if you align against the whole human genome and then extract reads mapping to the mtDNA.

edit: you may use GSNAP to specifically address mtDNA circularity:

Finally, you can provide information to gmap_build that certain
chromosomes are circular, with the -c or --circular flag.  The value
for these flags is a list of chromosomes, separated by commas.  The
gmap_build program will then allow GSNAP and GMAP to align reads
across the ends of the chromosome.  For example, the mitochondrial
genome in human beings is circular.
ADD COMMENTlink modified 3.8 years ago • written 3.8 years ago by h.mon27k

We have sequenced mitochondrial DNA samples. Initially the reads were mapped against the mitochondrial chromosome alone. There were some false positive hits. Inorder to avoid false positives hits. This time I am planning to map the reads against the whole human genome. 

Then I have been asked to take the reads that are failed to map in previous step should be remapped to mitochondrial chromosome. While remapping the reads, they stressed me to consider the fact of circular genome. Which will lead to a flat pileup of the reads across the whole mitochondrial chromosome, instead of the ends show a dropout in reads (pileup looks like inverted U-shape).

As you suggested, I am planning to use bwa-mem aligner. but I don't get the highlighted text. Can you explain if you are aware of it.

 

ADD REPLYlink written 3.8 years ago by nalandaatmi70
1
gravatar for Brian Bushnell
3.8 years ago by
Walnut Creek, USA
Brian Bushnell16k wrote:

If you want to align reads to a genome, I suggest BBMap.  If you want to align genomes to genomes, I suggest MUMmer - it's much better with very long sequences.  Circularity is not really all that important; it just affects the tips.

ADD COMMENTlink written 3.8 years ago by Brian Bushnell16k
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: 863 users visited in the last hour