I am trying to extract all chimeric and multi-map reads from either SAM/BAM file. Is there any simple command to do that? Can I use htslib for parsing sam/bam files and extracting based on any flags. Any help on how to start will be appreciated.
$\begingroup$
$\endgroup$
Add a comment
|
1 Answer
$\begingroup$
$\endgroup$
2
This is how I do it.
samtools view -@5 -f 0x800 -hb /path/sample.bam > sup.bam
Then if you want it as a fasta
samtools view -F 0x1 -hb sup.bam | samtools fasta -F 0x1 - > sup.fa
-
$\begingroup$ Thanks !! It works great. I will use samtools source code to write a small program to extract the reads based on flag. $\endgroup$– SBDK8219Feb 18, 2021 at 5:52
-