Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,231 A→C 19.4% G341G (GGT→GGG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,2310AC19.4% 134.5 / 35.9 72G341G (GGT→GGGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (8/6);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGG  >  minE/609162‑609300
                                                                     |                                                                     
cATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                                      <  1:1498445/71‑1 (MQ=255)
 aTGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                                      >  1:3725049/1‑70 (MQ=255)
  tGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATacca                                                                    >  1:1355781/1‑71 (MQ=255)
    gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgcc                                                                  >  1:3131886/1‑71 (MQ=255)
    gtgtGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATcccgc                                                                   >  1:4181715/1‑68 (MQ=255)
    gtgtACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATaccac                                                                   >  1:1664215/1‑70 (MQ=255)
       tACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGAt                                                               <  1:2213120/71‑1 (MQ=255)
         cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                                >  1:1558318/1‑68 (MQ=255)
         cAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATa                                                              <  1:4373535/70‑1 (MQ=255)
           gggTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATAACACCGATAGAt                                                           <  1:4219930/71‑1 (MQ=255)
                ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGt                                                      >  1:146097/1‑71 (MQ=255)
                  cAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                                      >  1:1681998/1‑53 (MQ=255)
                   aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                    <  1:1815458/70‑1 (MQ=255)
                   aGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGc                                                   <  1:3223435/71‑1 (MQ=255)
                       aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAg                                               <  1:5120824/71‑1 (MQ=255)
                        aCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGa                                              <  1:1618778/71‑1 (MQ=255)
                         cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCg                                                                 >  1:254371/1‑51 (MQ=255)
                         cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                                >  1:1655566/1‑52 (MQ=255)
                         cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGa                                                                >  1:3229867/1‑52 (MQ=255)
                         cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                  >  1:2858751/1‑66 (MQ=255)
                         cGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGAt                                             <  1:1468218/71‑1 (MQ=255)
                          ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAt                                                               >  1:2928662/1‑52 (MQ=16)
                          ggAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAt                                                               >  1:3848782/1‑52 (MQ=16)
                          ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                                       >  1:2470007/1‑44 (MQ=255)
                          ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGcc                                                  >  1:4980736/1‑65 (MQ=255)
                             aGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                          <  1:2476638/70‑1 (MQ=255)
                              ggagTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGt                                        <  1:2964678/68‑1 (MQ=255)
                                 gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                             >  1:1750265/1‑45 (MQ=16)
                                 gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                             >  1:2007762/1‑45 (MQ=16)
                                 gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                             >  1:417144/1‑45 (MQ=16)
                                 gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                                             >  1:5579429/1‑45 (MQ=16)
                                 gTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                                       >  1:4654558/1‑53 (MQ=255)
                                  tttttGGTCACCAGCGTACGACCGTCCAGCCCCATcccgc                                                                   <  1:2577542/40‑3 (MQ=12)
                                  tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtc                                                              <  1:3875433/45‑2 (MQ=12)
                                  tttttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTggg                                         >  1:4588029/1‑66 (MQ=255)
                                  tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                                       <  1:3431981/36‑1 (MQ=255)
                                  tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                                       <  1:3161205/36‑1 (MQ=255)
                                  tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                                       <  1:1202690/36‑1 (MQ=255)
                                  tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                                       <  1:1174576/36‑1 (MQ=255)
                                  tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                                       <  1:1735616/36‑1 (MQ=255)
                                  tttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                     <  1:151642/70‑1 (MQ=255)
                                   ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                     >  1:1515688/1‑69 (MQ=255)
                                   ttttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGc                                     >  1:2138134/1‑69 (MQ=255)
                                    tttGGTAACCAGGGTACGACCGTCGAGGCCCATacca                                                                    >  1:1993855/1‑37 (MQ=255)
                                    tttGGTAACCAGGGTACGACCGTCGAGGCCCATacc                                                                     >  1:266622/1‑36 (MQ=255)
                                    tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCg                                           >  1:2000489/1‑62 (MQ=255)
                                    tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                    >  1:2474570/1‑69 (MQ=255)
                                    tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                                    >  1:2712682/1‑69 (MQ=255)
                                     ttGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                  <  1:1269378/70‑3 (MQ=255)
                                       ggTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCagag                                >  1:1969865/1‑70 (MQ=255)
                                        gtcaCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGaa                              <  1:4104932/68‑1 (MQ=255)
                                        gTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAgaga                               <  1:2512797/70‑1 (MQ=255)
                                           aCCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                           >  1:845085/1‑71 (MQ=255)
                                           aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGcc                                    <  1:765651/62‑1 (MQ=255)
                                           aCCAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGCCGGAGAACAg                           <  1:3508966/71‑1 (MQ=255)
                                             cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGtt                         >  1:4221212/1‑71 (MQ=255)
                                             cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                          >  1:2647042/1‑70 (MQ=255)
                                             cAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGt                          >  1:3660695/1‑70 (MQ=255)
                                             cAGCGTACGACCGTCCAGCCCCATCCCGCCGATCACTTCCGTCGCCCAGATTGGGTCGccgg                                  >  1:2796722/1‑60 (MQ=255)
                                                ggTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAt                      <  1:1293650/71‑1 (MQ=255)
                                                  tACGACCGTCGAGGCCCATACCACCGATAGATTCgg                                                       >  1:4678409/1‑36 (MQ=255)
                                                  tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                    <  1:3889893/71‑1 (MQ=255)
                                                  tACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCg                    <  1:2079290/71‑1 (MQ=255)
                                                    cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                    <  1:3613668/37‑1 (MQ=255)
                                                    cGACCGTCGAGGCCCATACCACCGATAGATTCGGTTg                                                    <  1:3607041/37‑1 (MQ=255)
                                                      aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAtt                >  1:1698968/1‑71 (MQ=255)
                                                      aCCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTAt                 >  1:624818/1‑70 (MQ=255)
                                                       ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc               >  1:572776/1‑71 (MQ=255)
                                                       ccGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTc               >  1:4124624/1‑71 (MQ=255)
                                                         gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCgg                                          <  1:2518763/42‑1 (MQ=255)
                                                         gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg             <  1:850734/71‑1 (MQ=255)
                                                         gTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg             <  1:1691527/71‑1 (MQ=255)
                                                          tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                            <  1:4247526/55‑1 (MQ=255)
                                                          tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                            <  1:4986670/55‑1 (MQ=255)
                                                          tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                            <  1:3899602/55‑1 (MQ=255)
                                                          tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACa                            <  1:3307660/55‑1 (MQ=255)
                                                          tCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCg             >  1:4205076/1‑70 (MQ=255)
                                                             aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg          >  1:4514107/1‑70 (MQ=255)
                                                                 ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg      >  1:3484677/1‑70 (MQ=255)
                                                                  cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                             <  1:4291754/46‑1 (MQ=255)
                                                                  cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                             <  1:3358088/46‑1 (MQ=255)
                                                                  cATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                             <  1:2868013/46‑1 (MQ=255)
                                                                   aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                             <  1:4092524/45‑1 (MQ=255)
                                                                   aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg   >  1:4629968/1‑71 (MQ=255)
                                                                    tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAAc                             >  1:2804630/1‑44 (MQ=255)
                                                                    tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAg                           >  1:797551/1‑46 (MQ=255)
                                                                    tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg  <  1:1712763/71‑1 (MQ=255)
                                                                    tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg  <  1:3820930/71‑1 (MQ=255)
                                                                     |                                                                     
CATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGG  >  minE/609162‑609300

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: