Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,504 G→C 11.4% R107R (CGG→CGC)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,5040GC11.4% 94.2 / 5.0 35R107R (CGG→CGC) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base C (2/2);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CCTGGCAGGCAGCCAAGGTGCATCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGT  >  NC_002947/3528362‑3528629
                                                                                                                                               |                                                                                                                             
ccTGGCAGGCAGCCAAGGTGCATCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGTGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCa                                                                                                                         <  1:63152/149‑1 (MQ=255)
          agCCAAGGTGCATCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAggggg                                                                                                               >  5:132440/1‑149 (MQ=255)
                ggTGCATCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTAcc                                                                                                         <  4:100816/149‑1 (MQ=255)
                 gTGCATCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTAccc                                                                                                        >  4:53619/1‑148 (MQ=255)
                     aTCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTg                                                                                                    <  4:90617/149‑1 (MQ=255)
                          ttGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCg                                                                                               <  8:75402/149‑1 (MQ=255)
                            gCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTTGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGt                                                                                              <  2:130333/148‑1 (MQ=255)
                             cGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGtga                                                                                            >  5:11099/1‑149 (MQ=255)
                              ggCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGtgat                                                                                           >  8:128404/1‑149 (MQ=255)
                              ggCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGtgat                                                                                           >  1:165774/1‑149 (MQ=255)
                                   aTCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCCGAAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACt                                                                                       >  1:44529/1‑149 (MQ=255)
                                       ttGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTtgctg                                                                                  >  6:32272/1‑149 (MQ=255)
                                                      aaTTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGc                                                                   <  6:132440/149‑1 (MQ=255)
                                                               aCGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGCCGCCGGGCacgg                                                                                                                      <  6:190963/89‑4 (MQ=255)
                                                               aCGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGCCGCCGGGCacgg                                                                                                                      >  5:190963/1‑86 (MQ=255)
                                                                 gggCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCtt                                                        >  4:7270/1‑149 (MQ=255)
                                                                                  tgctgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCctgctg                                                             <  1:111088/128‑1 (MQ=255)
                                                                                  tgctgGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCctgctg                                                             >  2:111088/1‑128 (MQ=255)
                                                                                       gTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg                                   <  3:7270/149‑1 (MQ=255)
                                                                                       gTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg                                   <  6:11099/149‑1 (MQ=255)
                                                                                       gTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTgg                                   >  5:168656/1‑149 (MQ=255)
                                                                                          gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCct                                                                 >  2:244019/1‑116 (MQ=255)
                                                                                          gAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCct                                                                 <  1:244019/116‑1 (MQ=255)
                                                                                           aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAg                                                                                                       >  7:89397/1‑77 (MQ=255)
                                                                                           aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAg                                                                                                       <  8:89397/77‑1 (MQ=255)
                                                                                           aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTgg                                                            >  2:111198/1‑120 (MQ=255)
                                                                                           aaaaCCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTgg                                                            <  1:111198/120‑1 (MQ=255)
                                                                                                     gggTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACt                            <  2:229674/142‑1 (MQ=255)
                                                                                                     gggTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACt                            >  1:229674/1‑142 (MQ=255)
                                                                                                            cAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTg              <  2:44529/149‑1 (MQ=255)
                                                                                                                   tgtgCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACt                            <  4:78274/128‑1 (MQ=255)
                                                                                                                   tgtgCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACt                            >  3:78274/1‑128 (MQ=255)
                                                                                                                        gcgcTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAgtgt  <  2:165774/149‑1 (MQ=255)
                                                                                                                                     ctGCCAGCCGCCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTtgctgctgc                                                                              <  7:14999/59‑1 (MQ=255)
                                                                                                                                     ctGCCAGCCGCCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTtgctgctgc                                                                              >  8:14999/2‑60 (MQ=255)
                                                                                                                                               |                                                                                                                             
CCTGGCAGGCAGCCAAGGTGCATCAGTTGCGGCTGATCGTTGGCAGCGAGGTGCAATTGTGCGACGGGCCC‑AAGCTGGTGCTGCTGGTCGAAAACCTGACGGGTTACCAGAACCTGTGCGCGCTGATCACCCGTGCCCGGCGGCGGGCAGAGAAGGGGGCGTACCAGCTGTTTCGTGATGACTTGCTGCTGCACCACCAAGGCCTGCTGGCCTTGTGGGTGGCCGCTGACAGTGGCGATACTGCTACCGGCGCGTGGCTGCGCAGTGT  >  NC_002947/3528362‑3528629

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: ‑