metasv
metasv copied to clipboard
ValueError: dictionary update sequence element #0 has length 1; 2 is required
Hi,
I am trying a complete run of MetaSV 0.5.4 (installed from bioconda) using all 2 SV detectors (pindel and breakdancer), soft-clips based analysis, and local assembly. And got this error
INFO 2022-02-15 15:51:33,585 metasv.main Load native files
INFO 2022-02-15 15:51:33,586 metasv.pindel_reader File is LA3111t13-LA4330t13_D
INFO 2022-02-15 16:18:19,892 metasv.pindel_reader File is LA3111t13-LA4330t13_SI
INFO 2022-02-15 16:32:20,639 metasv.pindel_reader File is LA3111t13-LA4330t13_TD
INFO 2022-02-15 16:34:33,923 metasv.pindel_reader File is LA3111t13-LA4330t13_INV
INFO 2022-02-15 16:34:35,183 metasv.breakdancer_reader File is breakdancer.sv.vcf
Traceback (most recent call last):
File "/data/home/users/g.silvaarias/anaconda3/envs/metasv/bin/run_metasv.py", line 143, in <module>
sys.exit(run_metasv(args))
File "/data/home/users/g.silvaarias/anaconda3/envs/metasv/lib/python2.7/site-packages/metasv/main.py", line 106, in run_metasv
for record in svReader(native_file, svs_to_report=args.svs_to_report):
File "/data/home/users/g.silvaarias/anaconda3/envs/metasv/lib/python2.7/site-packages/metasv/breakdancer_reader.py", line 222, in next
self.header.parse_header_line(line)
File "/data/home/users/g.silvaarias/anaconda3/envs/metasv/lib/python2.7/site-packages/metasv/breakdancer_reader.py", line 74, in parse_header_line
self.header_dict[fields[0]] = dict(field.split(":") for field in fields[1:])
ValueError: dictionary update sequence element #0 has length 1; 2 is required
Here is the full command:
run_metasv.py --reference $ref \
--outdir $outdir \
--boost_sc \
--breakdancer_native breakdancer.sv.vcf \
--pindel_native LA3111t13-LA4330t13_D LA3111t13-LA4330t13_SI LA3111t13-LA4330t13_TD LA3111t13-LA4330t13_INV \
--sample LA3111t13 --sample LA4330t13 --bam LA3111t13_dedup_RG.bam LA4330t13_dedup_RG.bam --spades spades.py \
--age age_align --num_threads $threads \
--min_support_ins 10 --isize_mean 400 --isize_sd 100
I would appreciate any suggestion to fix that.
Best, Gustavo