Entering edit mode
5.8 years ago
MatthewP
★
1.4k
Hello, everyone! I use RED-ML to analyse RNA Editing
events in my RNA Seq
data. The result data is simple:
#P_edit: The probability of being a RNA editing site predicted by RED-ML
#Chromosome Position Read_depth Reference Reference_support_reads Alternative
Alternative_support_reads P_edit
chr1 14604 21 A 6 G 15 0.849387564524627
chr1 14610 21 T 5 C 16 0.516906327292332
chr1 16280 24 T 7 C 17 0.889600786631089
chr1 17606 27 A 21 G 6 0.602832411634867
chr1 19172 12 A 8 G 4 0.823503600263782
chr1 185125 6 A 2 G 4 0.752470897435155
chr1 185131 6 T 2 C 4 0.647400412860666
chr1 185428 7 A 1 G 6 0.885004885559825
chr1 186365 63 T 40 C 23 0.801998990609502
chr1 187378 22 A 9 G 13 0.785096322895598
chr1 188891 9 T 4 C 5 0.803351489912644
chr1 629218 978 A 0 G 978 0.991090135056072
chr1 631193 44 A 4 G 40 0.94447283158624
chr1 632346 339 T 0 C 339 0.845537964553492
chr1 633630 10 T 1 C 9 0.605907924765477
chr1 633714 16 A 2 G 14 0.539893855272783
chr1 634112 1211 T 4 C 1207 0.803272961009867
However I don't have experience at filtering such data. Shall I filter by comparing with data from RNA Editing
database(DARNED and RADAR)? Do I need to consider read deapth
between Reference and Alternate? P_edit
>= 0.9? Tell me your experience thanks!