Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 722,797 C→T 20.9% T680T (ACG→ACA 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
*W3110S.gb722,7970CT20.9% 88.2 / 18.4 43T680T (ACG→ACAkdpDfused sensory histidine kinase in two‑component regulatory system with KdpE
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base T (6/3);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

ATCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAGGGC  >  W3110S.gb/722745‑722841
                                                    |                                            
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:1261258/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:2307459/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:1423960/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:1491127/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:2622473/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:276374/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:2781481/1‑71 (MQ=255)
aTCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTaaa                            >  1:1558025/1‑71 (MQ=255)
 tCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAAt                           <  1:3332505/71‑1 (MQ=255)
 tCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAAt                           <  1:1510463/71‑1 (MQ=255)
 tCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAAt                           <  1:929515/71‑1 (MQ=255)
 tCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAAt                           <  1:3284325/71‑1 (MQ=255)
 tCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAAt                           <  1:2055392/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:843784/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:39894/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:985863/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:3287708/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:2966084/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:2820867/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:867889/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:910363/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:2428868/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:2269618/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:1280377/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:2078278/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:1524205/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGa                    <  1:2367229/71‑1 (MQ=255)
        cTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCGAACAGCACCGTAAGCGGCGTGAGTAAATCATGCGa                    <  1:1521249/71‑1 (MQ=255)
               gATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCTTAAATCATGCGAAAGcgcc             >  1:927674/1‑71 (MQ=255)
               gATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAATCATGCg                     >  1:260616/1‑63 (MQ=255)
               gATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAATCATGCGAAAgcg               >  1:524662/1‑69 (MQ=255)
               gATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAATCATGCGAAAGcgcc             >  1:157485/1‑71 (MQ=255)
               gATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAATCATGCGAAAGcgcc             >  1:460168/1‑71 (MQ=255)
               gATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACTGTAAGCGGCGTGCGTAAATCATGCGAAAGcgcc             >  1:640252/1‑71 (MQ=255)
                         tAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAggg   >  1:1665018/1‑71 (MQ=255)
                         tAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAggg   >  1:907118/1‑71 (MQ=255)
                         tAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAggg   >  1:1952620/1‑71 (MQ=255)
                          aaGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAgcgc              >  1:2097599/1‑59 (MQ=255)
                          aaGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAGGGc  >  1:2040232/1‑71 (MQ=255)
                          aaGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAGGGc  >  1:2043800/1‑71 (MQ=255)
                          aaGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAGGGc  >  1:209412/1‑71 (MQ=255)
                          aaGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAGGGc  >  1:1119311/1‑71 (MQ=255)
                          aaGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGAGAAAGCGCCGCCAGCAGGGc  >  1:1974938/1‑71 (MQ=255)
                                                    |                                            
ATCCTTCGCTTGCCAGATCGAGCGTTAAGATTTCTGCCTGACCAAACAGCACCGTAAGCGGCGTGCGTAAATCATGCGAAAGCGCCGCCAGCAGGGC  >  W3110S.gb/722745‑722841

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: