-
Notifications
You must be signed in to change notification settings - Fork 36
Description
Hi,
An "KeyError" happens when I try to run RAT. Here're my code and the last part of the output.
(cat_pack) [shennaijia@covid ~]$ CAT_pack reads -c /home/shennaijia/sewage_analyze/MH/1_siass/results/MH_220920_DNA_R1.contigs.fa
-t /home/shennaijia/database/CAT_pack/GTDB/GTDB_CAT_pack/tax
-m cr
-1 /home/shennaijia/sewage_analyze/MH/1_coass/data/MH_220920_DNA_R1_1.fastq
-2 /home/shennaijia/sewage_analyze/MH/1_coass/data/MH_220920_DNA_R1_2.fastq
--c2c /home/shennaijia/sewage_analyze/MH/3_classification/CAT_siass/0_rawresult/MH_220920_DNA_R1/MH_220920_DNA_R1.contig2classification.txt
-o /home/shennaijia/sewage_analyze/MH/3_classification/RAT_siass/0_rawresult/MH_220920_DNA_R1/MH_220920_DNA_R1
[2025-06-30 21:42:58] Chosen mode: cr. Classifying unclassified contigs and unmapped reads with diamond if no classification file is supplied.
[2025-06-30 21:42:58] No unmapped2classification file supplied. Grabbing unmapped and unclassified sequences...
[2025-06-30 21:42:59] Contigs written! Appending forward reads...
Traceback (most recent call last):
File "/home/shennaijia/software/CAT_pack/CAT_pack/CAT_pack", line 101, in
main()
File "/home/shennaijia/software/CAT_pack/CAT_pack/CAT_pack", line 85, in main
reads.run()
File "/home/shennaijia/software/CAT_pack/CAT_pack/reads.py", line 438, in run
make_unclassified_seq_fasta(reads_files[0], unmapped_reads['fw'],
File "/home/shennaijia/software/CAT_pack/CAT_pack/reads.py", line 1263, in make_unclassified_seq_fasta
fasta_dict[seq]))
KeyError: 'A00151:688:HLKY2DSX3:1:1628:31430:21230'
I check my .fq files, here're the details:
(cat_pack) [shennaijia@covid ~]$ grep 'A00151:688:HLKY2DSX3:1:1628:31430:21230' /home/shennaijia/sewage_analyze/MH/1_coass/data/MH_220920_DNA_R1_1.fastq
@A00151:688:HLKY2DSX3:1:1628:31430:21230/1
(cat_pack) [shennaijia@covid ~]$ grep 'A00151:688:HLKY2DSX3:1:1628:31430:21230' /home/shennaijia/sewage_analyze/MH/1_coass/data/MH_220920_DNA_R1_2.fastq
@A00151:688:HLKY2DSX3:1:1628:31430:21230/2
It seems that the program can't identify the sequence names correctly.
Sorry for my poor English.
Thanks,
Shen