Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,950,337 G→A 9.5% A295A (GCC→GCT)  oguA ← 8‑oxoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,950,3370GA 9.5% 112.4 / 2.7 42A295A (GCC→GCT) oguA8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAACTTCTCGCGGC  >  NC_002947/2950189‑2950476
                                                                                                                                                    |                                                                                                                                            
gcagGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTgg                                                                                                                                              >  3:264452/1‑149 (MQ=255)
     gCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGAc                                                                                                                                         >  5:228604/1‑149 (MQ=255)
                ttCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACgccggtgccgg                                                                                                                              >  7:141886/1‑149 (MQ=255)
                   aCCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCcggtgccggtgc                                                                                                                           >  6:337720/1‑149 (MQ=255)
                      gATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTgcgag                                                                                                                        <  6:228604/149‑1 (MQ=255)
                      gATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTgcgag                                                                                                                        >  3:361451/1‑149 (MQ=255)
                                gcgcCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTc                                                                                                             <  6:45869/149‑1 (MQ=255)
                                 cgcCGTCGTTGGAGGCAGAGCCATCTACCCCCAGGCCCACCGGTACGCCATGGTCGCGCATCTTGCGTACCGGGGCCATGCCAGAGGCCAGCCGCATGTTCAAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTc                                                                                                             >  7:234438/1‑148 (MQ=255)
                                            gAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTg                                                                                                 <  5:337720/149‑1 (MQ=255)
                                                    gCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCt                                                                                         <  8:105514/149‑1 (MQ=255)
                                                     ccATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTg                                                                                        >  7:95058/1‑149 (MQ=255)
                                                       aTCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCa                                                                                      <  8:335486/149‑1 (MQ=255)
                                                       aTCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCa                                                                                      >  2:66915/1‑149 (MQ=255)
                                                               ccAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTg                                                                               <  3:134467/148‑1 (MQ=255)
                                                               ccAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTg                                                                               >  4:134467/1‑148 (MQ=255)
                                                                cAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTggg                                                                             <  4:298074/149‑1 (MQ=255)
                                                                   gCCCACCGGTACGCCGTGGTCGCGCATCTTGCGAACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCAAGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGACGTCCTGGGCGAGCTGCACGCAGAGGGCAt                                                                          >  8:229863/1‑149 (MQ=255)
                                                                   gCCCACCGGTACGCCGTGGTCGCGCATCTCGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGc                                                                                                               <  6:179563/112‑1 (MQ=255)
                                                                   gCCCACCGGTACGCCGTGGTCGCGCATCTCGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGc                                                                                                               >  5:179563/1‑112 (MQ=255)
                                                                     ccACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGAGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGc                                                                        >  5:362982/1‑149 (MQ=255)
                                                                     ccACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGAGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGc                                                                        >  5:357764/1‑149 (MQ=255)
                                                                           gTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATAca                                                                  >  2:168898/1‑149 (MQ=255)
                                                                           gTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGAGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATAca                                                                  <  6:357764/149‑1 (MQ=255)
                                                                           gTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGAGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATAca                                                                  <  6:362982/149‑1 (MQ=255)
                                                                                ccGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATg                                                             <  1:31562/149‑1 (MQ=255)
                                                                                ccGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATg                                                             <  1:66915/149‑1 (MQ=255)
                                                                                ccGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGATGCGAGCGAGAGAGCAAGAAGCCGTTCAGGTAGAAGATCACCGAATGGGCATGCCATAcaaaat                                                              >  1:304653/1‑145 (MQ=255)
                                                                                cagTGGTCGCGCATGTTGCGTAGCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTTGTCGATCTGCACGCAGTGGGCATGCCATACATCATg                                                             <  8:141886/147‑1 (MQ=255)
                                                                                       cgcgCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGAcc                                                      >  3:252496/1‑149 (MQ=255)
                                                                                         cgcATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCa                                                    <  4:264452/149‑1 (MQ=255)
                                                                                                 gCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAg                                              <  8:95058/147‑1 (MQ=255)
                                                                                                       cGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCgg                                      >  3:329342/1‑149 (MQ=255)
                                                                                                              aTGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTc                               >  5:378649/1‑149 (MQ=255)
                                                                                                                 ccTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCgg                             <  8:211772/148‑1 (MQ=255)
                                                                                                                   tGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAg                          >  1:270416/1‑149 (MQ=255)
                                                                                                                   tGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGCCCGACCCAGCCAACGTCTTCGGCGTACTCGCCac                          >  2:320243/1‑148 (MQ=255)
                                                                                                                     aGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAggg                        <  2:270416/149‑1 (MQ=255)
                                                                                                                              gCATGTTCGAGCATGGGCAGTGGGCGCCGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGcc                >  1:242661/1‑148 (MQ=255)
                                                                                                                                 tGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAAc            <  1:168898/149‑1 (MQ=255)
                                                                                                                                   ttCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAACtt          >  2:247464/1‑149 (MQ=255)
                                                                                                                                     cGAGGATGGGCAGTGGGCGACGCCGGTTCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCTCGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCTAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAACTtc         <  7:229863/148‑1 (MQ=255)
                                                                                                                                      gAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAACTtct        <  1:247464/148‑1 (MQ=255)
                                                                                                                                           tGGGCAGGGGGCGACGCCGGTGCCGGTGCGAGCGA‑AGAGATCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAACTTCTCGCGGc  <  4:252496/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCCGTGGTCGCGCATCTTGCGTACCGGGGCGATGCCTGAGGCCAGCCGCATGTTCGAGCATGGGCAGTGGGCGACGCCGGTGCCGGTGCGAGCGA‑ACAGCTCGATGCCGTGCTGGTCGAGCTGCACGCAGTGGGCATGCCATACATCATGGCCGACCCAGCCAAGGTCTTCGGCGTACTCGGCAGGGGTCATGCCGAACTTCTCGCGGC  >  NC_002947/2950189‑2950476

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