Question about REDItools
Entering edit mode
4.1 years ago
tujuchuanli ▴ 100


I wanted to use REDItools to identify RNA editing events in my data. I had both RNA-seq and exome-seq data derived from the same sample.

I download the latest version of REDItools (, and test data (

When I run the scripts “” as “ -i rna.bam -j dna.bam -f reference.fa -o reditool-test -c 10,1 -q 25,25 -m 20,20 -s 2 -g 1 -u -a 6-0 -v 2 -n0.0 -N0.0 -V” basically following the README in test data folder to test the running. It returned err message as following:

Pysam version used: 0.15.3

Script time --> START: 02/09/2019 23:18:36

Analysis ID: 464424299

Analysis on 1 regions.

Started analysis on region: chr21:1-48129895

Process Process-1:

Traceback (most recent call last):

File "/usr/lib/python2.7/multiprocessing/", line 267, in _bootstrap

File "/usr/lib/python2.7/multiprocessing/", line 114, in run self._target(self._args, *self._kwargs)

File "", line 1177, in do_work exploreBAM(x)

File "", line 1121, in exploreBAM seq,qual,squal=normByStrand(seq,strand,squal,mystrand)

File "", line 621, in normByStrand if type(i)==int: qual+=ord(squal_[str(i)])-QVAL

TypeError: list indices must be integers, not str

Merging Tables.

Results saved on reditool-test/DnaRna_464424299/outTable_464424299

Script time --> END: 02/09/2019 23:18:38

Then I found that it could be no err message returned when I deleted “-s 2” from the command line. However, the output file could be empty this time.

Could you help me? Thanks.

REDItools • 1.6k views
Entering edit mode

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.

Entering edit mode

Excuse me, recently I've had the exactly same problem as you had 4 years ago......Could you please tell me how you dealt with it before, if you still remember? It would be very helpful to me! Thanks a lot!


Login before adding your answer.

Traffic: 1733 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6