Entering edit mode
2.9 years ago
Begonia_pavonina
▴
140
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?