Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,581,222	A→G	N29D (AAC→GAC)	yhhY →	N‑acetyltransferase YhhY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,581,222	0	A	G	100.0%	117.7 / NA	34	N29D (AAC→GAC)	yhhY	N‑acetyltransferase YhhY
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/16);  total (18/16)

GAGACAATAATAATCATTCTCATTCGCACTTGTCCAACACTTTTTGCAAAAAAATGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGT  >  NC_000913/3581017‑3581468                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       gagaCAATAATAATCATTCTCATTCGCACTTGTCCAACACTTTTTGCAAAAAAATGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGc                                                                                                                                                                                                           >  2:254617/1‑251 (MQ=255)           taatCATTCTCATTCGCACTTGTCCAACACTTTTTGCAAAAAAATGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCg                                                                                                                                                                                                  >  2:81666/1‑250 (MQ=255)                         cGCACTTGTCCAACACTTTTAGCAAAAAAATGTCTTTGACTCGCCTTTGAAAGCAAATGATGTTGAAAGGGAAATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAAAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGCGGTGTATTGCGACACACTCCAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCACCGATCGGCCCGGcctc                                                                                                                                                                                    >  2:64111/1‑247 (MQ=255)                           cACTTGTCCAACACTTTTTGAAAAAAACCGCATTTGACTCGCATTTGACAGCCAATGATGTTGAAAGGGACCTTTACCTCAAAGAGGACACAGGCATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTCCGAGGCCATCAGGCAGATTCACGCCCAGACGGAGGTGTATTGCGACACCCTACAGGTGCCTCATCCTTCCGATACTATGTGGCAGGAACACCTCGACGAACCTCCCGGCATcc                                                                                                                                                                                   >  2:200028/1‑248 (MQ=255)                                gTCCAACACTTGTTGCAATAAAATGCATGTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATATGAAGACGCCTCGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCACCCGGAGGTGTAGGGCGACACACTACATGTGCCTCATCCGCCCGCTCAAATGTTGCATGAGCGAGTCGCTGATCTTCCCgggttctagcaac                                                                                                                                                                             <  2:137479/250‑12 (MQ=255)                                             tGCAAAAAAATGCATTTGCCTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGCGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTAtt                                                                                                                                                               >  1:281755/1‑251 (MQ=255)                                              gCAAAAAAATGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTAtt                                                                                                                                                               <  1:158816/250‑1 (MQ=255)                                                 aaaaacTGCATTTGACTCGCATTTGAAAGTCCATGATGTTGAAAGGGAAATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTTGCCTGTATTGATg                                                                                                                                                           >  2:171079/1‑251 (MQ=255)                                                 aaaaaaTGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATg                                                                                                                                                           <  1:252739/251‑1 (MQ=255)                                                   aaaaTGCATTTAACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACCTTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATgg                                                                                                                                                          >  1:208319/1‑250 (MQ=255)                                                      aTGCATTTGACTCGCATTTAAACGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATCTGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGAc                                                                                                                                                      >  2:276216/1‑251 (MQ=255)                                                                         gAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCAcc                                                                                                                                    >  2:98917/1‑250 (MQ=255)                                                                                 atgatgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGGGGGCCATCTCACCATTGACGTg                                                                                                                           >  2:221810/1‑251 (MQ=255)                                                                                             gggACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCa                                                                                                               <  1:22218/251‑1 (MQ=255)                                                                                              ggACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAc                                                                                                              <  1:306015/251‑1 (MQ=255)                                                                                                            aaGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGtt                                                                                                >  1:155996/1‑251 (MQ=255)                                                                                                                   cAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTCCGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGAtt                                                                                         >  1:186938/1‑251 (MQ=255)                                                                                                                         aaTGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTa                                                                                   <  1:76375/251‑1 (MQ=255)                                                                                                                           tgagtgagATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGCTTTTGGTATc                                                                                 >  2:293363/1‑251 (MQ=255)                                                                                                                            gagtgagATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCt                                                                                <  1:266519/251‑1 (MQ=255)                                                                                                                                gaaaTAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCAACAGGCAGATTCACGCCCGGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCAGATGTGGCAGGAGGGCCTCGCCGATCGGCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCTATCTCACCATTGAGGTGCAACAGCGCGCACGCCGCAGTCATGTTGCCGATTTTGGTATCTgtgt                                                                            <  2:69614/248‑1 (MQ=255)                                                                                                                                  gaTAGTAATACGCCACGCAGAAACACGGGATTTCGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCg                                                                          <  1:221810/251‑1 (MQ=255)                                                                                                                                                 cGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGgaaga                                                           <  1:293363/251‑1 (MQ=255)                                                                                                                                                           cGGGATTACGAGGTCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGt                                                 <  1:276216/251‑1 (MQ=255)                                                                                                                                                            gggATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTc                                                >  2:256149/1‑251 (MQ=255)                                                                                                                                                               aTTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGcc                                             >  1:12446/1‑251 (MQ=255)                                                                                                                                                                  aCGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAgc                                          <  2:85522/251‑1 (MQ=255)                                                                                                                                                                    gAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTCTTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAgcgc                                        >  1:307919/1‑251 (MQ=255)                                                                                                                                                                                cagtttCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCACCCTTCCGATCATATGTGGCAGTAGCTACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                            <  2:12446/247‑1 (MQ=255)                                                                                                                                                                                cAGATTCACGCCCAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                            <  1:245637/251‑1 (MQ=255)                                                                                                                                                                                          cccAGCCGGAGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGaa                  <  1:255466/251‑1 (MQ=255)                                                                                                                                                                                                 ggaggaGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGTGATGATTGAAATGTGCg           >  2:141694/1‑251 (MQ=255)                                                                                                                                                                                                   aGGTGTATTGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCgtc         <  2:275520/251‑3 (MQ=255)                                                                                                                                                                                                          ttGCGACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt  >  2:254375/1‑251 (MQ=255)                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       GAGACAATAATAATCATTCTCATTCGCACTTGTCCAACACTTTTTGCAAAAAAATGCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGT  >  NC_000913/3581017‑3581468

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