Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,390 T→C 27.3% D632D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3900TC27.3% 29.5 / 10.2 22D632D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (12/4);  new base C (3/3);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAAC  >  NC_002947/196275‑196516
                                                                                                                   |                                                                                                                              
gcttGATCACTTCCACCGTTAGCGTTATAATTCCTGCCCAGACGCCGGTGACATTGACCCTGTCAACGGCCAAACCATTACCGTTTAAGCCGGCAAAACCCAGTGCAGCGTCGATTTCCAGACCCCGGCGAATGat                                                                                                           <  2:306038/133‑2 (MQ=11)
   tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGt                                                                                                         <  1:154978/136‑1 (MQ=255)
   tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGt                                                                                                         <  1:193195/136‑1 (MQ=255)
    gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGCCCCCGGCGAATGACGtg                                                                                                        >  1:87146/1‑135 (MQ=11)
              acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                              >  2:195743/1‑136 (MQ=255)
                             aacccTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgcg                                                                               >  1:163769/4‑134 (MQ=11)
                             aacccTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgcg                                                                               <  1:350023/133‑3 (MQ=11)
                             aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgtgg                                                                              >  2:232876/1‑135 (MQ=255)
                             aaTCCTGCCCAGACGACGGTGACAGAGACCCTGACCACGGCCAAACCATTAGCGTTGAAGCCGGCATAACCCCGGGCAGCGTCGATTTCCAGACGCCGGCGAATGACGTCTCCAACAACGGTGGGTCGGTCAgtgg                                                                              >  2:483139/1‑135 (MQ=11)
                                                 gACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcg                                                          >  1:435919/1‑136 (MQ=255)
                                                         cccTGTCCACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAAt                                                  >  2:455948/1‑136 (MQ=255)
                                                            tGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAAttt                                                >  1:11170/1‑136 (MQ=37)
                                                                  aCGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCa                                          >  1:499779/1‑136 (MQ=37)
                                                                       cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                   >  1:239871/1‑74 (MQ=34)
                                                                       cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                   <  2:239871/74‑1 (MQ=34)
                                                                                       aaGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATg                                                                                                             <  2:487669/48‑1 (MQ=38)
                                                                                       aaGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATg                                                                                                             >  1:487669/1‑48 (MQ=38)
                                                                                        aGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCg                    >  2:241438/1‑136 (MQ=40)
                                                                                         gCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACTGGCGGCAAATTCGAGCAGCTGACCCCGAAACCGACGCCgg                   >  1:231835/1‑136 (MQ=11)
                                                                                             gaaaaaCCCATGGCAGCGTCGATTTCCAGACCCCTGCGATTGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCa               <  1:290535/134‑1 (MQ=255)
                                                                                                          cAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGGAATTTCGAGCAACTGACCCCGAAACCGACGCCGGCTCAGACCACGATCAAc  >  1:90607/1‑136 (MQ=255)
                                                                                                           aGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGAc                                    >  1:51184/1‑101 (MQ=35)
                                                                                                                   |                                                                                                                              
GCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAAC  >  NC_002947/196275‑196516

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: