Hi all
I am trying to use the bwa for RNA-seq alignment, but looks like it reports all the reads without filtering the unmapped reads. I can't find any parameter related to filtering the unmapped reads. Does anyone try that before? Or we have to write a script to extract the aligned reads? Also, the format of the bwa generated sam file is bit different from the one I got from tophat. The output sam is shown:
D3B4KKQ1_0166:6:1101:2368:2227#NNNNAN 4 * 0 0 * * 0 0 ATTACTTGGAAACACCGGGATCATATGATTGTTCAATGAAGAAGCACGTCTGAT *
I am confused about
(4 * 0 0 * * 0 0)
Could you help me figure out the meaning of each item?
I am really grateful if anyone can provide me the related information!!
I am trying to use the bwa for RNA-seq alignment, but looks like it reports all the reads without filtering the unmapped reads. I can't find any parameter related to filtering the unmapped reads. Does anyone try that before? Or we have to write a script to extract the aligned reads? Also, the format of the bwa generated sam file is bit different from the one I got from tophat. The output sam is shown:
D3B4KKQ1_0166:6:1101:2368:2227#NNNNAN 4 * 0 0 * * 0 0 ATTACTTGGAAACACCGGGATCATATGATTGTTCAATGAAGAAGCACGTCTGAT *
I am confused about
(4 * 0 0 * * 0 0)
Could you help me figure out the meaning of each item?
I am really grateful if anyone can provide me the related information!!
Comment