Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 819,694 T→A 7.3% A58A (GCT→GCA)  ada → fused O(6)‑alkylguanine DNA alkyltransferase/transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947819,6940TA 7.3% 163.1 / 3.0 55A58A (GCT→GCA) adafused O(6)‑alkylguanine DNA alkyltransferase/transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base A (2/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AACGCTGGCAAGCCGTAGAGAGCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGC  >  NC_002947/819546‑819836
                                                                                                                                                    |                                                                                                                                              
aaCGCTGGCAAGCCGTAGAGAGCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCt                                                                                                                                                <  8:28095/149‑1 (MQ=255)
                 gagagCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTa                                                                                                                               >  2:155159/1‑149 (MQ=255)
                       cgcgACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGc                                                                                                                         >  8:76585/1‑149 (MQ=255)
                                  cGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCAGCCTCGGCAGCCgcggca                                                                                                                              >  4:116547/1‑128 (MQ=255)
                                  cGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCAGCCTCGGCAGCCgcggca                                                                                                                              <  3:116547/133‑6 (MQ=255)
                                    cTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGt                                                                                                            <  3:194566/149‑1 (MQ=255)
                                            cACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCaa                                                                                                    <  5:163852/149‑1 (MQ=255)
                                              cTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACCACACTGCCGCGGCTGCCGAGGCAGCCGGCTACCGTGCTTGCCAACGCTGTATCAGCAAgg                                                                                                  >  5:172032/1‑149 (MQ=255)
                                                 tGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGc                                                                                                 <  3:243002/147‑1 (MQ=255)
                                                                 accaccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCAcag                                                                               >  1:248397/1‑149 (MQ=255)
                                                                   caccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCaa                                                                                                    <  1:66232/126‑1 (MQ=255)
                                                                   caccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCaa                                                                                                    >  2:66232/1‑126 (MQ=255)
                                                                   caccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCAcagc                                                                              >  2:75865/1‑148 (MQ=255)
                                                                     ccGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCAcagccag                                                                           <  2:248397/149‑1 (MQ=255)
                                                                        gCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTg                                                                        <  8:82017/149‑1 (MQ=255)
                                                                             tACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGGTAc                                                                   >  3:47409/1‑149 (MQ=255)
                                                                               cTGCGACCCCGGCTGCAATTCACGGCAGGGCAAGCGTATGATCGTCGAGTTGTACGACACTGGCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACtg                                                                 <  4:47409/149‑2 (MQ=255)
                                                                                tGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCa                                                                >  2:12909/1‑149 (MQ=255)
                                                                                    aCCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCCCCCGCCCCAGCCCGCTGGTTACTCACGcc                                                            >  4:28374/1‑149 (MQ=255)
                                                                                    aCCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGcc                                                            >  4:171094/1‑149 (MQ=255)
                                                                                       ccGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCAc                                                               >  6:250975/1‑143 (MQ=255)
                                                                                       ccGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCAc                                                               <  5:250975/143‑1 (MQ=255)
                                                                                          gCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGc                                                      <  3:171094/149‑1 (MQ=255)
                                                                                          gCTGCAAATAACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGc                                                      >  8:207405/1‑149 (MQ=255)
                                                                                              cAACTCACGCCAGGCCACGCGTATGAAGGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTg                                                   <  6:172032/148‑1 (MQ=255)
                                                                                              cAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCa                                                                <  7:49903/135‑1 (MQ=255)
                                                                                              cAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCa                                                                >  8:49903/1‑135 (MQ=255)
                                                                                               aaaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                 <  1:155159/149‑1 (MQ=255)
                                                                                               aaaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                 <  5:132179/149‑1 (MQ=255)
                                                                                                aaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATc                                                >  6:84390/1‑149 (MQ=255)
                                                                                                aaTCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATc                                                <  3:163949/149‑1 (MQ=255)
                                                                                                  tCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCg                                               >  5:174420/1‑148 (MQ=255)
                                                                                                      gCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAg                                                                                              <  4:75299/97‑1 (MQ=255)
                                                                                                      gCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAg                                                                                              >  3:75299/1‑97 (MQ=255)
                                                                                                            ccAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGAc                                    >  2:267033/1‑149 (MQ=255)
                                                                                                             cAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGAcc                                   >  3:85601/1‑149 (MQ=255)
                                                                                                             cAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGAcc                                   <  5:84390/149‑1 (MQ=255)
                                                                                                                  gTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGtt                              <  4:52352/149‑1 (MQ=255)
                                                                                                                          gTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCt                      >  1:172078/1‑149 (MQ=255)
                                                                                                                           tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAgcg                                                                                            >  4:125417/1‑78 (MQ=255)
                                                                                                                           tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAgcg                                                                                            <  3:125417/78‑1 (MQ=255)
                                                                                                                           tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGcc                       >  1:898/1‑147 (MQ=255)
                                                                                                                           tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGcc                       <  2:898/147‑1 (MQ=255)
                                                                                                                           tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTc                     >  5:135724/1‑149 (MQ=255)
                                                                                                                           tCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCAGCCACAGCCCGCTGGTTACTAACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTGCCAAGCCTc                     >  8:240048/1‑149 (MQ=255)
                                                                                                                            cGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGc                        >  1:123019/1‑145 (MQ=255)
                                                                                                                            cGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGc                        <  2:123019/145‑1 (MQ=255)
                                                                                                                             gAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCgg                   <  1:75865/149‑1 (MQ=255)
                                                                                                                                  cTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCaa               <  6:174420/148‑1 (MQ=255)
                                                                                                                                  cTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCaa               >  8:233930/1‑148 (MQ=255)
                                                                                                                                  cTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCaa               <  2:172078/148‑1 (MQ=255)
                                                                                                                                      gACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAg          >  6:31636/1‑149 (MQ=255)
                                                                                                                                       acacTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGt         >  3:94051/1‑149 (MQ=255)
                                                                                                                                              cgcgGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGc  >  7:199198/1‑149 (MQ=255)
                                                                                                                                               ccgGCAGCCTCGGCAGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                 >  4:112129/2‑101 (MQ=255)
                                                                                                                                               ccgGCAGCCTCGGCAGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGAt                                                 <  3:112129/100‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AACGCTGGCAAGCCGTAGAGAGCCGCGACACTGCCGCTACAGGCCACTTTGTCTACGCCGTACGTACCACCGGCGTCTACTGCCACCCCGGCTGCAAATCACGCCAGGCCAAGCGTATGAACGTCGAGTTCTACGACACTGCCGCGGCTGCCGAGGCTGCCGGCTACCGTGCTTGCAAACGCTGTATCAGCAAGGCCAGCGCCACCCGCCACAGCCAGCTGGTTACTCACGCCTGCCGCTTGATCGAAACCAGCGACCCCGTTCCAAGCCTCGGCCAACTCAGTGCACGGC  >  NC_002947/819546‑819836

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