TRTools
TRTools copied to clipboard
dumpSTR is crashed to process the male chrX
dumpSTR --vcf output/gangstr/SC074198_chrX.vcf --out output/dumpstr_call/SC074201_chrX --gangstr-min-call-DP 20 --gangstr-max-call-DP 1000 --gangstr-filter-spanbound-only --gangstr-filter-badCI
Traceback (most recent call last):
File "/data/zhuw10/conda_envs/snakemake/bin/dumpSTR", line 10, in <module>
sys.exit(run())
File "/data/zhuw10/conda_envs/snakemake/lib/python3.10/site-packages/trtools/dumpSTR/dumpSTR.py", line 1272, in run
retcode = main(args)
File "/data/zhuw10/conda_envs/snakemake/lib/python3.10/site-packages/trtools/dumpSTR/dumpSTR.py", line 1231, in main
record.vcfrecord.INFO['HWEP'] = utils.GetHardyWeinbergBinomialTest(allele_freqs, genotype_counts)
File "/data/zhuw10/conda_envs/snakemake/lib/python3.10/site-packages/trtools/utils/utils.py", line 312, in GetHardyWeinbergBinomialTest
if gt[1] not in allele_freqs.keys():
IndexError: tuple index out of range
I think the cause of error is due to the genotype of TR in chrX is not diploidy as expected (for example 0 NOT 0/0; see below): chrX 156025977 . gagggagggagg . . . END=156025988;RU=gagg;PERIOD=4;REF=3;GRID=1,6;STUTTERUP=0.05;STUTTERDOWN=0.05;STUTTERP=0.9;EXPTHRESH=-1 GT:DP:Q:REPCN:REPCI:RC:ENCLREADS:FLNKREADS:ML:INS:STDERR:QEXP:GGL 0:31:0.333374:3:3-3:17,14,0,0:3,17:NULL:213.751:389.199,116.541:0,0:-1,-1,-1:-131.004,-131.004,-114.659,-131.004,-114.659,-92.8311,-131.004,-114.659,-92.8311,-115.651,-131.004,-114.659,-92.8311,-115.651,-132.975,-131.004,-114.659,-92.8311,-115.651,-132.975,-150.938
I would like to know whether there is a solution to fix it.
Thanks,
Wei Wei Zhu