Extract mapped SE and properly paired PE reads from a bam file
Hello, I have paired end (PE) sequencing data performed on primarily small fragments of DNA. During the trimming/filtering process, we merged overlapping reads such that the resulting bam file contains a mix of now single-end and PE reads. I'd like to filter the mapped output for properly mapped read pairs or mapped single reads, can anybody help me figure out how to do that rather than just extract all mapped reads?
Thanks!
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
small RNA-Seq analysis and doubts
written by Sergio •
-
Count reads per allele handling overlapping reads
written by pau.rodriguez.sodupe •
-
Flagstat vs stats Samtools
written by v.shapovalova1 •
-
samtools stats result different from samtools flagstat
written by snishtala03 •
-
How would you filter PCR duplicates for merged paired-end reads
written by hohoku •
-
Join mapped, overlapping, paired-end reads
written by Mr. Dave •
-
bowtie2 output all aligned reads from paired alignment (even unconcordant pairs or single reads)
written by JV
-
Counting Reads On Paired-End Strand-Specific Rnaseq Data
written by Arun
-
Convert Bam To Fastq With Pe & Se Reads
written by Pasta •
-
Fastq Quality Check
written by toshnam
one option is to use
samtools viewto pull out alignments with different alignment flags.For more info on alignment flags see https://broadinstitute.github.io/picard/explain-flags.html