Hi,
Is it possible to extract paired-end reads that map to a specific interval ( from a bam file ). I tried with intersectBed :
intersectBed -abam align.bam -b interval.gff3 -wa > result.bam
here's the result :
But I only want reads that map to the feature in bold blue (one of the paired reads is enough). For example, I don't want the reads that map either side of this feature (red arrow).
Is it possible with intersectbed or an other program ?
Thanks,
N.
•
written
by
Nicolas Rosewick
Use samtools view:
samtools view [options] in.bam "chr2:1000-078987"
Thanks, but I've also the same problem (reads that map either side of the feature (red arrow))
for paired-end bams, you should use pairToBed instead of intersectBed. the syntax would be almost the same:
pairToBed -abam align.bam -b interval.gff3 > result.bam
I tried pairToBed but I've an error with the bam file :
*****ERROR: BAM must be sorted by query name.
from the BEDTools' pairToBed docs: "pairToBed will fail if an aligner does not report alignments in pairs". I guess that if this is the case, you'll have to work with a compatible aligners or to sort your read pairs as requested.
[BEDOPS bedmap] could probably solve this by only reporting overlapping reads ("reference" elements), when the read is contained entirely within the set of map elements (the intervals of interest).
The bedmap --fraction-ref 1 option applies this 100% reference-overlap criteria, while --echo reports any qualifying reference elements (reads). We can use bash substitutions and convert2bed convenience scripts bam2bed and gff2bed to convert the reads and intervals to BED.
$ bedmap --echo --fraction-ref 1 <(bam2bed < align.bam) <(gff2bed < interval.gff3) > answer.bed
The use of bam2bed generates sorted BED-formatted reads, which makes the reads in answer.bed sorted and ready to use in downstream BED operations.
•
written
by
Alex Reynolds
Log in to answer this question.
The image link is broken, I cannot see it.
Could you please provide the content of the interval.gff3 file? I might want to do a very similar analysis, any help is welcome