Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,097 C→G 13.4% A719A (GCG→GCC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,0970CG13.4% 182.4 / 16.5 67A719A (GCG→GCC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (33/25);  new base G (2/7);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 7.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCG  >  minE/608034‑608166
                                                               |                                                                     
gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTc                                                                  >  1:2861031/1‑69 (MQ=255)
gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                                 <  1:2574956/70‑1 (MQ=255)
 aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCCTCGag                                                               <  1:589931/71‑1 (MQ=255)
  aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                                 >  1:1936877/1‑68 (MQ=255)
   aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                               >  1:566372/1‑69 (MQ=255)
   aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                               >  1:529629/1‑69 (MQ=255)
     gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                                                             >  1:3640/1‑69 (MQ=255)
       aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAt                                                         >  1:3378394/1‑71 (MQ=255)
        cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                        <  1:1265563/69‑1 (MQ=12)
        cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                        <  1:3645971/69‑1 (MQ=12)
          gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                       >  1:2362056/1‑70 (MQ=255)
           cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                                     >  1:2963865/1‑71 (MQ=255)
               gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                 <  1:3744660/71‑1 (MQ=255)
                 cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                                                  <  1:2729242/68‑2 (MQ=9)
                  agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg                                                   <  1:263141/66‑1 (MQ=9)
                  aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                              <  1:385533/71‑1 (MQ=255)
                  aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                              >  1:396873/1‑71 (MQ=255)
                   gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                             <  1:803100/68‑1 (MQ=9)
                   gATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                       <  1:303918/61‑1 (MQ=255)
                    aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                             >  1:3380607/1‑70 (MQ=255)
                      acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                          >  1:1149639/1‑71 (MQ=255)
                       cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                          <  1:251882/70‑1 (MQ=255)
                        acGGATGGTCAGCTGCGGAAATTTTTCCGGGTTTTCCATCGCGTCGagcagc                                                           >  1:25087/1‑52 (MQ=255)
                          ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                      <  1:1812077/71‑2 (MQ=255)
                          ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                       >  1:2612834/1‑70 (MQ=255)
                               gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                       <  1:731714/49‑1 (MQ=255)
                                tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                    <  1:160098/51‑1 (MQ=255)
                                tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                <  1:928150/71‑1 (MQ=255)
                                 cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                 >  1:254837/1‑69 (MQ=255)
                                  aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                              >  1:1010417/1‑71 (MQ=255)
                                  aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                              <  1:3241491/71‑1 (MQ=255)
                                  aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                              >  1:3155460/1‑71 (MQ=255)
                                   gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                            >  1:2916952/1‑56 (MQ=255)
                                   gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                               <  1:1912623/69‑1 (MQ=255)
                                    cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                            <  1:3583132/71‑1 (MQ=255)
                                    cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                            <  1:2538351/71‑1 (MQ=255)
                                    cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                            >  1:40621/1‑71 (MQ=255)
                                      taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                          >  1:1934025/3‑70 (MQ=255)
                                      gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                          >  1:3690568/1‑71 (MQ=255)
                                        ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                    >  1:2967272/1‑43 (MQ=255)
                                         gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                       <  1:1929460/71‑1 (MQ=255)
                                         gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                       >  1:616967/1‑71 (MQ=12)
                                            atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                        <  1:1328372/51‑1 (MQ=255)
                                            atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                    >  1:564393/1‑71 (MQ=255)
                                             tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                   <  1:2149654/71‑1 (MQ=255)
                                              ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                  >  1:855116/1‑71 (MQ=255)
                                              ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                  >  1:3251582/1‑71 (MQ=255)
                                              ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                  >  1:534188/1‑71 (MQ=255)
                                               tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                  <  1:3278536/70‑1 (MQ=255)
                                                   cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga             <  1:868636/71‑1 (MQ=255)
                                                   cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga             <  1:1278827/71‑1 (MQ=255)
                                                    gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                             <  1:558239/38‑1 (MQ=255)
                                                    gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat            <  1:1627687/71‑1 (MQ=255)
                                                      gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg           >  1:1553947/1‑70 (MQ=255)
                                                       ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                              <  1:1953432/50‑1 (MQ=255)
                                                       ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                              <  1:237535/50‑1 (MQ=255)
                                                       ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga         >  1:1659397/1‑71 (MQ=255)
                                                         ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat        >  1:1977977/1‑70 (MQ=255)
                                                          tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc      >  1:3689004/1‑71 (MQ=255)
                                                           ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc      <  1:3419300/70‑1 (MQ=255)
                                                           ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc      <  1:1213134/70‑1 (MQ=255)
                                                            cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                             <  1:904983/46‑1 (MQ=255)
                                                            cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt    >  1:3171065/1‑71 (MQ=255)
                                                             aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                >  1:1568133/1‑58 (MQ=255)
                                                             aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCtt                >  1:1577473/1‑58 (MQ=255)
                                                             aTCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTc   >  1:350886/1‑71 (MQ=255)
                                                              tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCg  >  1:136118/1‑71 (MQ=255)
                                                               cgcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                              >  1:1356558/1‑42 (MQ=255)
                                                               |                                                                     
GAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCG  >  minE/608034‑608166

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