SEQanswers (
-   Bioinformatics (
-   -   BWA MEM question (

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,

sarvidsson 02-24-2015 05:48 AM


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.