Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,204,395	T→A	intergenic (+106/‑201)	yehK → / → yehL	uncharacterized protein YehK/putative AAA(+) MoxR family ATPase YehL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,204,395	0	T	A	100.0%	104.6 / NA	35	intergenic (+106/‑201)	yehK/yehL	uncharacterized protein YehK/putative AAA(+) MoxR family ATPase YehL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (19/16);  total (19/16)

CCCGCTAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTG  >  NC_000913/2204151‑2204635                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 cccGCTAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAatcgttc                                                                                                                                                                                                                                            <  1:142902/251‑1 (MQ=255)          aaaTGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAaca                                                                                                                                                                                                                                   <  1:169740/251‑1 (MQ=255)          aaaTGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTAAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGTTTATAAATCGTTCGTTGGAac                                                                                                                                                                                                                                    >  2:420943/1‑250 (MQ=255)                      gATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGcac                                                                                                                                                                                                                        >  1:228029/1‑250 (MQ=255)                                          gtaatgatgaGCGGAAATTACGAGGTGAATCCCCTGGAACTGTAATAGCAATTGTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAACATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGTGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTg                                                                                                                                                                                                    <  2:13449/246‑1 (MQ=255)                                                 ggaGCGGAAATTACGAGGTGAATCAGCAGAAACTGTAATAGCAATTTTTAAGCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAAATCAAATTTGGCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAAGCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGt                                                                                                                                                                                            <  2:230522/250‑1 (MQ=255)                                                   aGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGAAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAAAACGATGTTGCAAAGGCTGGGGGGTTAGGCCTGAAAATAgaa                                                                                                                                                                                           >  2:4444/1‑248 (MQ=255)                                                       gAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGtaaataa                                                                                                                                                                                      <  1:65604/251‑1 (MQ=255)                                                             aCGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGAACCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAAc                                                                                                                                                                                <  1:317859/251‑1 (MQ=255)                                                                aGGTAAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAAATACAGAATAAATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGAGTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGGTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTg                                                                                                                                                                             <  2:184335/251‑1 (MQ=255)                                                                    gAATCACCTGGCACTGTAATAGCAATTTTTAATCGTGCACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGATCCGATCCGGCGCAACGTCCCACTGGCCTGGATTATAAATCTAATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACAAGATGTTGCACAGGCTGGGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTc                                                                                                                                                                         >  2:171573/1‑251 (MQ=255)                                                                           cTGAAACTGTAATAGCAATTTTTAATCGAGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTATTCCGATCCGGCGCAACGTCACACTGGCCTGGATTATAAATCTCATTATCTTACTGGCAGCGGGATCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGCTGATAACCAGTCTGTGGGTTTAGGCCTGTACGTCGTAAATAAAAGAACCTGCCTCACcaga                                                                                                                                                                   >  2:200550/1‑249 (MQ=255)                                                                                aCTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGttt                                                                                                                                                             <  1:138243/251‑1 (MQ=255)                                                                                       tAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTtttatt                                                                                                                                                      >  2:424246/1‑251 (MQ=255)                                                                                                           cTACAGAATATATTCAACCAAGAAATAGACTAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAgg                                                                                                                                                                                                        >  1:161829/1‑181 (MQ=255)                                                                                                           cTACAGAATATATCCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGGCCCAATGGCCTGGATTATAAATCTCATTATCGTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAgg                                                                                                                                                                                                        <  2:161829/181‑1 (MQ=255)                                                                                                                        tCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCTCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCt                                                                                                                      <  1:246509/250‑1 (MQ=255)                                                                                                                            ccAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGAGCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAaac                                                                                                                 <  1:211605/251‑1 (MQ=255)                                                                                                                                    tAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATg                                                                                                         >  1:421872/1‑251 (MQ=255)                                                                                                                                      gAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGCAAATAGTAAATAACAGAACCTGCCTAACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATgcg                                                                                                       >  1:262641/1‑251 (MQ=255)                                                                                                                                           atatACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCaaa                                                                                                  <  1:59605/251‑1 (MQ=255)                                                                                                                                             ataCTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAAc                                                                                                 <  1:425954/250‑1 (MQ=255)                                                                                                                                                       taatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGttt                                                                                                                                                             >  1:489033/1‑180 (MQ=255)                                                                                                                                                       taatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTgcgc                                                                                      >  1:363593/1‑251 (MQ=255)                                                                                                                                                       taatTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGGTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGttt                                                                                                                                                             <  2:489033/180‑1 (MQ=255)                                                                                                                                                                   gTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGt                                                                          <  1:424246/251‑1 (MQ=255)                                                                                                                                                                       gATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTAt                                                                      >  1:106332/1‑251 (MQ=255)                                                                                                                                                                        aTCCGGCGCAAAGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAAAACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGCTCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCACCAAACCTTTATTGCGCGTAACATATCGTTTAtt                                                                     >  2:333051/1‑251 (MQ=255)                                                                                                                                                                           cGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTc                                                                   >  2:441282/1‑250 (MQ=255)                                                                                                                                                                                    gTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAACATATCGTTTATTTCATTAATACAt                                                         >  1:121503/1‑251 (MQ=255)                                                                                                                                                                                                    ttATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAatat                                          >  1:133814/1‑250 (MQ=255)                                                                                                                                                                                                        aaaTCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTc                                     <  1:4712/251‑1 (MQ=255)                                                                                                                                                                                                                   tCTTAATTGCAACGGGGTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGaat                           >  1:175334/1‑250 (MQ=255)                                                                                                                                                                                                                                ggggTCCAGCCGTGGTTATAAATCGTTCGTTGTAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCTCCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAaaca                                                                                                                <  2:281833/152‑1 (MQ=255)                                                                                                                                                                                                                                ggggTCCAGCCGTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAaaca                                                                                                                >  1:281833/1‑152 (MQ=255)                                                                                                                                                                                                                                           gTGGTTATAAATCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGCACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTAAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCctgctg  >  2:79244/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 CCCGCTAATAAATGATGCGTTGATTTATATAAGAGATGAGTGTATTGATGAGCGGAAATTACGAGGTGAATCACCTGAAACTGTAATAGCAATTTTTAATCGTGAACTACAGAATATATTCAACCAAGAAATAGAATAATATACTCTAAATAATTCAAATTGGTCCGATCCGGCGCAACGTCCCAATGGCCTGGATTATAAATCTCATTATCTTAATTGCAACGGGGTCCAGCCGTGGTTATAATTCGTTCGTTGGAACACGATGTTGCACAGGCTGTGGTGTTAGGCCTGAAAATAGTAAATAAAAGAACCTGCCTCACCAGCAGGTTTTTTTATTTACTGTGATCTGCTTTCCAGATATTTTTCGCTCAAACAACTAATGCGCCAAACATTTATTGCGCGTAAAATATCGTTTATTTCATTAATACATTTCAGGGATGAATATATGTCACCGCAGAATAATCATCTGCAGCGTCCGCCTGCTG  >  NC_000913/2204151‑2204635

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