Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,464,830 C→A 16.0% P84P (CCC→CCA holA → DNA polymerase III subunit delta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,464,8300CA16.0% 69.1 / 2.6 25P84P (CCC→CCAholADNA polymerase III subunit delta
Reads supporting (aligned to +/- strand):  ref base C (11/10);  new base A (2/2);  total (13/12)
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.72e-01

CCGCGCGTCAGCAAGGTTTCGATGAACGCCAGGTATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGACCAAGTGGGGCAAGGCATT  >  NC_002947/5464697‑5464964
                                                                                                                                     |                                                                                                                                      
ccGCGCGTCAGCAAGGTTTCGATGAACGCCAGGTATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTc                                                                                                                                      >  1:72265/1‑136 (MQ=255)
      gTCAGCAAGGTTTCGATGAACGCCAGGTATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCaa                                                                                                                                >  2:295419/1‑136 (MQ=255)
      gTCAGCAAGGTTTCGATGAACGCCAGGTATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCATCGGGCaa                                                                                                                                >  1:318117/1‑136 (MQ=255)
      gTCAGCAAGGTTTCGATGAACGCCAGGTATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCATCGGGCaa                                                                                                                                >  1:41073/1‑136 (MQ=255)
                                 tATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCcgctgcgct                                                                                                     >  2:46152/1‑136 (MQ=255)
                                          ccGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATa                                                                                            >  1:384931/1‑136 (MQ=255)
                                             aCGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTg                                                                                         <  2:402590/136‑1 (MQ=255)
                                             aCGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTg                                                                                         >  2:287148/1‑136 (MQ=255)
                                                                 aCCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGaca                                                                     <  2:286574/136‑1 (MQ=255)
                                                                 aCCCTGCTCCAGGCCGGTGCAAGCCTGACGCTGTTCGCCCAGCGCCGCCTCCAGGAAATGCGCCTGCCCTCGGGAAAGCCAGGTGACAAGGGTCCCGCTCCGCTGATGGAAAAGTGCGCCACCCCCGCCGAAGacc                                                                     >  2:335832/1‑135 (MQ=255)
                                                                      gCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACgctg                                                                >  1:245975/1‑136 (MQ=255)
                                                                          cAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACgctgctgc                                                            >  2:458079/1‑136 (MQ=255)
                                                                               cGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATc                                                       <  1:313466/136‑1 (MQ=255)
                                                                                    cAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCt                                                  >  1:42897/1‑136 (MQ=255)
                                                                                          tGTCGCTGTTCGCCCAGCACCACCTGCTGGAACTCCGCCTACCCTCCCGCAACCCCGGTGACAAGAGTGCCGCTACGCTGATGGAATACTACACCAACCCCCCCGAAAACACACAGCAGCTGATCAACCTACCCaa                                            >  2:265699/1‑136 (MQ=255)
                                                                                               cTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTcc                                       >  2:266144/1‑135 (MQ=255)
                                                                                                    cGCCCAGCGCCGCCTGCTGGAACTGCGCCGGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGAcagc                                  <  2:336026/136‑1 (MQ=255)
                                                                                                         aGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTgc                             <  2:384931/136‑1 (MQ=255)
                                                                                                         aGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTgc                             <  1:287148/136‑1 (MQ=255)
                                                                                                                 ctgctgGAACTGCGCCTGCCATCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCagaaga                     <  2:318117/136‑1 (MQ=255)
                                                                                                                 ctgctgGAACTGCGCCTGCCATCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCagaaga                     <  2:41073/136‑1 (MQ=255)
                                                                                                                   gctgGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGAcc                   <  1:474683/136‑1 (MQ=255)
                                                                                                                   gctgGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGAcc                   <  1:46152/136‑1 (MQ=255)
                                                                                                                         aCTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGACCAAGTgt             <  1:458079/136‑2 (MQ=255)
                                                                                                                         aCTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGACCAAGTgg             <  2:355860/136‑1 (MQ=255)
                                                                                                                                    ccTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGACCAAGTGGGGCAAGGCAtt  >  1:32708/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                      
CCGCGCGTCAGCAAGGTTTCGATGAACGCCAGGTATTCAGTGCCGACGCCAACTTCGACTGGGGCACCCTGCTCCAGGCCGGTGCAAGCCTGTCGCTGTTCGCCCAGCGCCGCCTGCTGGAACTGCGCCTGCCCTCGGGCAAGCCCGGTGACAAGGGTGCCGCTGCGCTGATGGAATACTGCGCCAACCCCGCCGAAGACACGCTGCTGCTGATCAGCCTGCCCAAGCTCGACAGCAGTGCGCAGAAGACCAAGTGGGGCAAGGCATT  >  NC_002947/5464697‑5464964

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