Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,356,411 G→C 11.7% G94A (GGC→GCC)  PP_1181 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,356,4110GC11.7% 88.2 / 3.6 32G94A (GGC→GCC) PP_1181two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base C (2/2);  total (15/17)
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.91e-01

ATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATC  >  NC_002947/1356266‑1356550
                                                                                                                                                 |                                                                                                                                           
atcGACAGCCAGCCCGACCTGGTCATCCTCGACCTGATGCGGCCCCGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCATCCAGGACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTTCAGGCCCTg                                                                                                                                         <  2:124485/150‑1 (MQ=255)
   gACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCGTCATGCTGACCACGTGAAGTGACGAACCGGACCAGCGCCAGGGCCTGGAc                                                                                                                                      >  2:178063/1‑150 (MQ=255)
      agccagccCGACCTGGTGATCCTCGACCTTATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTg                                                                                                                                   >  1:91032/1‑150 (MQ=255)
        ccagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTggg                                                                                                                                 >  1:457018/1‑150 (MQ=255)
         cagccCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGgc                                                                                                                                <  1:439175/150‑1 (MQ=255)
                               aCCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAAc                                                                                                          <  1:307699/150‑1 (MQ=255)
                                     tGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgc                                                                                                    <  1:178063/150‑1 (MQ=255)
                                      gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTgcg                                                                                                   >  1:193097/1‑150 (MQ=255)
                                      gCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGACGGGACTCCCGAACGTTTTTCATGCTGATCCTCATCCGCTCCGCGCGCAGTGATGAACTGGACCAGGTCCATGGCCTGGACCTGGGTGCCGACGACTACGTCTGCAAACCCGTgcg                                                                                                   >  1:134280/1‑150 (MQ=255)
                                          cccGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCa                                                                                               <  2:91032/150‑1 (MQ=255)
                                              gCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCt                                                                                           <  2:296449/150‑1 (MQ=255)
                                                  ggACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtgct                                                                                       >  2:51035/1‑150 (MQ=255)
                                                      ggCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGctgctg                                                                                   >  1:250930/1‑150 (MQ=255)
                                                        ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtgct                                                                                       <  2:67273/144‑1 (MQ=255)
                                                        ccTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTtgct                                                                                       >  1:67273/1‑144 (MQ=255)
                                                          tGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCAc                                                                               <  1:68219/150‑1 (MQ=255)
                                                                               gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGcgcc                                                          >  2:257612/1‑150 (MQ=255)
                                                                               gcAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGcgcc                                                          >  1:16901/1‑150 (MQ=255)
                                                                                  gCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCa                                                       <  2:134280/150‑1 (MQ=255)
                                                                                    cAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGc                                                     <  1:361109/150‑1 (MQ=255)
                                                                                       tACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGATCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGaa                                                  >  2:69604/1‑150 (MQ=255)
                                                                                        aCCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGaaa                                                 <  2:457018/150‑1 (MQ=255)
                                                                                          cAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCACCCACGCTTGCTGCTGGCACGTATCTATGCCCTGCTGCGCCGCAGCGAAAcc                                               <  1:257612/150‑1 (MQ=255)
                                                                                                  gATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAg                                       >  2:316994/1‑150 (MQ=255)
                                                                                                                         gTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCg                <  1:416532/150‑1 (MQ=255)
                                                                                                                          tgatgaACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCGCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCgg               >  2:9380/1‑150 (MQ=255)
                                                                                                                                  tGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCGCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGc       <  1:9380/150‑1 (MQ=255)
                                                                                                                                       cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc  <  1:69604/150‑1 (MQ=255)
                                                                                                                                       cAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc  <  1:224221/150‑1 (MQ=255)
                                                                                                                                       cAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGctgctg                                                                                   >  1:47783/1‑69 (MQ=255)
                                                                                                                                       cAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGctgctg                                                                                   <  2:47783/69‑1 (MQ=255)
                                                                                                                                       cAGGTCCAGGCCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATc  >  2:150406/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
ATCGACAGCCAGCCCGACCTGGTGATCCTCGACCTGATGCTCCCCGGCGAGGACGGCCTGAGCATCTGCCGCCGGGTGCGCAGCCAGTACCAGGGGCCGATCCTCATGCTCACCGCGCGCAGTGATGAACTGGACCAGGTCCAGGGCCTGGACCTGGGCGCCGACGACTACGTCTGCAAACCCGTGCGCCCACGCTTGCTGCTGGCACGTATCCAGGCCCTGCTGCGCCGCAGCGAAACCGTGGACAGCAAGCGCCAGGACCTGGCCTTCGGTGCGCTGCGTATC  >  NC_002947/1356266‑1356550

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