超短reads(primer、barcode、UMI、index等)比對方法


二代reads最短都有50bp,所以大家常用的比對工具都是不支持50bp以下的reads的比對的。

但是,在實際中,我們確實又有比對super short reads的需求。

So,我找到了如下方法來比對super short reads:

1.msa.sh,bbmap里的,好用,但是太慢。

2.bbduk.sh,bbmap里的,不支持輸出sam文件,只能輸出map和unmap的reads。

3.bbmap.sh,好用,超快,但是不支持mismatch,可惜。

4.seal.sh,沒試過。

5.bowtie1,親測可用,輸出Sam文件。

 

參考:Good way to align very, very short reads?


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM