Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,953,465 G→T 19.0% P56P (CCC→CCA PP_3484 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,953,4650GT19.0% 55.3 / 2.4 21P56P (CCC→CCAPP_3484two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (8/9);  new base T (2/2);  total (10/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 = 6.99e-01

CTTGAACAGGACGTGACTCACTCCGAGTTCGGCGGCCCGCTCACGCACCTCGACGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGCCTGCAGATACGCTTGCAGGTTTTGCGCAAG  >  NC_002947/3953330‑3953596
                                                                                                                                       |                                                                                                                                   
cTTGAACAGGACGTGACTCACTCCGAGTTCGGCGGCCCGCTCACGCACCTCGACGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCg                                                                                                                                     <  1:21045/136‑1 (MQ=255)
                                               ccTCGACGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGc                                                                                      <  1:151290/136‑1 (MQ=255)
                                               ccTCGACGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGc                                                                                      <  1:324940/136‑1 (MQ=255)
                                               ccTCGACGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGc                                                                                      <  1:393754/136‑1 (MQ=255)
                                                     cGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTc                                                                                >  1:139153/1‑136 (MQ=255)
                                                             tgggCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCAc                                                                        >  2:48924/1‑136 (MQ=255)
                                                             tgggCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCAc                                                                        >  1:444580/1‑136 (MQ=255)
                                                             tgggCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGCCTTTCAGCCTGTTCGATACCAc                                                                        >  2:459067/1‑136 (MQ=255)
                                                               ggCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTg                                                                      <  2:444580/136‑1 (MQ=255)
                                                                       tCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACcgtc                                                              >  1:52692/1‑136 (MQ=255)
                                                                         aGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACcgtcgt                                                            <  1:48924/136‑1 (MQ=255)
                                                                                   ggCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATg                                                  >  1:173664/1‑136 (MQ=255)
                                                                                   ggCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATg                                                  >  2:257717/1‑136 (MQ=255)
                                                                                          ccTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTcc                                           <  1:148404/136‑1 (MQ=255)
                                                                                             gTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCATGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAgg                                        >  2:326554/1‑136 (MQ=255)
                                                                                                  tGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTg                                   <  1:95109/136‑1 (MQ=255)
                                                                                                      tACAGCCTCCAGAACCTGAAAACCTTCCATGTCTGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGcc                               >  2:404198/1‑136 (MQ=255)
                                                                                                      tACAGCCTCCAGAACCTGAAAACCTTCCATGTCTGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGcc                               >  2:502393/1‑136 (MQ=255)
                                                                                                                 caaCCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGCCTGCAGATACGc                    <  2:524544/135‑1 (MQ=255)
                                                                                                                  aaCCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGCCTGCAGATACGCt                   <  2:480190/136‑1 (MQ=255)
                                                                                                                                   tGTCTGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGCCTGCAGATACGCTTGCAGGTTTTGCGCAag  <  1:502393/136‑1 (MQ=255)
                                                                                                                                   tGTCTGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGCCTGCAGATACGCTTGCAGGTTTTGCGCAag  <  1:404198/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CTTGAACAGGACGTGACTCACTCCGAGTTCGGCGGCCCGCTCACGCACCTCGACGGTGGGGTGGGCGGTGATCAGCACGAAGTGGCATTGCCTGTTCTTGCGTACAGCCTCCAGAACCTGAAAACCTTCCATGTCGGGCAGGCGGTAATCCAGCACAACTACTTGCGGTGCCAGGCTTTCAGCCTGTTCGATACCACTGGCACCGTCGTGGGCAACATGAACTTCCAGGCCTTGCGCCTGCAGATACGCTTGCAGGTTTTGCGCAAG  >  NC_002947/3953330‑3953596

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