Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,258,403 G→T 14.4% intergenic (+74/+91) yecP → / ← torZ predicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,258,4030GT14.4% 229.7 / 18.0 95intergenic (+74/+91)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (44/37);  new base T (3/11);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 4.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCC  >  minE/1258333‑1258471
                                                                      |                                                                    
tCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                                      <  1:4844379/71‑1 (MQ=255)
tCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                                      <  1:1485836/71‑1 (MQ=255)
  tGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                                      >  1:3480889/1‑69 (MQ=255)
    gTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                                                  <  1:4004945/71‑1 (MQ=255)
       ttACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAc                                                               >  1:4767619/1‑71 (MQ=255)
             ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                           >  1:3734384/1‑69 (MQ=255)
             ccGGATGCGGCGTGAACGCCTTATCCGACCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                           <  1:4526893/69‑1 (MQ=255)
               ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAgg                                                       <  1:3870347/71‑1 (MQ=255)
                 aTGCGGCGTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGcc                                                     >  1:1312535/1‑71 (MQ=2)
                    cggcgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTg                                                   <  1:3147214/70‑1 (MQ=255)
                      gcgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATaa                                               <  1:3969092/71‑1 (MQ=2)
                       cgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                                              >  1:2912678/1‑71 (MQ=9)
                       cgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                                              <  1:163709/71‑1 (MQ=9)
                       cgTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGAt                                                 >  1:1196526/1‑69 (MQ=2)
                       cgTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGAt                                                 >  1:5601232/1‑69 (MQ=2)
                        gTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAgg                                                       >  1:4938926/1‑62 (MQ=2)
                        gTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAgg                                                       >  1:88908/1‑62 (MQ=2)
                         tGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCa                                                                    <  1:2533891/48‑1 (MQ=255)
                         tGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                           <  1:3963067/57‑1 (MQ=255)
                           aaCGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGcg                                            >  1:5304301/1‑68 (MQ=9)
                           aaCGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGcg                                            >  1:1237008/1‑68 (MQ=9)
                           aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                                                  >  1:5130665/1‑48 (MQ=255)
                           aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                                                  <  1:3962078/48‑1 (MQ=255)
                           aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                                                  >  1:4206099/1‑48 (MQ=255)
                           aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCAt                                           <  1:4433128/71‑1 (MQ=255)
                           aaCGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGATAAGcg                                            >  1:1345539/1‑69 (MQ=2)
                            aCGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAg                                         >  1:4359043/1‑71 (MQ=9)
                            aCGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGATAAGcg                                            >  1:4030252/1‑68 (MQ=2)
                             cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGa                                                                >  1:5654816/1‑48 (MQ=255)
                             cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGa                                                                >  1:106423/1‑48 (MQ=255)
                               ccTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                       >  1:1787019/1‑71 (MQ=255)
                               ccTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                       >  1:3335205/1‑71 (MQ=255)
                                 ttATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                       <  1:63454/69‑1 (MQ=255)
                                 ttATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                       <  1:955393/69‑1 (MQ=255)
                                 ttATCCGACCCACAAAGTCTTGCAAAATCAATAAATAGCAGGAACGATGTAGGCCTGATa                                                >  1:935178/1‑60 (MQ=2)
                                   aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                    >  1:1791272/1‑69 (MQ=9)
                                   aTCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCa                                     >  1:4474839/1‑69 (MQ=255)
                                   aTCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                    >  1:4121435/1‑70 (MQ=255)
                                   aTCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                    >  1:5182346/1‑70 (MQ=255)
                                    tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                                              <  1:144553/43‑3 (MQ=14)
                                    tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                                              <  1:4037079/43‑3 (MQ=14)
                                    tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                                              <  1:3417786/43‑3 (MQ=14)
                                    tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                                              <  1:3424206/43‑3 (MQ=14)
                                    tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCATCAg                                 <  1:4437046/71‑1 (MQ=9)
                                        gCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCATCa                                  >  1:3288905/1‑66 (MQ=9)
                                        gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                                              <  1:920956/55‑1 (MQ=255)
                                        gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                                              <  1:3345577/55‑1 (MQ=255)
                                           tACAAAACCTTGCTAATTCAATATATTGCAGGGACTAt                                                            <  1:4967192/38‑1 (MQ=255)
                                           tACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCa                                     <  1:2978919/61‑1 (MQ=255)
                                            aCAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCAt                                           <  1:1561115/54‑1 (MQ=255)
                                            aCAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                          >  1:4506173/1‑71 (MQ=255)
                                             cAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCt                                                    <  1:3539309/44‑1 (MQ=255)
                                             cAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCt                                                    <  1:2390062/44‑1 (MQ=255)
                                              aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                    <  1:3053980/59‑1 (MQ=255)
                                              aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                    <  1:959986/59‑1 (MQ=255)
                                              aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCttt                         <  1:1847400/70‑1 (MQ=255)
                                                aaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGc                                             >  1:1623060/1‑48 (MQ=255)
                                                 aCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGt                     <  1:2336628/71‑1 (MQ=255)
                                                 aCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGt                     <  1:4203399/71‑1 (MQ=255)
                                                  ccTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGtt                    >  1:1082953/1‑71 (MQ=255)
                                                     tGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                   >  1:1692566/1‑69 (MQ=255)
                                                      gCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAgg                                >  1:3310184/1‑55 (MQ=255)
                                                      gCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAgg                                >  1:2003473/1‑55 (MQ=255)
                                                       cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt               <  1:3921158/71‑1 (MQ=255)
                                                       cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt               <  1:3633435/71‑1 (MQ=255)
                                                        tAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt               >  1:3643659/1‑70 (MQ=255)
                                                        tAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt               >  1:2458677/1‑70 (MQ=255)
                                                        aagttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGttt                   >  1:2317224/4‑66 (MQ=9)
                                                         aaTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                          >  1:3091987/1‑58 (MQ=255)
                                                         aaTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                          >  1:831128/1‑58 (MQ=255)
                                                          aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATaa             >  1:4601760/1‑70 (MQ=255)
                                                           ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGctg                             <  1:4798301/53‑3 (MQ=9)
                                                           ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGctg                             <  1:5610084/53‑3 (MQ=9)
                                                           ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGctg                             <  1:1985576/53‑3 (MQ=9)
                                                           ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGctg                             <  1:4791809/53‑3 (MQ=9)
                                                           ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCtt                          <  1:1378907/56‑1 (MQ=9)
                                                           ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCtt                          <  1:5190587/56‑1 (MQ=9)
                                                           ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCa                                     <  1:2500841/45‑1 (MQ=255)
                                                            tCAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGgctgct                           <  1:2949309/54‑1 (MQ=9)
                                                             cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACg                      >  1:1236330/1‑58 (MQ=255)
                                                             cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt          <  1:943738/70‑1 (MQ=255)
                                                              aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAgcgc                                      <  1:5136185/41‑1 (MQ=12)
                                                              aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAgcgc                                      <  1:5605311/41‑1 (MQ=12)
                                                              aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAgcgc                                      <  1:5038193/41‑1 (MQ=12)
                                                              aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAgcgc                                      <  1:3025246/41‑1 (MQ=12)
                                                              aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                     >  1:3009687/1‑58 (MQ=9)
                                                              aaTATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                     >  1:1525217/1‑58 (MQ=9)
                                                               atatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAg       <  1:2701218/71‑1 (MQ=255)
                                                               atatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAg       <  1:357589/71‑1 (MQ=255)
                                                                tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                    <  1:785787/41‑1 (MQ=255)
                                                                tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCa        >  1:948591/1‑69 (MQ=255)
                                                                 atatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                  >  1:913594/1‑42 (MQ=255)
                                                                  tatTGGAGGGACTATGTAGGCCTGATAAGCATAGCGCATc                                   >  1:3183619/1‑40 (MQ=255)
                                                                  tatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATc                                   >  1:5316275/1‑40 (MQ=255)
                                                                  tatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATc                                   >  1:618650/1‑40 (MQ=255)
                                                                  tatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATc                                   >  1:2292898/1‑40 (MQ=255)
                                                                   atTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgcgc    <  1:3789462/70‑1 (MQ=255)
                                                                    ttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc  >  1:5067378/1‑71 (MQ=255)
                                                                      |                                                                    
TCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCC  >  minE/1258333‑1258471

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: