Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,293,682 C→T 37.4% H271H (CAC→CAT gldA → glycerol dehydrogenase, NAD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,293,6820CT37.4% 58.0 / 59.7 96H271H (CAC→CATgldAglycerol dehydrogenase, NAD
Reads supporting (aligned to +/- strand):  ref base C (39/21);  new base T (9/27);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 2.79e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TGGCTGCGGCGCACGCAGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTGG  >  minE/2293618‑2293740
                                                                |                                                          
tGGCTGCGGCGCACGCAGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACgg                                                          >  1:91612/1‑67 (MQ=255)
tGGCTGCGGCGCACGCAGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGt                                                         >  1:2805084/1‑68 (MQ=255)
           cacgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:1556233/69‑1 (MQ=255)
           cacgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:2053389/69‑1 (MQ=255)
           cacgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:2996549/69‑1 (MQ=255)
           cacgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:1561111/69‑1 (MQ=255)
           cacgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:1318076/69‑1 (MQ=255)
           cacgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:1122982/69‑1 (MQ=255)
            acgcaGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:2209509/68‑1 (MQ=255)
                aGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCgg                                        <  1:1381352/69‑1 (MQ=255)
                aGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCgg                                        <  1:610654/69‑1 (MQ=255)
                aGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCgg                                        <  1:395694/69‑1 (MQ=255)
                aGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCgg                                        <  1:235900/69‑1 (MQ=255)
                aGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCgg                                        <  1:329510/69‑1 (MQ=255)
                    cATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACg                                    >  1:1090398/1‑69 (MQ=255)
                    cATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACg                                    >  1:2783554/1‑69 (MQ=255)
                    cATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACg                                    >  1:1826994/1‑69 (MQ=255)
                    cATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACg                                    >  1:585030/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:2869352/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:704064/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:1296638/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:1490282/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:1875720/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:199064/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:2238986/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCt                                  >  1:622530/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATACGGTACGCt                                  >  1:973526/1‑69 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCt                                  <  1:3141775/69‑1 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCt                                  <  1:2071039/69‑1 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCt                                  <  1:2774866/69‑1 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCt                                  <  1:768121/69‑1 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCt                                  <  1:464163/69‑1 (MQ=255)
                      tAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCt                                  <  1:139776/69‑1 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAg                                                  <  1:1453860/51‑1 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAg                                                  <  1:2710267/51‑1 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTg                                 >  1:75110/1‑68 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:2976726/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:320091/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:3430618/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:729848/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:2340951/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:2326992/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:755466/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:1226146/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:939171/1‑69 (MQ=255)
                        aCGGCCTGACCGCTATCCCGGACGCGCATCACTATTACCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:2435295/1‑69 (MQ=255)
                          ggCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:187125/1‑67 (MQ=255)
                          ggCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:2355301/1‑67 (MQ=255)
                          ggCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:895298/1‑67 (MQ=255)
                          ggCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:3309599/1‑67 (MQ=255)
                          ggCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGa                                >  1:1017767/1‑67 (MQ=255)
                            ccTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCa                                             <  1:2012597/52‑1 (MQ=255)
                             cTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAg                           >  1:385594/1‑69 (MQ=255)
                             cTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAg                           >  1:2045247/1‑69 (MQ=255)
                             cTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAg                           >  1:2733429/1‑69 (MQ=255)
                             cTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAg                           >  1:1999765/1‑69 (MQ=255)
                             cTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAg                           >  1:2240493/1‑69 (MQ=255)
                             cTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAg                           >  1:1865279/1‑69 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGTTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:621326/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:1410643/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:400212/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:452437/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:1810043/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:47475/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:584187/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:1447358/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:2708133/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:2376826/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:232116/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:886213/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:213551/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:2030819/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:2765881/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:2961456/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:2864253/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:1737243/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3009385/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:301564/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3032030/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3054433/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3197494/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3261887/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:1609282/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3365359/58‑1 (MQ=255)
                                         ccGGACGCGCATCACTATTATCATGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGc                          <  1:3379600/58‑1 (MQ=255)
                                              cgcgCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATgc          >  1:1845985/1‑69 (MQ=255)
                                                 gcATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGcc       >  1:2575490/1‑69 (MQ=255)
                                                 gcATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGcc       >  1:1292590/1‑69 (MQ=255)
                                                 gcATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGcc       >  1:2490517/1‑69 (MQ=255)
                                                 gcATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGcc       >  1:568119/1‑69 (MQ=255)
                                                 gcATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGcc       >  1:2730744/1‑69 (MQ=255)
                                                      aCTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTgg  >  1:2824503/1‑69 (MQ=255)
                                                      aCTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTgg  >  1:622900/1‑69 (MQ=255)
                                                      aCTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTgg  >  1:802232/1‑69 (MQ=255)
                                                      aCTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTgg  >  1:2005151/1‑69 (MQ=255)
                                                      aCTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTgg  >  1:928383/1‑69 (MQ=255)
                                                                |                                                          
TGGCTGCGGCGCACGCAGTGCATAACGGCCTGACCGCTATCCCGGACGCGCATCACTATTATCACGGTGAAAAAGTGGCATTCGGTACGCTGACGCAGCTGGTTCTGGAAAATGCGCCGGTGG  >  minE/2293618‑2293740

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: