Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,210 G→C 27.2% T348T (ACC→ACG pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,2100GC27.2% 133.9 / 55.2 81T348T (ACC→ACGpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (34/25);  new base C (14/8);  total (48/33)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAT  >  minE/609141‑609280
                                                                     |                                                                      
gTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAggg                                                                       >  1:1284140/1‑71 (MQ=255)
 ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAggg                                                                       >  1:675934/1‑70 (MQ=255)
  tCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTa                                                                     >  1:251769/1‑71 (MQ=255)
  tCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTa                                                                     >  1:359129/1‑71 (MQ=255)
   cGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTAc                                                                    >  1:5628799/1‑71 (MQ=255)
        ccGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCg                                                               >  1:3917869/1‑71 (MQ=255)
         cGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCg                                                               <  1:4329675/70‑1 (MQ=255)
          ggAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTc                                                             <  1:5267283/71‑1 (MQ=255)
          ggAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTc                                                             <  1:4430501/71‑1 (MQ=255)
             gaCGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAg                                                          <  1:1688737/71‑1 (MQ=255)
                ggccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc                                                       >  1:88760/4‑71 (MQ=255)
                ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAgg                                                                        >  1:3280802/1‑54 (MQ=255)
                ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                       >  1:4277363/1‑71 (MQ=255)
                ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                       <  1:4683451/71‑1 (MQ=255)
                ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                       <  1:831753/71‑1 (MQ=255)
                 gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc                                                       >  1:5585024/3‑70 (MQ=255)
                 gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc                                                       >  1:1969608/3‑70 (MQ=255)
                   cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAt                                                    >  1:2876034/1‑71 (MQ=255)
                    ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa                                                           >  1:3538263/1‑63 (MQ=255)
                    ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa                                                           >  1:1727957/1‑63 (MQ=255)
                     cATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac                                                  <  1:1498445/71‑1 (MQ=255)
                     cATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAtcc                                                  <  1:28585/71‑3 (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)
                                     ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                       <  1:1232033/50‑1 (MQ=255)
                                     ttCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc                                                       <  1:2075826/50‑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)
                                            aaCGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCAt                                                    >  1:180679/1‑46 (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)
                                               ggAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCAt                                                    >  1:3627170/1‑43 (MQ=255)
                                               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:5579429/1‑45 (MQ=16)
                                                      gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                         >  1:417144/1‑45 (MQ=16)
                                                      gTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCATCCCGCCGAtca                                         >  1:2007762/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:1202690/36‑1 (MQ=255)
                                                       tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                   <  1:3431981/36‑1 (MQ=255)
                                                       tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                   <  1:3161205/36‑1 (MQ=255)
                                                       tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                   <  1:1735616/36‑1 (MQ=255)
                                                       tttttGGTAACCAGGGTACGACCGTCGAGGCCCATa                                                   <  1:1174576/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:2712682/1‑69 (MQ=255)
                                                         tttGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGcc                >  1:2474570/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)
                                                                     |                                                                      
GTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCAT  >  minE/609141‑609280

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: