I have a BED file with genomic regions. Is there a tool that, using this BED file, would help me to subset this BAM file so that all reads (and parts of reads) that fall outside of the regions in the BED file would be soft-clipped?
•
written
by
lechu •
1 answer
I looked at the PcrClipReads tool developed by @Pierre. It's great, but as far as I understand, it does not handle cases of reads that overlap two (or more) regions. I realized I was not precise in my question. In a situation like the one depicted below I would only like to clip bases marked x (read 1 should not be clipped). The goal is to use exonic regions and remove everything the protrudes outside of exonic ranges (this may require "internal' clips in the reads). Equally good alternative would be to somehow reset the qualities of the bases to be clipped to zero (instead of clipping). I tried to do it using sam2tsv (also from @Pierre), but ended up with file sizes that were impossible to handle (I need it for whole transcriptome data). Then I run out of ideas.
---------- --- read 1 (goes across two exons)
----------xxx --- read 2 (goes across two exons, some mapped to intron)
----------xxx read 3 (mapped to single regions + some protrusion)
|-REGION1-| |-REGION2-|
update 2021: there is now samtools ampliconclip samtools-ampliconclip http://www.htslib.org/doc/samtools-ampliconclip.html
Clip reads in a SAM compatible file based on data from a BED file.
Log in to answer this question.
-
trimming primers using iVar
written by binfo.sl •
-
Samtools ampliconclip hard clip not giving same results as soft clip
written by Kevin •
-
Clipping bases from non-overlapping regions of mate reads in BAM file
written by lechu •
-
Remove clipped reads
written by finswimmer
-
Double soft-clipped reads and how to romeve them
written by Carles Borreda •
-
Using reads aligned on the Negative strand of the genome to select soft clipped nucleotides
written by fusion.slope
-
Counting reads within regions of interest
written by EVR
-
tool for bam soft clipping reads within bed file regions
written by 14134125465346445
-
Samtools mpileup of soft-clipped alignment behaves non-deterministically?
written by iovercast •
-
Filter Exon Aligned Reads That Are Mapped Against Human Genome.
written by Coryza •
