agshumate/Liftoff

Error - polish.sam empty

KristinaGagalova opened this issue · 1 comments

Hi,
I am running liftoff with the polish option and I am getting the following error at the polishing step. Do you have any idea why that may happen? Thank you in advance

[M::main] Real time: 640.470 sec; CPU: 640.450 sec; Peak RSS: 9.110 GB
aligning features
lifting features
polishing annotations
Traceback (most recent call last):
  File "/projects/spruceup_scratch/dev/anaconda3/bin/liftoff", line 11, in <module>
    load_entry_point('Liftoff==1.6.1', 'console_scripts', 'liftoff')()
  File "/projects/spruceup_scratch/dev/anaconda3/lib/python3.7/site-packages/liftoff/run_liftoff.py", line 10, in main
    run_all_liftoff_steps(args)
  File "/projects/spruceup_scratch/dev/anaconda3/lib/python3.7/site-packages/liftoff/run_liftoff.py", line 42, in run_all_liftoff_steps
    unmapped_features, feature_db, ref_parent_order)
  File "/projects/spruceup_scratch/dev/anaconda3/lib/python3.7/site-packages/liftoff/run_liftoff.py", line 240, in find_and_polish_broken_cds
    if polish.polish_annotations(lifted_feature_list, ref_fa, target_fa, args, feature_hierarchy, target_feature):
  File "/projects/spruceup_scratch/dev/anaconda3/lib/python3.7/site-packages/liftoff/polish.py", line 24, in polish_annotations
    target_faidx,output_sam)
  File "/projects/spruceup_scratch/dev/anaconda3/lib/python3.7/site-packages/liftoff/polish.py", line 175, in polish_annotation
    target_interval = get_target_interval(exon_group_dict[i], target_exons, ref_interval)
  File "/projects/spruceup_scratch/dev/anaconda3/lib/python3.7/site-packages/liftoff/polish.py", line 236, in get_target_interval
    start = min([exon.start for exon in target_exons]) - (ref_interval[1]-ref_interval[0])
ValueError: min() arg is an empty sequence

hi,
i believe this should be fixed with the latest release (1.6.2). let me know if you still have troubles