Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,889,983 C→A 23.6% G111V (GGC→GTC)  PP_1695 ← sodium‑solute symporter/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_0029471,889,9830CA23.6% 32.2 / 6.5 17G111V (GGC→GTC) PP_1695sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base A (2/2);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TGCCGGTGGCGTCGGCGTTGGCGCCGATCAGCAGGTTCACGCCCAGTACGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACG  >  NC_002947/1889849‑1890112
                                                                                                                                      |                                                                                                                                 
tGCCGGTGGCGTCGGCGTTGGCGCCGATCAGCAGGTTCACGCCCAGTACGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGcc                                                                                                                                  >  1:75430/1‑136 (MQ=255)
                        cGATCAGCAGGTTCACGCCCAGTACGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGc                                                                                                          <  1:136394/136‑1 (MQ=255)
                        cGATCAGCAGGTTCACGCCCAGTACGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGc                                                                                                          <  1:212904/136‑1 (MQ=255)
                               cagGTTCACGCCCAGTACGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAg                                                                                                   <  1:80712/136‑1 (MQ=255)
                                               aCGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGACGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTgctgct                                                                                   >  2:61925/1‑136 (MQ=255)
                                               aCGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGACGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTgctgct                                                                                   >  2:60633/1‑136 (MQ=255)
                                                      ccTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGAt                                                                            >  1:19426/1‑136 (MQ=255)
                                                        ttGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCa                                                                          <  1:894/136‑1 (MQ=255)
                                                          gAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGc                                                                        <  2:48137/136‑1 (MQ=255)
                                                           aGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGcc                                                                       >  1:157406/1‑135 (MQ=255)
                                                           aGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCa                                                                       >  1:212556/1‑136 (MQ=255)
                                                                   ggCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCGGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCAtcttc                                                               >  1:173941/1‑136 (MQ=255)
                                                                                     acCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCg                                             >  2:130885/1‑136 (MQ=255)
                                                                                                                cacTGCCAGGGTCTGCGACTTGACGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAgg                  <  1:60633/136‑1 (MQ=255)
                                                                                                                cacTGCCAGGGTCTGCGACTTGACGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAgg                  <  1:61925/136‑1 (MQ=255)
                                                                                                                      cAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATg            >  2:153512/1‑136 (MQ=255)
                                                                                                                          gTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCa        <  2:5688/136‑1 (MQ=255)
                                                                                                                                gACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACg  >  1:61540/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                 
TGCCGGTGGCGTCGGCGTTGGCGCCGATCAGCAGGTTCACGCCCAGTACGATGCCCTTGAGCTGCAAGGCGATGTACGGCAGCACACCGACCAGGCAGATCAGTGCCACCACCACTGCCAGGGTCTGCGACTTGCCGTAGCGCGCGGCAATGAAGTCGGCGATCGAGGTGATGTTCTGCTGCTTGCTGATCAGCACCATCTTCTGCAGCACCCAAGGGGCGAAGATCAGCAGCAGTACCGGGCCCAGGTAGATGGGCAGGAACG  >  NC_002947/1889849‑1890112

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: