linsalrob / PhiSpy

Prediction of prophages from bacterial genomes
MIT License
70 stars 20 forks source link

bug report #5

Closed HaidYi closed 4 years ago

HaidYi commented 6 years ago

Hi maintainer,

when I use the phispy to detect the prophage of NZ_CP014862.1, I got the error like this.

Making Test Set... (need couple of minutes)
Start Classification Algorithm
Using training flag:  0
Done with classification Algorithm
As training flag is zero, considering unknown functions
Start evaluation...
Traceback (most recent call last):
  File "/home/ciillab/huangle/phiSpy.py", line 150, in <module>
    start_propgram(sys.argv)
  File "/home/ciillab/huangle/phiSpy.py", line 148, in start_propgram
    call_phiSpy(organismPath,output_dir,trainingFlag,INSTALLATION_DIR,args_parser.evaluate, args_parser.number, args_parser.window_size,args_parser.quiet)
  File "/home/ciillab/huangle/phiSpy.py", line 40, in call_phiSpy
    evaluation.call_start_end_fix(output_dir,organismPath,INSTALLATION_DIR,threshold_for_FN, phageWindowSize)
  File "/home/ciillab/huangle/source/evaluation.py", line 528, in call_start_end_fix
    fixing_start_end(output_dir+'initial_tbl.txt',output_dir+'prophage_tbl_temp.txt',organismPath,INSTALLATION_DIR)
  File "/home/ciillab/huangle/source/evaluation.py", line 352, in fixing_start_end
    repeat_list = find_repeat(dna[pp[i]['contig']][start:stop],start,INSTALLATION_DIR)
  File "/home/ciillab/huangle/source/evaluation.py", line 90, in find_repeat
    if math.fabs(int(temp[0]) - int(temp[3])) > 10000:
IndexError: list index out of range

Please address this problem, thank you.