Zhengdeng Lei, PhD

Zhengdeng Lei, PhD

2009 - Present Research Fellow at Duke-NUS, Singapore
2007 - 2009 High Throughput Computational Analyst, Memorial Sloan-Kettering Cancer Center, New York
2003 - 2007 PhD, Bioinformatics, University of Illinois at Chicago

Sunday, June 2, 2013

mapping rates from BWA

You can run bamtools stats -in path/to/your/alignments.bam, which will (quickly) report general alignment statistics.
You can also run bamtools filter -isMapped false -in path/to/your/alignments.bam > unmapped.bam (or something quite close to that) to get your unmapped reads.


you can tell from the sam file, by checking the flags, or make the file into a .bam with samtools, and use flagstat from samtools to count.

No comments:

Post a Comment