Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

dumpSTR is crashed to process the male chrX #185

Open
weizhu365 opened this issue Jul 26, 2023 · 0 comments
Open

dumpSTR is crashed to process the male chrX #185

weizhu365 opened this issue Jul 26, 2023 · 0 comments

Comments

@weizhu365
Copy link

weizhu365 commented Jul 26, 2023

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant