Quantcast
Channel: haplotypecaller — GATK-Forum
Viewing all articles
Browse latest Browse all 1335

Errors at the end of the SNP calling

$
0
0

Hi

I am following the next steps to call variants in RNA-seq data:

https://software.broadinstitute.org/gatk/guide/article?id=3891

I am runnig the step 6 of variant calling using this command line:

java -jar GenomeAnalysisTK.jar -T HaplotypeCaller -R /LUSTRE/usuario/egonzalez/referencias/Zea_mays.AGPv3.22.dna.genome.fasta -I PT_todo_splitN.bam -o PTxB73F1.vcf -dontUseSoftClippedBases -stand_call_conf 20.0 -stand_emit_conf 20.0

And that give me the next errors:

INFO  00:25:28,277 ProgressMeter -     1:301362893    1.766387238E9    33.1 h           67.0 s      100.0%   33.1 h       6.0 s 
INFO  00:25:29,258 VectorLoglessPairHMM - Time spent in setup for JNI call : 33.407765322 
INFO  00:25:29,259 PairHMM - Total compute time in PairHMM computeLikelihoods() : 87411.73185938601 
INFO  00:25:29,259 HaplotypeCaller - Ran local assembly on 221086 active regions 
INFO  00:25:29,756 ProgressMeter -            done    2.067864162E9    33.1 h           57.0 s      100.0%    33.1 h       0.0 s 
INFO  00:25:29,757 ProgressMeter - Total runtime 119221.45 secs, 1987.02 min, 33.12 hours 
INFO  00:25:29,758 MicroScheduler - 93331090 reads were filtered out during the traversal out of approximately 194330945 total reads (48.03%) 
INFO  00:25:29,758 MicroScheduler -   -> 4209 reads (0.00% of total) failing BadCigarFilter 
INFO  00:25:29,758 MicroScheduler -   -> 59464782 reads (30.60% of total) failing DuplicateReadFilter 
INFO  00:25:29,759 MicroScheduler -   -> 604373 reads (0.31% of total) failing FailsVendorQualityCheckFilter 
INFO  00:25:29,759 MicroScheduler -   -> 33257726 reads (17.11% of total) failing HCMappingQualityFilter 
INFO  00:25:29,760 MicroScheduler -   -> 0 reads (0.00% of total) failing MalformedReadFilter 
INFO  00:25:29,760 MicroScheduler -   -> 0 reads (0.00% of total) failing MappingQualityUnavailableFilter 
INFO  00:25:29,760 MicroScheduler -   -> 0 reads (0.00% of total) failing NotPrimaryAlignmentFilter 
INFO  00:25:29,761 MicroScheduler -   -> 0 reads (0.00% of total) failing UnmappedReadFilter 
INFO  00:25:32,255 GATKRunReport - Uploaded run statistics report to AWS S3 

I am using GATK/3.5.0. At the end I have the vcf file and everything looks normal in it, but I am not sure which is the problem with these errors. Do you know which is the problem?


Viewing all articles
Browse latest Browse all 1335