Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 44,126 C→A 32.3% I214I (ATC→ATA yaaU → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE44,1260CA32.3% 84.3 / 44.0 71I214I (ATC→ATAyaaUpredicted transporter
Reads supporting (aligned to +/- strand):  ref base C (29/19);  new base A (13/10);  total (42/29)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

ACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGT  >  minE/44060‑44193
                                                                  |                                                                   
aCTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCa                                                                    <  1:1971456/68‑1 (MQ=255)
aCTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCa                                                                    <  1:2065197/68‑1 (MQ=255)
aCTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCa                                                                    <  1:3237361/68‑1 (MQ=255)
   gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACt                                                                >  1:2771980/1‑69 (MQ=255)
   gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACt                                                                >  1:1738826/1‑69 (MQ=255)
   gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACt                                                                >  1:794693/1‑69 (MQ=255)
   gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACt                                                                >  1:1395598/1‑69 (MQ=255)
   gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACt                                                                >  1:1369381/1‑69 (MQ=255)
   gCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACt                                                                >  1:2547116/1‑69 (MQ=255)
      tGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTg                                                               <  1:2236350/67‑1 (MQ=255)
      tGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTg                                                               <  1:3208249/67‑1 (MQ=255)
      tGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTg                                                               <  1:2429800/67‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:2165958/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:2387830/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:1930755/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:1771297/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:2729242/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:945642/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:1498685/68‑1 (MQ=255)
                cGCTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:781007/68‑1 (MQ=255)
                 gCTGGTTGTTACGCAAAGGGCGAGTGAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:385394/67‑1 (MQ=255)
                  cTGGTTGTTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAAcc                                                    <  1:550296/66‑1 (MQ=255)
                    ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTgg                                               <  1:2427871/69‑1 (MQ=255)
                    ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTgg                                               <  1:1961541/69‑1 (MQ=255)
                    ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTgg                                               <  1:530320/69‑1 (MQ=255)
                    ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTgg                                               <  1:1277679/69‑1 (MQ=255)
                    ggTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTgg                                               <  1:2954631/69‑1 (MQ=255)
                       tattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCt                                             <  1:348129/68‑1 (MQ=255)
                       tattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCt                                             <  1:2504939/68‑1 (MQ=255)
                       tattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCt                                             <  1:239321/68‑1 (MQ=255)
                       tattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCt                                             <  1:2392656/68‑1 (MQ=255)
                       tattaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCt                                             <  1:525209/68‑1 (MQ=255)
                        attaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTg                                                               >  1:1088100/1‑49 (MQ=255)
                        attaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCtt                                            >  1:2555442/1‑68 (MQ=255)
                        attaCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCtt                                            >  1:3345573/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:1634084/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:51200/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2644557/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2383410/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2270142/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2888002/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2258641/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:848051/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2944654/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:304393/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:2878743/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:3201480/1‑68 (MQ=255)
                           aCGCAAAGGGCGAGTAAAAGAGTACGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCg                                         >  1:248597/1‑68 (MQ=255)
                                           aaaGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCgca                         >  1:1876595/1‑68 (MQ=255)
                                           aaaGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCgca                         >  1:899859/1‑68 (MQ=255)
                                                           tgatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTcgcg        <  1:2015392/69‑1 (MQ=255)
                                                           tgatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTcgcg        <  1:1971383/69‑1 (MQ=255)
                                                           tgatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTcgcg        <  1:818816/69‑1 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:962042/1‑69 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:3228189/1‑69 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:567558/1‑69 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:1129336/1‑69 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:1257794/1‑69 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:1397165/1‑69 (MQ=255)
                                                            gatgatCAAACTGTTTGGCGAACCGGTGCCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:3010907/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGa       >  1:1551375/1‑68 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:498584/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:501759/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:1442299/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:2146392/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:2149871/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:3061747/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:3020122/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:2195058/1‑69 (MQ=255)
                                                             atgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGAt      >  1:1340875/1‑69 (MQ=255)
                                                              tgatCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATc     >  1:2234580/1‑69 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:140003/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:804872/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:1429982/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:143879/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:1513573/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:2840531/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:3288828/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:2129516/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:3164944/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:2191036/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:2259873/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGagc                                <  1:2667764/38‑1 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:1606044/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:2116275/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:1495392/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:2146964/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:1414269/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:2994585/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:1134888/1‑68 (MQ=255)
                                                                  cAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGt  >  1:2716443/1‑68 (MQ=255)
                                                                  |                                                                   
ACTGCCTGAATCTCCCCGCTGGTTATTACGCAAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGT  >  minE/44060‑44193

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: