Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,128,106 C→T 35.0% S177S (TCC→TCT)  yhbV → predicted protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,128,1060CT35.0% 39.4 / 30.6 57S177S (TCC→TCT) yhbVpredicted protease
Reads supporting (aligned to +/- strand):  ref base C (8/29);  new base T (9/11);  total (17/40)
Fisher's exact test for biased strand distribution p-value = 7.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

GAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGAA  >  minE/2128038‑2128157
                                                                    |                                                   
gAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTcc                                                     >  1:2820032/1‑69 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2331572/69‑1 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2646008/69‑1 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:677123/69‑1 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:631438/69‑1 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1397317/69‑1 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2726345/69‑1 (MQ=255)
           tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:3165385/69‑1 (MQ=255)
            cGTAACCAGTTTGAAGTGGAAGTCCTGCGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:1244151/1‑69 (MQ=255)
            cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:2295240/1‑69 (MQ=255)
            cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:1578626/1‑69 (MQ=255)
            cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:870142/1‑69 (MQ=255)
            cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:890378/1‑69 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:3126623/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:3328465/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:493118/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2655582/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2232703/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1487865/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1909585/64‑1 (MQ=255)
                aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1415669/64‑1 (MQ=255)
                         aaGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCg                             <  1:2283154/68‑1 (MQ=255)
                          aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2408129/68‑1 (MQ=255)
                          aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2418542/68‑1 (MQ=255)
                          aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:399640/68‑1 (MQ=255)
                          aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:3375974/68‑1 (MQ=255)
                          aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2696484/68‑1 (MQ=255)
                          aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2892300/68‑1 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTctgc                                                   >  1:2051120/1‑43 (MQ=38)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGtt                               >  1:65518/1‑63 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTc                              >  1:30204/1‑64 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:437275/1‑69 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:3233521/1‑69 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:1936838/1‑69 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:872816/1‑69 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:401645/1‑69 (MQ=255)
                            tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:963389/1‑69 (MQ=255)
                            tGGAAGTCTTGAGCGACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:3199274/1‑69 (MQ=255)
                              gAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            >  1:1922316/1‑64 (MQ=255)
                                    cTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:1400269/1‑61 (MQ=255)
                                    cTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGccc                  >  1:1025287/1‑68 (MQ=255)
                                         caaCGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGaaa              <  1:2445872/65‑1 (MQ=255)
                                               tCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCg                             <  1:966459/46‑1 (MQ=255)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:671184/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:1333154/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:1366254/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:81154/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:1323364/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:3166482/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2635930/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2582118/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2549651/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2464271/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2392731/45‑1 (MQ=38)
                                                   cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:3354134/69‑1 (MQ=255)
                                                   cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:3247407/69‑1 (MQ=255)
                                                   cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:2986948/69‑1 (MQ=255)
                                                   cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:1290604/69‑1 (MQ=255)
                                                   cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:2274914/69‑1 (MQ=255)
                                                                    |                                                   
GAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGAA  >  minE/2128038‑2128157

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