Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,851 C→G 30.7% T1119T (ACC→ACG PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,8510CG30.7% 14.5 / 6.3 13T1119T (ACC→ACGPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base G (2/2);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAAC  >  NC_002947/197723‑197968
                                                                                                                                |                                                                                                                     
gTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCGCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAGCCATTACCg                                                                                                                >  2:297841/1‑136 (MQ=255)
                           gCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCg                                                                                     >  1:289778/1‑136 (MQ=255)
                                         gACCGAAGGCGGCGTGAACACCAACACCGTGACCCACACCAAACCTGCCCAGACCCCCGAAACCCAAACCATGACCAACGCCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGAcc                                                                       >  2:60518/1‑136 (MQ=11)
                                                                         cccTGAGCAACCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGAc                                       >  1:223671/1‑136 (MQ=21)
                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                           >  1:1219/1‑114 (MQ=17)
                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                           >  1:236672/1‑114 (MQ=17)
                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                           <  2:1219/114‑1 (MQ=17)
                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGAc                                                           <  2:236672/114‑1 (MQ=17)
                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGa                                                                                                            <  1:187185/65‑1 (MQ=35)
                                                                           cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCg                                     >  1:272715/1‑136 (MQ=40)
                                                                                       gCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCa                         >  1:171629/1‑136 (MQ=37)
                                                                                                       ccGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACggg         <  2:171629/136‑1 (MQ=40)
                                                                                                              cTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAAc  >  2:235799/1‑136 (MQ=37)
                                                                                                                                |                                                                                                                     
GTCGACACCACTACTGCGACCCTGACCGCCTCGCCAAGCGTGACCGAAGGCGGCGTGATCACCTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCTTGTCCAACGGCCAAACCATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGCTCGACCGTCAGCGTCACCATCGAGAGCGCCACGGGCGGCAAC  >  NC_002947/197723‑197968

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: