Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 57,596 C→T 82.3% E143E (GAG→GAA rluA ← pseudouridine synthase for 23S rRNA (position 746) and tRNAphe(position 32)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE57,5960CT82.3% 38.3 / 2.1 28E143E (GAG→GAArluApseudouridine synthase for 23S rRNA (position 746) and tRNAphe(position 32)
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (8/15);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTG  >  minE/57531‑57652
                                                                 |                                                        
cGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTcc                                                       <  1:737419/69‑1 (MQ=255)
  cAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTccac                                                     <  1:122535/69‑1 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:1637666/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:2691874/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:2525607/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:2311421/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:2766412/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:116477/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:1078174/2‑69 (MQ=255)
     ctGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATATTCcaccac                                                  >  1:1074272/2‑69 (MQ=255)
                        aTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCCCCACTTCATATTCCGTCTGCGCa                               >  1:1789054/1‑69 (MQ=255)
                        aTCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCa                               >  1:2076364/1‑69 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:505752/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:3291445/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:937558/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:3093956/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:2826914/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:2812999/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:1056061/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:2302993/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:2051557/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:1854066/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:1762494/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:1344803/68‑3 (MQ=255)
                                   caccacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCACAGGTTTAcctg                     <  1:1381018/68‑3 (MQ=255)
                                    accacTCTTGCCGTGTTATCCGCCGCATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:2601480/67‑3 (MQ=255)
                                                     aTCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTg  <  1:3002424/69‑1 (MQ=255)
                                                          ccgcATATTCCACCACTTCATACTCCGTCTGCGCAGGTTTAcctg                     <  1:1344346/45‑3 (MQ=38)
                                                                 cTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTaa             >  1:2937057/1‑46 (MQ=255)
                                                                 cTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTaa             >  1:1323651/1‑46 (MQ=255)
                                                                 cTCAACCACTTCATATTCCGTCTGCGCAGGTTTACCAGTTTCGTaa             >  1:2559132/1‑46 (MQ=255)
                                                                 |                                                        
CGCAATTGATGCGAACGCCCGGTAATCGGTTTTAACACCACTCTTGCCGTGTTATCCGCCGCATACTCCACCACTTCATATTCCGTCTGCGCAGGTTTACCCGTTTCGTAACAGACTTTCTG  >  minE/57531‑57652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: