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

reditools.py error #10

Open
hzongyao opened this issue Sep 30, 2021 · 0 comments
Open

reditools.py error #10

hzongyao opened this issue Sep 30, 2021 · 0 comments

Comments

@hzongyao
Copy link

Dear @tizianoflati @tflati

I am very keen to use REDItools2, which seems a fantastic tool for analysing RNA editing with very useful instructions. But I am getting this error :
('[SYSTEM]', 'PYSAM VERSION', '0.16.0.1')
('[SYSTEM]', 'PYSAM PATH', ['/home/tmp/tmp/anaconda3/envs/redtools/lib/python2.7/site-packages/pysam'])
[INFO] [None] START=2021-10-01 00:28:37.051437
[INFO] Opening BAM file=/home/tmp/tmp/project/tmp/GATGTAACGAACCACC.bam
[INFO] Fetching data from bam /home/tmp/tmp/project/tmp/GATGTAACGAACCACC.bam
[INFO] Narrowing REDItools to region None
Opening bamfile within region=None
[INFO] REFERENCE NAME=chr1 (2021-10-01 00:28:38.684963) [00:00:01]
Traceback (most recent call last):
File "/mnt/raid/tmp/software/biosoft/REDItools2/src/cineca/reditools.py", line 1374, in
analyze(options)
File "/mnt/raid/tmp/software/biosoft/REDItools2/src/cineca/reditools.py", line 1131, in analyze
pos_based_read_dictionary = update_reads(reads, i)
File "/mnt/raid/tmp/software/biosoft/REDItools2/src/cineca/reditools.py", line 149, in update_reads
read["qual"] = read["query_qualities"][read["alignment_index"]]
IndexError: list index out of range

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