Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,394 G→C 15.9% T620T (ACC→ACG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,3940GC15.9% 156.2 / 21.8 69T620T (ACC→ACGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (41/17);  new base C (6/5);  total (47/22)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGT  >  minE/608324‑608464
                                                                      |                                                                      
cTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                                                        >  1:161795/1‑71 (MQ=255)
cTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                                                        >  1:669612/1‑71 (MQ=255)
  ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                                      >  1:1575891/1‑71 (MQ=255)
   tgcggTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgt                                                                      <  1:3566356/67‑1 (MQ=255)
   tCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGTCGCGCCAGCACGACGACCGTCTgggg                                                                       <  1:3913036/69‑1 (MQ=255)
   tCTGGTCACGACCGCGCATCGGCTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTgggg                                                                       <  1:3622849/69‑1 (MQ=255)
       gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTa                                                                  <  1:3237536/70‑1 (MQ=255)
       gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTa                                                                  <  1:319054/70‑1 (MQ=255)
       gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                                                 >  1:1344474/1‑71 (MQ=255)
       gTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAc                                                                 <  1:802030/71‑1 (MQ=255)
        tCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTa                                                                  >  1:3712209/1‑69 (MQ=255)
          aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccgg                                                              >  1:2130713/1‑69 (MQ=255)
            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)
                     tCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                                      >  1:1597242/1‑52 (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)
                           tAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                                                        >  1:2723754/1‑44 (MQ=255)
                           tAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                              >  1:3024232/1‑70 (MQ=255)
                             gCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAc                                           >  1:3832078/1‑71 (MQ=255)
                             gCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAAc                                           <  1:3354176/71‑1 (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:1373502/1‑71 (MQ=255)
                                          aCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                              >  1:2005067/1‑71 (MQ=255)
                                          aCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                              >  1:2077016/1‑71 (MQ=255)
                                           cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATg                              >  1:2907238/1‑70 (MQ=255)
                                            ggCGCGCCAGCACGACGACCGTCTGGGGTCTTACCCGTTTTCTTACCATacac                                              <  1:273582/53‑1 (MQ=255)
                                              cgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                            >  1:3077494/1‑69 (MQ=255)
                                              cgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                            >  1:2267719/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:3952751/1‑49 (MQ=255)
                                                              ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                >  1:3139450/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:152984/48‑1 (MQ=255)
                                                                      ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                         <  1:1055417/48‑1 (MQ=255)
                                                                      ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGt  <  1:2153441/71‑1 (MQ=255)
                                                                      |                                                                      
CTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGT  >  minE/608324‑608464

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: