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

QcConpair exits with error #924

Open
anoronh4 opened this issue Sep 15, 2021 · 1 comment
Open

QcConpair exits with error #924

anoronh4 opened this issue Sep 15, 2021 · 1 comment
Assignees
Labels
backburner probably won't address in a near future bug Something isn't working qc related to qc modules

Comments

@anoronh4
Copy link
Collaborator

anoronh4 commented Sep 15, 2021

A small minority of samples (so far 1 in wes_repo) repeatedly error out with the following:

Traceback (most recent call last):
  File "/usr/bin/conpair/scripts/verify_concordances.py", line 315, in <module>
    main()
  File "/usr/bin/conpair/scripts/verify_concordances.py", line 183, in main
    str_return = run_verify_concordance(N_pileup, T_pileup, MARKER_FILE, MMQ, MBQ, COVERAGE_THRESHOLD, AA_BB_only)
  File "/usr/bin/conpair/scripts/verify_concordances.py", line 220, in run_verify_concordance
    Tumor_genotype_likelihoods = genotype_likelihoods_for_markers(Markers, T_pileup, min_map_quality=MMQ, min_base_quality=MBQ)
  File "/usr/bin/conpair/modules/ContaminationMarker.py", line 106, in genotype_likelihoods_for_markers
    pileup = parse_mpileup_line(line, min_map_quality=min_map_quality, min_base_quality=min_base_quality)
  File "/usr/bin/conpair/modules/ContaminationMarker.py", line 76, in parse_mpileup_line
    mapqs_above_threshold = set([i for i in range(0, len(mapqs)) if int(mapqs[i]) >= min_map_quality])
ValueError: invalid literal for int() with base 10: 'A00227:20:H3CFFDMXX:1:1440:7500:17'

I will try to remove the cached directories for the previous step (QcPileup) to see if that gets rid of the error. If it does, then we should try to see if we can validate the Pileup file before proceeding to QcConpair.

@anoronh4 anoronh4 self-assigned this Sep 20, 2021
@anoronh4
Copy link
Collaborator Author

anoronh4 commented Sep 20, 2021

after forcing nextflow to repeat QcPileup for the normal sample, the QcConpair completed with success. So far this occurred in 1 sample, so this is definitely not high priority

@anoronh4 anoronh4 added bug Something isn't working qc related to qc modules backburner probably won't address in a near future labels Oct 5, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
backburner probably won't address in a near future bug Something isn't working qc related to qc modules
Projects
None yet
Development

No branches or pull requests

1 participant