Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,968,171 C→G 8.5% M153I (ATG→ATC) ‡ tktA ← transketolase 1, thiamin‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,968,1710CG 8.5% 245.4 / 10.6 82M153I (ATG→ATC) ‡tktAtransketolase 1, thiamin‑binding
Reads supporting (aligned to +/- strand):  ref base C (33/42);  new base G (4/3);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

CCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTT  >  minE/1968104‑1968241
                                                                   |                                                                      
ccAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     <  1:3230136/71‑1 (MQ=255)
 cAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     >  1:247989/1‑70 (MQ=255)
 cAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     >  1:845661/1‑70 (MQ=255)
  aGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     >  1:2606726/1‑69 (MQ=255)
  aGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     >  1:2610651/1‑69 (MQ=255)
  aGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     >  1:2691993/1‑69 (MQ=255)
  aGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGaa                                                                   <  1:427668/71‑1 (MQ=255)
   gCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAg                                                                  >  1:2594630/1‑71 (MQ=255)
   gCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAg                                                                  >  1:3072947/1‑71 (MQ=255)
       cAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGc                                                                <  1:3096455/69‑1 (MQ=255)
       cAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGc                                                                <  1:1178532/69‑1 (MQ=255)
        atcGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTa                                                             <  1:2030367/69‑1 (MQ=255)
        aGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGc                                                                <  1:2654315/68‑1 (MQ=255)
        aGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTa                                                             <  1:1958560/71‑1 (MQ=255)
            tACCCGCCAGAGAGCAAACTTCGTGGGAGATTCCTTCCATCATGCAGCCGTCGCCCATg                                                                     <  1:2020731/59‑1 (MQ=255)
            tACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     <  1:3140370/59‑1 (MQ=255)
            tACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     <  1:3254777/59‑1 (MQ=255)
            tACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGaa                                                                   <  1:3078283/61‑1 (MQ=255)
            tACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGaa                                                                   <  1:389248/61‑1 (MQ=255)
              cccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                        <  1:1348888/70‑1 (MQ=255)
              cccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                        <  1:3218328/70‑1 (MQ=255)
              cccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                        <  1:2516158/70‑1 (MQ=255)
              cccGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                        <  1:351745/70‑1 (MQ=255)
                    agagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                                 <  1:1766970/71‑1 (MQ=255)
                     gagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCg                                                <  1:1318057/71‑1 (MQ=255)
                     gagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCg                                                <  1:591497/71‑1 (MQ=255)
                     gagagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCg                                                <  1:1864581/71‑1 (MQ=255)
                      agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAgg                                                           >  1:2664212/1‑59 (MQ=255)
                      agagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAgg                                                           >  1:987709/1‑59 (MQ=255)
                       gagCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGAc                                              <  1:2442227/71‑1 (MQ=255)
                        agCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAg                                                            >  1:47163/1‑56 (MQ=255)
                        agCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACa                                             <  1:2500579/71‑1 (MQ=255)
                        agCAAACTTCGTGGGAGATGCCTTCCATCATCCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACa                                             <  1:588497/71‑1 (MQ=255)
                         gCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtg                                                         <  1:584748/58‑1 (MQ=255)
                            aaCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGaa                                                                   <  1:2456592/45‑1 (MQ=255)
                             aCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTc                                        <  1:105465/71‑1 (MQ=255)
                              cTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                     <  1:3076438/41‑1 (MQ=255)
                                tCGTGGGAGATGCCTTCCATCATGCAGCCGTCGccc                                                                        >  1:942992/1‑36 (MQ=255)
                                tCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTg                                     >  1:3228926/1‑71 (MQ=255)
                                tCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTg                                     >  1:939055/1‑71 (MQ=255)
                                       agaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtg                                                         <  1:2713901/44‑1 (MQ=255)
                                       agaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCgg                               >  1:2201190/1‑70 (MQ=255)
                                         aTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtg                                                         <  1:862404/42‑1 (MQ=255)
                                         aTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtg                                                         <  1:2363072/42‑1 (MQ=255)
                                           gCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACa                                             >  1:266404/1‑52 (MQ=255)
                                           gCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGt                          >  1:3164596/1‑71 (MQ=255)
                                            ccTTCCATCATGCAGCCGTCGCCCATGAAGGCGTa                                                             <  1:2864262/35‑1 (MQ=255)
                                            ccTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                        <  1:2898709/40‑1 (MQ=255)
                                            ccTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGAc                                              <  1:365292/50‑1 (MQ=255)
                                             cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTa                        <  1:2532176/71‑1 (MQ=255)
                                             cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTa                        <  1:2998872/71‑1 (MQ=255)
                                             cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTa                        <  1:1314880/71‑1 (MQ=255)
                                             cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTa                        <  1:651863/71‑1 (MQ=255)
                                              ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAgg                                                           <  1:3110940/35‑1 (MQ=255)
                                              ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                 >  1:192236/1‑61 (MQ=255)
                                              ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                 >  1:1553926/1‑61 (MQ=255)
                                              ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                 >  1:1354808/1‑61 (MQ=255)
                                              ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                 >  1:1870855/1‑61 (MQ=255)
                                                ccATCATGCAGCCGTCGCCCATGAAGGCGTAGgtg                                                         >  1:1096484/1‑35 (MQ=255)
                                                ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTgg                                    >  1:3126093/1‑56 (MQ=255)
                                                ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTgg                                    >  1:1925630/1‑56 (MQ=255)
                                                ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTgg                                    >  1:2222264/1‑56 (MQ=255)
                                                ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTaa                       <  1:2390893/69‑1 (MQ=255)
                                                ccATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAAc                     >  1:837664/1‑71 (MQ=255)
                                                   tcatGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTc                                                 <  1:3151942/40‑1 (MQ=255)
                                                   tcatGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTc                                                 <  1:2912398/40‑1 (MQ=255)
                                                   tcatGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTc                                                 <  1:3155845/40‑1 (MQ=255)
                                                   tcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAg                                                      <  1:1659893/35‑1 (MQ=255)
                                                   tcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAg                                                      <  1:3015124/35‑1 (MQ=255)
                                                    catGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCg                                       >  1:1623396/1‑49 (MQ=255)
                                                     atGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTCGa                                               >  1:151091/1‑40 (MQ=38)
                                                     atGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTCGa                                               >  1:3036865/1‑40 (MQ=38)
                                                     atGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTCGa                                               >  1:2901893/1‑40 (MQ=38)
                                                     atGCAGCCGTCGCCGAGGAAGGCGTAGGTGTAGTGGTCGa                                               >  1:1090580/1‑40 (MQ=38)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCg                                       >  1:895493/1‑48 (MQ=255)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCg                                       >  1:2106534/1‑48 (MQ=255)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCg                                       >  1:2457194/1‑48 (MQ=255)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                 >  1:1338810/1‑54 (MQ=255)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                 >  1:2238623/1‑54 (MQ=255)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTgcgc                >  1:1063905/1‑71 (MQ=255)
                                                     atGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTgcgc                >  1:2358366/1‑71 (MQ=255)
                                                            cGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGc         >  1:2413018/1‑71 (MQ=255)
                                                                gCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGtt      <  1:657213/70‑1 (MQ=255)
                                                                 cccATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGtttt    <  1:61412/71‑1 (MQ=255)
                                                                 cccATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGtttt    <  1:1298228/71‑1 (MQ=255)
                                                                   cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTa                        <  1:801916/49‑1 (MQ=255)
                                                                   cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTa                        <  1:2632289/49‑1 (MQ=255)
                                                                   cATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGtttttt  <  1:470905/71‑1 (MQ=255)
                                                                   |                                                                      
CCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTT  >  minE/1968104‑1968241

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: