extract chimeric and multimap reads from bam file

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.

One Answer

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

Add your own answers!

Related Questions

Generating 3D coordinates error

1  Asked on January 15, 2021 by shahbaaz


BAM file filteing to remain best isoform

0  Asked on January 10, 2021 by user977828


Somatic mutations for normal WES samples

0  Asked on January 6, 2021 by lot_to_learn


Get list of urls of GSM data set of a GSE set

1  Asked on January 6, 2021 by user432797


Biohackers Netflix – DNA to binary and video

1  Asked on January 3, 2021 by xamax


DNA sequence error annotation

0  Asked on December 30, 2020 by matthew-jones


samtools / bamUtil | Meaning of as Reference Name

1  Asked on December 25, 2020 by paul-endymion


How to remove batch effect from TCGA and GTEx data

2  Asked on December 22, 2020 by kai-he


Ask a Question

Get help from others!

© 2023 All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir