metasv icon indicating copy to clipboard operation
metasv copied to clipboard

KeyError: 'Assuming'

Open Qijie0615 opened this issue 2 years ago • 0 comments

Error after running Metasv;can you give me some suggestion?

$METASV --pindel_native /bak01/yangqj/pindel/20220116/vcf_format/ps0001/ps0001_* --cnvnator_native /bak01/yangqj/cnvnator/ps0001.call --reference /bak01/yangqj/Metasv/hg19.fa --outdir out --sample ps0001 --filter_gaps --minsvlen 500 --maxsvlen 500000 --disable_assembly --keep_standard_contigs

INFO 2022-01-25 16:34:31,082 metasv.main Running MetaSV 0.5.4 INFO 2022-01-25 16:34:31,082 metasv.main Command-line /lustre/yangqj/software/miniconda3/envs/metasv/bin/run_metasv.py --pindel_native /bak01/yangqj/pindel/20220116/vcf_format/ps0001/ps0001_DEL.vcf /bak01/yangqj/pindel/20220116/vcf_format/ps0001/ps0001_TD.vcf --cnvnator_native /bak01/yangqj/cnvnator/ps0001.call --reference /bak01/yangqj/Metasv/hg19.fa --outdir out --sample ps0001 --filter_gaps --minsvlen 500 --maxsvlen 500000 --disable_assembly --keep_standard_contigs INFO 2022-01-25 16:34:31,083 metasv.main Arguments are Namespace(age=None, age_timeout=300, age_window=20, assembly_max_tools=1, assembly_pad=500, bams=[], boost_sc=False, breakdancer_native=[], breakdancer_vcf=[], breakseq_native=[], breakseq_vcf=[], chromosomes=[], cnvkit_vcf=[], cnvnator_native=['/bak01/yangqj/cnvnator/ps0001.call'], cnvnator_vcf=[], disable_assembly=True, enable_per_tool_output=False, extraction_max_read_pairs=10000, filter_gaps=True, gaps=None, gatk_vcf=[], gt_normal_frac=0.05, gt_window=100, inswiggle=100, isize_mean=350.0, isize_sd=50.0, keep_standard_contigs=True, lumpy_vcf=[], manta_vcf=[], max_ins_cov_frac=1.5, max_ins_intervals=10000, max_nm=10, maxsvlen=500000, mean_read_coverage=50, mean_read_length=100, min_avg_base_qual=20, min_del_subalign_len=50, min_ins_cov_frac=0.5, min_inv_subalign_len=50, min_mapq=5, min_matches=50, min_soft_clip=20, min_support_frac_ins=0.05, min_support_ins=15, minsvlen=500, num_threads=1, outdir='out', overlap_ratio=0.5, pindel_native=['/bak01/yangqj/pindel/20220116/vcf_format/ps0001/ps0001_DEL.vcf', '/bak01/yangqj/pindel/20220116/vcf_format/ps0001/ps0001_TD.vcf'], pindel_vcf=[], reference='/bak01/yangqj/Metasv/hg19.fa', sample='ps0001', sc_other_scale=5, spades=None, spades_max_interval_size=50000, spades_options='', spades_timeout=300, stop_spades_on_fail=False, svs_to_assemble=set(['DUP', 'INV', 'INS']), svs_to_report=set(['INV', 'CTX', 'INS', 'DEL', 'ITX', 'DUP']), svs_to_softclip=set(['DUP', 'INV', 'INS']), wham_vcf=[], wiggle=100, workdir='work') INFO 2022-01-25 16:34:31,084 metasv.main Only SVs on the following contigs will be reported: ['chr1', 'chr10', 'chr11', 'chr12', 'chr13', 'chr14', 'chr15', 'chr16', 'chr17', 'chr18', 'chr19', 'chr2', 'chr20', 'chr21', 'chr22', 'chr3', 'chr4', 'chr5', 'chr6', 'chr7', 'chr8', 'chr9', 'chrM', 'chrX', 'chrY'] INFO 2022-01-25 16:34:31,084 metasv.sv_interval Loading the gaps in the genome from /lustre/yangqj/software/miniconda3/envs/metasv/lib/python2.7/site-packages/metasv/resources/hg19.gaps.bed INFO 2022-01-25 16:34:31,117 metasv.main Load native files INFO 2022-01-25 16:34:31,117 metasv.cnvnator_reader File is /bak01/yangqj/cnvnator/ps0001.call Traceback (most recent call last): File "/lustre/yangqj/software/miniconda3/envs/metasv/bin/run_metasv.py", line 143, in sys.exit(run_metasv(args)) File "/lustre/yangqj/software/miniconda3/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 "/lustre/yangqj/software/miniconda3/envs/metasv/lib/python2.7/site-packages/metasv/cnvnator_reader.py", line 123, in next record = CNVnatorRecord(line.strip()) File "/lustre/yangqj/software/miniconda3/envs/metasv/lib/python2.7/site-packages/metasv/cnvnator_reader.py", line 38, in init self.sv_type = sv_type_dict[fields[0]] KeyError: 'Assuming'

Qijie0615 avatar Jan 25 '22 08:01 Qijie0615