Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,283 A→G 8.5% F657F (TTT→TTC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,2830AG 8.5% 194.0 / 7.4 70F657F (TTT→TTC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base G (2/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTT  >  minE/608213‑608351
                                                                      |                                                                    
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                      <  1:1712302/71‑1 (MQ=255)
cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                      <  1:2255876/71‑1 (MQ=255)
   cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                                   >  1:1025265/1‑71 (MQ=255)
   cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                                   >  1:3407785/1‑71 (MQ=255)
    gtcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAAc                                                                   <  1:3206845/70‑1 (MQ=255)
       gtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                                 >  1:2550577/1‑69 (MQ=255)
        tcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACAATCTTTAGCGTAAGCAAACGGCAg                                                              <  1:1804639/71‑1 (MQ=255)
         cTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGt                                                             >  1:1149312/1‑71 (MQ=255)
            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGttt                                                           <  1:3813684/70‑1 (MQ=255)
            tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGttt                                                           <  1:821148/70‑1 (MQ=255)
              cccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                          <  1:1176120/69‑1 (MQ=255)
                cAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                          >  1:1048478/1‑67 (MQ=255)
                 aGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                     >  1:7101/1‑71 (MQ=255)
                  gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                    <  1:3339879/71‑1 (MQ=255)
                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                   <  1:269442/71‑1 (MQ=255)
                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                   <  1:1722830/71‑1 (MQ=255)
                   tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTATAGc                                                        >  1:2105193/1‑66 (MQ=255)
                     cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                   >  1:3867587/1‑69 (MQ=255)
                     cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                   >  1:3834838/1‑69 (MQ=255)
                      gTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAg                                                >  1:369184/1‑71 (MQ=255)
                        tCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                               >  1:361587/1‑70 (MQ=255)
                         cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACg                                                                  >  1:101027/1‑50 (MQ=255)
                         cGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                               >  1:3140270/1‑69 (MQ=255)
                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                      >  1:1670850/1‑61 (MQ=255)
                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                            <  1:3333443/71‑1 (MQ=255)
                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                            <  1:2044591/71‑1 (MQ=255)
                          ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                            <  1:694873/71‑1 (MQ=255)
                           gAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                            <  1:3827985/70‑1 (MQ=255)
                            aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                             >  1:1982606/1‑68 (MQ=255)
                            aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                             >  1:3352267/1‑68 (MQ=255)
                            aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                             >  1:1412528/1‑68 (MQ=255)
                             aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                           >  1:2992609/1‑67 (MQ=255)
                             aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                         <  1:324297/71‑1 (MQ=255)
                             aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                           >  1:3302705/1‑69 (MQ=255)
                             aCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                           >  1:3283272/1‑69 (MQ=255)
                                attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                      <  1:2742794/68‑1 (MQ=255)
                                attgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                      <  1:267180/68‑1 (MQ=255)
                                aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                      <  1:2445083/71‑1 (MQ=255)
                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                   >  1:1607299/1‑57 (MQ=255)
                                 tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                      >  1:428323/1‑70 (MQ=255)
                                  agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagag                                          >  1:2221826/1‑65 (MQ=255)
                                    agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                      <  1:2212037/67‑1 (MQ=255)
                                    agaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACt                                  >  1:2131970/1‑71 (MQ=255)
                                      aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                 >  1:523529/1‑54 (MQ=255)
                                      aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                 >  1:1123783/1‑70 (MQ=255)
                                          tgtACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGac                                                 <  1:2837786/50‑2 (MQ=255)
                                           gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                           >  1:953217/1‑71 (MQ=255)
                                            tacgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCttt                          <  1:1753908/68‑1 (MQ=255)
                                              ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                         <  1:2401660/70‑1 (MQ=255)
                                                  aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                    <  1:1654394/71‑1 (MQ=255)
                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                               <  1:2029132/42‑1 (MQ=255)
                                                    aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTc                    >  1:1331851/1‑69 (MQ=255)
                                                      cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGa                >  1:720036/1‑71 (MQ=255)
                                                       aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                      <  1:3942664/48‑1 (MQ=255)
                                                       aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACATTTCTGGTCACGAc               >  1:3260619/1‑71 (MQ=255)
                                                        tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAc               >  1:1051499/1‑70 (MQ=255)
                                                          tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt            >  1:1980820/1‑71 (MQ=255)
                                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGTACCTTTCTGGTCACGACCGTg           >  1:87529/1‑70 (MQ=255)
                                                            tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc              <  1:2037412/67‑1 (MQ=255)
                                                               cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc       <  1:3250035/71‑1 (MQ=255)
                                                               cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTAACGACCGTGc          <  1:913851/68‑1 (MQ=255)
                                                               cGTAAGAAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc       <  1:3164314/71‑1 (MQ=255)
                                                                 tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTc                         <  1:3559758/51‑1 (MQ=255)
                                                                 tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg     >  1:1019022/1‑71 (MQ=255)
                                                                 tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTAATg                                 <  1:355408/43‑1 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                                >  1:2251034/1‑42 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCATCTTTCTGGTCACGACCGTGCATCGGGt   <  1:3696793/71‑1 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt   >  1:850934/1‑71 (MQ=255)
                                                                   aGCAAACGGCAGTTTAGCAACGGAAGGCAGCGCGGCTACTGCACCTTTCTGGTCACGACCGTGCATCggg    <  1:3344432/70‑1 (MQ=255)
                                                                    gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                     <  1:1104843/52‑1 (MQ=255)
                                                                    gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCgg     >  1:1086037/1‑68 (MQ=255)
                                                                    gCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGtt  <  1:3646045/71‑1 (MQ=255)
                                                                      |                                                                    
CTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTT  >  minE/608213‑608351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑