Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,716,861 A→T 39.4% intergenic (+119/+79) yphH → / ← glyA predicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,716,8610AT39.4% 99.3 / 75.2 66intergenic (+119/+79)yphH/glyApredicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base T (10/16);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAC  >  minE/1716793‑1716931
                                                                     |                                                                      
aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGt                                                                      >  1:2554609/1‑71 (MQ=255)
 cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtt                                                                     >  1:1456999/1‑71 (MQ=255)
 cGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGtt                                                                     >  1:3195851/1‑71 (MQ=255)
    caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTg                                                                  <  1:2660934/70‑1 (MQ=11)
    aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                                                  <  1:2893220/71‑1 (MQ=255)
    aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                                                  >  1:1727867/1‑71 (MQ=255)
     agcttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg                                                                  <  1:3602840/67‑1 (MQ=11)
     aaaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTg                                                                  <  1:1597560/71‑1 (MQ=2)
     aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGt                                                                 <  1:1290110/71‑1 (MQ=255)
      gcttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTg                                                                  >  1:3071000/3‑69 (MQ=11)
      aaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTGt                                                                 <  1:3293701/71‑1 (MQ=0)
      aaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTa                                                                <  1:2306066/71‑1 (MQ=255)
      aaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTa                                                                <  1:2052954/71‑1 (MQ=255)
       cttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAg                                                               <  1:1647948/70‑1 (MQ=11)
       cttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAg                                                               <  1:1199022/70‑1 (MQ=11)
       cttGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAg                                                               <  1:198412/70‑1 (MQ=11)
        ttGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAgg                                                              <  1:2345146/71‑1 (MQ=255)
        ttGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAgg                                                              <  1:1777636/71‑1 (MQ=11)
          gCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc                                                            >  1:533538/1‑71 (MQ=255)
          gCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc                                                            >  1:318986/1‑71 (MQ=255)
            cTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg                                                          <  1:2518958/71‑1 (MQ=255)
            cTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                                                           >  1:1166864/1‑70 (MQ=255)
             tGATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa                                                         >  1:896482/1‑71 (MQ=11)
              gATGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa                                                         >  1:2193709/1‑70 (MQ=11)
               aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                                                        >  1:1539611/1‑70 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                                                        <  1:1708868/70‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                                       <  1:3098876/71‑1 (MQ=255)
               aTGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                                                       <  1:3632001/71‑1 (MQ=255)
               aTGCGCGATGCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAt                                                        <  1:2435653/70‑1 (MQ=11)
                tGCGCTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGAt                                                        >  1:83255/1‑69 (MQ=11)
                 gcgcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTg                                                                  <  1:1786627/58‑1 (MQ=255)
                  cgcTACGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAg                                                     <  1:920009/70‑1 (MQ=11)
                   gcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                                   <  1:2621248/71‑1 (MQ=255)
                   gcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                                   >  1:2081729/1‑71 (MQ=255)
                   gcTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGc                                                   <  1:3736786/71‑1 (MQ=255)
                     tACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                                                     >  1:3572999/1‑67 (MQ=255)
                     gatgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                                                 <  1:3627719/68‑1 (MQ=11)
                      atgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt                                                <  1:76576/69‑1 (MQ=11)
                      atgCTTATCAGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGtt                                                <  1:3319468/69‑1 (MQ=11)
                       cGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                                                <  1:650953/70‑1 (MQ=255)
                       cGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTc                                               <  1:3132762/71‑1 (MQ=255)
                       cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATaa                                                      >  1:1928570/1‑64 (MQ=2)
                       cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATaa                                                      >  1:2901644/1‑64 (MQ=2)
                       cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgc                                                 <  1:152775/69‑2 (MQ=11)
                       cGCTTATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTc                                               >  1:1089706/1‑71 (MQ=2)
                          ttATCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg                                            >  1:3753857/1‑71 (MQ=2)
                            aTCAGGCCTACAAGTGC‑ATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCgct                                                <  1:1792260/65‑3 (MQ=2)
                             tCAGGCCTACACGGAG‑ATTGCAATATATTTAATTTGCAAAGTTTTGt                                                                 >  1:108144/1‑47 (MQ=255)
                             tCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTa                                                                >  1:3556636/1‑48 (MQ=255)
                               aGGCCTACCAGAAG‑ATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGt                                                 >  1:2792880/1‑61 (MQ=11)
                                         ggAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                      >  1:3395276/1‑62 (MQ=255)
                                         ggAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                      >  1:2561440/1‑62 (MQ=255)
                                          aag‑aTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCg                                                           <  1:2694764/39‑1 (MQ=11)
                                           ag‑aTTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaac                           <  1:1855622/71‑1 (MQ=11)
                                            g‑aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcg                                            >  1:2952213/1‑53 (MQ=255)
                                             caTTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTCGCGCCGCATCCGGCATCaaca                          <  1:1112107/70‑1 (MQ=11)
                                               ttGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATGaacaat                        >  1:345471/1‑70 (MQ=11)
                                                      atatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca                 >  1:736985/1‑71 (MQ=255)
                                                      atatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca                 >  1:2907564/1‑71 (MQ=255)
                                                       tatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGcaca                 <  1:2241551/70‑1 (MQ=255)
                                                        atTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAt                <  1:3313617/70‑1 (MQ=255)
                                                         ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAtt               <  1:376076/70‑1 (MQ=255)
                                                           gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACAtt               >  1:235304/1‑68 (MQ=255)
                                                           gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa             <  1:3779635/70‑1 (MQ=255)
                                                               ttGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa           >  1:2126778/1‑68 (MQ=255)
                                                                  cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCaaa      <  1:2845940/70‑1 (MQ=255)
                                                                   aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTg              <  1:588361/61‑1 (MQ=255)
                                                                     aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc  <  1:2498404/71‑1 (MQ=255)
                                                                     aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAc  >  1:692247/1‑71 (MQ=255)
                                                                     |                                                                      
ACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAG‑ATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAC  >  minE/1716793‑1716931

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: ‑