Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,350,711 G→A 16.0% G389S (GGC→AGC)  PP_2945 → 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_0029473,350,7110GA16.0% 54.0 / 4.1 25G389S (GGC→AGC) PP_2945two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (11/10);  new base A (2/2);  total (13/12)
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.85e-01

AAGCCGACCTGGTAAAACGGTACCAGCCGTTACCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGA  >  NC_002947/3350562‑3350858
                                                                                                                                                     |                                                                                                                                                   
aaGCCGACCTGGTAAAACGGTACCAGCCGTTACCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAg                                                                                                                                                     <  2:344027/150‑1 (MQ=255)
             aaaaCGGTACCAGCCGTTACCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGtgtg                                                                                                                                        <  1:272112/150‑1 (MQ=255)
                    tACCAGCCGTTACCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGaa                                                                                                                                 >  1:103759/1‑150 (MQ=255)
                         gCCGTTACCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGc                                                                                                                            >  2:444328/1‑150 (MQ=255)
                               aCCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATAc                                                                                                                      >  1:44087/1‑150 (MQ=255)
                                  cccTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGgt                                                                                                                                           >  1:215802/1‑126 (MQ=255)
                                  cccTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGgt                                                                                                                                           <  2:215802/126‑1 (MQ=255)
                                                    gTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGga                                                                                                 <  1:444328/150‑1 (MQ=255)
                                                       gCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAAGCGCCt                                                                                                                                               <  1:309252/101‑1 (MQ=255)
                                                       gCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAAGCGCCt                                                                                                                                               <  1:309655/101‑1 (MQ=255)
                                                       gCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAAGCGCCt                                                                                                                                               >  2:309655/1‑101 (MQ=255)
                                                       gCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAAGCGCCt                                                                                                                                               >  2:309252/1‑101 (MQ=255)
                                                                                                accAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCg                                                     >  2:477129/2‑150 (MQ=255)
                                                                                                  cAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTc                                                   >  2:316897/1‑150 (MQ=255)
                                                                                                        aTGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCaa                                             <  2:304192/150‑1 (MQ=255)
                                                                                                         tGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCaa                                                          >  2:223322/1‑136 (MQ=255)
                                                                                                         tGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCaa                                                          <  1:223322/136‑1 (MQ=255)
                                                                                                          gggCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCtttttt                                                                   <  2:294470/126‑1 (MQ=255)
                                                                                                          gggCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCtttttt                                                                   >  1:294470/1‑126 (MQ=255)
                                                                                                          gggCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCAAgg                                           >  1:55861/1‑150 (MQ=255)
                                                                                                                            cATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGgaga                                                                                               <  2:227809/80‑1 (MQ=255)
                                                                                                                            cATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGgaga                                                                                               >  1:227809/1‑80 (MQ=255)
                                                                                                                            cATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTGTTATTATCATTTCCGTCGGCCAAGGTACCGGGCTAGGCCTGTc                         >  2:424403/1‑150 (MQ=255)
                                                                                                                                acaTTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCtt                     <  2:292184/150‑1 (MQ=255)
                                                                                                                                      ggCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTAt               <  1:316897/150‑1 (MQ=255)
                                                                                                                                                   aaGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGa  <  2:55861/150‑1 (MQ=255)
                                                                                                                                                     ggCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTcac                                                                >  1:60628/1‑86 (MQ=255)
                                                                                                                                                     ggCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTcac                                                                <  2:60628/86‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAGCCGACCTGGTAAAACGGTACCAGCCGTTACCCCCTGTCGAATGCCTGCCGTCGCAGATCAACCAAGTGATAATGAACCTGGTGCTCAACGCCGCCCAAGCCATGGGCCCCGAGCGGGGCACCATCACATTGGGCACTGCTCATCAAGGCGCCTGGGTGTGCATCGAAGTTGCAGATACAGGCGCGGGGATAGCACCGGAGAACATCAAGAAAATCTTCGATCCTTTTTTCACCACCAAGCCCGTCGGCCAAGGTACCGGGCTAGGCCTGTCCCTTTCCTATGGCATTGTAAAGA  >  NC_002947/3350562‑3350858

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