看报错信息,可能输入文件格式有问题吧,建议对比例子的输入文件格式准备文件;
建议 截图说明一下,输入输出文件还有命令行都截图,才能更能好的给解决问题
这是我的报错:
rm: cannot remove ‘/work/out/tmp_in/*’: No such file or directory
Options done
Output selected CAPS file name is: CAPS_output/selected_CAPS_primers.txt
ids is []
SNP_A, SNP_B G A
CAPS found with enzyme BstNI,PspGI,19
One nt can be changed to fit enzyme NlaIII,CviAII,FatI,132
change position and primer end postions are 298 [299]
CAPS found with enzyme ScrFI,StyD4I,66
CAPS found with enzyme FspEI,530
One nt can be changed to fit enzyme MspJI,530
change position and primer end postions are 297 [298, 299]
CAPS found with enzyme AvrII,680
One nt can be changed to fit enzyme BseYI,660
change position and primer end postions are 297 [298, 299]
One nt can be changed to fit enzyme AlwNI,132
change position and primer end postions are 294 [295, 296, 297, 298, 299]
One nt can be changed to fit enzyme DdeI,66
change position and primer end postions are 298 [299]
CAPS found with enzyme LpnPI,530
One nt can be changed to fit enzyme DraIII,67
change position and primer end postions are 298 [299]
One nt can be changed to fit enzyme BsaJI,66
change position and primer end postions are 296 [297, 298, 299]
CAPS found with enzyme BfaI,142
One nt can be changed to fit enzyme BslI,63
change position and primer end postions are 291 [292, 293, 294, 295, 296, 297, 298, 299]
CAPS found with enzyme StyI,23
One nt can be changed to fit enzyme BpmI,660
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme BsrI,63
change position and primer end postions are 297 [298, 299]
One nt can be changed to fit enzyme HphI,63
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme BmrI,660
change position and primer end postions are 297 [298, 299]
One nt can be changed to fit enzyme PflMI,63
change position and primer end postions are 291 [292, 293, 294, 295, 296, 297, 298, 299]
One nt can be changed to fit enzyme Bpu10I,330
change position and primer end postions are 298 [299]
caps_list is ['BstNI,PspGI,19', 'ScrFI,StyD4I,66', 'FspEI,530', 'AvrII,680', 'LpnPI,530', 'BfaI,142', 'StyI,23']
dcaps_list is ['NlaIII,CviAII,FatI,132', 'MspJI,530', 'BseYI,660', 'AlwNI,132', 'DdeI,66', 'DraIII,67', 'BsaJI,66', 'BslI,63', 'BpmI,660', 'BsrI,63', 'HphI,63', 'BmrI,660', 'PflMI,63', 'Bpu10I,330']
Alignment command: /work/script/SNP_Primer_Pipeline-master-old/bin/muscle -in /work/out/tmp_in/chrA1:2590051.fa -out alignment_raw.fa -quiet
Primer3 command 1st time: /biosoft/miniconda/envs/reseq/bin/primer3_core -default_version=2 -output=CAPS_output/primer3.output.SNP CAPS_output/primer3.input.SNP
/bin/sh: /biosoft/miniconda/envs/reseq/bin/primer3_core: No such file or directory
Traceback (most recent call last):
File "/work/script/SNP_Primer_Pipeline-master-old/bin/getCAPS-with-user-input.py", line 1027, in <module>
caps(seqfile, target, SNP_A, SNP_B, snp_pos, max_price)
File "/work/script/SNP_Primer_Pipeline-master-old/bin/getCAPS-with-user-input.py", line 776, in caps
primerpairs = parse_primer3output(primer3output, 5)
File "/work/script/SNP_Primer_Pipeline-master-old/bin/getCAPS-with-user-input.py", line 508, in parse_primer3output
with open(primer3output) as infile:
IOError: [Errno 2] No such file or directory: 'CAPS_output/primer3.output.SNP'
open /work/out/CAPS_output/selected_CAPS_primers.txt failed
mv: cannot stat ‘/work/out/CAPS_output/selected_CAPS_primers.txt’: No such file or directory
Options done
Output selected CAPS file name is: CAPS_output/selected_CAPS_primers.txt
ids is []
SNP_A, SNP_B C T
One nt can be changed to fit enzyme Hpy99I,660
change position and primer end postions are 304 [301, 302, 303]
CAPS found with enzyme FspEI,530
One nt can be changed to fit enzyme NmeAIII,272
change position and primer end postions are 303 [301, 302]
One nt can be changed to fit enzyme BsiEI,66
change position and primer end postions are 296 [297, 298, 299]
One nt can be changed to fit enzyme ApeKI,TseI,272
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme MmeI,660
change position and primer end postions are 298 [299]
One nt can be changed to fit enzyme AciI,330
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme DdeI,66
change position and primer end postions are 303 [301, 302]
One nt can be changed to fit enzyme LpnPI,530
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme BspCNI,630
change position and primer end postions are 303 [301, 302]
CAPS found with enzyme HaeIII,20
One nt can be changed to fit enzyme MspI,HpaII,13
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme BsrFαI,63
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme BslI,63
change position and primer end postions are 308 [301, 302, 303, 304, 305, 306, 307]
One nt can be changed to fit enzyme Sau96I,61
change position and primer end postions are 296 [297, 298, 299]
One nt can be changed to fit enzyme AluI,66
change position and primer end postions are 297 [298, 299]
One nt can be changed to fit enzyme NgoMIV,NaeI,63
change position and primer end postions are 302 [301]
CAPS found with enzyme EaeI,315
One nt can be changed to fit enzyme BbvI,220
change position and primer end postions are 302 [301]
One nt can be changed to fit enzyme EagI,126
change position and primer end postions are 296 [297, 298, 299]
caps_list is ['FspEI,530', 'HaeIII,20', 'EaeI,315']
dcaps_list is ['Hpy99I,660', 'NmeAIII,272', 'BsiEI,66', 'ApeKI,TseI,272', 'MmeI,660', 'AciI,330', 'DdeI,66', 'LpnPI,530', 'BspCNI,630', 'MspI,HpaII,13', 'BsrF\xce\xb1I,63', 'BslI,63', 'Sau96I,61', 'AluI,66', 'NgoMIV,NaeI,63', 'BbvI,220', 'EagI,126']
Alignment command: /work/script/SNP_Primer_Pipeline-master-old/bin/muscle -in /work/out/tmp_in/chrA1:27536272.fa -out alignment_raw.fa -quiet
Primer3 command 1st time: /biosoft/miniconda/envs/reseq/bin/primer3_core -default_version=2 -output=CAPS_output/primer3.output.SNP CAPS_output/primer3.input.SNP
/bin/sh: /biosoft/miniconda/envs/reseq/bin/primer3_core: No such file or directory
Traceback (most recent call last):
File "/work/script/SNP_Primer_Pipeline-master-old/bin/getCAPS-with-user-input.py", line 1027, in <module>
caps(seqfile, target, SNP_A, SNP_B, snp_pos, max_price)
File "/work/script/SNP_Primer_Pipeline-master-old/bin/getCAPS-with-user-input.py", line 776, in caps
primerpairs = parse_primer3output(primer3output, 5)
File "/work/script/SNP_Primer_Pipeline-master-old/bin/getCAPS-with-user-input.py", line 508, in parse_primer3output
with open(primer3output) as infile:
IOError: [Errno 2] No such file or directory: 'CAPS_output/primer3.output.SNP'
open /work/out/CAPS_output/selected_CAPS_primers.txt failed
mv: cannot stat ‘/work/out/CAPS_output/selected_CAPS_primers.txt’: No such file or directory
perl /work/script/dcaps_primer_result.pl -id /work/out/CAPS_output -out /work/out/primers.result.xls