Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 195,539 G→A 13.3% V349I (GTC→ATC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947195,5390GA13.3% 89.9 / 8.7 45V349I (GTC→ATC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (17/22);  new base A (2/4);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA  >  NC_002947/195397‑195687
                                                                                                                                              |                                                                                                                                                    
ccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAc                                                                                                                                                  >  3:318482/1‑147 (MQ=34)
ccTGAGCAACCCTGCCCAGACGCAGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAAACCAGGGCAGCGTCGATTTCCAGACACCGGCCAATGACGTCTACAACAACCGTTCGACCGTCATCGTCACCa                                                                                                                                                >  1:180663/1‑149 (MQ=11)
  tGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCATCACCATc                                                                                                                                              >  7:216053/1‑149 (MQ=21)
    aGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATc                                                                                                                                              <  8:258531/147‑1 (MQ=34)
          ccTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                                                                                      <  2:98428/149‑1 (MQ=32)
                         gtgaccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGAGAAtt                                                                                                                        >  5:117657/1‑148 (MQ=17)
                             ccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                                                                                   <  7:97899/149‑1 (MQ=35)
                             ccgtgaccCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACAGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                                                                                   <  6:122772/149‑1 (MQ=25)
                                   cccTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAAACCATGGCAGCGTCGATTTCCAGACACAGGCCAATGACGACTACAACAACGGTTAGACCGTCAGCGTCACAATCTAGAACGCAAACGGCGGGTATTTCGAGCAAcg                                                                                                              >  8:262953/1‑147 (MQ=11)
                                        tCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAcccc                                                                                                        <  7:31965/149‑1 (MQ=39)
                                                 gaGGTCTTCACTTTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACGCCGGCCAATGACTTCTACAAGAACGGTTCGACCGTCAGCATGACCATCGAGAACTCCACCGGCGGTAATTTCGAGCAACTGACCGCGAATCCGAc                                                                                               <  5:10322/148‑1 (MQ=14)
                                                         caCTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCg                                                                                                                     <  3:152626/119‑1 (MQ=35)
                                                         caCTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCg                                                                                                                     >  4:152626/1‑119 (MQ=35)
                                                            tGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                                    <  3:42783/149‑1 (MQ=38)
                                                              ttGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGa                                                                                    <  6:71584/147‑1 (MQ=37)
                                                                      ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCTGCCAATGACGTCTACAACAACTGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCaa                                                                          <  7:262953/149‑1 (MQ=255)
                                                                            aCCCAGGGCAGCGTCGAGTTCCAGACCCCGGACAATGACGTCTACAACAAAGGTTCGACCGTCAGCGTCACCATCGAGAACGCCCCCGGCGGTAATTTCGAGAACCTGCCCCCGAATACGACCCCGGCTCAGACCACGATCAACCACTc                                                                    >  1:248537/1‑149 (MQ=11)
                                                                              ccAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCgg                                                                  <  4:318482/149‑1 (MQ=34)
                                                                              ccAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCgg                                                                  <  7:130/149‑1 (MQ=34)
                                                                                aGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAGCAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                                                                <  2:35156/149‑1 (MQ=21)
                                                                                aGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTc                                                                <  4:21260/149‑1 (MQ=34)
                                                                                  ggCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGCCCCCGGCTCCGACCACGATCAACGACTCGGTCGa                                                              >  6:24726/1‑149 (MQ=21)
                                                                                   gccgcggcgATTTCCAGACCCCGGCCTATGACGTCTACAACAACGCTGCGACCTTCAGCATCCCCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGac                                                             <  8:327196/146‑1 (MQ=11)
                                                                                   gCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGac                                                             <  4:150676/149‑1 (MQ=35)
                                                                                             attCCAGCCCCCGGCCAATGACGTCTACAACAACGCATCGACCGTCAGCGTCACCATCCAGAACGCCACCGGCGTTAATTTCGAGCACCTGACCCCCAATCCCCCCCCGGCTCAGACCACGATCAACGACTCTGTCGACACCACCACcg                                                   <  2:180663/148‑1 (MQ=14)
                                                                                                        cGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACg                                        >  8:108173/1‑149 (MQ=39)
                                                                                                              aTGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGTCGGAAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACGACCACCGCGACCCTGACGGCGAGc                                  <  6:117657/149‑1 (MQ=255)
                                                                                                                        caacaaCGGTTCGACCGTCAGCATCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGt                          <  8:216053/147‑1 (MQ=37)
                                                                                                                         aacaacGGTTCAACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAc                       <  5:185992/149‑1 (MQ=12)
                                                                                                                         aacaacGGTGCGACCGGCAGCGGCACCATCGAGAACGCCACAGTGGGTAATTTCGAGCAAATGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCAc                       <  2:3809/149‑1 (MQ=11)
                                                                                                                           caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGAcc                                              <  5:104542/124‑1 (MQ=38)
                                                                                                                           caacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGAcc                                              >  6:104542/1‑124 (MQ=38)
                                                                                                                            aacGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                     >  2:96222/1‑148 (MQ=17)
                                                                                                                            aacGGTTCGACCGTCAGCATCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCg                     >  7:173724/1‑148 (MQ=37)
                                                                                                                             acGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGaa                   <  1:226939/149‑1 (MQ=17)
                                                                                                                                  tCGACCGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAggcgg              >  5:325954/1‑149 (MQ=17)
                                                                                                                                     aCCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt           >  2:351891/1‑149 (MQ=37)
                                                                                                                                     aCCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGt           >  8:50763/1‑149 (MQ=37)
                                                                                                                                      ccGTCAGCGTCACCATCGAGAACGCCACAGGCGTTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGAGTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTg          <  3:110012/149‑1 (MQ=255)
                                                                                                                                      ccGTCAGCATCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTg          <  8:173724/149‑1 (MQ=38)
                                                                                                                                       cGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGa         >  6:44850/1‑149 (MQ=17)
                                                                                                                                        gTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGa         >  4:212254/1‑148 (MQ=17)
                                                                                                                                         tCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTCAt        >  6:112787/1‑148 (MQ=12)
                                                                                                                                             cGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACtt   <  3:139177/149‑1 (MQ=17)
                                                                                                                                             cGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACt    >  4:139177/1‑148 (MQ=17)
                                                                                                                                              gTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGCCTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTa  >  8:329241/1‑149 (MQ=255)
                                                                                                                                              gTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAATGCGGCGTGATCACTTa  <  7:339586/149‑1 (MQ=12)
                                                                                                                                              |                                                                                                                                                    
CCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTA  >  NC_002947/195397‑195687

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