Hi everybody
I am Juan and I am very very new in this community. I am starting doing NGS data analysis and I have some doubts.
I using BWA to map my reads against a reference sequence. Afterwards I convert the sam file into a bam file that I can use to get a consensus sequence.My doubt is whether the sam/bam file I get after the mapping includes paired reads as well as the unpaired reads matching the reference sequence or only the paired reads and removes the reads as they don't have a mate.
Thanks in advance
Juan
•
written
by
Juan Ledesma •
1 answer
It will have both reads for a pair irrespective of whether one of them remained unmapped or mapped at too many positions. The second column in the BAM file will be properly flagged as explained here: http://broadinstitute.github.io/picard/explain-flags.html
•
written
by
Ashutosh Pandey
Log in to answer this question.
-
Reconstruct reference sequence from SAM file?
written by eholleman •
-
HTseq error with paired end alignments
written by bnayer26 •
-
Should I keep or remove the alignment reads whose mate reads mapped to a different chr ?
written by dingailuma •
-
Bowtie fails to align reads that appear in the reference sequence
written by AM •
-
Mapping using multiple reference sequences
written by Juan Ledesma •
-
Where to find PCR duplicate reads in bam file?
written by bioinforesearchquestions
-
Mapping BAM/SAM file back to reference?!?!
written by ab.tsubaki •
-
How to generate contigs from bam alignment
written by goubert.clement •
-
flagstats results show 0.0 % properly paired reads
written by SOHAIL
-
What Better Way To Get The Paired Reads Aligned Against The Reference Genome?
written by thiagomafra
Thanks for your help