Remove clipped reads
Hello,
is there a tool or script which can remove reads from my bam files which contain soft clipped bases? It would be nice if I could also define to just remove reads that have at least n soft clipped bases, or the fraction of soft clipped bases is more than x.
Thanks for any advice.
fin swimmer
• 6,475 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
finswimmer
•
written
by
finswimmer
1 answer
up-to-date answer using samjdk: http://lindenb.github.io/jvarkit/SamJdk.html , removing clipped reads where more than 33% of the read is clipped.
java -jar dist/samjdk.jar -e 'return record.getReadUnmappedFlag() || record.getCigar()==null || record.getCigar().getCigarElements().stream().filter(C->C.getOperator().isClipping()).mapToInt(C->C.getLength()).sum() / (double)record.getCigar().getReadLength() < 0.33;'
• 0 views
•
link
written
by
Pierre Lindenbaum
replace : C.getOperator().isClipping()
with : C.getOperator().name().equals("H")
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
paired-end bam clip bases outside insert range
written by pau.rodriguez.sodupe •
-
How to selectively turn soft-clipped bases to hard-clipped bases based on part of the sequence?
written by goldberry88 •
-
bamtools filtering syntax - reads containing more than X mismatches and soft-clipped reads
written by joelepaul •
-
Calculate Depth with and without Soft-Clipping
written by ccnn •
-
Counting soft clipped bases and reads
written by finswimmer
-
Excluding reads for Coverage Calculation in given regions
written by finswimmer
-
Double soft-clipped reads and how to romeve them
written by Carles Borreda •
-
Soft-clipping of reads in Amplicon-sequenced data
written by jsneaththompson
-
bcftools consensus: include the OVERHANG
written by Marvin
-
Disable Soft-Clipping Bwa Aln
written by Irsan
@Pierre to rescue: Remove Soft Clipped Bases
How to remove reads with hard/soft clipping along with its mate?
This seems to remove only the soft clipped bases. But I'd like to get rid of the whole reads.
fin swimmer
There are two links in the post. Did you look at the second?
Oh, sorry haven't seen your edit.
But in your first link I also found a link to bamutils. A combination of removeclipping and filter seems to fit my needs. It's only a pity that bamutils seems not to support streaming.
Thanks for showing me the right direction.
fin swimmer
It's probably a bit late as you already have a BAM but if you really don't want any read clipped, your probably best off forcing your aligner to align entire reads (e.g. bowtie2 end-to-end).