Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,632 T→C 20.2% T1046T (ACT→ACC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,6320TC20.2% 125.5 / 25.3 64T1046T (ACT→ACCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (26/25);  new base C (4/9);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGAC  >  NC_002947/197491‑197766
                                                                                                                                             |                                                                                                                                      
cGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGt                                                                                                                                >  1:54230/1‑150 (MQ=25)
 gTGACCCTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACGCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTc                                                                                                                               <  2:225344/150‑1 (MQ=11)
      ccTGAGCAACCCTGCCCACACGCCGGAGACCGTGACCCTGTCCAACGGCAAAACGATTACCGTTGAAGCCGGCCAAACCCACGGCAGCGTCGACTTCCAGACCCCGGCGAATAACGTCTACAACAACGAGCCGGCTGTCAGCGTCACCAt                                                                                                                          >  2:296972/1‑150 (MQ=11)
       cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAGCAACGGTTCGACCGTCAGCGTCACCAtg                                                                                                                         <  2:514642/150‑2 (MQ=9)
            cAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTc                                                                                                                                     >  2:465777/1‑133 (MQ=255)
            cAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTc                                                                                                                                     <  1:465777/133‑1 (MQ=255)
            cAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGaa                                                                                                                    <  2:377545/150‑1 (MQ=14)
               tCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAAAGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                                                                 >  1:105116/1‑150 (MQ=12)
                  tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAATCCATTTTCTTTTGAGACGGCAAGAACCGAGGCCTCCTCGTTCTCTGGCCCTCTGTCACTACAGGCCGCAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCcac                                                                                                              <  2:668695/150‑1 (MQ=2)
                  tGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCcac                                                                                                              >  1:89131/1‑150 (MQ=39)
                    cccAGACACCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAg                                                                                                            <  2:211317/150‑1 (MQ=17)
                      cAGACGCCGGTAACCGTGACCCTGTCCAAGGGCCAAACCATTACCGTTGAAGCCGGCAAGAACCAGGGCAGCGTCGATTTACAGACACCGGCGAATGACGACTACAACGATGGTTCGACCGTCAGCCTGGCCATCGGGATCGCCACAGGc                                                                                                          >  1:287306/1‑150 (MQ=11)
                       aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                                                                         >  1:458332/1‑150 (MQ=17)
                       aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                                                                         <  1:742933/150‑1 (MQ=17)
                       aGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCg                                                                                                         <  1:742936/150‑1 (MQ=17)
                        gACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                                                                        <  2:361564/150‑1 (MQ=17)
                              ggTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAttt                                                                                                  >  1:711776/1‑150 (MQ=11)
                               gTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                                                                 <  1:323258/150‑1 (MQ=14)
                               gTAACCGTGACCCTGTCCAACGCCCAAACCATTACCGTTGAAGCCGGCAATATCCATGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                                                                 <  2:351145/150‑1 (MQ=11)
                                  aCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                               >  2:404760/1‑149 (MQ=255)
                                  aCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                               <  1:404760/149‑1 (MQ=255)
                                  aCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                                                                              <  2:557936/150‑1 (MQ=255)
                                  aCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                                                                              <  2:711776/150‑1 (MQ=255)
                                    cGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCa                                                                                            <  2:157901/150‑1 (MQ=33)
                                           cTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcc                                                                                     <  2:508012/150‑1 (MQ=12)
                                            tGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAccc                                                                                    <  2:100113/150‑1 (MQ=9)
                                            tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGGAATTTCGAGCAACTGAccc                                                                                    >  2:742675/1‑150 (MQ=17)
                                                cAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGaa                                                                                >  1:657035/1‑150 (MQ=32)
                                                      ccAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                               <  1:291351/129‑1 (MQ=255)
                                                      ccAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                                                               >  2:291351/1‑129 (MQ=255)
                                                       cAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACg                                                                         >  2:525732/1‑150 (MQ=35)
                                                          aCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAAACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTACTTTCGAGCAACTGACCCCGCACCCGACGCCg                                                                      >  2:116869/1‑150 (MQ=255)
                                                              ttACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTc                                                                  >  1:526785/1‑150 (MQ=255)
                                                                 ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGt                                                                                                                                <  1:620339/85‑1 (MQ=35)
                                                                 ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGt                                                                                                                                >  2:620339/1‑85 (MQ=35)
                                                                 ccGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                                                               >  1:493347/1‑150 (MQ=255)
                                                                    ttGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGc                                                                    >  1:576057/1‑142 (MQ=255)
                                                                    ttGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGc                                                                    <  2:576057/142‑1 (MQ=12)
                                                                        aGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGAt                                                        <  2:526785/150‑1 (MQ=255)
                                                                         gCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATc                                                       <  2:493347/150‑1 (MQ=255)
                                                                               aaaaCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGAc                                                 >  2:148179/4‑150 (MQ=255)
                                                                                    ccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGt                                            >  1:83496/1‑150 (MQ=40)
                                                                                    ccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGCAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGt                                            >  2:185167/1‑150 (MQ=255)
                                                                                         gCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGc                                                                    >  2:262640/1‑121 (MQ=9)
                                                                                         gCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGc                                                                    <  1:262640/121‑1 (MQ=16)
                                                                                          cAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGacac                                      <  2:657035/150‑1 (MQ=255)
                                                                                            gCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGAcacca                                    <  2:83496/150‑1 (MQ=40)
                                                                                                        aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCt                                                                   >  1:672481/1‑107 (MQ=255)
                                                                                                        aGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCt                                                                   <  2:672481/107‑1 (MQ=255)
                                                                                                               ggCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCt                                                                   <  2:169940/100‑1 (MQ=21)
                                                                                                               ggCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCt                                                                   >  1:169940/1‑100 (MQ=21)
                                                                                                                     tGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAtt                                                                                                   <  1:160638/62‑1 (MQ=37)
                                                                                                                     tGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAAtt                                                                                                   >  2:160638/1‑62 (MQ=37)
                                                                                                                       acttCTACAACTACGGTTCGACTGTCAGCGTCACCATCAAGAACGCCACAGGCGATAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCaa         <  1:742675/147‑1 (MQ=255)
                                                                                                                          tCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGa                                                         >  2:507590/1‑99 (MQ=21)
                                                                                                                          tCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGa                                                         <  1:507590/99‑1 (MQ=21)
                                                                                                                            tACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAAc                                                                                          <  2:324742/64‑1 (MQ=37)
                                                                                                                            tACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAAc                                                                                          >  1:324742/1‑64 (MQ=37)
                                                                                                                            tACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCAAACCCGACGCCGGCTCAGACCACGAACACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCGCGCCAAGCGTg    >  1:733380/1‑150 (MQ=255)
                                                                                                                             acaacaACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                                                                        <  1:637996/49‑1 (MQ=40)
                                                                                                                             acaacaACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                                                                        >  2:637996/1‑49 (MQ=40)
                                                                                                                              caacaaCGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGAc  <  1:650607/150‑1 (MQ=255)
                                                                                                                                 caacGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcccc                                                                                   >  1:323097/1‑66 (MQ=35)
                                                                                                                                 caacGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGAcccc                                                                                   <  2:323097/66‑1 (MQ=37)
                                                                                                                                             |                                                                                                                                      
CGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGAC  >  NC_002947/197491‑197766

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