Bioinformatics Asked by SBDK8219 on March 26, 2021
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.
This is how I do it.
samtools view [email protected] -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
Answered by Liam McIntyre on March 26, 2021
1 Asked on January 15, 2021
1 Asked on January 13, 2021
0 Asked on January 10, 2021 by user977828
0 Asked on January 6, 2021 by lot_to_learn
1 Asked on January 6, 2021 by user432797
1 Asked on January 4, 2021 by manuel-milla
covid 19 interactions protein protein interaction protein structure sars cov 2
1 Asked on January 2, 2021 by marilu
0 Asked on December 30, 2020 by matthew-jones
1 Asked on December 30, 2020 by ryan-fahy
haplotypes networks phylogenetics phylogeny population genetics
0 Asked on December 29, 2020
1 Asked on December 29, 2020 by anamaria
1 Asked on December 26, 2020
1 Asked on December 25, 2020 by paul-endymion
1 Asked on December 23, 2020
2 Asked on December 22, 2020 by kai-he
0 Asked on December 22, 2020 by poccia
Get help from others!
Recent Questions
Recent Answers
© 2023 AnswerBun.com. All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir