Does A Bwa Bam File Contain Entries For All Unaligned Paired-End Reads?
i have a basic question: I am using default parameter to run the bwa to align to a a reference geneome. I would like to know if bam represents all the reads used in the paired end fastq file? (including reads from pairs where one mate maps, and where no mates map). Why Iam asking this is, I need to rebuild the fastq using bam?
• 5,911 views
•
link
updated
by
Istvan Albert
•
written
by
Rm
•
written
by
Rm
1 answer
Yes the file contains all reads from your input file. You may use samtools to filter for reads that have proper pairs, only one end etc.
• 1 views
•
link
written
by
Istvan Albert
yes, here is a way to look at the flags: http://picard.sourceforge.net/explain-flags.html
• 0 views
•
link
updated
by
h.mon
•
written
by
Istvan Albert
•
written
by
Istvan Albert
Log in to answer this question.
More posts like this
-
bwa mem skip orientation when aligning a small subset of reads
written by Yihang •
-
paired reads have different names BWA MEM after samtools bam > fastq
written by robert.murphy •
-
featureCounts paired end mode
written by henry-keen •
-
Best SAM flags to evaluate host contamination
written by ctseto
-
To extract unmapped pair end reads
written by Bioinfonext
-
Uniquely mapped PAIRS of reads
written by abascalfederico
-
Mapping BAM/SAM file back to reference?!?!
written by ab.tsubaki •
-
Bwa Mapping Of Solid Paired End Reads
written by William
-
Bwa Mismatches And Seed Length
written by GPR
-
Tophat Parameters Misunderstanding
written by Nicolas Rosewick
