Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,920,579 G→C 9.1% D389H (GAC→CAC)  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,5790GC 9.1% 154.8 / 7.3 55D389H (GAC→CAC) PP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base G (20/28);  new base C (2/3);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

ACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCTGCATACCAA  >  NC_002947/5920432‑5920727
                                                                                                                                                   |                                                                                                                                                     
aCCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGa                                                                                                                                                      <  2:368994/149‑1 (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                                                                                                                                                   <  8:297510/132‑1 (MQ=255)
                    cAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCt                                                                                                                                                   >  7:297510/1‑132 (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                                                                                                                                                   <  6:387646/106‑1 (MQ=255)
                                              tCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTCCAGCGCGGGCTGCACCt                                                                                                                                                   >  5:387646/1‑106 (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)
                                                                                                                                                cTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCTGCATAc     <  8:144961/149‑1 (MQ=255)
                                                                                                                                                   gACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCTGCATAcc    <  1:383768/147‑1 (MQ=255)
                                                                                                                                                   gACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCTGCATAcc    <  5:85340/147‑1 (MQ=255)
                                                                                                                                                   gACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCTGCATACCaa  <  1:361396/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCGGCCCCACCGGCTCCGGCAAGACCAGCACGCTGTACACCAGCCTCAAACTGCTGGCCACACCACAGGTCAACCTGTGCACCATCGAAGACCCGATCGAGCGCCTGGAGCCAGCCTTCAACCAGTTGCAGGTGCAGCCCGCGCTGGACCTGGGGTTTGCCAACGCGGTCAGGGCCATGCTGCGCCAAGACCCCGACATCATCATGATCGGCGAAATCCGCGACCGCGAAACCGCTCTGGTGGCGGTGCAGGCAGCGCTGACCGGGCACCTGGTGCTATCGACCCTGCATACCAA  >  NC_002947/5920432‑5920727

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: