Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,137,325 G→A 49.2% T182T (ACC→ACT truB ← tRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,137,3250GA49.2% ‑2.0 / 59.0 67T182T (ACC→ACTtruBtRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base G (3/31);  new base A (18/15);  total (21/46)
Fisher's exact test for biased strand distribution p-value = 5.95e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCG  >  minE/2137264‑2137387
                                                             |                                                              
gACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGa                                                         <  1:2423797/69‑1 (MQ=255)
gACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGa                                                         <  1:968483/69‑1 (MQ=255)
   ggCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                                                      <  1:1082734/69‑1 (MQ=255)
    gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                                                      <  1:2261311/68‑1 (MQ=255)
    gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGatat                                                      <  1:1746818/68‑1 (MQ=255)
    gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                                     <  1:2297205/69‑1 (MQ=255)
    gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                                     <  1:2173170/69‑1 (MQ=255)
    gcgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                                     <  1:166751/69‑1 (MQ=255)
     cgcAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGAtata                                                     <  1:1760907/68‑1 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:355264/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:2788442/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:3134899/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1099491/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:3160136/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1292105/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:3308389/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1338693/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:194312/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:582229/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1886243/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1876516/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1363783/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:1453502/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:658112/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTGCGGATATAAGTGCCTTTTg                                          >  1:749617/1‑69 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCGTCAATGATAGTCCGGATATAAGTCCCtttt                                           >  1:611783/1‑68 (MQ=255)
               tAACATGCGCGCCACAACCGAGTTTTTCGCCGAGATCCTCAATGATAGTGCGGATAAAAGTGCCTTTTg                                          >  1:3351491/1‑69 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   <  1:882267/69‑1 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   >  1:991624/1‑69 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   <  1:315293/69‑1 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   >  1:1065793/1‑69 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   >  1:3088524/1‑69 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   <  1:793578/69‑1 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   <  1:1095463/69‑1 (MQ=255)
                      cgcgcCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTg                                   <  1:2206041/69‑1 (MQ=255)
                             caGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc                            <  1:2040328/69‑1 (MQ=255)
                             caGCCGAGTTTTTCACCCAGGGCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc                            <  1:2767693/69‑1 (MQ=255)
                                ccGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtccagc                         <  1:890786/69‑1 (MQ=255)
                                ccGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtccagc                         <  1:3264520/69‑1 (MQ=255)
                                           aCCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCAt               <  1:188888/68‑1 (MQ=255)
                                               aGGTCATCAATGATGGTGCGGTTATAAGTGCCTTTTGAGCAGTGTTTTtcc                            <  1:437214/51‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCTGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:2847361/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:1700469/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:2747634/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:954193/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:860468/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:1597362/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:3075196/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:2666655/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:2704888/69‑1 (MQ=255)
                                                ggTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCtt         <  1:2294807/69‑1 (MQ=255)
                                                 gTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTtcc                            <  1:2369379/49‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:855015/66‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:1374887/66‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:237570/66‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:1551478/66‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:755335/66‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:1636474/66‑1 (MQ=255)
                                                  tcgtcaATGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCa       <  1:652050/66‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:3187561/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:444869/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:861462/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:1525028/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:884066/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:611364/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:989011/69‑1 (MQ=255)
                                                       aaTGATAGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:540902/69‑1 (MQ=255)
                                                          gataGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCg  <  1:2005173/66‑1 (MQ=255)
                                                             |                                                              
GACGGCGCAGGTAAATAACATGCGCGCCACAGCCGAGTTTTTCACCCAGGTCATCAATGATGGTGCGGATATAAGTGCCTTTTGAGCAGTGAATTTCCAGCTCCAGCTCATTGCCTTCATGGCG  >  minE/2137264‑2137387

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: