The Frontline Support team will be offline February 18 for President's Day but will be back February 19th. Thank you for your patience as we get to all of your questions!
Error wit BaseRecalibrator: The stop position ... is less than start ...
I am using GATK version 2.4.9. When trying to execute BaseRecalibrator I get this message:
ERROR MESSAGE: Badly formed genome loc: Parameters to GenomeLocParser are incorrect:The stop position 4580 is less than start 4581 in contig chr1
I'm not specifying any of the advanced options, just an input BAM, the reference genome and a BED file for the -knownSites argument. The BAM file has been correctly validated with Picard and I've checked the BED and all stop positions are greater or equal to the start positions.
I also tried with an older version of GATK so that I could use CountCovariates, but it gave a similar error (same error with different start and stop positions).
The exact call to the method is:
java -Xmx4g -jar GenomeAnalysisTK.jar -T BaseRecalibrator -I XXXX.marked.realigned.fixed.bam -R hg18KO.fa -knownSites dbsnp130_hg18_KO.bed -o XXXX_recal_data.grp
Any help would be greatly appreciated.