Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,403 C→G 13.0% T617T (ACG→ACC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,4030CG13.0% 182.1 / 15.7 69T617T (ACG→ACC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (40/20);  new base G (6/3);  total (46/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGT  >  minE/608336‑608470
                                                                   |                                                                   
gACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGtt                                                                  >  1:2110747/1‑71 (MQ=255)
    atGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCt                                                              <  1:2165478/70‑1 (MQ=255)
     tGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCtt                                                             <  1:2506416/71‑1 (MQ=255)
       cATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCtggc                                                           <  1:536931/71‑4 (MQ=255)
       cATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAc                                                           >  1:186636/1‑71 (MQ=255)
        aTCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                                                            >  1:1333899/1‑69 (MQ=255)
         tCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         <  1:2027753/71‑1 (MQ=255)
         tCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         <  1:1570096/71‑1 (MQ=255)
          cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         >  1:3414294/1‑70 (MQ=255)
          cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         >  1:477167/1‑70 (MQ=255)
              ttAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTacac                                                    >  1:1069453/1‑71 (MQ=255)
              ttAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAccc                                                                     <  1:2251354/54‑1 (MQ=255)
               tAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                    >  1:3024232/1‑70 (MQ=255)
                 gCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAc                                                 <  1:3354176/71‑1 (MQ=255)
                 gCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAc                                                 >  1:3832078/1‑71 (MQ=255)
                    cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGtt                                              <  1:1391101/71‑1 (MQ=255)
                    cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                               >  1:2929608/1‑70 (MQ=255)
                     ccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                                                            <  1:2315160/56‑1 (MQ=255)
                     ccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                               >  1:3540482/1‑69 (MQ=255)
                        gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAc                                                 <  1:435448/64‑1 (MQ=255)
                        gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                            >  1:3297967/1‑69 (MQ=255)
                         tCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACAATACACaa                                                  >  1:2691225/1‑62 (MQ=255)
                           cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACAGGTTTTCTGGCCGTACACCACGTTAGAAGt                                        <  1:1287419/70‑1 (MQ=255)
                           cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGt                                        >  1:1936727/1‑70 (MQ=255)
                            gAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                      >  1:3005560/1‑71 (MQ=255)
                            gAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         <  1:303295/52‑1 (MQ=255)
                            gAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                      <  1:2975766/71‑1 (MQ=255)
                              aCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                                    >  1:2005067/1‑71 (MQ=255)
                              aCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                                    >  1:2077016/1‑71 (MQ=255)
                              aCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                                    >  1:1373502/1‑71 (MQ=255)
                               cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                                    >  1:2907238/1‑70 (MQ=255)
                                ggCGCGCCAGCACGACGACCGTCTGGGGTCTTACCCGTTTTCTTACCATacac                                                    <  1:273582/53‑1 (MQ=255)
                                  cgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                  >  1:2267719/1‑69 (MQ=255)
                                  cgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                  >  1:3077494/1‑69 (MQ=255)
                                   gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGttt                                                                 <  1:2517088/37‑1 (MQ=255)
                                   gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         >  1:1721901/1‑45 (MQ=255)
                                     gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                         >  1:1526192/1‑43 (MQ=255)
                                     gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                             >  1:30672/1‑71 (MQ=255)
                                     gcCAGCACGACGACCGTCTGGGGCGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                             <  1:1787674/71‑1 (MQ=255)
                                      ccAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                               <  1:3632670/68‑1 (MQ=255)
                                        aGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGa                              >  1:1684856/1‑67 (MQ=255)
                                         gCACGTCGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                         >  1:1463806/1‑71 (MQ=255)
                                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCtt                                                             >  1:3749690/1‑35 (MQ=255)
                                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                         >  1:554475/1‑71 (MQ=255)
                                         gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                         >  1:3029703/1‑71 (MQ=255)
                                           gcgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt                       >  1:318/2‑71 (MQ=255)
                                           acgacgACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                             >  1:1709023/1‑65 (MQ=255)
                                           acgacgACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                       >  1:310005/1‑71 (MQ=255)
                                             gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                       >  1:3646235/1‑69 (MQ=255)
                                               cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                   <  1:3425335/71‑1 (MQ=255)
                                                gacCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTggg                    >  1:948367/1‑67 (MQ=255)
                                                 acCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                 >  1:3583383/1‑71 (MQ=255)
                                                  ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                      >  1:3139450/1‑49 (MQ=255)
                                                  ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                      >  1:3952751/1‑49 (MQ=255)
                                                  ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                >  1:3254601/1‑71 (MQ=255)
                                                  ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                >  1:857464/1‑71 (MQ=255)
                                                   cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGaa                                          >  1:693657/1‑44 (MQ=255)
                                                   cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAg                                         >  1:1246570/1‑45 (MQ=255)
                                                     tCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa                                >  1:2094052/1‑52 (MQ=255)
                                                     tCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGTGa              >  1:2063821/1‑70 (MQ=255)
                                                         gggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATa            <  1:2808616/68‑1 (MQ=255)
                                                          ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                               <  1:1055417/48‑1 (MQ=255)
                                                          ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                               <  1:152984/48‑1 (MQ=255)
                                                          ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGt        <  1:2153441/71‑1 (MQ=255)
                                                           gtgtTACCCGTTTTCTTACCATACAAAACGTTAGAAGTGATGGt                                  <  1:2053358/44‑1 (MQ=255)
                                                              ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGaat                            >  1:3065960/1‑46 (MQ=19)
                                                              ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg    >  1:1948356/1‑71 (MQ=255)
                                                              ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg    <  1:2197671/71‑1 (MQ=255)
                                                                aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGacg            >  1:2458262/1‑59 (MQ=255)
                                                                aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg   >  1:2542331/1‑70 (MQ=255)
                                                                 cccGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt  <  1:1640801/70‑1 (MQ=255)
                                                                 cccGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt  <  1:3708927/70‑1 (MQ=255)
                                                                 cccGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt  <  1:2457293/70‑1 (MQ=255)
                                                                   |                                                                   
GACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGT  >  minE/608336‑608470

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: ‑