Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,030 T→G 7.4% intergenic (+85/‑95) lpcA → / → prfH D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE223,0300TG 7.4% 242.2 / 13.2 94intergenic (+85/‑95)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base T (40/47);  new base G (3/4);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGAAAATTCTT  >  minE/222960‑223099
                                                                      |                                                                     
ggACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGc                                                                     <  1:2257632/71‑1 (MQ=2)
  aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGc                                                                     >  1:368178/1‑71 (MQ=255)
     aggcGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                                                                  >  1:3539000/4‑71 (MQ=11)
     aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                   <  1:330049/67‑1 (MQ=11)
     aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                   <  1:1592322/67‑1 (MQ=11)
     aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                   <  1:1566096/67‑1 (MQ=11)
     aggcGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTa                                                                  >  1:2336579/4‑71 (MQ=11)
      gACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                                                                  >  1:1374841/1‑70 (MQ=255)
       aCGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCTTGCCTACa                                                                <  1:1527663/71‑1 (MQ=255)
        cGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaa                                                               <  1:2412485/71‑1 (MQ=14)
        cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa                                                                >  1:1412261/1‑70 (MQ=255)
         gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaaa                                                              >  1:770164/1‑71 (MQ=255)
          tttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt                                                             <  1:2860960/71‑1 (MQ=255)
           ttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                                                            >  1:2275068/1‑71 (MQ=255)
             aCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttcc                                                          >  1:2556266/1‑68 (MQ=11)
             aCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttct                                                          <  1:3662595/71‑4 (MQ=11)
               gccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTa                                                        >  1:323003/1‑71 (MQ=255)
                ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc                                                       <  1:3931188/71‑1 (MQ=255)
                 tgcATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                                  <  1:442082/57‑1 (MQ=14)
                  gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc                                                     >  1:1371132/1‑71 (MQ=255)
                  gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc                                                     >  1:2146681/1‑71 (MQ=255)
                    aTCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGa                                                   <  1:483843/71‑1 (MQ=255)
                     tCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAt                                                             >  1:747356/1‑60 (MQ=11)
                      ccGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt                                                             >  1:457429/1‑59 (MQ=255)
                         gCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCAtcc                                                                     >  1:2705439/1‑46 (MQ=255)
                          cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg                                               >  1:1868962/1‑69 (MQ=255)
                          cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg                                               >  1:2067074/1‑69 (MQ=255)
                          cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTa                                             <  1:2612615/71‑1 (MQ=255)
                           aTTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGt                                              <  1:381638/69‑1 (MQ=255)
                           aTTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAg                                            >  1:1276661/1‑71 (MQ=255)
                            actgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                                                   <  1:3002515/45‑1 (MQ=14)
                            actgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                                                   <  1:225542/45‑1 (MQ=14)
                             ttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGc                                          >  1:768215/1‑71 (MQ=255)
                             tcgtgCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACa                                                                <  1:905556/47‑1 (MQ=11)
                              cgtgCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACa                                                                <  1:430042/47‑1 (MQ=11)
                                tgCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCg                                        <  1:2390605/70‑1 (MQ=255)
                                 gcgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaa                                      >  1:2012727/1‑71 (MQ=255)
                                  cgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaat                                     <  1:3570070/71‑1 (MQ=255)
                                   gcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaata                                    <  1:2282270/71‑1 (MQ=255)
                                    ctgatgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaataa                                   >  1:3428121/1‑71 (MQ=255)
                                     ctatgcctgatgcgaCGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatc                                  <  1:2371030‑M1/69‑31 (MQ=255)
                                     ctatgcctgatgcgaCGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatc                                  <  1:2375456‑M1/69‑31 (MQ=255)
                                     ctatgcctgatgcgaCGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatc                                  <  1:3491278‑M1/69‑31 (MQ=255)
                                      tatgcctgatgcgACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatcc                                 >  1:790998‑M1/2‑40 (MQ=255)
                                      gatgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATg                                 >  1:1990384/1‑71 (MQ=255)
                                       atgcctgatgcGACGCTGACGCGTTTTATCATGCCTAc                                                                 <  1:2269318/38‑1 (MQ=255)
                                       atgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaa                                      <  1:2156567/65‑1 (MQ=255)
                                       atgcctgatgcGACGCTGAAGCGTCTTATCATGCCTAc                                                                 >  1:74542/1‑38 (MQ=255)
                                       atgcctgatgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccg                                >  1:243104‑M1/1‑39 (MQ=255)
                                        tgcctgatgcGACGCTGACGCGTTTTATCATGCCt                                                                   <  1:973100/35‑1 (MQ=255)
                                           ctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgcg                                                    <  1:1743280/47‑1 (MQ=255)
                                           ctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgcg                                                    <  1:3070658/47‑1 (MQ=255)
                                           ctgatgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcat                            <  1:3283024‑M1/71‑37 (MQ=255)
                                           ctgatgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggca                             >  1:1676041‑M1/1‑35 (MQ=255)
                                            tgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaa                                      >  1:2771634/1‑60 (MQ=255)
                                            tgatgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtt                                             <  1:1167118‑M1/53‑20 (MQ=255)
                                             gatgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgt                                                  >  1:731175‑M1/1‑33 (MQ=255)
                                              atgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcg                         <  1:618101/71‑1 (MQ=255)
                                              atgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaa                         <  1:2519465‑M1/71‑40 (MQ=255)
                                              atgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaa                         <  1:1089040‑M1/71‑40 (MQ=255)
                                              atgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaa                         <  1:1079391‑M1/71‑40 (MQ=255)
                                              atgcGACGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataa                          <  1:2603307‑M1/70‑39 (MQ=255)
                                               tgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcgc                        >  1:1892864/1‑71 (MQ=255)
                                               tgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcg                         >  1:1441439/1‑70 (MQ=255)
                                                gcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAtgcc                       >  1:3115940/1‑71 (MQ=255)
                                                  gACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcgccg                      >  1:596923/1‑70 (MQ=255)
                                                  gACGCTGAACCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaaaaca                     <  1:2384417‑M1/71‑44 (MQ=255)
                                                   aCGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCa                    <  1:72554/71‑1 (MQ=255)
                                                   aCGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCa                    >  1:3291355/1‑71 (MQ=255)
                                                   aCGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaaaacac                    <  1:1965463‑M1/71‑45 (MQ=255)
                                                   aCGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcata                           <  1:297767‑M1/64‑38 (MQ=255)
                                                   aCGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatcaggcata                           <  1:697423‑M1/64‑38 (MQ=255)
                                                   aCGCTGAAGCGTCTTATCATGCCTACAtagcactgccac                                                    <  1:1691388‑M1/39‑13 (MQ=255)
                                                    cGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccg                                >  1:1409815‑M1/1‑26 (MQ=255)
                                                    cGCTGAAGCGTCTTATCATGCCTACAtagcactgccacgtat                                                >  1:360096‑M1/1‑26 (MQ=255)
                                                      cTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATcc                 <  1:2250511/71‑1 (MQ=255)
                                                       tGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATcc                 >  1:2119144/1‑70 (MQ=255)
                                                       tGAAGCGTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaaaacacgcg                 >  1:3902961‑M1/1‑23 (MQ=255)
                                                         aagcgtcTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaaaacacgcgcac              >  1:1828998‑M1/3‑21 (MQ=255)
                                                         aagcgtcTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaaaacacgcgca               <  1:526000‑M1/68‑50 (MQ=255)
                                                          cgcgTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCAcgc                        >  1:2417102/1‑60 (MQ=255)
                                                            cgTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTg           <  1:2875335/71‑1 (MQ=255)
                                                             gTCTTATCATGCCTACAtagcactgccacgtatgtttacaccgcatccggcataaaaacacgcgcactttg          <  1:2519438‑M1/71‑55 (MQ=255)
                                                             gTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTg           >  1:1636809/1‑70 (MQ=255)
                                                             gTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTg           >  1:340662/1‑70 (MQ=255)
                                                              tCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaat                                     <  1:1711661/43‑1 (MQ=255)
                                                              tCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaat                                     <  1:1802677/43‑1 (MQ=255)
                                                              tCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGa          <  1:1035164/70‑1 (MQ=255)
                                                                ttATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaataa                                   >  1:3178201/1‑43 (MQ=255)
                                                                 tATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGa               >  1:1527205/1‑62 (MQ=255)
                                                                  aTCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCa                           >  1:428038/1‑49 (MQ=255)
                                                                  aTCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGaaaa       <  1:3833477/69‑1 (MQ=255)
                                                                  aTCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGaaaa       <  1:3731410/69‑1 (MQ=255)
                                                                  aTCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGaaaa       <  1:309547/69‑1 (MQ=255)
                                                                    cATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGAAAATTCt   >  1:1575060/1‑71 (MQ=255)
                                                                      tGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGAAAATTCtt  <  1:1095527/70‑1 (MQ=255)
                                                                      |                                                                     
GGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAATGCGTTCACGCCGCATCCGACCTGAAAATTCTT  >  minE/222960‑223099

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: