Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,920,571 C→G 10.7% P386R (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,920,5710CG10.7% 165.1 / 9.8 56P386R (CCC→CGC) PP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (19/31);  new base G (3/3);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ATCCTGGTCACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCT  >  NC_002947/5920423‑5920718
                                                                                                                                                    |                                                                                                                                                    
atcCTGGTCACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGcc                                                                                                                                                      >  5:297360/1‑149 (MQ=255)
      gTCACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGGGCAGCCCGCGCt                                                                                                                                                <  5:237600/149‑1 (MQ=255)
       tCACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTg                                                                                                                                               <  2:123611/149‑1 (MQ=255)
         aCCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGa                                                                                                                                             <  2:368994/149‑1 (MQ=255)
         aCCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTgcg                                                                                                                                             >  5:70098/1‑147 (MQ=255)
             gCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCTg                                                                                                                                         <  6:70098/149‑1 (MQ=255)
              ccccACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTgg                                                                                                                                        <  8:86175/149‑1 (MQ=255)
                ccACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTgggg                                                                                                                                      <  2:264582/149‑1 (MQ=255)
                 cACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGt                                                                                                                                     <  5:248847/149‑1 (MQ=255)
                     ggCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGc                                                                                                                                 <  8:8525/149‑1 (MQ=255)
                       cTCCGGCAAGACCAGCACGCTGTACTCCAGCCTCAAACTGCTGGCGACACGACAGGTCAAGCTGTGCACCATCGTAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTggct                                                                                                                               <  4:82769/149‑4 (MQ=255)
                          cGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAAcg                                                                                                                            >  2:361396/1‑149 (MQ=255)
                          cGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAAc                                                                                                                             <  1:351985/148‑1 (MQ=255)
                             cAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCt                                                                                                                                          >  7:297510/1‑132 (MQ=255)
                             cAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCt                                                                                                                                          <  8:297510/132‑1 (MQ=255)
                                gACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCa                                                                                                                     <  5:350132/149‑1 (MQ=255)
                                             tACACCAGCCTCAACCTGCTGGCCCCACCACAGGTCAACCTGTGCACCATCGAAGACCCGAGCGCGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTgcg                                                                                                        <  2:313436/149‑1 (MQ=255)
                                               cacCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTgcgc                                                                                                       >  7:144961/1‑148 (MQ=255)
                                               cacCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGcc                                                                                                      <  5:77227/149‑1 (MQ=255)
                                                 ccAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCaa                                                                                                    <  4:176360/149‑1 (MQ=255)
                                                      cTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGAccc                                                                                               >  6:85340/1‑149 (MQ=255)
                                                       tCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGACGCCCGCGCTGGCACTGGGGTTTGCCAACGCGGTCAGGGCCATGCTTCGCCAAGAcacc                                                                                              >  8:159940/1‑149 (MQ=255)
                                                       tCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCt                                                                                                                                          >  5:387646/1‑106 (MQ=255)
                                                       tCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCt                                                                                                                                          <  6:387646/106‑1 (MQ=255)
                                                                  gCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGAcatcatc                                                                                     <  4:215596/147‑1 (MQ=255)
                                                                    cacaccacaGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGa                                                                                 >  5:180557/1‑149 (MQ=255)
                                                                         cacaGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGc                                                                            >  2:277726/1‑149 (MQ=255)
                                                                            aGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGa                                                                          >  4:36026/1‑148 (MQ=255)
                                                                               tCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATc                                                                      >  2:383768/1‑149 (MQ=255)
                                                                                          aCCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGc                                                                            >  8:112678/1‑132 (MQ=255)
                                                                                          aCCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGc                                                                            <  7:112678/132‑1 (MQ=255)
                                                                                          aCCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATccgcgaccgcga                                                           <  1:277726/149‑1 (MQ=255)
                                                                                          aCCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATccgcgaccgcg                                                            >  6:109930/1‑148 (MQ=255)
                                                                                             aTCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAAc                                                        <  1:298003/149‑1 (MQ=255)
                                                                                             aTCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAAc                                                        <  4:120804/149‑1 (MQ=255)
                                                                                             aTCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAAc                                                        >  7:164090/1‑149 (MQ=255)
                                                                                                 aaGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGct                                                    >  4:157107/1‑149 (MQ=255)
                                                                                                      ccGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCtggt                                               >  5:182349/1‑149 (MQ=255)
                                                                                                      ccGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCAGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCtggt                                               >  7:331458/1‑149 (MQ=255)
                                                                                                                       gccagccTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAgcgc                              <  6:171264/149‑1 (MQ=255)
                                                                                                                           gccTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGAc                          <  3:36026/149‑1 (MQ=255)
                                                                                                                           gaattCAACCAGTTGCAGGTGCAGCCCTCGCTGTACCTGGGGTTTGCAAACGCGGTCAGGGCCATGCTGCGCCAAGACGCCGACATCATCATGATCGGCGAAATCCGCGACCTCGAAAGCGCTCTGGTGGCGGTGCAGGCAGCGCTGAc                          <  4:85214/146‑1 (MQ=255)
                                                                                                                             cgtCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTGTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCg                        <  8:104103/147‑1 (MQ=255)
                                                                                                                             cTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCg                        <  8:314321/149‑1 (MQ=255)
                                                                                                                                cAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCggg                      <  8:164090/148‑1 (MQ=255)
                                                                                                                                 aaCCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGc                     <  7:231228/148‑1 (MQ=255)
                                                                                                                                  aCCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGAcatca                                                                                       >  7:318810/1‑81 (MQ=255)
                                                                                                                                  aCCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGAcatca                                                                                       <  8:318810/81‑1 (MQ=255)
                                                                                                                                   ccAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCAcc                  >  1:243166/1‑149 (MQ=255)
                                                                                                                                    cAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCt                 >  5:252038/1‑149 (MQ=255)
                                                                                                                                        ttctggTGCTCCCCGCGCTGTCCCGTGGGGTTGCCAACCGGTGCATGGCCTTGTTTCGCCAATTCCCCTTAATCTTCATTTTCTGCGAAATCCGCGACCGCGTAACCGCTCTGGTGTCTGTGCAGGCAGCGCTGACCGGGCACCTGGt              <  2:349324/144‑1 (MQ=255)
                                                                                                                                          cAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCt           <  6:252038/149‑1 (MQ=255)
                                                                                                                                           aGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTa          >  3:247144/1‑149 (MQ=255)
                                                                                                                                             gTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATc        <  6:182349/149‑1 (MQ=255)
                                                                                                                                             ctGCGGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGGCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATCATCGGCGAAATCCGCGACCGCGAAAc                                                        <  2:176615/100‑1 (MQ=255)
                                                                                                                                               gCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCg                                                                           <  7:166765/80‑1 (MQ=255)
                                                                                                                                               gCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCg                                                                           >  8:166765/1‑80 (MQ=255)
                                                                                                                                                   cccGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCt  >  6:157444/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCCTGGTCACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCT  >  NC_002947/5920423‑5920718

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