Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,950,489 C→T 9.3% D245N (GAC→AAC)  oguA ← 8‑oxoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,950,4890CT 9.3% 118.3 / 2.5 43D245N (GAC→AAC) oguA8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (2/2);  total (24/19)
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

GCCGGTGCCGGTGCGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACGGCGCATCGAGCCGTGGC  >  NC_002947/2950343‑2950626
                                                                                                                                                   |                                                                                                                                         
gccggtgccggtGCGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGt                                                                                                                                         >  6:80719/1‑149 (MQ=255)
gccggtgccggtGCGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGt                                                                                                                                         <  1:8808/149‑1 (MQ=255)
 ccggtgccggtgCGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGtt                                                                                                                                        <  5:114263/149‑1 (MQ=255)
     tgccggtgcGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACg                                                                                                                                    <  1:125051/149‑1 (MQ=255)
          gtgcGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTtctc                                                                                                                               >  4:259278/1‑149 (MQ=255)
               agcgaACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGc                                                                                                                            >  4:66346/1‑147 (MQ=255)
                  gaACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTTGTTGACGTTCTCGGCCAggt                                                                                                                       >  2:143852/1‑149 (MQ=255)
                  gaACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATCGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTTGTTGACGTTCTCGGCCAggt                                                                                                                       >  2:148694/1‑149 (MQ=255)
                      aGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAggtgggt                                                                                                                   <  3:66346/149‑1 (MQ=255)
                        cTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGgtgt                                                                                                                 >  2:271103/1‑149 (MQ=255)
                                cGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTTGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGac                                                                                                         <  1:148694/149‑1 (MQ=255)
                                cGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTTGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGac                                                                                                         <  1:143852/149‑1 (MQ=255)
                                 gTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGaca                                                                                                        >  6:147891/1‑149 (MQ=255)
                                 gTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGaca                                                                                                        <  5:252512/149‑1 (MQ=255)
                                     tGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAggt                                                                                                                       >  6:151995/1‑130 (MQ=255)
                                     tGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAggt                                                                                                                       <  5:151995/130‑1 (MQ=255)
                                            gCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTgg                                                                                             >  5:167015/1‑149 (MQ=255)
                                                               ccATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCt                                                                           <  1:271103/148‑1 (MQ=255)
                                                                cATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTc                                                                         >  6:138440/1‑149 (MQ=255)
                                                                  tACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTcgc                                                                       >  8:262679/1‑149 (MQ=255)
                                                                      tcatGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCAt                                                                   >  5:288445/1‑149 (MQ=255)
                                                                          ggCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAGACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAg                                                                >  1:226189/1‑149 (MQ=255)
                                                                                       aaGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCg                                                  <  6:149436/149‑1 (MQ=255)
                                                                                         ggTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGAcc                                                   >  1:136232/1‑146 (MQ=255)
                                                                                                 gCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGa                                        <  5:138440/149‑1 (MQ=255)
                                                                                                  cGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAg                                       <  6:167015/149‑1 (MQ=255)
                                                                                                    tACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCa                                     >  6:285576/1‑149 (MQ=255)
                                                                                                      cTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATg                                   >  3:147801/1‑149 (MQ=255)
                                                                                                          gCAGGGGTCATACCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGgcg                                                                                           <  1:109964/89‑1 (MQ=255)
                                                                                                          gCAGGGGTCATACCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGgcg                                                                                           >  2:109964/1‑89 (MQ=255)
                                                                                                                gTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCAc                         <  6:288445/149‑1 (MQ=255)
                                                                                                                gTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCAc                         >  6:125345/1‑149 (MQ=255)
                                                                                                                      ccGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGc                                          <  5:298603/126‑1 (MQ=255)
                                                                                                                      ccGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGc                                          >  6:298603/1‑126 (MQ=255)
                                                                                                                       cGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCa                          <  5:8059/141‑1 (MQ=255)
                                                                                                                       cGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCa                          >  6:8059/1‑141 (MQ=255)
                                                                                                                       cGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACgg                  >  7:150504/1‑149 (MQ=255)
                                                                                                                       cGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACgg                  >  7:154731/1‑149 (MQ=255)
                                                                                                                                  gcgGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACGGCGCATCGAg         <  2:226189/148‑1 (MQ=255)
                                                                                                                                   cgGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACGGCGCATCGAGcc       >  1:51302/1‑149 (MQ=255)
                                                                                                                                        gTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACGGCGCATCGAGCCGtggc  <  7:262679/149‑1 (MQ=255)
                                                                                                                                               aTGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACgg                  <  6:9426/126‑1 (MQ=255)
                                                                                                                                               aTGTCGTAGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACgg                  >  5:9426/1‑126 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GCCGGTGCCGGTGCGAGCGAACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGA‑ACTTCTCGCGGCTGTAGGCGATGTCGTTGACGTTCTCGGCCAGGTGGGTGTGCAGCGACACCCCATACTGGCGCGCCAGCACGGCAGCTTCGCGCATCAGGTCGCGGCTGACCGAGAATGGCGAGCATGGCGCGACCACGATACGGCGCATCGAGCCGTGGC  >  NC_002947/2950343‑2950626

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