Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,579 G→A 14.0% V695V (GTG→GTA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,5790GA14.0% 84.3 / 7.6 43V695V (GTG→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (3/3);  total (20/23)
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.34e-01

CTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACG  >  NC_002947/196431‑196724
                                                                                                                                                    |                                                                                                                                                 
cagTCAGTGTCACCAACGAGAACGCCACAGGCGGTCATTTCGAGCAACTGACCTCGACTCCGACTACGGCTCAGGCCACGATCAACGACTCCACAGGCGCCACCACCCTTTCCCTGCCGGCGGGTCTTTCGGTCACCGAAGGCGGCGTGa                                                                                                                                                  <  2:1692211/148‑1 (MQ=2)
 tGTCAGTGTCACCATCGAGAACGGCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                                                                 >  2:1187583/1‑150 (MQ=255)
 tGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                                                                 <  2:2138808/150‑1 (MQ=255)
 tGTCAGTGTCACCATCGAGAACGACACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGGCCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGAt                                                                                                                                                 >  2:1187993/1‑150 (MQ=255)
    cAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                                                              <  1:2198130/150‑1 (MQ=39)
     aGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAcc                                                                                                                                             >  1:1460182/1‑150 (MQ=255)
      gtgtCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCt                                                                                                                                            >  1:1286079/1‑150 (MQ=39)
         ccACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTaca                                                                                                                                         <  2:1348995/149‑1 (MQ=39)
           cccACCGAGATCGCCACCGGCGGTAATTCCGAGCATCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTCCAcc                                                                                                                                       <  2:601184/149‑1 (MQ=11)
                    aaCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTg                                                                                                                              >  1:1137204/1‑150 (MQ=38)
                     aCGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGa                                                                                                                             >  1:317985/1‑150 (MQ=39)
                           ccGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                                                                                                                       <  2:317985/150‑1 (MQ=37)
                           ccGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                                                                                                                       >  1:375507/1‑150 (MQ=37)
                           ccGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                                                                                                                       <  1:594155/150‑1 (MQ=38)
                              gcggTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGCGGACGAGCAACGACTCGGGCGATGCCGCCGCCGCGTGTGTAGCTCCGAGTCCTTGCCGCACCATTGGCGGCGAGATCACCTACACCGTGACCCTGAGCAATCCTg                                                                                                                    >  2:72310/1‑150 (MQ=2)
                                ggTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGcc                                                                                                                  <  2:1460182/150‑1 (MQ=35)
                                     tttCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGCTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGAc                                                                                                             <  2:1286079/150‑1 (MQ=25)
                                           gCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAg                                                                                                                            >  2:1311317/1‑129 (MQ=37)
                                           gCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAg                                                                                                                            <  1:1311317/129‑1 (MQ=35)
                                                   ccccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACCCCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa                                                                                               <  2:2301536/150‑1 (MQ=17)
                                                   ccccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa                                                                                               <  2:552462/150‑1 (MQ=32)
                                                    cccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAAt                                                                                                                        >  1:2031261/1‑124 (MQ=32)
                                                    cccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAAt                                                                                                                        <  2:2031261/124‑1 (MQ=33)
                                                     ccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTa                                                                                                                                           >  1:1851476/1‑104 (MQ=33)
                                                     ccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTa                                                                                                                                           <  2:1851476/104‑1 (MQ=32)
                                                     ccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGa                                                                                                                             <  2:936825/118‑1 (MQ=33)
                                                     ccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGa                                                                                                                             >  1:936825/1‑118 (MQ=33)
                                                     ccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc                                                                                             <  2:1913290/150‑1 (MQ=32)
                                                     ccGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc                                                                                             >  1:2263081/1‑150 (MQ=32)
                                                         aTCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt                                                                                         <  2:2163236/150‑1 (MQ=31)
                                                                                                 cgccgccttcggtgaccgacgggctcgccgtcagggtcgcggtggtGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCa                                                 >  2:830259‑M2/47‑150 (MQ=255)
                                                                                                  gccgccttcggtgaccgacgggctcgccgtcagggtcgcggtggtGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAg                                                >  1:1423112‑M2/46‑150 (MQ=255)
                                                                                                  gccgccttcggtgaccgacgggctcgccgtcagggtcgcggtggtGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTa                                                                       >  1:1133488‑M2/46‑127 (MQ=255)
                                                                                                  gccgccttcggtgaccgacgggctcgccgtcagggtcgcggtggtGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTa                                                                       <  2:1133488‑M2/82‑1 (MQ=255)
                                                                                                                                               GGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac   <  2:1145342‑M2/150‑1 (MQ=255)
                                                                                                                                               GGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac   <  1:830259‑M2/150‑1 (MQ=255)
                                                                                                                                               GGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaa                                                       <  1:1415435‑M2/98‑1 (MQ=255)
                                                                                                                                               GGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaa                                                       >  2:1415435‑M2/1‑98 (MQ=255)
                                                                                                                                               GGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACaacaac   >  1:1038473‑M2/1‑150 (MQ=255)
                                                                                                                                                GCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                   >  1:201934‑M2/1‑69 (MQ=255)
                                                                                                                                                GCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACg                                                                                   <  2:201934‑M2/69‑1 (MQ=255)
                                                                                                                                                GCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGCCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg  >  1:1738258‑M2/1‑150 (MQ=255)
                                                                                                                                                GCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGCCCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg  <  2:1652334‑M2/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACG  >  NC_002947/196431‑196724

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