AnswerBun.com

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 AnswerBun.com. All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir