How to extract unmapped reads if you have paired end reads (fastq) and mapped BAM file?
I have the following files:
- sampleA_1.fastq.gz - Forward reads
- sampleA_2.fastq.gz - Reverse reads
- mapped.sorted.bam - Sorted BAM file of mapped reads (doesn't contain unmapped reads)
How can I use these 3 files to get the following files:
- sampleA_1.unmapped.fastq.gz
- sampleA_2.unmapped.fastq.gz
Can I use anything from the following tools I have installed in my current environment?
- bbmap (and all the software in this suite)
- samtools
- seqkit
1 answer
If your BAM file contains the unmapped reads (not all aligners include them) then you can use bbmap (reformat.sh) and samtools to get the unmapped reads (search for samtools threads here or look at inline help for reformat.sh).
If your BAM file does not contain unmapped reads then you will need to get the read headers from your BAM (mapped reads) and then use filterbyname.sh to get the unmapped reads from your original data files.
Log in to answer this question.
More posts like this
-
Understanding Discrepancies in Deduplication Rates Between Forward and Reverse Reads in Metagenomic…
written by Stiofan •
-
filterBam (Augustus) results explanation
written by Sanam •
-
How to get unaligned reads and aligned reads into separate files from SAM/BAM?
written by O.rka
-
How can I determine the number of detected genes and detected transcripts/isoforms in a BAM file + …
written by O.rka
-
Trimmomatic - can`t find a path for to file
written by vvs.hazia •
-
Bamtools stats and samtools flagstat produce different results and confusing my sambamba output.
written by DNAngel
-
Total number of trimmed reads less than number of mapped reads
written by kspata
-
Trimmomatic reverse reads not improve
written by saleh.sereshki •
-
scRNA-seq: Kallisto processing of bioproject data (fastq-dump)
written by bsmith030465
-
Question About Samtools Flag Number 117
written by Alby •
