Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	254,249	A→T	intergenic (+1088/+10)	yafP → / ← pepD	putative N‑acetyltransferase YafP/peptidase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	254,249	0	A	T	100.0%	112.2 / NA	33	intergenic (+1088/+10)	yafP/pepD	putative N‑acetyltransferase YafP/peptidase D
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/16);  total (17/16)

GTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACC  >  NC_000913/254011‑254499                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                           gTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGACTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAGCGCAATCCGGCAGCATATCAAATGAtt                                                                                                                                                                                                                                                 >  2:390285/1‑250 (MQ=255)         tcagtcagAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTGCTTCGcc                                                                                                                                                                                                                                         <  2:220495/250‑1 (MQ=255)            gtcagAGCGTATGCAGCATACTAACAAGCGGCTGGCACGATTGCTGATGGCTTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCTGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCAGCTACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAATTGCCGGAAGGCGATTAAACGCCAGCCGGCTGCATATCAAATAATTCCTGCGCCggga                                                                                                                                                                                                                                     <  2:381309/251‑3 (MQ=255)             tcagAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAt                                                                                                                                                                                                                                    >  1:262440/1‑251 (MQ=255)                   cGTAGTCAGCATGCTAACAAGCGGCTGGCCCGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCTCGGCGCTGAAATCGCAGCGGCGAATGTTCCATAACCAGATTGAACGTGGCAAACCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttcttt                                                                                                                                                                                                                              >  1:172298/1‑251 (MQ=255)                              tGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAACTAATTACTACGCCGGAATTTCTTTCAGCAGTtc                                                                                                                                                                                                                     >  1:301782/1‑249 (MQ=255)                                  aaCAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtca                                                                                                                                                                                                                <  2:439902/250‑1 (MQ=255)                                                    ttGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGACAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc                                                                                                                                                                                             >  1:70803/1‑251 (MQ=255)                                                             gCCTGGAAGATGGAGCAGCAACAACAGGAAAAAAGCGCGGCGCTGCAATAGCCGCGGCGAATGTTCCATCACCCGACTGAACGTGGCAACCCGCGACGGACCTTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCCGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc                                                                                                                                                                                    >  2:238346/1‑251 (MQ=255)                                                                             aacagtaactgGAACATAGCGCGGCGCTGACATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAAAACGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATCATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCCGTAATGACCTACGCTTTCGATGTGAACTTGCTc                                                                                                                                                                     >  2:242699/1‑250 (MQ=255)                                                                              acagcagcagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATTGCTTTTATCGAATGATAATGAAGGATGAAACTGCCGGAATGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTGCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCAt                                                                                                                                                                   <  1:490006/251‑1 (MQ=255)                                                                                              aGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGAcc                                                                                                                                                   <  1:154187/251‑1 (MQ=255)                                                                                                  cggcgCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTg                                                                                                                                               <  1:127584/251‑1 (MQ=255)                                                                                                              cGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAAACCGCGACGGACATTTACAGGGATGGCTTTTATAGACGGATAATGCAGGATGAAACAGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGcacg                                                                                                                                   >  1:524934/1‑248 (MQ=255)                                                                                                                cAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTg                                                                                                                                               >  1:142951/1‑237 (MQ=255)                                                                                                                cAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTg                                                                                                                                               <  2:142951/237‑1 (MQ=255)                                                                                                                cAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGAt                                                                                                                                 <  2:262440/251‑1 (MQ=255)                                                                                                                       gAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAcc                                                                                                                          >  1:11837/1‑251 (MQ=255)                                                                                                                        aaTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTCCTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTGCGATGTGAACTTGCTCATCGGGAGAGTGTGGACCGGTGAGAGTTGGCCCGATAGAAACCa                                                                                                                         <  2:329116/251‑1 (MQ=255)                                                                                                                         aTGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGTCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                                                                                                        <  1:176274/251‑1 (MQ=255)                                                                                                                                aTCACCAGATTGAAAGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCAtt                                                                                                                 >  2:88603/1‑251 (MQ=255)                                                                                                                                  cACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTc                                                                                                               <  1:88603/251‑1 (MQ=255)                                                                                                                                      aGATTGACCGTGGCAACCGGCGGCGGACTTATCTAGGGATGGCTTTTATCGAAGGATACTGATGGATGCAACTGCCGGAAGGCGATTAACCGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTATGCAGTGCCCAGTCATGACATACGCTTTCGCTGTGCACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAACATGTCCCTTTCCgg                                                                                                            >  2:423009/1‑250 (MQ=255)                                                                                                                                           gAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACgg                                                                                                      >  1:208101/1‑251 (MQ=255)                                                                                                                                                  gCAACCCGCGCAGGCCATTTACAGGGATGGCTTTTATCGCAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCAATCAGGCAGCATATCACATAATTAATTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTAAGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGCGATAGTTGGACCGATAGCAACCATGTCCATTTCCGGTTACGGTTTTTTg                                                                                               >  2:358324/1‑251 (MQ=255)                                                                                                                                                     aCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTTTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAAc                                                                                            <  1:280173/251‑1 (MQ=255)                                                                                                                                                       ccGCGAAGGTCATCTACCGGCATGGCTTTTATCGCAGGATCATGCACGCCGAAACTGCCGGCAGCCGATTCAACCCCAACCGGCAGAATATCAAATAATTACTTCGCCGGCATTTCTTTCAGCAGTCCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGAGATAGTTGGCCCGATAGAACCCATGTCCATTTCCGGATACGGTTTTTTGAAccg                                                                                          >  2:487759/1‑249 (MQ=255)                                                                                                                                                           gACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACca                                                                                      <  1:64062/251‑1 (MQ=255)                                                                                                                                                                                    tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCGGCGTGGATAATc                                                              <  1:390285/250‑1 (MQ=255)                                                                                                                                                                                         aaGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATg                                                        >  1:424665/1‑251 (MQ=255)                                                                                                                                                                                                                  ggAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCt                               >  1:258574/1‑251 (MQ=255)                                                                                                                                                                                                                   gAAGTCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACTGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTg                              <  2:98041/251‑1 (MQ=255)                                                                                                                                                                                                                             aaaCGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTGTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTc                     <  2:208101/250‑1 (MQ=255)                                                                                                                                                                                                                                               aatcaaatAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCGGTAATGACCTACGGTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGATCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATGCCAGGCCCTCGTGGATAATCTGGATGTACGGCGTCTTGTTGGACGAGCGCTGAAAGGTTGGACGTACCGGAGGTATCAcc  <  2:59131/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                           GTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACC  >  NC_000913/254011‑254499

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