Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,299,959	G→A	intergenic (‑44/‑7)	yhbP ← / → yhbQ	putative FMN binding protein YhbP/DNA damage response nuclease YhbQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,299,959	0	G	A	100.0%	194.9 / NA	55	intergenic (‑44/‑7)	yhbP/yhbQ	putative FMN binding protein YhbP/DNA damage response nuclease YhbQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/23);  total (32/23)

CACCAGAGTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCG  >  NC_000913/3299830‑3300087                                                                                                                                  |                                                                                                                                 cACCAGAGTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGa                                                                                                                        >  1:652322/1‑140 (MQ=255)  aCCAGAGTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGac                                                                                                                       >  2:21644/1‑140 (MQ=255)     agagTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACAcc                                                                                                                    <  2:421862/140‑1 (MQ=255)     agagTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACAcc                                                                                                                    >  1:421862/1‑140 (MQ=255)     agagTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACAcc                                                                                                                    <  1:386022/140‑1 (MQ=255)      gagTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGa                                                                                                                              >  1:194061/1‑129 (MQ=255)      gagTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGa                                                                                                                              <  2:194061/129‑1 (MQ=255)         ttCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTgg                                                                                                                >  1:1149565/1‑140 (MQ=255)            cccTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGttt                                                                                                             >  2:760489/1‑140 (MQ=255)            cccTTCCTGCTGCACACACCAAGTGACAACATGTTGTTGTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGttt                                                                                                             >  2:233475/1‑140 (MQ=255)              cTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCt                                                                                                           >  1:91616/1‑140 (MQ=255)                  ctgctgCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTAc                                                                                                       <  2:1149565/140‑1 (MQ=255)                  ctgctgCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTAc                                                                                                       <  1:625227/140‑1 (MQ=255)                   tgctgcACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGAGGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACt                                                                                                      >  1:180244/1‑140 (MQ=255)                         acacacCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATcc                                                                                                >  1:110947/1‑140 (MQ=255)                          cacacCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCg                                                                                               >  1:947917/1‑140 (MQ=255)                             acCAAGTGACAACATGTTGTTTTGCTAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTAc                                                                                            <  2:1049753/140‑1 (MQ=255)                                  gTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTAccgccg                                                                                       >  2:13242/1‑140 (MQ=255)                                       aacaTGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGt                                                                                              <  2:1087934/128‑1 (MQ=255)                                       aacaTGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGt                                                                                              >  1:1087934/1‑128 (MQ=255)                                       aacaTGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAAt                                                                                  >  1:222535/1‑140 (MQ=255)                                            gttgttTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCt                                                                             <  2:222535/140‑1 (MQ=255)                                                   tGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATAcc                                                                      <  2:91616/140‑1 (MQ=255)                                                      cagccaACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCggg                                                                   <  2:652322/140‑1 (MQ=255)                                                        gccaACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGAt                                                                 <  1:13242/140‑1 (MQ=255)                                                             cGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATcacca                                                            <  2:947917/140‑1 (MQ=255)                                                                   aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATg                                                      >  2:991937/1‑140 (MQ=255)                                                                   aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATg                                                      >  2:941051/1‑140 (MQ=255)                                                                   aTGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATg                                                      >  1:116167/1‑140 (MQ=255)                                                                       cGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGa                                                  <  1:941051/140‑1 (MQ=255)                                                                          tGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                               <  2:116167/140‑1 (MQ=255)                                                                              tgtTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                               <  2:1001559/136‑1 (MQ=255)                                                                              tgtTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                               >  1:1001559/1‑136 (MQ=255)                                                                                tttCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTa                                         >  2:925007/1‑140 (MQ=255)                                                                                      ttGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                               <  1:876319/128‑1 (MQ=255)                                                                                      ttGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAAcg                                               >  2:876319/1‑128 (MQ=255)                                                                                      ttGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATcagca                                   >  2:497716/1‑140 (MQ=255)                                                                                           gTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCAcc                              <  1:760489/140‑1 (MQ=255)                                                                                            ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCa                             >  1:619484/1‑140 (MQ=255)                                                                                            ttCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCa                             >  2:907945/1‑140 (MQ=255)                                                                                                    tataCTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCa                     <  1:991937/140‑1 (MQ=255)                                                                                                      taCTGTGGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCaaa                   >  2:211983/1‑140 (MQ=255)                                                                                                           tgGGCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGc              >  1:147470/1‑140 (MQ=255)                                                                                                             ggCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCaa                            <  1:911610/124‑1 (MQ=255)                                                                                                             ggCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCaa                            >  2:911610/1‑124 (MQ=255)                                                                                                             ggCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCaaa                   >  1:93783/1‑133 (MQ=255)                                                                                                             ggCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCaaa                   <  2:93783/133‑1 (MQ=255)                                                                                                             ggCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGa            >  2:864744/1‑140 (MQ=255)                                                                                                             ggCCTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGa            >  2:815836/1‑140 (MQ=255)                                                                                                                cTGACCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAg         >  2:271306/1‑140 (MQ=255)                                                                                                                   aCCTTAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCAc      >  2:707007/1‑140 (MQ=255)                                                                                                                       tAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg  <  1:864744/140‑1 (MQ=255)                                                                                                                       tAACATGCGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg  <  1:211983/140‑1 (MQ=255)                                                                                                                              cGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg  >  2:850573/1‑133 (MQ=255)                                                                                                                              cGACACAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTgcg  <  1:850573/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                                 CACCAGAGTTCCCCTTCCTGCTGCACACACCAAGTGACAACATGTTGTTTTGCCAGCCAACGGCTGATGGCGATGAGTGTTTCCATTGCTGTTCTCTTTTATACTGTGGGCCTGACCTTAACATGCGACGCAGACGATGACACCCTGGTTTCTTTACTTGATCCGTACCGCCGACAATAAGCTTTATACCGGGATCACCACGGATGTCGAACGCCGCTATCAGCAGCACCAAAGCGGCAAAGGGGCGAAAGCACTGCG  >  NC_000913/3299830‑3300087

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