Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

Somehow the requested coordinate is not covered by the read?

$
0
0

I was running some samtools calls through UG with genotyping_mode GENOTYPE_GIVEN_ALLELES after reading this thread:

http://gatkforums.broadinstitute.org/discussion/1842/variant-annotator-not-annotating-mappingqualityranksumtest-and-readposranksumtest-for-indels

I got this error message:

##### ERROR MESSAGE: Somehow the requested coordinate is not covered by the read. Alignment 440969 | 36S26M1D6M1I2M1D1M3D27M1S

What did I do wrong? It works for most other fragments (chrom20:pos1-pos2). Is it possible to avoid this error message and just do ./. calls at coordinates not covered by reads? I'm not sure, why samtools called at a coordinate not covered by a read. I don't know at which coordinate this happened. I was at this coordinate, when the error happened:

INFO  10:21:06,036 ProgressMeter -       20:479014   2329564.0    55.0 m      23.6 m        0.6%     6.7 d       6.7 d 

I was using this command:

vcf=samtools.vcf.gz
$java -Djava.io.tmpdir=tmp -Xmx63900m \
 -jar $jar \
 --analysis_type UnifiedGenotyper \
 --reference_sequence $ref \
 --input_file bam.list \
 --dbsnp $dbsnp138 \
 --num_cpu_threads_per_data_thread 3 \
 --num_threads 4 \
 --out $out \
 --intervals $vcf \
 -stand_call_conf 10 \
 -stand_emit_conf 10 \
 --genotype_likelihoods_model BOTH \
 --output_mode EMIT_ALL_SITES \
 --annotation Coverage -A FisherStrand -A HaplotypeScore -A MappingQualityRankSumTest -A QualByDepth -A RMSMappingQuality -A ReadPosRankSumTest \
 --alleles $vcf \
 --genotyping_mode GENOTYPE_GIVEN_ALLELES \
 -L 20:$pos1-$pos2 \

Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>