Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,258,382 A→T 19.4% 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,3820AT19.4% 141.1 / 27.3 63intergenic (+53/+112)yecP/torZpredicted methyltransferase/trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base T (5/7);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

ATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTA  >  minE/1258312‑1258449
                                                                      |                                                                   
aTTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                     >  1:3021810/1‑71 (MQ=255)
 ttGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAc                                                                    >  1:1390246/1‑71 (MQ=255)
  tGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCTGGCCTACAAAAcc                                                                   >  1:308100/1‑71 (MQ=255)
   gcgcgcAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCt                                                                  >  1:2389357/1‑71 (MQ=255)
    cgcgcAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCt                                                                  <  1:3302930/70‑1 (MQ=255)
    cgcgcAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCt                                                                  <  1:3626820/70‑1 (MQ=255)
       gcAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGc                                                               >  1:3292417/1‑70 (MQ=255)
                 aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGGCTACAAAACCTTGCTAATTCAata                                                     >  1:2532500/1‑70 (MQ=255)
                 aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAAt                                                           >  1:1865022/1‑64 (MQ=255)
                 aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatat                                                    <  1:3380771/71‑1 (MQ=255)
                 aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                     >  1:745899/1‑70 (MQ=255)
                 aaGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAata                                                     >  1:2257478/1‑70 (MQ=255)
                   ggTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatatat                                                  <  1:1030772/71‑1 (MQ=255)
                         gTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                            >  1:3930305/1‑71 (MQ=255)
                               cTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTa                                       <  1:3881990/70‑1 (MQ=255)
                                  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:2502827/1‑71 (MQ=255)
                                      aTGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGcc                               >  1:1344960/1‑71 (MQ=255)
                                        gcggcgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCa                                                        <  1:40830/44‑1 (MQ=11)
                                        gcggcgTGAACGCCTTATCCGGCCTACAAATCCTTGCGAATTCa                                                        <  1:2575894/44‑1 (MQ=11)
                                        gcggcgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                        >  1:810565/1‑44 (MQ=255)
                                        gcggcgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                        >  1:1357186/1‑44 (MQ=255)
                                        ccggcgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCa                                                        <  1:1813584/43‑1 (MQ=255)
                                            cgTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAc                                         >  1:2948248/1‑55 (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)
                                                aaCGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTgt                                               <  1:2042568/45‑2 (MQ=11)
                                                aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAatatat                                                  <  1:138719/42‑1 (MQ=255)
                                                aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAt                                      >  1:3255345/1‑54 (MQ=255)
                                                aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTAt                                      >  1:3302282/1‑54 (MQ=255)
                                                aaCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCa                      <  1:1560218/70‑1 (MQ=255)
                                                 aCGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTg                                                <  1:1621388/43‑1 (MQ=14)
                                                 aCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCa                                              <  1:440427/45‑1 (MQ=255)
                                                 aCGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCa                                              <  1:3812904/45‑1 (MQ=255)
                                                  cGCCTTATCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAgc                  <  1:3361755/71‑1 (MQ=9)
                                                  cGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAg                   >  1:3384363/1‑71 (MQ=255)
                                                    ccTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAgcg                 >  1:2996392/1‑71 (MQ=255)
                                                     cTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                            <  1:774345/43‑1 (MQ=255)
                                                      ttATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                     <  1:2185908/49‑1 (MQ=255)
                                                      ttATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                     <  1:2190403/49‑1 (MQ=255)
                                                        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt              >  1:3007148/1‑69 (MQ=9)
                                                        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt              >  1:2934480/1‑69 (MQ=9)
                                                        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt              <  1:2422172/69‑1 (MQ=9)
                                                        aTCCGGCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGCGTAGCGCAt              <  1:1781149/69‑1 (MQ=9)
                                                        aTCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAgg                                            <  1:2323138/40‑1 (MQ=255)
                                                         tCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCAt              >  1:714290/1‑69 (MQ=255)
                                                            ggCCTACAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAGAGTAGCGCATCAGGct        >  1:3661923/1‑70 (MQ=9)
                                                             gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                     <  1:2533442/42‑1 (MQ=255)
                                                             gCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATg                                     <  1:527739/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:3090375/1‑49 (MQ=9)
                                                                  cAAATCCTTGCGAATTCAATATATTGTATGGAAATTCAGGCCTGATAAg                        >  1:58022/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)
                                                                      |                                                                   
ATTGCGCGCAAGCCGTAAAGGTCTGGTATTACTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTA  >  minE/1258312‑1258449

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