Extract unmapped reads OR reads covering chromosome 6 from a bam file?
Is there a one-liner approach to extracting reads from a .bam file that are unmapped OR aligned to chromosome 6?
Edit: ChatGPT give this response to the question: samtools view -f 4 -b input.bam 6 > chr6_and_unmapped.bam
Edit: Changed AND to OR
• 2,232 views
•
link
written
by
bioinf_sci •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
bam file: unmapped but matching reads (samtools view -f 4)
written by Tommy •
-
Samtools mapped and unmapped flags, what is the difference between the flags?
written by DNAngel
-
extract unmapped reads from WGS data - which samtools parameters to choose
written by Palgrave •
-
Get unmapped reads from cram/bam file
written by Palgrave •
-
Samtools SORT or bedtools bamToFastq giving error when converting reads to fastq.
written by DNAngel
-
FASTQ exctract ID's
written by User000
-
Samtools view Exec format error Help
written by aman-coasab •
-
Extract reads where both pairs maps properly and none of the pairs are split from a bam file
written by rimjhim.roy.ch
-
Extracting mapped paired end FastQ reads from BAM file using Samtools - trying to understand -f 0x0…
written by acxcrew007 •
-
Sam Reads Not Aligned To Chromosome
written by bic.uwaterloo.ca •
That would seem to be the answer to the question you posed. You are sure you do not want reads that aligned to Chr 6 OR unmapped reads?
I'm doing HLA typing and need reads that aligned to chromosome 6 or that were unmapped.
If a read was aligned to a chromosome how would it be unmapped at the same time?
I mean or as opposed to and. I need reads that are mapped to chromosome 6. I also need unmapped reads.
I see. So get the reads that were unmapped in one operation and the get the ones mapped to Chr6 in other to keep things simple.
That's right, but after some time looking into it I don't think it can be done as a one liner, with samtools at least.
You may be able to use
teeprogram to split the input into two streams and operate on both (doing two operations in parallel) but it may be better to keep things simple.I'll look into that, thanks! I'm leaning towards two separate samtools commands, as you've suggested.
Is that a logical OR? Unless I am not missing something obvious, I don't see how a read can both map to chromosome 6 and be unmapped at the same time?
You are correct, I need OR