Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 447,521 T→C 54.8% R65R (CGA→CGG kdpD ← fused sensory histidine kinase in two‑component regulatory system with KdpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE447,5210TC54.8% 5.0 / 30.9 42R65R (CGA→CGGkdpDfused sensory histidine kinase in two‑component regulatory system with KdpE
Reads supporting (aligned to +/- strand):  ref base T (6/13);  new base C (16/7);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAACAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCTTC  >  minE/447449‑447586
                                                                        |                                                                 
atGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAA‑‑‑‑CAGCCCTTCCAGCATGGCGGCGGTATCtttt                                                                   >  1:2066667/1‑69 (MQ=255)
atGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAA‑‑‑‑CAGCCCTTCCAGCATGGCGGCGGTATCtttt                                                                   >  1:270207/1‑69 (MQ=255)
atGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAA‑‑‑‑CAGCCCTTCCAGCATGGCGGCGGTATCtttt                                                                   >  1:1987021/1‑69 (MQ=255)
atGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAA‑‑‑‑CAGCCCTTCCAGCATGGCGGCGGTATCtttt                                                                   >  1:3106340/1‑69 (MQ=255)
atGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAA‑‑‑‑CAGCCCTTCCAGCATGGCGGCGGTATCtttt                                                                   >  1:524938/1‑69 (MQ=255)
               cGCCTGGCGTTTTAACGGCAGAACAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGt                                                        >  1:1361620/1‑69 (MQ=255)
                          ttAACGGCAGAACAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCTACCACg                                              >  1:1474324/1‑68 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:3065060/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2049308/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:3008477/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2959331/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2402362/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:3170829/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2244194/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2934089/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2105531/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2881784/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:2063331/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:1931892/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:414163/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:613646/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:1097519/1‑57 (MQ=255)
                           tAACGGCAGAACAGCCAGACCTTTCAGCATGGCGGCGGTATCTTTCCGACCGTGGGt                                                        >  1:1056969/1‑57 (MQ=255)
                                     aCAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCt                                                    >  1:2343130/1‑51 (MQ=255)
                                     aCAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAata                                  >  1:1926024/1‑69 (MQ=255)
                                     aCAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAata                                  >  1:937411/1‑69 (MQ=255)
                                                  ccAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTgcg                     >  1:755811/1‑69 (MQ=255)
                                                            ggcggTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTc            <  1:2045275/68‑1 (MQ=255)
                                                            ggcggTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTc            <  1:2229375/68‑1 (MQ=255)
                                                            ggcggTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTc            <  1:470540/68‑1 (MQ=255)
                                                            ggcggTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTc            <  1:593848/68‑1 (MQ=255)
                                                            ggcggTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTc            <  1:2599199/68‑1 (MQ=255)
                                                            ggcggTATCTTTCCTACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:1011434/69‑1 (MQ=255)
                                                            ggcggTATCTTTCCGACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:2683148/69‑1 (MQ=255)
                                                            ggcggTATCTTTCCGACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:105167/69‑1 (MQ=255)
                                                            ggcggTATCTTTCCGACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:3044911/69‑1 (MQ=255)
                                                            ggcggTATCTTTCCGACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:2413642/69‑1 (MQ=255)
                                                            ggcggTATCTTTCCGACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:2263285/69‑1 (MQ=255)
                                                            ggcggTATCTTTCCGACCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCg           <  1:433436/69‑1 (MQ=255)
                                                             gcggTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTc            <  1:2707101/67‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:2703020/69‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:1127855/69‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:1667805/69‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:1705432/69‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:2417711/69‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:246374/69‑1 (MQ=255)
                                                                   tCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCt    <  1:2615796/69‑1 (MQ=255)
                                                                     ttttCTCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCtt   >  1:932343/1‑68 (MQ=255)
                                                                     ttttCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCtt   >  1:1010416/1‑68 (MQ=255)
                                                                     ttttCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCTTc  >  1:154776/1‑69 (MQ=255)
                                                                     ttttCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCTTc  >  1:1333548/1‑69 (MQ=255)
                                                                        |                                                                 
ATGCCGCCCACGGTACGCCTGGCGTTTTAACGGCAGAACAGCCAGCCCTTCCAGCATGGCGGCGGTATCTTTTCGCCCGTGGGTTTCTACCACGCCAACCACAATATCCAGCCCTTGCGCCCGCAGTCGCTGGGCTTC  >  minE/447449‑447586

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: