BackofenLab / CRISPRtracrRNA

7 stars 3 forks source link

Error with missing library but output still generated? #6

Open mhyleung opened 1 year ago

mhyleung commented 1 year ago

Dear all

I tried to run CRISPRtracrRNA with a contig containing a Class V Cas operon:

python CRISPRtracrRNA.py \
--input_folder [input_directory] \
--output_folder [output_directory] \
--run_type complete_run --model_type V --perform_type_v_anti_repeat_analysis True

I got the following output

Creating folders if necessary...
Running CRISPR array search, it may take a while...
Running CRISPR cas search, it may take a while...
        Working with 1 files

        Working with file 1 - contig01.fasta
1. Searching for anti repeats
*** error [comp_lib9.c:642] cannot open library [input_directory]_x
*** error [comp_lib9.c:642] cannot open library [input_directory]_x
2. Filtering Anti-repeats
3. Forming initial TracrRNA candidates
4. Improving the results using CRISPRRNA-TracrRNA interaction prediction
5. Searching for the terminator sequences
6. Searching the hits with the provided tail model
7. Writing results

        Creating the final_summary

        Adding consistency scores to the final summary file
Traceback (most recent call last):
  File "CRISPRtracrRNA.py", line 46, in <module>
    main()
  File "CRISPRtracrRNA.py", line 36, in main
    dict_weights=dict_weights)
  File "[CRISPRtracrRNA_directory]/modules/pipelines.py", line 35, in __init__
    self._pipeline_run()
  File "[CRISPRtracrRNA_directory]/modules/pipelines.py", line 316, in _pipeline_run
    consistency_score_maker(self.output_file_name, ",", self.output_file_name)
  File "[CRISPRtracrRNA_directory]/modules/consistency_score_maker.py", line 21, in consistency_score_maker
    interval_anti_repeat = get_interval_anti_repeat(index_anti_repeat, line, separator)
  File "[CRISPRtracrRNA_directory]/modules/consistency_score_maker.py", line 51, in get_interval_anti_repeat
    anti_repeat_interval = line_elements[anti_repeat_index]
IndexError: list index out of range

However, I am still able to get the output .csv file. I thought it was not able to predict the anti-repeat region, but the anti-repeat sequences are there in the final .csv file.

I just want to know what the basis of the error is and how I would go about removing it. Thanks

Marc