Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,397 G→A 14.3% N619N (AAC→AAT)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,3970GA14.3% 174.6 / 17.6 70N619N (AAC→AAT) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (39/20);  new base A (6/4);  total (45/25)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGT  >  minE/608331‑608464
                                                                  |                                                                   
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)
              tCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                   <  1:2027753/71‑1 (MQ=255)
              tCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                   <  1:1570096/71‑1 (MQ=255)
               cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                   >  1:477167/1‑70 (MQ=255)
               cGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCa                                                   >  1:3414294/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: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:3029703/1‑71 (MQ=255)
                                              gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                   >  1:554475/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:152984/48‑1 (MQ=255)
                                                               ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                         <  1:1055417/48‑1 (MQ=255)
                                                               ggTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGt  <  1:2153441/71‑1 (MQ=255)
                                                                gtgtTACCCGTTTTCTTACCATACAAAACGTTAGAAGTGATGGt                            <  1:2053358/44‑1 (MQ=255)
                                                                  |                                                                   
GTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGT  >  minE/608331‑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: ‑