-
Notifications
You must be signed in to change notification settings - Fork 18
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
genomic position #76
Comments
hi @acarmas1 , sorry for the delay in my reply. I believe you have replied to the issue in [here](#54 (comment) and you encountered an error which I am not quite sure about unless you share with me the function that you used and the full error log? Thanks! |
Hi Chris, no problem Also sorry for the late response. On this folder I added the code I run (gen_pos.sbatch), the error I got (out_genomic_position) and the files I used as inputs, which are in .fasta .gtf and .txt format. https://1drv.ms/u/s!AokqkR3muxL0jvIA82csbixN-ANRcQ?e=KWig3S I'll appreciate your help and thank you. |
hi @acarmas1 seems like a syntax error on gen_pos.py, specifically the last line, your dictionary is missing one curly braces at the end |
hi @acarmas1, may I know which gtf file you are using? It might be formatted differently |
Sure, on this link you can find the gtf file I'm using https://1drv.ms/u/s!AokqkR3muxL0jvIA82csbixN-ANRcQ?e=KWig3S |
Please can someone help me with this |
hi @acarmas1, sorry for the delay I looked at your GTF file and it seems that the problem arises because your gtf file format is a bit different than what we used to deal with. I modified your gen_pos.py, try this out, t2g_dict will give transcript_position -> genomic_position for every transcript in your all_transcripts.txt file `#!/usr/bin/env python
` |
hi @acarmas1 , yeah you encountered that error because the t2g dict does not have a closing bracket. Try
|
By running it with the closing brackets I didn't get any errors, but I also didn't get any new files or results. |
hi @acarmas1 , the script gives you a Python dictionary that provides a mapping between each transcriptomic position to a genomic position in t2g_dict variable. Will you be able to use this to map the position in your data yourself? I'm a bit packed with other projects recently so I might take a while to help with code implementation |
The problem is that after running the script I didn't get any t2g_dict file. Therefore, I was wondering if you could help me please. After getting the dictionary, yes run the mapping by myself. |
yeah the dictionary is not saved as a file, you can use it directly to map your position. If you want to save it, then you have to save the t2g_dict within that script into a file |
hi @acarmas1 if you're more familiar with R, you can also check this discussion thread |
Yes, I'm more familiar with R. I'll give it a try. Thanks for sharing the discussion thread. |
Hello,
I have used m6anet to identified the m6a modification in direct RNA reads of honey bees, and I got I table like this:
I added the column id_pos.
I want to see the distribution of my m6A sites across 5' UTR CDS and 3' UTR, I found the code used on your article, which is in the codeocean clapsule. The table used to create the figure needs a transcript annotations file and a .csv file with genomic position. However, when running m6anet the .csv file only has transcript_positions.
Is there anyway to obtain the genomic positions as well?, and also how can I create the transcript annotation file that contains the end of the 5' utr end of CDS and end tx.?
Thank you and I'm looking forward to hearing from you.
Best,
Camila
The text was updated successfully, but these errors were encountered: