Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,879,390 C→T 18.7% L212L (CTG→CTA ytfE ← predicted regulator of cell morphogenesis and cell wall metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,879,3900CT18.7% 175.0 / 25.0 111L212L (CTG→CTAytfEpredicted regulator of cell morphogenesis and cell wall metabolism
Reads supporting (aligned to +/- strand):  ref base C (32/58);  new base T (12/9);  total (44/67)
Fisher's exact test for biased strand distribution p-value = 8.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 111 total aligned reads displayed.

ATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATT  >  minE/2879326‑2879443
                                                                |                                                     
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:841510/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:36363/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:475715/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:3304868/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:3265808/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:3200889/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:1551885/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:915685/1‑69 (MQ=255)
aTCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTa                                                   >  1:2911271/1‑69 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCTCGCGCGGGAATAGTAc                                                  <  1:3131044/64‑1 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGt                                             >  1:2253824/1‑69 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGt                                             >  1:2236991/1‑69 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTAc                                                  <  1:1281040/64‑1 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTAc                                                  <  1:2366529/64‑1 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTAc                                                  <  1:2639575/64‑1 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTAc                                                  <  1:1532088/64‑1 (MQ=255)
      cTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTAc                                                  <  1:2103278/64‑1 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:978444/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:2820051/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:552790/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:1205638/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:860433/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:2337038/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:867696/1‑69 (MQ=255)
        ggATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGttt                                           >  1:2486267/1‑69 (MQ=255)
            ggCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGTGAACAGTACATTGTTTTCCa                                       <  1:1751567/69‑1 (MQ=255)
            ggCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                       <  1:2781767/69‑1 (MQ=255)
            ggCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                       <  1:1979121/69‑1 (MQ=255)
            ggCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                       <  1:3212421/69‑1 (MQ=255)
            ggCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                       <  1:159193/69‑1 (MQ=255)
            ggCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCa                                       <  1:2439673/69‑1 (MQ=255)
               gcAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGa                                     >  1:1993383/1‑68 (MQ=255)
               gcAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGa                                     >  1:669422/1‑68 (MQ=255)
               gcAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGa                                     >  1:1220018/1‑68 (MQ=255)
                cAATGCCATCCAGCTTTTAGATAACTCACCCGCCAGCGCGCGCGGGAATAGTAc                                                  <  1:876220/54‑1 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCCGCGCGCGCGGGAATAGTACATTGt                                             >  1:167874/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCCGCGCGCGCGGGAATAGTACATTGt                                             >  1:94400/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:2712921/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:1163450/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:1311631/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:3033353/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:3076915/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:1548665/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:2162992/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:737114/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGCGGGAATAGTACATTGt                                             >  1:656970/1‑57 (MQ=255)
                  aTGCCATCCAGCTTTTAGATCACTCACCCGCAAGCGCGCGCGGGAATAGTACATTGt                                             >  1:1444208/1‑57 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:762438/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:320997/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:3290534/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:2731143/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:1059265/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:1056003/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:1636119/1‑69 (MQ=255)
                     ccATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgt                              >  1:2212209/1‑69 (MQ=255)
                       aTCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgtg                             <  1:1089570/68‑1 (MQ=255)
                       aTCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgtg                             <  1:1672241/68‑1 (MQ=255)
                       aTCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgtg                             <  1:1669579/68‑1 (MQ=255)
                       aTCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGAtgtg                             <  1:1527074/68‑1 (MQ=255)
                          cAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACt                                   >  1:2708597/1‑59 (MQ=255)
                          cAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATc                          <  1:1525886/68‑1 (MQ=255)
                                    aTCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtca                <  1:880771/68‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:498819/69‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:641653/69‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:701253/69‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:1992695/69‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:2130975/69‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:2639699/69‑1 (MQ=255)
                                     tcactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtcatc              <  1:2795345/69‑1 (MQ=255)
                                      cactcacCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGtc                 <  1:3238180/65‑1 (MQ=255)
                                                 cAGCGCGCGTTGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1383145/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTGAtt  <  1:1899181/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:785079/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:841170/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1485373/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1637289/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1358230/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:583369/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1727464/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2773708/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:36056/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:265121/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1917851/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:197724/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:3213891/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2597424/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2352352/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:235529/69‑1 (MQ=255)
                                                 cAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:294465/69‑1 (MQ=255)
                                                  aGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:628070/68‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2045658/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1146516/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2800397/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2969020/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:2988694/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:1685123/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:3380855/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:3426286/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:610170/56‑1 (MQ=255)
                                                              aaCAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCAtt  <  1:378249/56‑1 (MQ=255)
                                                                |                                                     
ATCATGCTGGATGGCGCAATGCCATCCAGCTTTTAGATCACTCACCCGCCAGCGCGCGTGGGAACAGTACATTGTTTTCCAGACTGATGTGATCCATCAGGTCATCAATCAGTTCATT  >  minE/2879326‑2879443

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