Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 212,011 C→A 23.5% V5839V (GTC→GTA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947212,0110CA23.5% 29.6 / 2.9 17V5839V (GTC→GTAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base A (2/2);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

AGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAG  >  NC_002947/211878‑212143
                                                                                                                                     |                                                                                                                                    
aGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCaa                                                                                                                                    >  1:79576/1‑136 (MQ=34)
            gCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGc                                                                                                                        <  2:500479/136‑1 (MQ=34)
            ccACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGc                                                                                                                        <  2:437027/135‑1 (MQ=21)
                          gacgaGCCAGGTACGCCGGGCAACGAACGCGCACTGATCAAGGTCACCATCCCCCCCGTCCAGACTTCCGTAGCCGAGAACGTCACACCGCCCTTCAATGTTCACTTCAACCAGCCGCTGGCTCGCGCGCTGGTCg                                                                                                          >  1:28778/1‑136 (MQ=11)
                             gaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTgg                                                                                                       >  2:7016/1‑135 (MQ=21)
                               gCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAcc                                                                                                     <  2:194480/136‑1 (MQ=34)
                               gCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAcc                                                                                                     <  2:382394/136‑1 (MQ=34)
                                                                                                      agaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAggg                              >  1:335295/1‑136 (MQ=32)
                                                                                                      agaACGTCAAACCGACCTTCACTGTTCACGTAAACCAGCCGCTGGCTCACGACCTGGTAGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAggg                              >  1:333439/1‑136 (MQ=17)
                                                                                                      agaACGTCAAACCGACCTTCACTGTTCACGTAAACCAGCCGCTGGCTCACGACCTGGTAGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAggg                              >  1:129894/1‑136 (MQ=17)
                                                                                                       gaACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGt                             <  2:259210/136‑1 (MQ=32)
                                                                                                       gaACGTCAAACCGACCTTCACTGTTCACGTAAACCAGCCGCTGGCTCACGACCTGGTAGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGt                             <  2:129894/136‑1 (MQ=17)
                                                                                                       gaACGTCAAACCGACCTTCACTGTTCACGTAAACCAGCCGCTGGCTCACGACCTGGTAGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGt                             <  2:333439/136‑1 (MQ=17)
                                                                                                        aaCGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCCGCGCGCCGTACACCCACGCTGCGCAGGGtg                            >  2:425748/1‑136 (MQ=17)
                                                                                                        aaCGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATACCCAGGTCACCATCAAGGCCCGTGAAACCAGCGCGCCGTACACCCACGCAGCGCACGGtg                            >  2:379480/1‑136 (MQ=255)
                                                                                                             cAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGAc                       <  1:464115/136‑1 (MQ=32)
                                                                                                                                  cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAg  <  1:379480/136‑1 (MQ=34)
                                                                                                                                     |                                                                                                                                    
AGACCCCGGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAG  >  NC_002947/211878‑212143

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