Yes I believe this will solve my problem! Thank you! GenoMax
• 0 views
•
link
written
by
kaufman.eli.j •
How to filter bam file for only reads that span certain regions of the genome?
Hello, is there a tool that can find only the reads that span the entire section(s) that I am interested in. Ideally I would just need a sorted and indexed bam file and a bed file with the regions that I am interested in and it would output a new bam with just the reads of interest. I do not think that bedtools intersect would get me what I want, though I would be happy if someone could explain what arguments I might need if it can do what I want.
I want only the reads that span the entire region(s) of interest, as I am looking at highly repetitive regions using long read sequencing technologies.
Just checking in before I go ahead and make such a tool...
Thanks!
edit: it is possible that the bedtools intersect flag -f can let me do this, reading more now... "-f Minimum overlap required as a fraction of A. Default is 1E-9 (i.e. 1bp)." if I set f to a value of 1, it should only include reads that capture the entire region. I will continue to report back
• 1,926 views
•
link
updated
by
Ram
•
written
by
kaufman.eli.j •
•
written
by
kaufman.eli.j •
1 answer
Log in to answer this question.
More posts like this
-
Efficiently merge two BAM files while retaining reads from only one file in overlapping regions
written by kiranchari •
-
Find count of reads in BAM file that fall within BED region
written by KH •
-
Get nonoverlapping regions of two bed files
written by vctrm67 •
-
Paired-end reads that span a position in an assembly
written by saladi •
-
Extraction of read subsequences based on a region from a sam file
written by gordon.freeman •
-
Clarification on Bedtools Intersect Usage
written by Noam Teyssier
-
Does samtools view only show reads that completely span the region?
written by Noam Teyssier
-
extract reads from BAM whose fragments overlap a specified region
written by Richard
-
Counting reads within regions of interest
written by EVR
-
bedtools multicov and fraction overlap flag
written by fitziano •
How to Extract Specific Region on Bam file
If you need to extract the reads then a similar strategy works. Looks for similar threads here.
Thank you for the quick reply, but I believe that solution will include partial reads that are at that location. I want only reads that span the entire region(s) to build a consensus from. I have the read depth for it.