Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 449,111 A→G 100% G409G (GGT→GGC kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE449,1110AG83.5% 56.5 / 3.1 36G409G (GGT→GGCkdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base G (9/21);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 2.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTC  >  minE/449046‑449173
                                                                 |                                                              
aGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGaccacc                                                              >  1:2248874/1‑68 (MQ=255)
aGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGaccacc                                                              >  1:915084/1‑68 (MQ=255)
aGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGaccacc                                                              >  1:2631354/1‑68 (MQ=255)
         gCTGCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTc                                                     >  1:1646830/1‑68 (MQ=255)
         gCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTATGGAAGTGACCACCGTTAGCCTc                                                     >  1:182133/1‑68 (MQ=255)
         gCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGtt                                                           >  1:2476725/1‑62 (MQ=255)
         gCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTc                                                     >  1:2760592/1‑68 (MQ=255)
         gCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTc                                                     >  1:2196470/1‑68 (MQ=255)
         gCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTc                                                     >  1:1849517/1‑68 (MQ=255)
                                                  ggcaGGGAAGTGCCCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:443781/4‑69 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCATCATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:2225445/4‑69 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGCCAGAACCTTTACg                 >  1:637880/4‑63 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCa            >  1:1623003/4‑68 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:1754553/4‑69 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:2540300/4‑69 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:1759392/4‑69 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:1281331/4‑69 (MQ=255)
                                                  ggcaGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCAt           >  1:3212948/4‑69 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCATCATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1050649/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:977065/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:3147274/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:3268597/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:330032/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:3341385/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:3376470/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:409155/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:435481/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:908938/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:2832415/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1059267/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1317115/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1416924/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1547845/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1751327/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1922011/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:2354441/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:2761724/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGGCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:3025353/67‑1 (MQ=255)
                                                    caGGGAAGTGACCGCCGTTAGCCTCAACATGGCGCCGAATGGCATCGACAGAACCTTTACGGATCATg          <  1:1529007/67‑1 (MQ=255)
                                                              accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc  >  1:1807764/1‑66 (MQ=255)
                                                              accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc  >  1:162942/1‑66 (MQ=255)
                                                              accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc  >  1:3258227/1‑66 (MQ=255)
                                                              accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc  >  1:78362/1‑66 (MQ=255)
                                                              accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc  >  1:882696/1‑66 (MQ=255)
                                                              accaccGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTc  >  1:112406/1‑66 (MQ=255)
                                                                 |                                                              
AGCGGCGTGGCTCCCTGACGCGCAACCTGATCGACTTTTTGATCAACATCGGTAGGGAAGTGACCACCGTTAGCCTCAACATGGCGACGAATGGCATCGACAGAACCTTTACGGATCATGCGGTTGTC  >  minE/449046‑449173

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: