Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,012,553 Δ1 bp intergenic (+160/+61) yqeA → / ← yqeB putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,5520G.96.8% 253.3 / 1.6 63intergenic (+159/+62)yqeA/yqeBputative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (30/31);  total (31/32)
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.91e-01

TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGT  >  NC_000913/3012419‑3012683
                                                                                                                                     |                                                                                                                                   
tCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAAtt                                                                                                                               >  2:1270284/1‑139 (MQ=255)
   ggCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTgg                                                                                                                            <  2:630696/139‑1 (MQ=255)
             tGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:1326426/139‑1 (MQ=255)
              gTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcatc                                                                                                                 >  1:469733/1‑139 (MQ=255)
                   gtaagtGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGtct                                                                                                            <  2:469733/139‑1 (MQ=255)
                           gcTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:99971/139‑1 (MQ=255)
                                  tgATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTa                                                                                             >  2:969397/1‑139 (MQ=255)
                                     tGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATa                                                                                          <  2:318240/139‑1 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGttat                                                                                                                  >  2:87776/1‑86 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  <  1:1040012/86‑1 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtcat                                                                                                                  >  2:1040012/1‑86 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACtt                                                             >  1:654104/1‑139 (MQ=255)
                                                                  tCCTACAAAACATCGCAAATAAAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGTAATTTGGCGTTTGttat                                                                                                                  <  1:87776/86‑1 (MQ=21)
                                                                   ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtc                                                                                                                    >  2:548532/1‑83 (MQ=255)
                                                                   ccTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGtc                                                                                                                    <  1:548532/83‑1 (MQ=255)
                                                                    cTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTcc                                                           >  2:74672/1‑139 (MQ=255)
                                                                      aCAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAg                                                         <  1:969397/139‑1 (MQ=255)
                                                                       cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                                                                                 <  1:1201018/114‑1 (MQ=255)
                                                                       cAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGc                                                                                 >  2:1201018/1‑114 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                     >  2:550939/1‑58 (MQ=255)
                                                                            cATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                                                     <  1:550939/58‑1 (MQ=255)
                                                                              tCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtct                                                  >  2:1493261/1‑139 (MQ=255)
                                                                               cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                                       >  2:1453377/1‑84 (MQ=255)
                                                                               cGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGc                                                                                                       <  1:1453377/84‑1 (MQ=255)
                                                                                           aGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCAc                                    <  1:638286/139‑1 (MQ=255)
                                                                                            gATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                   <  1:74672/139‑1 (MQ=255)
                                                                                              ttGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  2:1532465/1‑72 (MQ=255)
                                                                                              ttGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:1532465/72‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  1:1195255/1‑71 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  2:1195255/71‑1 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAttt                                                                            >  1:977982/1‑95 (MQ=255)
                                                                                               tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCAttt                                                                            <  2:977982/95‑1 (MQ=255)
                                                                                                gCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCt                               <  2:654104/139‑1 (MQ=255)
                                                                                                  aaTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGa                             >  1:1113075/1‑139 (MQ=255)
                                                                                                   aTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGa                             <  1:1493261/139‑1 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  2:495860/64‑1 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  1:495860/1‑64 (MQ=255)
                                                                                                      cACTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         >  1:1350993/1‑139 (MQ=255)
                                                                                                       aCTTGTAGGCCTGATAAGTGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                                     >  1:461074/1‑126 (MQ=255)
                                                                                                       aCTTGTAGGCCTGATAAGTGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCa                                     <  2:461074/126‑1 (MQ=255)
                                                                                                       aCTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt                        >  1:522331/1‑139 (MQ=255)
                                                                                                       aCTTGTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt                        <  2:1350993/139‑1 (MQ=255)
                                                                                                           gTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAttt                                                                      >  2:1057811/1‑89 (MQ=255)
                                                                                                           gTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAttt                                                                      <  1:1057811/89‑1 (MQ=255)
                                                                                                           gTAGGCCTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc                    <  1:286449/139‑1 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTGTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAcatcat               >  2:1084319/1‑139 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGa                             <  1:670897/125‑1 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGa                             >  2:670897/1‑125 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         >  1:1511812/1‑129 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         <  2:1511812/129‑1 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         >  1:1144948/1‑129 (MQ=255)
                                                                                                                ccTGATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCa                         <  2:1144948/129‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                               <  2:263715/56‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgggtggg                                                                                               >  1:263715/1‑56 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    <  1:760833/51‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCgg                                                                                                    >  2:760833/1‑51 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGTGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAc                    <  2:1044760/131‑1 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTaa            >  2:490258/1‑139 (MQ=255)
                                                                                                                   gATAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGAGGCATCAAc                    >  1:1044760/1‑131 (MQ=255)
                                                                                                                     tAAGCGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACg          >  1:1372129/1‑139 (MQ=255)
                                                                                                                         cGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg                                                                   <  2:605760/78‑1 (MQ=255)
                                                                                                                         cGTAGCGCATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAg                                                                   >  1:605760/1‑78 (MQ=255)
                                                                                                                             gcgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                                                                                  >  2:1018528/1‑59 (MQ=255)
                                                                                                                             gcgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCgg                                                                                  <  1:1018528/59‑1 (MQ=255)
                                                                                                                             gcgcATCAGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGt  <  2:1113075/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
TCCGGCCAAAGATTGTCAAGTAAGTGCGCTTTGTTGATGCCGGATGCGGCGCGAGCGCCTTATCCGTCCTACAAAACATCGCAAATTCAATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGT  >  NC_000913/3012419‑3012683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: