Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,879 G→A 19.5% V795V (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,8790GA19.5% 76.4 / 14.1 41V795V (GTG→GTA) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (18/15);  new base A (6/2);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGT  >  NC_002947/196736‑196937
                                                                                                                                                |                                                          
aGCGTCACCATCGCTAACGCCACTGGGGGCAACTTCTAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCTGTCACCGAAGGCGGCTTGATCAc                                                       <  8:145383/149‑1 (MQ=11)
aGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCATCGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                       <  5:351045/149‑1 (MQ=17)
aGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCTTCGGTCACCGAAGGCGGCGTGATCAc                                                       <  1:218405/149‑1 (MQ=17)
aGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                       <  4:224829/149‑1 (MQ=32)
        cATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCg                                               >  1:70357/1‑149 (MQ=32)
         aTCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGCG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGt                                              >  5:252618/1‑149 (MQ=17)
          tCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGt                                              >  2:195155/1‑148 (MQ=32)
            gagaACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGAc                                           >  1:145/1‑149 (MQ=32)
             agaACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGAcc                                          >  1:271273/1‑149 (MQ=33)
                  gCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACCCCGTGACCCTGAg                                     >  2:273730/1‑149 (MQ=21)
                  gCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAg                                     <  3:377593/149‑1 (MQ=34)
                  gCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAg                                     <  1:19148/149‑1 (MQ=34)
                   ccACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCCCTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAg                                     <  1:287822/149‑1 (MQ=255)
                    cACTGGTGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGc                                    >  4:325387/1‑148 (MQ=21)
                       tGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                       <  6:322075/126‑1 (MQ=11)
                       tGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                                >  7:271007/1‑149 (MQ=32)
                       tGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                                >  6:181810/1‑149 (MQ=32)
                         gcggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCCCCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCt                              >  6:8700/1‑149 (MQ=17)
                         gcggcAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATc                                <  8:271007/147‑1 (MQ=32)
                                cTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAATGCGGCCTGATCACCTACACCGTGACCCTGAGCAATCCTCCCCAGa                       >  8:243672/1‑149 (MQ=12)
                                 ttCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGAc                      <  1:381046/149‑1 (MQ=32)
                                    gAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                   >  8:92415/1‑149 (MQ=32)
                                     aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                  >  3:250105/1‑149 (MQ=32)
                                     aaCAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCg                  <  7:92415/149‑1 (MQ=32)
                                      accgCTGACCCCGAACCCGACGCCGGCTCAGACCACTATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCgg                 <  4:356045/146‑1 (MQ=17)
                                             aCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAg           <  2:70357/148‑1 (MQ=32)
                                             aCCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAg           <  7:243672/148‑1 (MQ=32)
                                             aCCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGACACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGAAATCCTGCCCAGACGCCGGTGACAGt          >  8:354810/1‑149 (MQ=255)
                                              ccccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGt          >  6:384842/1‑148 (MQ=32)
                                              ccccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGt          <  7:354810/148‑1 (MQ=17)
                                               cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa        >  1:167746/1‑149 (MQ=32)
                                               cccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa        >  7:205466/1‑149 (MQ=32)
                                               cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa        >  6:381818/1‑149 (MQ=17)
                                               cccGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGa        >  6:334323/1‑149 (MQ=17)
                                                ccGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAc       >  6:40036/1‑149 (MQ=32)
                                                 cGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc      <  6:252618/149‑1 (MQ=32)
                                                 cGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc      >  3:402755/1‑149 (MQ=32)
                                                 cGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAcc      >  2:303119/1‑149 (MQ=17)
                                                     cccGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt  <  5:334323/149‑1 (MQ=11)
                                                     cccGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt  >  5:308013/1‑149 (MQ=11)
                                                     cccGACACCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGt  >  8:345878/1‑149 (MQ=11)
                                                                                                                                                |                                                          
AGCGTCACCATCGAGAACGCCACTGGCGGCAACTTCGAACAGCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAG‑CCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGT  >  NC_002947/196736‑196937

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