Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 505,424 T→C 8.3% G128G (GGT→GGC ybhE → 6‑phosphogluconolactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE505,4240TC 8.3% 159.5 / 6.5 60G128G (GGT→GGCybhE6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base C (2/3);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

TGAGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCGCATTT  >  minE/505357‑505494
                                                                   |                                                                      
tGAGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGc                                                                     <  1:1186661/71‑1 (MQ=255)
 gAGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGc                                                                     >  1:91908/1‑70 (MQ=255)
 gAGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGc                                                                     >  1:31801/1‑70 (MQ=255)
  aGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGTGGCTGGACGGTTGCCa                                                                   >  1:4034089/1‑71 (MQ=255)
  aGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGcc                                                                    <  1:902417/70‑1 (MQ=255)
  aGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCa                                                                   >  1:3219566/1‑71 (MQ=255)
  aGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCa                                                                   <  1:134789/71‑1 (MQ=255)
   gCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTTGACGGCTGCCAt                                                                  <  1:2310320/71‑1 (MQ=255)
   gCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCAt                                                                  <  1:2866200/71‑1 (MQ=255)
    cGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCAt                                                                  <  1:1932866/70‑1 (MQ=255)
      tAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTcc                                                               >  1:4134700/1‑71 (MQ=255)
      tAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTcc                                                               >  1:2707845/1‑71 (MQ=255)
         cgcgTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGcc                                                            <  1:1088997/71‑1 (MQ=255)
         cgcgTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGcc                                                            <  1:1415101/71‑1 (MQ=255)
             tCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGtt                                                                       <  1:1745232/56‑1 (MQ=255)
             tCTGGAAGATGGCCTGCCAGTGGGCGTCGGCGATGTGGGCGAGGGGCTGGACGGTTGc                                                                     <  1:3700325/58‑1 (MQ=255)
               tGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATAtc                                                      >  1:495321/1‑71 (MQ=255)
                 gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATAtct                                                     <  1:3935198/70‑1 (MQ=255)
                 gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATAtct                                                     <  1:1343170/70‑1 (MQ=255)
                 gAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATAtct                                                     >  1:2935586/1‑70 (MQ=255)
                     aTGGCCTGCCAGTGGGCGTCGTCGATGTGTTCGAGGTGCTGGACGGTTGCCATTCCGCCAATATCTCACCg                                                >  1:1371080/1‑71 (MQ=255)
                        gCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGa                                              >  1:3304726/1‑70 (MQ=255)
                        gCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGa                                              >  1:1121240/1‑70 (MQ=255)
                        gCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGa                                              >  1:939776/1‑70 (MQ=255)
                             ccAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCg                                        <  1:1731934/71‑1 (MQ=255)
                                gTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTAc                                     <  1:919116/71‑1 (MQ=255)
                                 tGGGCGGCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTAc                                     <  1:1032755/70‑1 (MQ=255)
                                   ggCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGc                                   <  1:2436825/70‑1 (MQ=255)
                                    gCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtg                                 <  1:3844535/71‑1 (MQ=255)
                                    gCGTCGTCGATGTGGGCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtg                                 <  1:1855220/71‑1 (MQ=255)
                                     cgtcgtCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgt                                <  1:3590873/71‑1 (MQ=255)
                                      gtcgtcGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgtg                               <  1:2463787/71‑1 (MQ=255)
                                      gtcgtcGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgtg                               <  1:3814532/71‑1 (MQ=255)
                                       tcgtcgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATAtctc                                                    <  1:2035392/49‑1 (MQ=255)
                                         gtcgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAata                                                        >  1:3710745/1‑43 (MQ=255)
                                         gtcgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATAtc                                                      >  1:789500/1‑45 (MQ=255)
                                         gtcgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGt                            <  1:1020992/71‑1 (MQ=255)
                                         gtcgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGt                            >  1:801258/1‑71 (MQ=255)
                                         gtcgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGt                            >  1:4056240/1‑71 (MQ=255)
                                           cgATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTggggtt                          >  1:1233959/1‑68 (MQ=255)
                                            gATGTGGTCTAGGGGCTGGACGGTTGCCATTCCGCCaa                                                          >  1:991752/1‑38 (MQ=255)
                                            gATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCaa                                                          >  1:2506401/1‑38 (MQ=255)
                                            gATGTGGTCGAGGGGCTGGACGGCTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCtgt                                >  1:2696492/1‑64 (MQ=255)
                                              tgtgGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCgg                       <  1:1922203/71‑1 (MQ=255)
                                               gtgGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCgg                       >  1:1325666/1‑70 (MQ=255)
                                                tgGTCGAGGGGCTGGACGGTTGCCATTCCGCGAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCa                     <  1:3330334/71‑1 (MQ=255)
                                                tgGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCa                     <  1:3843048/71‑1 (MQ=255)
                                                  gTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTa                                      >  1:2491241/1‑52 (MQ=255)
                                                  gTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTa                                      >  1:2365700/1‑52 (MQ=255)
                                                  gTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCAtt                   <  1:1819138/71‑1 (MQ=255)
                                                  gTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCAtt                   <  1:698910/71‑1 (MQ=255)
                                                  gTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCAt                    <  1:2595541/70‑1 (MQ=255)
                                                  gTCGAGGGGCTGGACGGTTGCCATTACGCCAATATCTCACCGGACAACCGTa                                      >  1:3527613/1‑52 (MQ=255)
                                                        gggCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGc              <  1:1815065/70‑1 (MQ=255)
                                                         ggCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAg            <  1:3637621/71‑1 (MQ=255)
                                                            tGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTa                  <  1:948061/62‑1 (MQ=255)
                                                              gACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGAtgg       >  1:1218677/1‑69 (MQ=255)
                                                              gACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCg       <  1:3081932/71‑1 (MQ=255)
                                                              gACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCg       >  1:3392173/1‑71 (MQ=255)
                                                                cGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCGCa     <  1:1743799/71‑1 (MQ=255)
                                                                cGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCGCa     <  1:2396405/71‑1 (MQ=255)
                                                                   ttGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCGCAttt  <  1:2464217/71‑1 (MQ=255)
                                                                   |                                                                      
TGAGCGTAACGCGTCTGGAAGATGGCCTGCCAGTGGGCGTCGTCGATGTGGTCGAGGGGCTGGACGGTTGCCATTCCGCCAATATCTCACCGGACAACCGTACGCTGTGGGTTCCGGCATTAAAGCAGGATCGCATTT  >  minE/505357‑505494

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: