The frontline support team will be unavailable to answer questions on April 15th and 17th 2019. We will be back soon after. Thank you for your patience and we apologize for any inconvenience!
Different coverage metrics between GATK and PICARD
I want to evaluate the panel performance by calculating the target coverage.
Both GATK and PICARD offered great convenience for the calculation. However, slightly difference was observed.
And it seems that PICARD can not calculate the percentage of target base larger than pre-defined coverage as GATK, by using -ct option. Is there any option?
java -jar /bioinfo/software/bin/GenomeAnalysisTK.jar -T DepthOfCoverage -R /bioinfo/data/iGenomes/genome.fa -o gatk -I ../4D_S1_dedup.bam -ct 1 -ct 2 -ct 10 -ct 30 -ct 50 -ct 100 -L ../../../Agilent_ClearSeq_Inherited.v1.bed java -jar /bioinfo/software/packages/picard-2.9.2/picard.jar CollectHsMetrics I=../4D_S1_dedup.bam O=picard.metrics R=/bioinfo/data/iGenomes/genome.fa BI=../agilent.intervals TI=../agilent.intervals MQ=0 Q=0
agilent.intervals was obtained by PICARD
java -jar /bioinfo/software/packages/picard-2.9.2/picard.jar BedToIntervalList I=Agilent_ClearSeq_Inherited.v1.bed O=Agilent_ClearSeq_Inherited.v1.interval SD=/bioinfo/data/iGenomes/genome.dict