SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   How Do Mappers Handle Chromosome Ends? (http://seqanswers.com/forums/showthread.php?t=70548)

nw328 07-28-2016 07:44 AM

How Do Mappers Handle Chromosome Ends?
 
I am trying to figure out how BW-based mappers like BWA and Bowtie handle reads that map at the end of the reference and extend off the end. One example would be reads from the origin of a circular chromosome, where part maps to the beginning of the sequence used to build the index, and part of it maps to the very end. How are such reads scored? Is each extending base considered a mismatch, gap extension, or something else?

Thanks in advance!


NOTE:Also posted on Biostars

tristan dubos 07-28-2016 11:59 PM

Hi,
I was with the same problematic when I studies repeats region to determine number copy. I saw in BWA(mem) no mapping to the extremity of my reference (500 bases). To solve the problem by put 2 times my references aside. I tried before this to put N to the extremity of the references and low stringency alignment parameters but i was i low complexity sequences the i had everything mapping on the position that is why i used the first solution.
I hope it could help

Tristan

maubp 07-30-2016 01:18 PM

Very few mappers have support for circular genomes (e.g. CLC mapper), and worse, the SAM/BAM format does not handle this well. As suggested by Tristan, mapping to a doubled linear reference can be a useful alternative (but brings with it a lot of complications).

If you are mainly interested in mapping to the origin region of a circular chromosome, it would be simpler to make a linear chromosome split elsewhere in the circle, and map to that.


All times are GMT -8. The time now is 07:21 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.