Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,247,796 T→C 51.0% G503G (GGT→GGC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,247,7960TC51.0% 5.1 / 92.5 106G503G (GGT→GGCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base C (31/23);  total (58/48)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGA  >  minE/2247728‑2247861
                                                                    |                                                                 
ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGt                                                                   >  1:2283559/1‑69 (MQ=255)
ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGt                                                                   >  1:263650/1‑69 (MQ=255)
ggTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGt                                                                   >  1:3188762/1‑69 (MQ=255)
ggTGCGGCGATGATTGGCTGCTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGt                                                                   >  1:1261034/1‑69 (MQ=255)
 gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTc                                                                  >  1:600658/1‑69 (MQ=255)
 gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTc                                                                  >  1:872281/1‑69 (MQ=255)
 gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTc                                                                  >  1:2752563/1‑69 (MQ=255)
 gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTc                                                                  >  1:2952108/1‑69 (MQ=255)
 gTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTc                                                                  >  1:687796/1‑69 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAGGAGCATCTGGGCCTGc                                                               >  1:3009780/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:816384/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:708739/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:2849470/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:3320165/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:3434621/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:261194/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:102224/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:2216756/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:2114421/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:2015281/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:875263/1‑68 (MQ=255)
     ggcgATGATTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGc                                                               >  1:1767631/1‑68 (MQ=255)
      gcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                                             >  1:2228674/1‑69 (MQ=255)
      gcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                                             >  1:2644679/1‑69 (MQ=255)
      gcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                                             >  1:2776149/1‑69 (MQ=255)
      gcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                                             >  1:2151587/1‑69 (MQ=255)
      gcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                                                             >  1:148201/1‑69 (MQ=255)
      gcgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGcc                                                              >  1:314709/1‑68 (MQ=255)
       cgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                                                            <  1:832852/69‑1 (MQ=255)
       cgATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                                                            <  1:592876/69‑1 (MQ=255)
        gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                                                            <  1:183725/68‑1 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAGCGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:87443/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2794868/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:817453/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:848249/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2930824/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:864697/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:275909/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2709129/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:889290/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:3183121/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2585097/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:1252773/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2218583/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2102769/1‑67 (MQ=255)
          tgatTGGCTGGTTTGGTTGTGCAATGCTCTGCTACGTAACGCCGAAAGAGCATCTGGGCCTGCCaaa                                                           >  1:2074205/1‑67 (MQ=255)
               ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:3137992/69‑1 (MQ=255)
               ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:1751914/69‑1 (MQ=255)
               ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:874174/69‑1 (MQ=255)
                  tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAATAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                                 <  1:2570713/69‑1 (MQ=255)
                  tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                                 <  1:298364/69‑1 (MQ=255)
                  tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                                 <  1:1185692/69‑1 (MQ=255)
                  tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                                 <  1:455748/69‑1 (MQ=255)
                  tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                                 <  1:542468/69‑1 (MQ=255)
                  tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt                                                 <  1:848231/69‑1 (MQ=255)
                       tGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                                             >  1:845570/1‑68 (MQ=255)
                       tGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                                             >  1:1002472/1‑68 (MQ=255)
                       tGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                                             >  1:1233943/1‑68 (MQ=255)
                       tGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                                             >  1:2274385/1‑68 (MQ=255)
                       tGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa                                             >  1:1438205/1‑68 (MQ=255)
                                    ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:1602954/48‑1 (MQ=255)
                                    ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:786592/48‑1 (MQ=255)
                                    ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:670474/48‑1 (MQ=255)
                                    ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:3377719/48‑1 (MQ=255)
                                    ctctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:3424629/48‑1 (MQ=255)
                                      ctGTTACGTAACGCCAAAAGAGGATCTGGGTCTGCCCAATAaagaa                                                    <  1:1177201/46‑1 (MQ=255)
                                      ctGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa                                                    <  1:2616169/46‑1 (MQ=255)
                                                tcGCCGAAAGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2651775/67‑1 (MQ=255)
                                                aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt                   >  1:2027847/1‑69 (MQ=255)
                                                aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt                   >  1:2075410/1‑69 (MQ=255)
                                                aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt                   >  1:1923813/1‑69 (MQ=255)
                                                aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt                   >  1:3424539/1‑69 (MQ=255)
                                                aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt                   >  1:189717/1‑69 (MQ=255)
                                                aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt                   >  1:1731092/1‑69 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:98066/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:898183/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:821246/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:401369/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:1055898/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:1187109/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:1736096/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2392021/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2647695/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2783251/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:27962/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2845445/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2860474/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2917107/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:319057/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:3295808/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:329633/62‑1 (MQ=255)
                                                    cgaaaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:538773/62‑1 (MQ=255)
                                                       aaGAGCATCTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:2779373/61‑1 (MQ=255)
                                                         gagCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAc          <  1:2846067/69‑1 (MQ=255)
                                                         gagCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAc          <  1:696927/69‑1 (MQ=255)
                                                         gagCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAc          <  1:211815/69‑1 (MQ=255)
                                                         gagCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAc          <  1:159443/69‑1 (MQ=255)
                                                         gagCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCAc          <  1:635565/69‑1 (MQ=255)
                                                               cTGGGCCTGCCAAATAAAGAAGATGTTAAACAGGGGCTTATCACCTATAAGAt                    <  1:3189139/53‑1 (MQ=255)
                                                                    tCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGa  <  1:3359972/66‑1 (MQ=255)
                                                                    |                                                                 
GGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATTGCTGCCCACGCCGCTGA  >  minE/2247728‑2247861

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: