15 老师好,我在利用在SNP设计dcasp标记时得到的primers.result.xls文件为空。

这是我的报错:

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

请先 登录 后评论

1 个回答

omicsgene - 生物信息
擅长:重测序,遗传进化,转录组,GWAS

看报错信息,可能输入文件格式有问题吧,建议对比例子的输入文件格式准备文件;

建议 截图说明一下,输入输出文件还有命令行都截图,才能更能好的给解决问题

请先 登录 后评论
  • 1 关注
  • 0 收藏,492 浏览
  • 马杰宇 提出于 2024-08-09 21:19

相似问题