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

KeyError: 'Assuming' #136

Open
Qijie0615 opened this issue Jan 25, 2022 · 0 comments
Open

KeyError: 'Assuming' #136

Qijie0615 opened this issue Jan 25, 2022 · 0 comments

Comments

@Qijie0615
Copy link

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 Qijie0615 reopened this Jan 25, 2022
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