Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,590 T→C 14.5% D1032D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,5900TC14.5% 108.2 / 9.9 48D1032D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base C (2/5);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCA  >  NC_002947/197442‑197732
                                                                                                                                                    |                                                                                                                                              
cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAtt                                                                                                                                               <  2:661489/150‑1 (MQ=39)
 ccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAttt                                                                                                                                              <  1:551858/150‑1 (MQ=39)
   tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACGTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACAGCCAAAACATTACCGTTGAATCCGGCAAGACCCAGGGAAACGTCGATTTca                                                                                                                                            >  2:263638/1‑149 (MQ=11)
   tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCATGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAACCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTcc                                                                                                                                            <  2:620624/150‑1 (MQ=255)
        gCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAcc                                                                                                                                       <  2:25227/150‑1 (MQ=39)
         cGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                                      >  1:23619/1‑150 (MQ=39)
             cccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                                  <  1:424801/150‑1 (MQ=39)
                gTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGCGTAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGaa                                                                                                                               <  1:86252/150‑1 (MQ=255)
                 tCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCTTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                              <  1:498719/150‑1 (MQ=38)
                    gTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                                                                                                           <  2:702692/150‑1 (MQ=39)
                     tCACCGAAGGCGGTGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGTTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACg                                                                                                                          >  1:496609/1‑150 (MQ=255)
                        ccGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGATCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCt                                                                                                                       <  1:659231/150‑1 (MQ=255)
                         cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTa                                                                                                                      >  2:794521/1‑150 (MQ=39)
                               cggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                                <  1:703451/150‑1 (MQ=39)
                                   gTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTa                                                                                                                      >  1:200286/1‑140 (MQ=40)
                                   gTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTa                                                                                                                      <  2:200286/140‑1 (MQ=40)
                                    ttatcatcTCCAACATGCTCCTAAGAAATCCTGACCAACAGCCGCTAACCCTGACCCTGCACAAAGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGCTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                                           <  2:365925/148‑1 (MQ=11)
                                    tGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                                           <  2:641997/150‑1 (MQ=38)
                                     gATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                                           <  2:646687/117‑1 (MQ=255)
                                     gATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                                           >  1:646687/1‑117 (MQ=255)
                                       tCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                        >  1:173806/1‑150 (MQ=38)
                                         aCCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                      >  1:287871/1‑150 (MQ=38)
                                         aCCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACAGCGACACCGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGAcc                                                                                                      >  1:614358/1‑149 (MQ=255)
                                          ccTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                                           <  2:263095/112‑1 (MQ=255)
                                          ccTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                                           >  1:263095/1‑112 (MQ=255)
                                               acCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCNGCGAATGACGTCTACAACAACGGTTCGACCGTCAGc                                                                                                >  2:647826/1‑150 (MQ=25)
                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCg                                                                                      >  1:242684/1‑150 (MQ=18)
                                                          gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                        >  1:589032/1‑131 (MQ=35)
                                                          gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCga                                                                                     >  2:277188/1‑150 (MQ=38)
                                                          gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                     <  2:260554/150‑1 (MQ=17)
                                                          gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTGCGACTGTCAGCGTCACCATCga                                                                                     <  2:669628/150‑1 (MQ=38)
                                                          gAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                        <  2:589032/131‑1 (MQ=18)
                                                           aGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCgag                                                                                    <  2:310421/150‑1 (MQ=39)
                                                              aaTCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAAc                                                                                 <  1:794521/150‑1 (MQ=17)
                                                                       cAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGc                                                                        <  1:229821/150‑1 (MQ=11)
                                                                             ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                                       <  2:199142/129‑1 (MQ=40)
                                                                             ccGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                                       >  1:199142/1‑129 (MQ=40)
                                                                               ggTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                     <  1:153796/129‑1 (MQ=17)
                                                                               ggTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                     >  2:153796/1‑129 (MQ=18)
                                                                                gtgaCCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTc                                                               <  2:365480/147‑1 (MQ=255)
                                                                                          cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGCCTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGa                                                     >  1:617331/1‑150 (MQ=17)
                                                                                                      gCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGa                                         >  1:778161/1‑150 (MQ=255)
                                                                                                        cAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                             >  2:550813/1‑128 (MQ=255)
                                                                                                        cAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                             <  1:550813/128‑1 (MQ=255)
                                                                                                                tACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGTGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTcg                               >  2:313024/1‑149 (MQ=37)
                                                                                                                tACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCa                               <  1:591679/150‑1 (MQ=255)
                                                                                                                          gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATCTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATc                     >  2:575864/1‑150 (MQ=38)
                                                                                                                              gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCg                 >  1:387270/1‑150 (MQ=40)
                                                                                                                               cAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGa                >  1:298412/1‑150 (MQ=255)
                                                                                                                                             gCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGAcacca  <  2:298412/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCA  >  NC_002947/197442‑197732

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: