Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,709,494 C→A 20.1% C10* (TGC→TGA PP_2372 → sugar‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,709,4940CA20.1% 49.5 / 4.6 24C10* (TGC→TGAPP_2372sugar‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (13/6);  new base A (3/2);  total (16/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGT  >  NC_002947/2709366‑2709634
                                                                                                                                |                                                                                                                                            
aaaaGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCa                                                                                                                         >  1:271901/1‑150 (MQ=255)
        cAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCa                                                                                                                 >  1:200910/1‑150 (MQ=255)
          aaTGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATCCGCAAc                                                                                                               >  1:73443/1‑150 (MQ=255)
                         actgccTGCCACTTGTGCGCTGAGCATCGCCCGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGgtgt                                                                                                >  2:213060/6‑150 (MQ=255)
                                 ccATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTtc                                                                                             <  1:165953/145‑1 (MQ=255)
                                 ccATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTtc                                                                                             >  2:165953/1‑145 (MQ=255)
                                             tGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGGTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTcc                                                                            >  1:446219/1‑150 (MQ=255)
                                                            gTTTATAATCAGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTgg                                                             >  2:581184/1‑150 (MQ=255)
                                                                ataATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                    >  1:166455/1‑107 (MQ=255)
                                                                ataATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                    <  2:166455/107‑1 (MQ=255)
                                                                    tCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                    <  1:226093/103‑1 (MQ=255)
                                                                    tCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCggtg                                                                                                    >  2:226093/1‑103 (MQ=255)
                                                                                     agatATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGACTTCTGGTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAgg                                    >  1:475749/1‑150 (MQ=255)
                                                                                              aCCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCgg                           >  2:494847/1‑150 (MQ=255)
                                                                                               ccccATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACAAGCTATACCCGCTTCATGCAGGCCGcggcggg                          >  2:296816/1‑149 (MQ=255)
                                                                                                      cTCAAGGCCCTGTGCCAAAGCTTGTGACTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGc                                                       <  2:73843/114‑1 (MQ=255)
                                                                                                      cTCAAGGCCCTGTGCCAAAGCTTGTGACTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGc                                                       <  2:73907/114‑1 (MQ=255)
                                                                                                      cTCAAGGCCCTGTGCCAAAGCTTGTGACTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGc                                                       >  1:73907/1‑114 (MQ=255)
                                                                                                      cTCAAGGCCCTGTGCCAAAGCTTGTGACTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGc                                                       >  1:73843/1‑114 (MQ=255)
                                                                                                         aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTgg                  <  2:3988/148‑1 (MQ=255)
                                                                                                         aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTgg                  <  2:4027/148‑1 (MQ=255)
                                                                                                         aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTgg                  >  1:4027/1‑148 (MQ=255)
                                                                                                         aaGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTgg                  >  1:3988/1‑148 (MQ=255)
                                                                                                                       aaGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGt  <  1:213060/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                            
AAAAGAACCAAATGGACTATGTGAGATTGACTGCCATTTGTGCGCTGAGCATCGCCAGGCGTTTATAATCCGGTGTTTGCTTCGCAGAGATGTTACCCCATGCTCAAGGCCCTGTGCCAAAGCTTGTGCCTTGTCTTGCCACTGGTGGCAAATGCGCAACCGGCTTCGGTGGTGTTTCTCAACCCGGGGCTGTCCACCGAGACATTCTGGACCAGCTATACCCGCTTCATGCAGGCCGCCGCGGACGAGCTGGGCATGACCCTGCAGGT  >  NC_002947/2709366‑2709634

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