Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,258,382 A→T 25.0% intergenic (+53/+112) 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,3820AT25.0% 145.8 / 46.3 77intergenic (+53/+112)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base T (8/10);  total (37/40)
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.05e-01

TGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGT  >  minE/1258314‑1258452
                                                                    |                                                                      
tGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAcc                                                                      >  1:3975321/1‑71 (MQ=255)
         gCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                             >  1:4526559/1‑71 (MQ=255)
          ccGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTc                                                            <  1:1501724/71‑1 (MQ=255)
            gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                             <  1:927460/68‑1 (MQ=255)
            gTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCaa                                                          <  1:5475036/71‑1 (MQ=255)
             tAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCaa                                                          >  1:412870/1‑70 (MQ=255)
               aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatat                                                       <  1:5024919/71‑1 (MQ=255)
               aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                        >  1:2229955/1‑70 (MQ=255)
               aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                        >  1:2076606/1‑70 (MQ=255)
                 ggTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatatat                                                     >  1:3478701/1‑71 (MQ=255)
                   tCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                   <  1:1485836/71‑1 (MQ=255)
                   tCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                   <  1:4844379/71‑1 (MQ=255)
                     tGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTg                                                   >  1:3480889/1‑69 (MQ=255)
                       gTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                               <  1:4004945/71‑1 (MQ=255)
                          ttACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                             <  1:2476714/54‑1 (MQ=255)
                          ttACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                        <  1:3528017/59‑1 (MQ=255)
                          ttACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAc                                            >  1:4767619/1‑71 (MQ=255)
                            aCTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                        <  1:1564739/57‑1 (MQ=255)
                              tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                        <  1:5606571/55‑1 (MQ=255)
                                ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                        >  1:3734384/1‑69 (MQ=255)
                                ccGGATGCGGCGTGAACGCCTTATCCGACCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                        <  1:4526893/69‑1 (MQ=255)
                                 cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                           <  1:1022153/49‑1 (MQ=255)
                                 cGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                           <  1:857020/49‑1 (MQ=255)
                                  ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAgg                                    <  1:3870347/71‑1 (MQ=255)
                                   gATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAtt                                                             <  1:2538819/45‑1 (MQ=255)
                                    aTGCGGCGTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGcc                                  >  1:1312535/1‑71 (MQ=2)
                                     tGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                        >  1:887837/1‑48 (MQ=255)
                                       cggcgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTg                                <  1:3147214/70‑1 (MQ=255)
                                         gcgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAata                                                        <  1:5157068/44‑1 (MQ=11)
                                         gcgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATaa                            <  1:3969092/71‑1 (MQ=2)
                                          cgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAata                                                        <  1:2385923/43‑1 (MQ=14)
                                          cgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                           >  1:2912678/1‑71 (MQ=9)
                                          cgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                           <  1:163709/71‑1 (MQ=9)
                                          cgTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGAt                              >  1:5601232/1‑69 (MQ=2)
                                          cgTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGAt                              >  1:1196526/1‑69 (MQ=2)
                                           gTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCAatata                                                      >  1:2725210/1‑44 (MQ=11)
                                           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)
                                              aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCaa                                                          >  1:4045185/1‑37 (MQ=255)
                                              aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                               >  1:4206099/1‑48 (MQ=255)
                                              aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                               <  1:3962078/48‑1 (MQ=255)
                                              aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                               >  1:5130665/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)
                                                cGCCTTATCCGGCCTACAAATCCTTGCGAATTCAata                                                        <  1:5378064/37‑1 (MQ=14)
                                                cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGa                                             >  1:5654816/1‑48 (MQ=255)
                                                cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGa                                             >  1:106423/1‑48 (MQ=255)
                                                  ccTTATCCGGCCTACAAATCCTTGCGAATTCAata                                                        >  1:3253930/1‑35 (MQ=14)
                                                  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:4037079/43‑3 (MQ=14)
                                                       tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                           <  1:3424206/43‑3 (MQ=14)
                                                       tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                           <  1:3417786/43‑3 (MQ=14)
                                                       tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGaaa                                           <  1:144553/43‑3 (MQ=14)
                                                       tCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCATCAg              <  1:4437046/71‑1 (MQ=9)
                                                           gCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCATCa               >  1:3288905/1‑66 (MQ=9)
                                                           gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                           <  1:3345577/55‑1 (MQ=255)
                                                           gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                           <  1:920956/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:2390062/44‑1 (MQ=255)
                                                                cAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCt                                 <  1:3539309/44‑1 (MQ=255)
                                                                 aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                 <  1:959986/59‑1 (MQ=255)
                                                                 aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                 <  1:3053980/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)
                                                                    |                                                                      
TGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGT  >  minE/1258314‑1258452

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: