SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   BWA MEM question (http://seqanswers.com/forums/showthread.php?t=50532)

Liy 02-24-2015 05:38 AM

BWA MEM question
 
Dear all,

I just have a question about bwa mem alignment. I ran the following simple command:

bwa mem <fa> <fastq> -t 5

Before the alignment I checked the number of reads with fastqc and then after alignment I counted the IDs present in the sam file with grep '<common ID> | wc -l and I got different results.
What causes the difference? Is the multi-mapper reads present in the sam file?

Total reads: 24 961 320
Aligned reads: 24 961 535

Thank you for the answere,
Regards,
Liy

sarvidsson 02-24-2015 05:48 AM

Run
Code:

samtools view -S -f 0x100 <file.sam>
to get the multimappers.


All times are GMT -8. The time now is 02:46 AM.

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