Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,919,686 C→G 8.6% P91R (CCC→CGC) ‡ PP_5190 → type II secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,919,6860CG 8.6% 162.6 / 4.2 53P91R (CCC→CGC) ‡PP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base G (2/2);  total (26/27)
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.72e-01

ACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTC  >  NC_002947/5919542‑5919833
                                                                                                                                                |                                                                                                                                                   
aCCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAg                                                                                                                                                 <  3:214687/149‑1 (MQ=255)
  ccGCTGGAGCACATCGCCGCCAGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGACCACCTTTGCCAGTGGCTCCCCGACAAGGTCGGCCAACCCTCCCTCCGCCTCGAACCCAAGCAACACGACCTGCCCCAGGt                                                                                                                                               >  3:11258/1‑149 (MQ=255)
         aGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGc                                                                                                                                        <  2:925/149‑1 (MQ=255)
          gCACATCGCGGGCTGCAGCCTGGAAAGCCGGCAGCACACCGGCCAGCCCCTGGACTTGGAGCACCTTTGGCAGTGGCTGGCCGACAAGGTCGGCCAACCCTTCCTGCGCATCGAGCACATTCAACTCGACCTGCTCCAGTTCACCAGc                                                                                                                                        <  8:213433/148‑1 (MQ=255)
               tCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGAtc                                                                                                                                  <  7:155807/149‑1 (MQ=255)
                   cgccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTccc                                                                                                                              >  5:158348/1‑149 (MQ=255)
                     ccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCa                                                                                                                            <  6:158348/149‑1 (MQ=255)
                      cTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAg                                                                                                                           >  8:28523/1‑149 (MQ=255)
                        gCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                           <  2:122631/131‑4 (MQ=255)
                        gCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                           >  1:122631/1‑128 (MQ=255)
                           gCCTGGAACGCCGCCAGCACCCCGGCCTGCCCCTTGTGTTTGAGCACCTTTGCCAGTGGCTGGCGGACAAGGGCTGCCAACCCTACCTGCGCCTCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTc                                                                                                                      <  6:62299/149‑1 (MQ=255)
                                   aGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAgcg                                                                                                              <  4:226573/149‑1 (MQ=255)
                                    gccgccAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                           <  7:115281/119‑4 (MQ=255)
                                    gccgccAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGGGGCAGGTcgag                                                                                                                                           >  8:115281/1‑116 (MQ=255)
                                      cgccAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCa                                                                                                           >  8:216287/1‑149 (MQ=255)
                                             aCCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATc                                                                                                    >  5:74475/1‑149 (MQ=255)
                                             aCCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATc                                                                                                    <  5:273182/149‑1 (MQ=255)
                                              ccccGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCTTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATcc                                                                                                   >  1:91428/1‑149 (MQ=255)
                                                 cGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCa                                                                                                <  4:181031/149‑1 (MQ=255)
                                                           cTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCAAGCGCCACGGCATCCTCATCGTCGCGGc                                                                                      >  4:187785/1‑149 (MQ=255)
                                                               aCTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGAt                                                                                  <  7:28523/149‑1 (MQ=255)
                                                                 ttGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATcc                                                                                                   >  3:182559/1‑130 (MQ=255)
                                                                 ttGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATcc                                                                                                   <  4:182559/130‑1 (MQ=255)
                                                                    gAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCTTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCt                                                                              <  2:91428/148‑1 (MQ=255)
                                                                     aGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTcc                                                                            <  7:216287/149‑1 (MQ=255)
                                                                      gCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTc                                                                             >  1:115524/1‑147 (MQ=255)
                                                                                 aGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTc                                                                <  2:115524/149‑1 (MQ=255)
                                                                                        ggcCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgcg                                                         >  7:93428/1‑149 (MQ=255)
                                                                                         gcCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgcg                                                         <  6:74475/148‑1 (MQ=255)
                                                                                           cGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGcc                                                       >  7:124873/1‑148 (MQ=255)
                                                                                                     ggCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCa                                            <  3:163992/149‑1 (MQ=255)
                                                                                                       ccAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGcc                                                                                                                   <  6:44653/76‑1 (MQ=255)
                                                                                                       ccAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGcc                                                                                                                   >  5:44653/1‑76 (MQ=255)
                                                                                                         aaCCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgc                                                          <  4:21524/131‑1 (MQ=255)
                                                                                                         aaCCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgc                                                          <  4:234186/131‑1 (MQ=255)
                                                                                                         aaCCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgc                                                          >  3:21524/1‑131 (MQ=255)
                                                                                                         aaCCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAgc                                                          >  3:234186/1‑131 (MQ=255)
                                                                                                           cccTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATc                                                                                                    >  4:221654/1‑87 (MQ=255)
                                                                                                           cccTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATc                                                                                                    <  3:221654/87‑1 (MQ=255)
                                                                                                             cTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACt                                    <  8:93428/149‑1 (MQ=255)
                                                                                                                   gcgcATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAg                                           >  7:192147/1‑136 (MQ=255)
                                                                                                                   gcgcATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAg                                           <  8:192147/136‑1 (MQ=255)
                                                                                                                   gcgcATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAgg                              >  6:257700/1‑149 (MQ=255)
                                                                                                                           agcccATTCAACTCGAGCTGGCCACGGTCACGAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGGGGCCGATGCCTCGAGCGTCACCGTCGCCAGCGCCCAGTCTTTCCAGGATGACTGGCAGGCCGACGTg                      <  4:11258/147‑1 (MQ=255)
                                                                                                                              ccATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGcc                   >  3:105960/1‑149 (MQ=255)
                                                                                                                               cATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGccc                  >  4:198983/1‑149 (MQ=255)
                                                                                                                                       tCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCg                 <  2:238048/142‑1 (MQ=255)
                                                                                                                                       tCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCg                 >  1:238048/1‑142 (MQ=255)
                                                                                                                                       tCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTc          >  8:10745/1‑149 (MQ=255)
                                                                                                                                       tCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTc          >  1:58709/1‑149 (MQ=255)
                                                                                                                                         gACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCgg        >  5:227494/1‑149 (MQ=255)
                                                                                                                                              gCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGcc                            >  2:199912/1‑124 (MQ=255)
                                                                                                                                              gCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGcc                            <  1:199912/124‑1 (MQ=255)
                                                                                                                                               ccccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGt   <  2:58709/148‑1 (MQ=255)
                                                                                                                                                cccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTc          <  2:127276/140‑1 (MQ=255)
                                                                                                                                                cccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTc          >  1:127276/1‑140 (MQ=255)
                                                                                                                                                cccAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTc  >  8:195413/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGACAAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTTTCCAGGATGACTGGCAGGCCGACCTGGCCCGCAGCCTCGGCCGGTC  >  NC_002947/5919542‑5919833

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