Hello there,
I have been running the RILseq package on M. smegmatis RIL samples but have been running into issues on the calling S-chimera step. I downloaded the Biocyc files and stored the genes.dat, promoters.dat, terminators.dat, the genes.gff & transcripts.gff files that were made, and all *.fsa files into my ecocycdata directory. I ensured the chromosome name would be the same from bam to these files (which is NC_018289).
When I run: RILseq_significant_regions.py --bc_dir /Users/yfan/Desktop/RILRawFilesMsmeg/Attempt2/ecocycdata --ribozero --min_odds_ratio 1 HA2_all_fragments.txt > HA2_significant_interactions.txt
I get this error:
Traceback (most recent call last):
File "/Users/yfan/opt/anaconda3/bin/RILseq_significant_regions.py", line 305, in
status = main()
File "/Users/yfan/opt/anaconda3/bin/RILseq_significant_regions.py", line 177, in main
rr_pos.append([chr_dict[uid_pos[rrgene][0]]] +
KeyError: 'NC_018289'
I am unsure how to fix this issue and was hoping to get your input. Thank you!
Hello there,
I have been running the RILseq package on M. smegmatis RIL samples but have been running into issues on the calling S-chimera step. I downloaded the Biocyc files and stored the genes.dat, promoters.dat, terminators.dat, the genes.gff & transcripts.gff files that were made, and all *.fsa files into my ecocycdata directory. I ensured the chromosome name would be the same from bam to these files (which is NC_018289).
When I run: RILseq_significant_regions.py --bc_dir /Users/yfan/Desktop/RILRawFilesMsmeg/Attempt2/ecocycdata --ribozero --min_odds_ratio 1 HA2_all_fragments.txt > HA2_significant_interactions.txt
I get this error:
Traceback (most recent call last):
File "/Users/yfan/opt/anaconda3/bin/RILseq_significant_regions.py", line 305, in
status = main()
File "/Users/yfan/opt/anaconda3/bin/RILseq_significant_regions.py", line 177, in main
rr_pos.append([chr_dict[uid_pos[rrgene][0]]] +
KeyError: 'NC_018289'
I am unsure how to fix this issue and was hoping to get your input. Thank you!