Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,951,136 T→C 8.5% K29R (AAG→AGG)  oguA ← 8‑oxoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,951,1360TC 8.5% 131.0 / 2.8 47K29R (AAG→AGG) oguA8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (25/18);  new base C (2/2);  total (27/20)
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.39e-01

GCTCTGGTACATGTGGTGGTGGGTGTTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAG  >  NC_002947/2950988‑2951283
                                                                                                                                                    |                                                                                                                                                   
gCTATGGTACATGTGGTGGTGGGTGTTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTAGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTAGATGAACAGACCGCCGTTCt                                                                                                                                                     >  4:311317/1‑149 (MQ=255)
     ggTACATGTGGTGGTGGGTGTTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGAtt                                                                                                                                                >  2:260509/1‑149 (MQ=255)
         cATGTGGTGGTGGGTGTTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTcgc                                                                                                                                            >  8:4547/1‑149 (MQ=255)
          aTGTGGTGGTGGGTGTTGACCAGGCCTGGGATGCCCACCTTGGCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGTGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTcgcg                                                                                                                                           <  6:319979/149‑1 (MQ=255)
               gtggtggGTGTTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGtt                                                                                                                                      <  8:25123/149‑1 (MQ=255)
                      gtgtTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTcgcg                                                                                                                                           <  8:41494/137‑1 (MQ=255)
                      gtgtTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTcgcg                                                                                                                                           >  7:41494/1‑137 (MQ=255)
                        gtTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCa                                                                                                                             >  2:280403/1‑149 (MQ=255)
                         ttGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCa                                                                                                                             <  3:52254/148‑1 (MQ=255)
                              cAGGCCTGGGATGACCATCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGa                                                                                                                                                  >  4:137041/1‑122 (MQ=255)
                              cAGGCCTGGGATGACCATCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGa                                                                                                                                                  <  3:137041/122‑1 (MQ=255)
                               aGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGAc                                                                                                                      <  7:4547/149‑1 (MQ=255)
                                   cTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACcag                                                                                                                   >  2:20195/1‑148 (MQ=255)
                                   cTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACcag                                                                                                                   >  7:247029/1‑148 (MQ=255)
                                     gggATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATg                                                                                                                           <  1:10414/138‑1 (MQ=255)
                                     gggATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATg                                                                                                                           >  2:10414/1‑138 (MQ=255)
                                      ggATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGt                                                                                                               >  6:294223/1‑149 (MQ=255)
                                               ccTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGttt                                                                                                      <  3:312588/149‑1 (MQ=255)
                                                 ttGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGttttt                                                                                                    >  2:33138/1‑149 (MQ=255)
                                                   gCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGttttt                                                                                                    <  3:301136/147‑1 (MQ=255)
                                                   gCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGttttt                                                                                                    <  1:260509/147‑1 (MQ=255)
                                                                       cGTCGGCGTGCTGCGTGAGGGTTTCGCTGGGGACGACTTGCTTGATCAGCATGCCTTCGATGAAAAGACCGCCGTTCTTGAATTCGCTGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAggg                                                                              <  6:250271/149‑1 (MQ=255)
                                                                       cGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGt                                                                                                                                >  5:294618/1‑99 (MQ=255)
                                                                       cGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGt                                                                                                                                <  6:294618/99‑1 (MQ=255)
                                                                           ggCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGttt                                                                           <  1:280403/148‑1 (MQ=255)
                                                                           ggCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCCTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTg                                                                          >  6:81289/1‑149 (MQ=255)
                                                                           ggCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCCTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTg                                                                          >  6:76711/1‑149 (MQ=255)
                                                                            gCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGt                                                                         <  3:311317/149‑1 (MQ=255)
                                                                            gCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGt                                                                         <  1:20195/149‑1 (MQ=255)
                                                                                 ctgcGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGACGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGAt                                                                     >  3:265924/1‑148 (MQ=255)
                                                                                  tgcGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGa                                                                   >  8:26845/1‑149 (MQ=255)
                                                                                           gTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGc                                                           >  1:44984/1‑148 (MQ=255)
                                                                                            tttCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGc                                                           >  3:224406/1‑147 (MQ=255)
                                                                                             ttCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAg                                                        >  3:139232/1‑149 (MQ=255)
                                                                                                 cTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATg                                                                    >  2:270691/1‑133 (MQ=255)
                                                                                                 cTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATg                                                                    <  1:270691/133‑1 (MQ=255)
                                                                                                  tGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTATGCAGGGTTTGTAGATGACATGGTGCAAGCCTgc                                                   >  4:98189/1‑149 (MQ=255)
                                                                                                                       ggTTGCCTTCGACGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGACGTCCATGGTGACCAGCAGTGCGGCGTTTTTTACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTAc                              >  5:2534/1‑149 (MQ=255)
                                                                                                                             cTTCGATGAACAGACCGCCGTTCCTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGt                        <  5:81289/149‑1 (MQ=255)
                                                                                                                             cTTCGATGAACAGACCGCCGTTCCTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGt                        <  5:76711/149‑1 (MQ=255)
                                                                                                                              ttCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGt                        >  5:268741/1‑148 (MQ=255)
                                                                                                                                       cAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCt              >  1:134048/1‑149 (MQ=255)
                                                                                                                                        aGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAAt                                              <  7:112412/116‑1 (MQ=255)
                                                                                                                                        aGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAAt                                              >  8:112412/1‑116 (MQ=255)
                                                                                                                                           ccgccgTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCgg           >  2:124303/1‑148 (MQ=255)
                                                                                                                                            cgccgTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTa         >  2:253718/1‑149 (MQ=255)
                                                                                                                                                gTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGt     <  5:294223/149‑1 (MQ=255)
                                                                                                                                                   cTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTaag  >  8:226892/1‑149 (MQ=255)
                                                                                                                                                    ttGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCt                                                     >  7:239116/1‑97 (MQ=255)
                                                                                                                                                    ttGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCt                                                     <  8:239116/97‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTCTGGTACATGTGGTGGTGGGTGTTGACCAGGCCTGGGATGACCACCTTGCCGGCCATGTCCAGAACCACGTCGGCGTGCTGCGGGAGGGTTTCGCTGGGGCCGACTTGCTTGATCAGGTTGCCTTCGATGAACAGACCGCCGTTCTTGATTTCGCGGCGTTCGCCGTCCATGGTGACCAGCAGTGCGGCGTTTTTGACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAG  >  NC_002947/2950988‑2951283

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