Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 572,826 A→G 61.1% intergenic (‑80/‑346) aqpZ ← / → ybjD aquaporin/conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE572,8260AG61.1% 18.7 / 23.5 41intergenic (‑80/‑346)aqpZ/ybjDaquaporin/conserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base A (12/4);  new base G (18/7);  total (30/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

TGAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGATAAATAAATCCAGTAAATTTGATTTA  >  minE/572774‑572894
                                                    |                                                                    
tgAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGc                                                             >  1:1286780/1‑62 (MQ=255)
tgAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGc                                                             >  1:2689318/1‑62 (MQ=255)
tgAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGc                                                             >  1:3169817/1‑62 (MQ=255)
tgAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGc                                                             >  1:2244298/1‑62 (MQ=255)
tgAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGc                                                             >  1:3030139/1‑62 (MQ=255)
tgAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGAGGc                                                             >  1:1932982/1‑62 (MQ=255)
        atGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGCTATTTTAAAAACGaa                                              <  1:859167/69‑1 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAaataat                                                                    >  1:3344582/1‑46 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATgc                                                                  >  1:2678306/1‑48 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCgg                                                              >  1:802801/1‑52 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:3062615/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:2518011/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:428458/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:2864734/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:245910/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:319831/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:3276441/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:3382717/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:590379/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:2283792/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:707268/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:192984/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:187141/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:1373508/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:903631/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATGATGCGCGGCTa                                                           >  1:974858/1‑55 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATg                                                                   >  1:1000151/1‑47 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATg                                                                   >  1:2017615/1‑47 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATg                                                                   >  1:712349/1‑47 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATg                                                                   >  1:420749/1‑47 (MQ=255)
         tGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATg                                                                   >  1:1341407/1‑47 (MQ=255)
                                 aaaaaGTGACGACGGAAATAATGCGCGGCTATTTTAAGAACGAAGGCGAGTCa                                     >  1:2840427/1‑53 (MQ=255)
                                 aaaaaGTGACGACGGAAATAATGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGataaataa                     <  1:3380732/69‑1 (MQ=255)
                                 aaaaaGTGACGACGGAAATAATGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGataaataa                     <  1:455087/69‑1 (MQ=255)
                                 aaaaaGTGACGACGGAAATAATGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGataaataa                     <  1:1123188/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:3098008/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:2486698/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:2242126/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:1909399/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:1513776/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:876929/69‑1 (MQ=255)
                                              ggAAATGATGCGCGGCTATTTTAAAAACGAAGGCAAGTCATTCACCAGATAAATAAATCCAGTAAAttt        <  1:1247694/69‑1 (MQ=255)
                                                  ataatGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGATAAATAAATCCAGTAAATTTGAtt    <  1:1651359/69‑1 (MQ=255)
                                                    aatGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGATAAATAAATCCAGTAAATTTGATTTa  >  1:2801147/1‑69 (MQ=255)
                                                    |                                                                    
TGAAAAATATGGTCGTTTTATAGGGTCGTCGTAAAAAAGTGACGACGGAAATAATGCGCGGCTATTTTAAAAACGAAGGCGAGTCATTCACCAGATAAATAAATCCAGTAAATTTGATTTA  >  minE/572774‑572894

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