The error message shows can't find 1.fa in my outputdir. but I did not find anything wrong in install and script.
INFO 2019-11-01 09:32:46,750 /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/scripts/run_breakseq2.py Command-line: /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/scripts/run_breakseq2.py --reference /home/bpvast/01.test/1031_SV_test/breakseq_sv --bams /home/bpvast/01.test/1031_SV_test/PHY19092499.deduped.bam --work /home/bpvast/01.test/1031_SV_test/breakseq_sv --bwa bwa --samtools samtools --bplib_gff /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/breakseq2_bplib_20150129.gff --nthreads 12 --sample PHY19092499 INFO 2019-11-01 09:32:46,750 breakseq2_workflow Generating breakpoint-library using /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/breakseq2_bplib_20150129.gff Traceback (most recent call last): File "/home/bpvast/01.test/1031_SV_test/breakseq2-2.2/scripts/run_breakseq2.py", line 28, in <module> args.keep_temp, args.window, args.junction_length)) File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/breakseq_top.py", line 96, in breakseq2_workflow breakseq_index.generate_bplib(bplib_gff, reference, bplib, junction_length) File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/breakseq_index.py", line 47, in generate_bplib flanks = sv.get_flanks() File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/biopy/io/SV.py", line 120, in get_flanks seqA=seqs.get_sequence(*seqs.get_window(self.name, self.start)) File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/biopy/io/Fasta.py", line 16, in get_sequence self.seqs[name] = parse(self.base+"/"+name+".fa").next() File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/biopy/io/Fasta.py", line 31, in parse handle = open(file, "rU") IOError: [Errno 2] No such file or directory: '/home/bpvast/01.test/1031_SV_test/breakseq_sv/1.fa'
The error message shows can't find 1.fa in my outputdir. but I did not find anything wrong in install and script.
INFO 2019-11-01 09:32:46,750 /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/scripts/run_breakseq2.py Command-line: /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/scripts/run_breakseq2.py --reference /home/bpvast/01.test/1031_SV_test/breakseq_sv --bams /home/bpvast/01.test/1031_SV_test/PHY19092499.deduped.bam --work /home/bpvast/01.test/1031_SV_test/breakseq_sv --bwa bwa --samtools samtools --bplib_gff /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/breakseq2_bplib_20150129.gff --nthreads 12 --sample PHY19092499 INFO 2019-11-01 09:32:46,750 breakseq2_workflow Generating breakpoint-library using /home/bpvast/01.test/1031_SV_test/breakseq2-2.2/breakseq2_bplib_20150129.gff Traceback (most recent call last): File "/home/bpvast/01.test/1031_SV_test/breakseq2-2.2/scripts/run_breakseq2.py", line 28, in <module> args.keep_temp, args.window, args.junction_length)) File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/breakseq_top.py", line 96, in breakseq2_workflow breakseq_index.generate_bplib(bplib_gff, reference, bplib, junction_length) File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/breakseq_index.py", line 47, in generate_bplib flanks = sv.get_flanks() File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/biopy/io/SV.py", line 120, in get_flanks seqA=seqs.get_sequence(*seqs.get_window(self.name, self.start)) File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/biopy/io/Fasta.py", line 16, in get_sequence self.seqs[name] = parse(self.base+"/"+name+".fa").next() File "/home/bpvast/.local/lib/python2.7/site-packages/breakseq2/biopy/io/Fasta.py", line 31, in parse handle = open(file, "rU") IOError: [Errno 2] No such file or directory: '/home/bpvast/01.test/1031_SV_test/breakseq_sv/1.fa'