Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,238	G→A	intergenic (‑173/‑85)	yjjI ← / → deoC	DUF3029 domain‑containing protein YjjI/deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,238	0	G	A	100.0%	105.6 / NA	34	intergenic (‑173/‑85)	yjjI/deoC	DUF3029 domain‑containing protein YjjI/deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/15);  total (19/15)

AAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCG  >  NC_000913/4617000‑4617484                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       aaGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAgtgt                                                                                                                                                                                                                                             <  1:249630/250‑1 (MQ=255) aaGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAgtgt                                                                                                                                                                                                                                             <  1:290702/250‑1 (MQ=255)   gCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGtt                                                                                                                                                                                                                                          >  1:153668/1‑251 (MQ=255)                                              aTTTTCATCAGAAGTGGGAATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGACATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATATTAACAAACTAGAAAGGGGAAtt                                                                                                                                                                                               >  1:169379/1‑251 (MQ=255)                                                tttCATCAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGAATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAAttt                                                                                                                                                                                              <  1:52836/250‑1 (MQ=255)                                                  tcatcaGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCAATTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAGTACTAACAAACTCGCAAGGTGAATTTTAt                                                                                                                                                                                           >  1:196023/1‑251 (MQ=255)                                                   catcagaAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTAtt                                                                                                                                                                                          >  1:412100/5‑251 (MQ=255)                                                      cagaAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGc                                                                                                                                                                                       >  2:378491/2‑251 (MQ=255)                                                               gCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCa                                                                                                                                                                              <  1:11637/251‑1 (MQ=255)                                                                 aTCTTCTTTTCCTTTTATGCCGCAGGAGATGCGCCATTGTAAGAAGTTTTGTGTTGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAAATTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGAAGATGTTAGAATACTAACAAACTAGCAAGGTGAATTTTATTGGCGACAAGcatgg                                                                                                                                                                            >  2:330186/1‑247 (MQ=255)                                                                                         ggTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAaag                                                                                                                                                    <  1:50844/251‑1 (MQ=255)                                                                                          gTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAaagc                                                                                                                                                   <  1:126843/251‑1 (MQ=255)                                                                                                   cATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCt                                                                                                                                          >  1:205971/1‑251 (MQ=255)                                                                                                   cATTGTAAGAAGTTTCGTGATGTTAACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGGGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCt                                                                                                                                          >  1:417018/1‑251 (MQ=255)                                                                                                                     gATGTTCACTTTGATCCTGATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGAt                                                                                                                        >  1:500770/1‑251 (MQ=255)                                                                                                                           cACTTTGATCCTGATGAGTTTGCCACCCCTTACGCATTCATTTGAAAGTGAATTATTTTAACCTGATCGCATTACAGTGATGCTAACTTGTAAGTAGATTTCCTTAATTGGGATGTATATCGAAGTGTGTTGCGGAGTCGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCGGCGTGCACGGAAATTGATGTACCt                                                                                                                  <  2:405606/251‑1 (MQ=255)                                                                                                                               ttGATCCTGGTGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGac                                                                                                               <  2:334153/250‑1 (MQ=255)                                                                                                                                       gATGAGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGa                                                                                                       >  2:9279/1‑250 (MQ=255)                                                                                                                                                  caccacTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAACCTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGACTTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGCCCACCCTGAATgccgacgac                                                                                            >  2:225161/1‑250 (MQ=255)                                                                                                                                                  caccacTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCt                                                                                                                                          >  1:300290/1‑204 (MQ=255)                                                                                                                                                  caccacTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCt                                                                                                                                          <  2:300290/204‑1 (MQ=255)                                                                                                                                                         gACGCATTCATTTGAAAGGGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGTCGACTCcgacg                                                                                     <  2:324659/250‑1 (MQ=255)                                                                                                                                                               ttcattTGAAAGCGAATTATTTGAAACAGATCGCATTACAGTGATGCAAACTTGTGAGTAGATTTCCTTAATTGTGATGGATATAGAAGTGTGTTGCGGCGTAGATGTTAAAATACTAACAAAATCGAAAGGTGAATTTTATTGGGGACAAGCAAGGAGAATGAAATGAATGATATGACAGCAAGCAGCCTGCGTGCACTGAAATTGATGGTCCTGACCAACATGAACGCTGACAACACCGACGAGAAAGt                                                                              >  1:20794/1‑251 (MQ=255)                                                                                                                                                                 cattTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGtaa                                                                            <  2:237734/251‑3 (MQ=255)                                                                                                                                                                              attattTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGtca                                                               <  2:284023/251‑1 (MQ=255)                                                                                                                                                                                                   tACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGCCGAGAAAGTAATTGCTCTGTGTCCTCAGGCCAAACCCCCGGTCgg                                          >  2:152945/1‑251 (MQ=255)                                                                                                                                                                                                       gTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCCCCCTGAATGACGACGACACCGACGATAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCaag                                      <  1:174703/251‑2 (MQ=255)                                                                                                                                                                                                          aTGCAAACTTGTAAGTAGATTTCCTTAATTTTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATAcc                                   <  2:385390/251‑1 (MQ=255)                                                                                                                                                                                                             cAAACTTGTAAGTAGATTTCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGTAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACAACCCTGATTGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTAATAAGGCCAAAACCCCGGTCGGCAATAccgcc                                >  2:281002/1‑251 (MQ=255)                                                                                                                                                                                                                              ttCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATATATcc               >  2:455910/1‑251 (MQ=255)                                                                                                                                                                                                                               tCCTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGAACTGACCACCCTGAATGACTACGACACCGACGCGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTTTATCTATCCt              >  1:435031/1‑251 (MQ=255)                                                                                                                                                                                                                                ccTTAATTGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAAAATACTAACAAACTCGAAGGGTGATTTTTATTGGCGACAAGCCAGGAGAATGAAATGATGGATCGGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGCCCGCACCACCCTGAATGACGACGACCCCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCAAAACCCCCGGTCGTCAATACCGCCGCTATCTGTATCTATCCt              <  1:239559/250‑1 (MQ=255)                                                                                                                                                                                                                                      ttGTGATGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCttt        >  1:70633/1‑250 (MQ=255)                                                                                                                                                                                                                                           aTGTATATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCg  >  1:191879/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       AAGCACTGGGCTGGTGACAATTTGCTGGCAACGTTGTTGCAGTGCATTTTCATGAGAAGTGGGCATCTTCTTTTCCTTTTATGCCGAAGGTGATGCGCCATTGTAAGAAGTTTCGTGATGTTCACTTTGATCCTGATGCGTTTGCCACCACTGACGCATTCATTTGAAAGTGAATTATTTGAACCAGATCGCATTACAGTGATGCAAACTTGTAAGTAGATTTCCTTAATTGTGATGTGTATCGAAGTGTGTTGCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAATTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCACCCTGAATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGGTCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCG  >  NC_000913/4617000‑4617484

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: ‑