Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,495,409 A→G 8.0% *206W (TGA→TGG)  PP_2188 → tRNA 2‑methylthio‑N6‑isopentenyl adenosine(37) hydroxylase MiaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,495,4090AG 8.0% 136.3 / 3.1 50*206W (TGA→TGG) PP_2188tRNA 2‑methylthio‑N6‑isopentenyl adenosine(37) hydroxylase MiaE
Reads supporting (aligned to +/- strand):  ref base A (27/19);  new base G (2/2);  total (29/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 = 7.54e-01

CTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTGCT  >  NC_002947/2495266‑2495554
                                                                                                                                                |                                                                                                                                                 
cTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATggg                                                                                                                                               >  1:164557/1‑148 (MQ=255)
       tGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTAGGCCGCCTGATGGGCACCAggc                                                                                                                                       >  6:167334/1‑149 (MQ=255)
          cgcACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAggcggc                                                                                                                                    <  4:21699/149‑1 (MQ=255)
            cacaACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGGTGCCCATCAggcggcca                                                                                                                                  <  6:141409/149‑3 (MQ=255)
             acaACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGgcggcg                                                                                                                                   >  7:178953/1‑147 (MQ=255)
              ctaCTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGGACGCGCGGCGGAAATGGAGCTG‑ATCCAGGCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCGCGCAGGCTTTTGCCGCCTGATGGGCACCAGGCGGCGCAg                                                                                                                                <  5:167334/147‑1 (MQ=255)
                 cTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGTCCGCCTGATGGGCACCAGGCGGCGCAGGcc                                                                                                                             <  2:41645/149‑1 (MQ=255)
                  tACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCt                                                                                                                            >  1:167597/1‑149 (MQ=255)
                           gaAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGGTGCCCATCAggcggccaaag                                                                                                                               >  8:223566/1‑132 (MQ=255)
                           gaAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGGTGCCCATCAggcggccaaag                                                                                                                               <  7:223566/137‑6 (MQ=255)
                              gCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCgg                                                                                                                >  1:110118/1‑149 (MQ=255)
                                    aTCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTGCATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAg                                                                                                           >  1:224725/1‑149 (MQ=255)
                                            ctGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGcc                                                                                                  >  7:263458/1‑149 (MQ=255)
                                                   gAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGc                                                                                            >  8:93827/1‑148 (MQ=255)
                                                         gTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGccc                                                                                      >  4:10099/1‑148 (MQ=255)
                                                          tACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGa                                                                                    <  7:126661/149‑1 (MQ=255)
                                                              cgcgGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGcc                                                                                <  8:178953/149‑1 (MQ=255)
                                                               gcggcgGAAATGGAGCTG‑ATCCAGTCACCAGACAAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGAGCCGGCCCAGTAGACGGCCTACCGCCACCCCCAGGCCg                                                                               >  8:126938/1‑149 (MQ=255)
                                                                   cggAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCGTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAg                                                                           <  4:15366/149‑1 (MQ=255)
                                                                         tGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAg                                                                     <  2:164557/149‑1 (MQ=255)
                                                                            aGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCaa                                                                  <  7:93827/149‑1 (MQ=255)
                                                                             gCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCaa                                                                  >  7:160673/1‑148 (MQ=255)
                                                                                g‑ATCCAGTCACCAGACCAGGAGCTGCGCTTACACAGCGGCATACCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAg                                                                                                           >  1:140628/1‑104 (MQ=255)
                                                                                 gatggaGTCACGAGAGCAGGAGCTGCGCTTCCACAGCTGCATGCCGCAGGGTTTGGCCGCCTGATGGTCACCAGGCGGCGCAGGGCTTAACTGCGCCGTCCCAg                                                                                                           <  2:140628/99‑1 (MQ=255)
                                                                                    ccAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTg                                                           >  4:168077/1‑149 (MQ=255)
                                                                                      aGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCATTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGAAACAGTCTGc                                                          >  2:87555/1‑148 (MQ=255)
                                                                                      aGTCACCAGACCAGGAGCTGCGATTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGcc                                                         >  8:15899/1‑149 (MQ=255)
                                                                                       gTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGc                                                          >  2:291723/1‑147 (MQ=255)
                                                                                          accacaccacGCGCCGCGCCGCCACAGCGGCATCCCCCAGGCTTTGGCCGCCTGATGGGCACCACGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAggg                                                     <  4:286696/149‑1 (MQ=255)
                                                                                           ccagaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAgggg                                                    <  8:165630/149‑1 (MQ=255)
                                                                                             agaccagGAGCTGCGCTTCCTCAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTg                                                  <  7:126938/149‑1 (MQ=255)
                                                                                             agaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTg                                                  >  6:40963/1‑149 (MQ=255)
                                                                                             agaccagGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTg                                                  >  3:236346/1‑149 (MQ=255)
                                                                                                       cTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGTCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGACGCCGCCAGAAATAGCAACAGACTGCCAGCGCTGgcacatcgg                                         >  6:74673/1‑141 (MQ=255)
                                                                                                          cgcTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTa                                      <  2:167597/148‑1 (MQ=255)
                                                                                                             ttCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTcc                                  >  6:191426/1‑149 (MQ=255)
                                                                                                             ttCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTcc                                  >  6:44498/1‑149 (MQ=255)
                                                                                                               ccACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGt                                       >  1:136547/1‑142 (MQ=255)
                                                                                                               ccACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGt                                       <  2:136547/142‑1 (MQ=255)
                                                                                                                cacaGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTgg                               <  4:68767/149‑1 (MQ=255)
                                                                                                                cacaGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTgg                               >  3:68767/1‑149 (MQ=255)
                                                                                                                   aGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTg                            <  3:10099/149‑1 (MQ=255)
                                                                                                                     cGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGc                          >  1:146697/1‑149 (MQ=255)
                                                                                                                     cGGCATCCCGCAGGCTTTGGCCGCCTGAGGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGc                          <  2:115338/149‑1 (MQ=255)
                                                                                                                       gCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGAAGTAGTCCTGGGAGGcct                        >  7:239486/1‑147 (MQ=255)
                                                                                                                                aGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGTCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCg               >  2:67624/1‑149 (MQ=255)
                                                                                                                                    tttGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTgg                               >  6:313395/1‑129 (MQ=255)
                                                                                                                                    tttGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTgg                               <  5:313395/129‑1 (MQ=255)
                                                                                                                                    ctgcgccgccTGGTGCCCATCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTg           >  2:6533/5‑149 (MQ=255)
                                                                                                                                             cTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTtgct  <  8:160673/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CTGAAACTGGCGCACAACTACGGTGATGAAGCGGACATCGCTCGCTGCGTCGAGCTGGTACGCGCGGCGGAAATGGAGCTG‑ATCCAGTCACCAGACCAGGAGCTGCGCTTCCACAGCGGCATCCCGCAGGCTTTGGCCGCCTGATGGGCACCAGGCGGCGCAGGCCTTAACTGCGCCGGCCCAGTAGCAGGCCTACCGCCAGCCCGAGGCCGGCAGAAATAGCAACAGTCTGCCAGGGGTGGCCGCCGATGTAGTCCTGGGTGGCGTCCACCACCGGTTGCGCCTTGCT  >  NC_002947/2495266‑2495554

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