cnvkit icon indicating copy to clipboard operation
cnvkit copied to clipboard

Unable to use fix without antitarget file for WGS samples

Open GACGAMA opened this issue 1 year ago • 0 comments

Hello

This is related to #242

It seems like the newest version contains some old bugs for some reason

When trying to generate .cnr files with fix command without antitarget for whole-genomes:

cnvkit.py fix /scratch4/nsobrei2/ggama1/somatic_SVs/cnvkit/call_genomes/BH10027_1_TUMOR.targetcoverage.cnn /scratch4/nsobrei2/ggama1/somatic_SVs/cnvkit/needed_files/genomes_reference.cnn -o /scratch4/nsobrei2/ggama1/somatic_SVs/cnvkit/call_genomes/BH10027_1_TUMOR.cnr

I get: cnvkit.py fix: error: the following arguments are required: reference

If I create a decoy empty file antitarget, as proposed on older version, I get:

Blank tab file?: /scratch4/nsobrei2/ggama1/somatic_SVs/cnvkit/call_genomes/BH10027_1_TUMOR.antitargetcoverage.cnn
Processing target: BH10027_1_TUMOR
Traceback (most recent call last):
  File "/home/ggama1/.conda/envs/cnvkit/bin/cnvkit.py", line 10, in <module>
    sys.exit(main())
  File "/home/ggama1/.conda/envs/cnvkit/lib/python3.10/site-packages/cnvlib/cnvkit.py", line 10, in main
    args.func(args)
  File "/home/ggama1/.conda/envs/cnvkit/lib/python3.10/site-packages/cnvlib/commands.py", line 939, in _cmd_fix
    target_table = fix.do_fix(
  File "/home/ggama1/.conda/envs/cnvkit/lib/python3.10/site-packages/cnvlib/fix.py", line 24, in do_fix
    cnarr, ref_matched = load_adjust_coverages(
  File "/home/ggama1/.conda/envs/cnvkit/lib/python3.10/site-packages/cnvlib/fix.py", line 108, in load_adjust_coverages
    ref_matched = match_ref_to_sample(ref_cnarr, cnarr)
  File "/home/ggama1/.conda/envs/cnvkit/lib/python3.10/site-packages/cnvlib/fix.py", line 203, in match_ref_to_sample
    raise ValueError(
ValueError: Reference is missing 199 bins found in BH10027_1_TUMOR

When processing with batch, using the same reference, everything works fine.

GACGAMA avatar Jul 02 '24 22:07 GACGAMA