SplitNCigarReads: Hard clipping of overhangs

Hi,

Running SplitNCigarReads doesn't seem to hard-clip intronic overhangs despite using aggressive clipping parameters
I ran gatk (version: 4.1.0.0) 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,
A

Best Answer

Answers

Sign In or Register to comment.