Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,5590CA31.5% 89.3 / 38.6 54T355T (ACC→ACA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (14/23);  new base A (10/7);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAG  >  NC_002947/195411‑195702
                                                                                                                                                    |                                                                                                                                               
cccAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCCTCACCATCGAGAACGCCACcg                                                                                                                                                >  1:365038/1‑150 (MQ=17)
           gtgaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCCTCACCATCGAGAACGCCACCGGCGGTAATTTc                                                                                                                                     <  2:365038/150‑1 (MQ=38)
             gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                                                                                                   <  2:269120/150‑1 (MQ=37)
             gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                                                                                                   <  2:141562/150‑1 (MQ=37)
             gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                                                                                                   <  2:664059/150‑1 (MQ=37)
             gaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGACCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                                                                                                   <  1:646592/150‑1 (MQ=37)
               ccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGTGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGc                                                                                                                                 <  1:537789/150‑1 (MQ=37)
                cgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAGCGCCACCGGCGGCAATTTCGAGCa                                                                                                                                >  2:1106868/1‑150 (MQ=17)
                 gtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCaa                                                                                                                               >  2:213001/1‑150 (MQ=38)
                 gtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCACACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCaa                                                                                                                               <  1:465872/150‑1 (MQ=255)
                  tgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTACAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAAc                                                                                                                              <  2:418607/150‑1 (MQ=255)
                      ccTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCg                                                                                                                                    <  2:5325/140‑1 (MQ=37)
                      ccTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCg                                                                                                                                    >  1:5325/1‑140 (MQ=37)
                      ccTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGCAATTTCGAGCAACTGAc                                                                                                                          >  2:359904/1‑150 (MQ=21)
                      ccTGTCCAACGGCCAGGTCATCACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGCCGGTAATTTCGAGCAACTGAc                                                                                                                          >  2:286865/1‑150 (MQ=255)
                                 gCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGa                                                                                                               <  1:228144/150‑1 (MQ=255)
                                      gTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCg                                                                                                          <  2:267356/150‑1 (MQ=255)
                                           caCTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCa                                                                                                     <  2:293975/150‑1 (MQ=40)
                                                        ggCAAAACCCAGGGCAGCGTCGATTTCCAGCCCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCCACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAAc                                                                                        <  1:1403789/150‑1 (MQ=255)
                                                              aCCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCg                                                                                  <  1:213001/150‑1 (MQ=35)
                                                               cccAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCgg                                                                                 >  1:339242/1‑150 (MQ=35)
                                                                ccAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                <  1:286865/150‑1 (MQ=34)
                                                                ccAGGGCAGCGTCGATTTCCAGACCCCGGCCAATCACGTCTTCCTCGCCTGCCCGTGCGGCAGCTTCACCATAAAGCGCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGt                                                                                <  1:1435769/150‑1 (MQ=11)
                                                                 cAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACTGTTCGACCGTCAGCGTCACCATAGAGAACGCAACCGGCGGTCCATACGAGCAATTGACTCCGAGTCCGACCCCGGCTCAGACCTCGATCAACGACTCTGTc                                                                               <  2:754508/150‑1 (MQ=11)
                                                                      cAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGacac                                                                          <  1:101332/150‑1 (MQ=38)
                                                                      cAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGTCCCCGGCTCAGACCACGATCAACGACTCGGTCGacac                                                                          <  1:1196776/150‑1 (MQ=14)
                                                                           tCGATTTCCAGACCCCGGCCCATGACGTCTACAACAACGGTTTTACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGAtaccacca                                                                     <  1:1348966/150‑1 (MQ=255)
                                                                                    aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGaca                                                                           <  2:151420/135‑1 (MQ=38)
                                                                                    aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGaca                                                                           >  1:151420/1‑135 (MQ=37)
                                                                                    aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGa                                                               <  1:1364264/147‑1 (MQ=40)
                                                                                    aGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGa                                                               >  2:1364264/1‑147 (MQ=40)
                                                                                             cgaaTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGa                                                   <  1:1154111/148‑1 (MQ=2)
                                                                                                    tgTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACTCCGAATCCGACCCCGGCGCAGTCCACGATCAACGTCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTc                                            <  1:1321373/149‑1 (MQ=11)
                                                                                                         acaacaACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                       <  2:439238/150‑1 (MQ=38)
                                                                                                         acaacaACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCa                                       >  1:316752/1‑150 (MQ=18)
                                                                                                           aacaacTGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGTCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                     <  1:1130450/150‑1 (MQ=255)
                                                                                                           aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                     >  1:531475/1‑150 (MQ=14)
                                                                                                           aacaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAcc                                     >  1:377455/1‑150 (MQ=17)
                                                                                                            acaacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                                    >  1:111601/1‑150 (MQ=12)
                                                                                                                  gTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                              >  1:102458/1‑150 (MQ=12)
                                                                                                                  gTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcg                              <  1:947101/150‑1 (MQ=12)
                                                                                                                   ttCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcgg                             >  1:1469951/1‑150 (MQ=39)
                                                                                                                   ttCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcgg                             <  2:1297288/150‑1 (MQ=12)
                                                                                                                        ccGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGa                        >  1:553351/1‑150 (MQ=39)
                                                                                                                              gCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAAt                                                                                                                   >  2:621993/1‑53 (MQ=33)
                                                                                                                              gCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAAt                                                                                                                   <  1:621993/53‑1 (MQ=33)
                                                                                                                                 tCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTaca               >  2:839792/1‑150 (MQ=12)
                                                                                                                                  cACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTacac              >  1:1384606/1‑150 (MQ=255)
                                                                                                                                    ccATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGGCGGTCACCGAAGGCGGCGTGATCACTTACACCg            >  2:809838/1‑150 (MQ=255)
                                                                                                                                     cATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTCCACCGt           >  2:652089/1‑150 (MQ=255)
                                                                                                                                     cATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGt           >  1:1026234/1‑150 (MQ=18)
                                                                                                                                       tCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGa         >  1:1455556/1‑150 (MQ=12)
                                                                                                                                            aaCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTg    >  2:1092503/1‑150 (MQ=12)
                                                                                                                                             aCGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGa   <  2:1034375/150‑1 (MQ=12)
                                                                                                                                              cGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAg  >  2:300024/1‑150 (MQ=12)
                                                                                                                                                 cacaGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGa                                                               >  1:907761/1‑86 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAG  >  NC_002947/195411‑195702

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