Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,258,414 G→T 12.1% intergenic (+85/+80) 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,4140GT12.1% 189.8 / 9.7 66intergenic (+85/+80)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base T (5/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

CCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCC  >  minE/1258346‑1258477
                                                                    |                                                               
ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGt                                                                <  1:3215171/70‑1 (MQ=255)
ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTa                                                               <  1:2321149/71‑1 (MQ=255)
  ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAg                                                              >  1:2153568/1‑70 (MQ=255)
   gATGCGGCGTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGc                                                            <  1:1830028/71‑1 (MQ=2)
    aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                           >  1:1344960/1‑71 (MQ=255)
    aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                                                           >  1:2502827/1‑71 (MQ=255)
          cgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGAt                                                       >  1:1902686/1‑69 (MQ=255)
          cgTGAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAAGCCTGAt                                                       >  1:1839011/1‑69 (MQ=2)
           gTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                                                    <  1:2098586/71‑1 (MQ=255)
            tGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                 <  1:28413/57‑1 (MQ=255)
            tGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGt                                                                <  1:2089834/58‑1 (MQ=255)
             gAACGCCTTATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGATAAGcg                                                  >  1:846600/1‑70 (MQ=2)
              aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCa                                                  <  1:1560218/70‑1 (MQ=255)
                cGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAgc                                              <  1:3361755/71‑1 (MQ=9)
                cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAg                                               >  1:3384363/1‑71 (MQ=255)
                  ccTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                             >  1:2996392/1‑71 (MQ=255)
                    ttATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                 <  1:2190403/49‑1 (MQ=255)
                    ttATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                 <  1:2185908/49‑1 (MQ=255)
                      aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                          >  1:3007148/1‑69 (MQ=9)
                      aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                          >  1:2934480/1‑69 (MQ=9)
                      aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                          <  1:1781149/69‑1 (MQ=9)
                      aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt                                          <  1:2422172/69‑1 (MQ=9)
                       tCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt                                          >  1:714290/1‑69 (MQ=255)
                          ggCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGAGTAGCGCATCAGGct                                    >  1:3661923/1‑70 (MQ=9)
                           gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                 <  1:527739/42‑1 (MQ=255)
                           gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                                                 <  1:2533442/42‑1 (MQ=255)
                           gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAgg                                      >  1:3797587/1‑69 (MQ=255)
                            ccTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAg                                   >  1:2191031/1‑71 (MQ=255)
                              tACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                             <  1:2401853/59‑1 (MQ=255)
                              tACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                             <  1:1502178/59‑1 (MQ=255)
                               aCAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                                                    <  1:1940867/51‑1 (MQ=255)
                                cAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                                                    >  1:58022/1‑49 (MQ=9)
                                cAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                                                    >  1:3090375/1‑49 (MQ=9)
                                cAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCa                                        >  1:2874553/1‑62 (MQ=255)
                                 aaaaCCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTa                              >  1:3594583/1‑71 (MQ=255)
                                  aaaCCTTGCTAATTCAATATATTGCAGGGACTATg                                                                 <  1:1095669/35‑1 (MQ=255)
                                     ccTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGtt                          >  1:2534001/1‑71 (MQ=255)
                                        tGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGAt                                                       >  1:1734014/1‑39 (MQ=255)
                                         gCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGttt                         >  1:3749742/1‑68 (MQ=255)
                                         gCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGc                       <  1:2375998/70‑1 (MQ=255)
                                          cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcgc                                            >  1:1255422/1‑48 (MQ=255)
                                          cTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCAt                     <  1:3446304/71‑1 (MQ=255)
                                            aaTTCAATATATTGCAGGGACTATGTAGGCCTGATAAg                                                    <  1:1387976/38‑1 (MQ=255)
                                             aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                    <  1:54899/69‑1 (MQ=255)
                                             aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                    <  1:2988513/69‑1 (MQ=255)
                                             aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                  <  1:3152461/71‑1 (MQ=255)
                                             aTTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                  >  1:791146/1‑71 (MQ=255)
                                              ttCAATATATTTAATGGGTTATGTAGACCTGATAAGCAt                                                 <  1:287292/39‑1 (MQ=9)
                                              ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATa                    >  1:594385/1‑68 (MQ=255)
                                              ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAcc                 <  1:1051038/71‑1 (MQ=255)
                                              ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAcc                 >  1:3091160/1‑71 (MQ=255)
                                               tCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCtt                                >  1:3411357/1‑55 (MQ=255)
                                                cAATATATTTAATGGGTTATGTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGt                           >  1:1087146/1‑59 (MQ=9)
                                                cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                >  1:1664035/1‑70 (MQ=255)
                                                cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                <  1:3869756/70‑1 (MQ=255)
                                                cAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCt                <  1:3874387/70‑1 (MQ=255)
                                                  atatatTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                                             <  1:2482192/39‑1 (MQ=255)
                                                   tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCttt                               >  1:199104/1‑52 (MQ=255)
                                                   tatatTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCt                                 >  1:1520544/1‑50 (MQ=255)
                                                       ttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc        >  1:874143/1‑71 (MQ=255)
                                                       ttGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc        <  1:453880/71‑1 (MQ=255)
                                                        tGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGccc        >  1:2095191/1‑70 (MQ=255)
                                                          cAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGt      <  1:944924/70‑1 (MQ=255)
                                                          cAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGt      <  1:3854815/70‑1 (MQ=255)
                                                           aGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGttt    <  1:3293059/71‑1 (MQ=255)
                                                             ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACg                            <  1:3300766/45‑1 (MQ=255)
                                                             ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTcc  <  1:1931014/71‑1 (MQ=255)
                                                             ggACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTcc  <  1:3560809/71‑1 (MQ=255)
                                                                  aTGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                  >  1:2049194/1‑50 (MQ=255)
                                                                  aTGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAAc                  >  1:1884088/1‑50 (MQ=255)
                                                                   tGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgcgc          >  1:791239/1‑57 (MQ=255)
                                                                   tGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAgcgc          >  1:1514191/1‑57 (MQ=255)
                                                                    gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGTTTgt                       <  1:409011/43‑2 (MQ=11)
                                                                    gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGTTTgt                       <  1:2320350/43‑2 (MQ=11)
                                                                    gTAGACCTGATAAGCATAGCGCATCAGGCTGCTTTGCGTTTgt                       <  1:3617267/43‑2 (MQ=11)
                                                                    |                                                               
CCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCC  >  minE/1258346‑1258477

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: