We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

SplitNCigarReads: Hard clipping of overhangs


Running SplitNCigarReads doesn't seem to hard-clip intronic overhangs despite using aggressive clipping parameters
I ran gatk (version: with the following command

gatk SplitNCigarReads -R /path/to/genome/hg19.fa -I /path/to/$inputbam --max-mismatches-in-overhang 0 --max-bases-in-overhang 5 -O $outbam

I'm assuming the exon boundaries are inferred. if not, is there an option to provide transcript annotations to SplitNCigarReads.


Best Answer


Sign In or Register to comment.