Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,785,373 C→A 41.2% intergenic (‑33/‑179) pepE ← / → yjbC (alpha)‑aspartyl dipeptidase/23S rRNA pseudouridine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,785,3730CA41.2% 49.9 / 70.4 80intergenic (‑33/‑179)pepE/yjbC(alpha)‑aspartyl dipeptidase/23S rRNA pseudouridine synthase
Reads supporting (aligned to +/- strand):  ref base C (35/12);  new base A (7/26);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 3.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 109 total aligned reads displayed.

CGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCGCACATTGCTTCTGGAAAGT  >  minE/2785314‑2785441
                                                           |                                                                    
cGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCa                                                             >  1:1821415/1‑69 (MQ=255)
cGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCa                                                             >  1:754565/1‑69 (MQ=255)
cGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCa                                                             >  1:1596166/1‑69 (MQ=255)
cGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCa                                                             >  1:2881242/1‑69 (MQ=255)
cGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCa                                                             >  1:1404784/1‑69 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:2899697/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:2025733/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:2569793/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:2863566/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:3189562/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:547865/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:1153526/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAAGCa                                                                     >  1:256545/1‑60 (MQ=255)
 gAGTTACTCAATAAAAGCAGTTCCATTACTTCTACGGTTAGTAAGAAATCACGGGAGGCa                                                                     >  1:2854424/1‑60 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:342175/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:435346/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:902913/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:536579/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:2623314/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:2588043/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:712364/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:738501/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:2428312/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:2082435/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:898188/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:1173568/68‑1 (MQ=255)
     tACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:1596985/68‑1 (MQ=255)
      aCTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGc                                                         <  1:1309382/67‑1 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAAttt                                                   <  1:58811/67‑1 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAAttt                                                   <  1:3364338/67‑1 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAAttt                                                   <  1:2146274/67‑1 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:2939806/1‑68 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:2980364/1‑68 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:2476939/1‑68 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:3197912/1‑68 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:222204/1‑68 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:2140679/1‑68 (MQ=255)
            tAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTa                                                  >  1:396744/1‑68 (MQ=255)
             aaaaGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAAc                                                <  1:1334038/69‑1 (MQ=255)
             aaaaGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAAc                                                <  1:2744240/69‑1 (MQ=255)
                 gCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACt                                            <  1:1938312/69‑1 (MQ=255)
                 gCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACt                                            <  1:181585/69‑1 (MQ=255)
                       ccATTACTTCTCCGGTTATTAAGAAATCACGGCAGGc                                                                      <  1:1587684/37‑1 (MQ=38)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:2245571/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:2833082/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:1933908/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:606509/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:617713/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:2582591/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:1944875/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:3110762/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:2416821/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGc                                                                      <  1:2110433/37‑1 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGa                                      >  1:1248816/1‑69 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGa                                      >  1:279287/1‑69 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGa                                      >  1:283254/1‑69 (MQ=255)
                       ccATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGa                                      >  1:1457854/1‑69 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:525080/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:1291568/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:19261/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:570192/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:57587/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:2769243/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:879560/69‑1 (MQ=255)
                           tACTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:181432/69‑1 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACt                                            >  1:828423/1‑58 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACt                                            >  1:3286308/1‑58 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACt                                            >  1:1956776/1‑58 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATgg                                 >  1:2654337/1‑69 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATgg                                 >  1:1530167/1‑69 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATgg                                 >  1:829811/1‑69 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATgg                                 >  1:1987999/1‑69 (MQ=255)
                            aCTTCTCCGGTTAGTAAGAAAACACGGAAGGAATGCTACCACTGCGAATTTAACCACTTACTGACATg                                  <  1:1505747/68‑1 (MQ=255)
                                       tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            <  1:110267/63‑1 (MQ=255)
                                       tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            <  1:3155152/63‑1 (MQ=255)
                                       tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            <  1:3033807/63‑1 (MQ=255)
                                       tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            <  1:2153913/63‑1 (MQ=255)
                                       tAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGt                       >  1:2868257/1‑68 (MQ=255)
                                       tAGTAAGAAATCACGGGAGGCATGCTACCACCGCGAATTTAACCACTTACTGACATGGCTCAc                            <  1:476644/63‑1 (MQ=255)
                                        aGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg                     >  1:3162397/1‑69 (MQ=255)
                                        aGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg                     >  1:1228954/1‑69 (MQ=255)
                                        aGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg                     >  1:1429317/1‑69 (MQ=255)
                                        aGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg                     >  1:1281604/1‑69 (MQ=255)
                                        aGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg                     >  1:1982243/1‑69 (MQ=255)
                                        aGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCg                     >  1:3172122/1‑69 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCa                             >  1:2585670/1‑42 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1652500/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1219811/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1939771/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:718120/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1239977/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:670219/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1502961/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:2646192/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1486041/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:1456128/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:301454/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:3086224/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:390791/1‑43 (MQ=255)
                                                           cATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCAc                            >  1:3105440/1‑43 (MQ=255)
                                                           |                                                                    
CGAGTTACTCAATAAAAGCAGTTCCATTACTTCTCCGGTTAGTAAGAAATCACGGGAGGCATGCTACCACTGCGAATTTAACCACTTACTGACATGGCTCACGATGTCGCACATTGCTTCTGGAAAGT  >  minE/2785314‑2785441

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