Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,716,833 C→A 13.3% intergenic (+91/+107) 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,8330CA13.3% 226.7 / 14.9 75intergenic (+91/+107)yphH/glyApredicted DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (29/36);  new base A (2/8);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA  >  minE/1716764‑1716901
                                                                     |                                                                     
cTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACg                                                                      >  1:306367/1‑71 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:3790146/70‑1 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:1661947/70‑1 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:3585983/70‑1 (MQ=255)
   cccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:1290254/70‑1 (MQ=255)
    ccGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgag                                                                   <  1:3879586/70‑1 (MQ=255)
     cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:3261819/68‑1 (MQ=255)
     cGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGAt                                                                <  1:1044909/71‑1 (MQ=255)
            agCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAata                                                         >  1:2977423/1‑71 (MQ=255)
             gCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:4078590/60‑1 (MQ=255)
              ccTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAatata                                                       >  1:3051359/1‑71 (MQ=255)
                ttGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGga                                                                    <  1:2090217/57‑1 (MQ=255)
                 tgatgaTCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATAtt                                                     >  1:1846672/1‑70 (MQ=255)
                  gatgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTg                                                    <  1:3964185/70‑1 (MQ=255)
                   atgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCaa                                                           >  1:1685762/1‑62 (MQ=255)
                   atgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCaa                                                           >  1:3674266/1‑62 (MQ=255)
                    tgatCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAAt                                                 >  1:1039740/1‑71 (MQ=255)
                     gatCATAAACGAAAAATTGCCTGATGCGCTGCGCTTATCAGGCCTACACGgag                                                                   <  1:301944/53‑1 (MQ=255)
                       tCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                              <  1:1006996/71‑1 (MQ=255)
                        cATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                              >  1:2701873/1‑70 (MQ=255)
                        cATAAAAGAAAAATTGCCTGATGCTCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTg                                              >  1:1212670/1‑70 (MQ=255)
                         aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACTga                                                                    <  1:3571552/48‑1 (MQ=255)
                         aTAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAata                                                         >  1:1481500/1‑58 (MQ=255)
                          tAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAatata                                                       <  1:3940519/59‑1 (MQ=255)
                           aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCa                                            >  1:2118127/1‑69 (MQ=255)
                           aaaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCa                                            >  1:2461888/1‑69 (MQ=255)
                            aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGgaga                                                                 >  1:194831/1‑47 (MQ=255)
                            aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAg                                         >  1:3265183/1‑71 (MQ=255)
                            aaCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAg                                         >  1:1661596/1‑71 (MQ=255)
                             aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGt                                        <  1:1174803/71‑1 (MQ=255)
                             aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGt                                        <  1:2285852/71‑1 (MQ=255)
                             aCGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGt                                        <  1:1419646/71‑1 (MQ=255)
                               gAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtt                                       <  1:1172683/70‑1 (MQ=255)
                                 caaaTTGCCTGATGCGCTACGCTTATCGGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                    <  1:585976/70‑1 (MQ=11)
                                 caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                    <  1:1180270/70‑1 (MQ=11)
                                 caaaTTGCCTGATGCGCTACGCTTATCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTg                                    <  1:2980459/70‑1 (MQ=11)
                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTg                                    <  1:3799106/71‑1 (MQ=255)
                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTg                                    <  1:2767673/71‑1 (MQ=255)
                                 aaaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTg                                    <  1:2501861/71‑1 (MQ=255)
                                  aaaTTGCCTGATGCGCTGCGCTTATCAGGCCTACAAGTTGAATTGCAATTTATTGAATTTGCACATTTTTg                                    <  1:1206041/71‑1 (MQ=2)
                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTTCAATATATTGAATTTGCAAAGTTTTGt                                   <  1:3970829/71‑1 (MQ=255)
                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTg                                    >  1:1864276/1‑70 (MQ=255)
                                  aaaTTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGt                                   <  1:811366/71‑1 (MQ=255)
                                       gCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGcc                              >  1:911080/1‑71 (MQ=255)
                                        ccTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                             <  1:61138/71‑1 (MQ=255)
                                         cTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg                            <  1:1305034/71‑1 (MQ=255)
                                           gATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                           >  1:1815686/1‑70 (MQ=255)
                                           gATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGa                           >  1:3293141/1‑70 (MQ=255)
                                           gATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGa                           >  1:834745/1‑70 (MQ=11)
                                            aTGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATa                         <  1:1412157/71‑1 (MQ=255)
                                            aTGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCGGAt                          <  1:1877370/70‑1 (MQ=11)
                                                 cTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCg                    <  1:346244/71‑1 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCaaa                                          >  1:1564827/1‑47 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGtttt                                     >  1:1861805/1‑52 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                  >  1:3359635/1‑71 (MQ=255)
                                                   aCGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                  >  1:2876386/1‑71 (MQ=255)
                                                    tgCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                    <  1:3181163/51‑1 (MQ=11)
                                                    tgCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTg                                    <  1:1662940/51‑1 (MQ=11)
                                                    cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTc                 <  1:843926/71‑1 (MQ=255)
                                                      cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGttt                                      >  1:2490392/1‑48 (MQ=255)
                                                      cTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc               >  1:2571736/1‑71 (MQ=255)
                                                        tATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCact                                          >  1:3949887/1‑40 (MQ=11)
                                                         aTCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgc             <  1:2465711/70‑1 (MQ=255)
                                                         aTCAGGCCTACAAGTGCATTGCAATCTATTGAATTTGCGATGTTTTGTAGGACGGATAAGGCGCTcgcg              <  1:2009239/69‑1 (MQ=2)
                                                            aGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccg           >  1:181541/1‑69 (MQ=255)
                                                              gCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc       >  1:528223/1‑71 (MQ=255)
                                                              gCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATc       <  1:2128009/71‑1 (MQ=255)
                                                                cTACACGGAGATTGCAATATATTGAATTTGCAAAGttt                                      <  1:1282849/38‑1 (MQ=255)
                                                                cTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGtt                  >  1:939185/1‑58 (MQ=255)
                                                                cTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg     >  1:1275493/1‑71 (MQ=255)
                                                                 tACACGGAGATTGCAATATATTGAATTTGCAAAGttt                                      <  1:3633273/37‑1 (MQ=255)
                                                                 tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCgg                            <  1:1101214/47‑1 (MQ=255)
                                                                 tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCg                             <  1:3525913/46‑1 (MQ=255)
                                                                 tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                          <  1:1665877/49‑1 (MQ=255)
                                                                 tACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAg                       <  1:798626/52‑1 (MQ=255)
                                                                  acacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGAt                          <  1:575918/48‑1 (MQ=255)
                                                                  acacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg    >  1:2866669/1‑70 (MQ=255)
                                                                   cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCa                >  1:2585801/1‑57 (MQ=255)
                                                                   cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAcgccgc          >  1:3393738/1‑63 (MQ=255)
                                                                   cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg    <  1:307968/69‑1 (MQ=255)
                                                                   cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa  <  1:413119/71‑1 (MQ=255)
                                                                   cacGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa  >  1:829181/1‑71 (MQ=255)
                                                                     |                                                                     
CTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCA  >  minE/1716764‑1716901

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: