im trying to calculate the coverage of a exome sequencing project.
My Code is
java -jar $gatk_jar \
-T DepthOfCoverage \
-R $ref_genome \
-I "$sample_name".recal_reads.bam \
-o $result_dir/$sample_name/coverage/"$sample_name" \
-geneList $ref_refseq \
-L test.bed \
-ct 10 -ct 20 -ct 30
but i get an ERROR that i cant solve...
Could not parse location from line: chr1 66999814 67000061 NM_032291_exon_0_10_chr1_66999825_f 0 +
This error relates to my refseq_file but i cant find the problem with it. This is the first line from my geneList so something is messed up.
I created it along http://www.broadinstitute.org/gatk/guide/article?id=1329 and used the script sortByRef.pl to sort it by the fai file from your bundle
Can you give me any advise to check on ?
Thank you very much for helping!