KeyError in the pipeline SECAPR
0
0
Entering edit mode
3.0 years ago

Hello everyone,

I am using the Hyb-Seq analysis pipeline SECAPR, and I have an error that I don't understand. I am using the tool for locus selection:

secapr locus_selection \                               
     --input remapped_reads \
     --output locus_selection/ \
     --n 2500 \
     --read_cov 3 \
     --reference remapped_reads/reference_seqs/joined_fasta_library.fasta

And I get a Python KeyError in the error report file:

Reading read-depth info for B.........
Generating locus database.........
Traceback (most recent call last):
  File "/mnt/shared/scratch/usr/apps/conda/envs/secapr_env/bin/secapr", line 10, in <module>
    sys.exit(main())
  File "/mnt/shared/scratch/usr/apps/conda/envs/secapr_env/lib/python3.8/site-packages/secapr/__main__.py", line 55, in main
    args.func(args)
  File "/mnt/shared/scratch/usr/apps/conda/envs/secapr_env/lib/python3.8/site-packages/secapr/locus_selection.py", line 312, in main
    reference_pickle = reference_file_dict[sample]
KeyError: 'B'

Would anyone have an idea about what could be the cause of the error?

SECAPR Hyb-Seq locus_selection Python Key_error • 436 views
ADD COMMENT

Login before adding your answer.

Traffic: 2980 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6