Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,719,906 G→T 36.5% G48W (GGG→TGG)  PP_4173 → two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,719,9060GT36.5% 9.4 / 4.5 11G48W (GGG→TGG) PP_4173two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (2/2);  total (5/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

GCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGT  >  NC_002947/4719775‑4720015
                                                                                                                                   |                                                                                                             
gCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCa                                                                                                           <  1:400204/136‑1 (MQ=255)
gCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCa                                                                                                           >  2:363890/1‑136 (MQ=255)
      gCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCTGGCAACTGAt                                                                                                     >  2:118034/1‑136 (MQ=255)
      gCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCTGGCAACTGAt                                                                                                     >  2:478919/1‑136 (MQ=255)
                                     caGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCATACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAatat                                                                      <  2:191666/136‑1 (MQ=255)
                                                                 ttGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCTGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGc                                          <  1:118034/136‑1 (MQ=255)
                                                                 ttGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCTGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGc                                          <  1:478919/136‑1 (MQ=255)
                                                                 ttGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGc                                          >  1:51000/1‑136 (MQ=255)
                                                                                                         aCCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGt  <  1:363890/136‑1 (MQ=255)
                                                                                                                 cAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCAc    >  1:275816/1‑126 (MQ=255)
                                                                                                                 cAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCAc    <  2:275816/126‑1 (MQ=255)
                                                                                                                                   |                                                                                                             
GCTTGAGCTGGGACATTCATACCCGCACGCAGATCATCAGTCTCGGCCCCGCCCTGCTACTGACCTTGCTCTTGATCAGCTTCTTTACCTTCGTACGCATCCAGGACCTGCGCCAGGAACTCAACCACACCGGGCAACTGATCGCCAACCAGTTGGCACCGGCGTCGGAATATGGGGTGATCTCGGGCAACAACGAAGTGCTTGAAGGGCTGATGCGTGCCACCCTCAGCATCCCGCACGT  >  NC_002947/4719775‑4720015

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: