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?