Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Question on RILseq_significant_regions.py step #23

Open
fany8 opened this issue Feb 13, 2024 · 0 comments
Open

Question on RILseq_significant_regions.py step #23

fany8 opened this issue Feb 13, 2024 · 0 comments

Comments

@fany8
Copy link

fany8 commented Feb 13, 2024

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!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant