Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,179 A→G 15.9% V895V (GTA→GTG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,1790AG15.9% 45.2 / 2.5 25V895V (GTA→GTGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (7/14);  new base G (2/2);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGT  >  NC_002947/197047‑197313
                                                                                                                                    |                                                                                                                                      
cGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATc                                                                                                                                     <  2:25701/136‑1 (MQ=37)
 gagaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCa                                                                                                                                    <  2:333160/136‑1 (MQ=37)
  ccaccGCCACTGGCGTTAATTTCTAGCAACTCACCCCGAACCCGCCCCCGGCTCAGACCACGTTCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                                                                                                                   <  2:233116/132‑1 (MQ=11)
  agaACGCCACAGGCGTTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCAc                                                                                                                                   <  2:231112/136‑1 (MQ=37)
  agaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAc                                                                                                                                   <  2:458417/136‑1 (MQ=17)
   gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAcc                                                                                                                                  >  1:167417/1‑136 (MQ=17)
   gaACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCAcc                                                                                                                                  >  1:169676/1‑136 (MQ=17)
       gCCACAGGGGGTAATTTCGAGAAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTaca                                                                                                                              <  2:60225/136‑1 (MQ=255)
       gCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTaca                                                                                                                              <  2:436581/136‑1 (MQ=37)
            aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                                                                                                         <  2:241954/136‑1 (MQ=35)
            aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                                                                                                         >  1:352207/1‑136 (MQ=35)
            aGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                                                                                                         >  2:319210/1‑136 (MQ=35)
                    aTTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAg                                                                                                                 >  2:102127/1‑136 (MQ=34)
                      ttCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCa                                                                                                               >  2:426916/1‑136 (MQ=34)
                      ttCGAGAAATTTACCCCTAATCGCACCCCGGGTCAGACAACCTTAAACGTCTCGGTCGTTGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCa                                                                                                               <  1:462618/136‑1 (MQ=11)
                                      ccaagccGACCCCGGCCCAGCCCACGATCAACGAGTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGcc                                                                                               <  2:357753/131‑1 (MQ=11)
                                           cccaCCCCGCCTCCCACCAGTATCAAGGACGCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACGGAAGGCGGCGTAATCACCTACACCGTGACGCTGAGCAATCCTGCCCAGACGCCGGTaa                                                                                          <  2:8029/133‑1 (MQ=11)
                                                                    cccaCTCGGTCGTTGCCACCCCCGCGACCCTGACGTCGTGCCCGTCGGTCACGGAAGGCGTCGTATTCACCTTCACCGTGAGCCTGAGCATTCGTGCCCTGACGCCGGTATCCGTGACCCTGTCCAACGGCCAGGt                                                                 <  2:467905/133‑1 (MQ=11)
                                                                    aCGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGt                                                                 <  1:102127/136‑1 (MQ=37)
                                                                      caCTCGGTCCATGCCACCACCGCGACCCTGCCTGGGTGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGtca                                                               <  1:281569/135‑1 (MQ=14)
                                                                       aCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                                                                                                         <  2:302147/77‑1 (MQ=34)
                                                                       aCTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTg                                                                                                                         >  1:302147/1‑77 (MQ=34)
                                                                                     accaccGCGACCCTGACGGCGAGCCGGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCGTGCCCACACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCg                                                <  2:187296/136‑1 (MQ=255)
                                                                                                             ccGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAtt                        >  1:289991/1‑136 (MQ=255)
                                                                                                                                   tAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGt  >  1:443419/1‑136 (MQ=37)
                                                                                                                                    |                                                                                                                                      
CGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACCCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCACCACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTAATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAGGTCATCACTGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGT  >  NC_002947/197047‑197313

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: