Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,078,954	A→G	intergenic (+23/+78)	yoeF → / ← tmaR	putative uncharacterized protein YoeF/PTS system regulator TmaR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,078,954	0	A	G	100.0%	141.5 / NA	37	intergenic (+23/+78)	yoeF/tmaR	putative uncharacterized protein YoeF/PTS system regulator TmaR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/16);  total (21/16)

GGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCT  >  NC_000913/2078718‑2079185                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        ggttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGAt                                                                                                                                                                                                                           <  2:266884‑M2/39‑1 (MQ=255)     aacgtccgacaaaatgcgtcattgacagcggtttcggtggattgacacatttcatacagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATacac                                                                                                                                                                                                                       <  2:294391‑M2/43‑1 (MQ=255)               aaaatgcgtcattgacagaggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAAt                                                                                                                                                                                                             >  2:415986‑M2/199‑251 (MQ=255)                                     ttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGcc                                                                                                                                                                                        <  2:60933‑M2/74‑1 (MQ=255)                                              ttgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCt                                                                                                                                                                              >  1:209088‑M2/168‑251 (MQ=255)                                               tgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGCTGATGACGGTGATACACAGTTAAGAATCAATGAGCACTGAGGGTTGCCGGGCAACCCTc                                                                                                                                                                             >  2:481786‑M2/167‑251 (MQ=255)                                                      tttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggcatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAAt                                                                                                                                                                      >  2:404491‑M2/160‑251 (MQ=255)                                                         cattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATaaa                                                                                                                                                                   <  1:290473‑M2/95‑1 (MQ=255)                                                               gataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGa                                                                                                                                                             <  1:303642‑M2/101‑1 (MQ=255)                                                                  aatgaattaatgctactgccggaacaacaagcaaacgggcattgtgttctgaaatccaggctattgatgcaacgtccgcgacagctctgcgctgtaaaagaacttgtaaaacaacgaaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGAAACCCTCATAGAATAAATCGGGAAtt                                                                                                                                                          >  1:412464‑M2/148‑251 (MQ=255)                                                                         taatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgc                                                                                                                                                   >  1:284750‑M2/141‑251 (MQ=255)                                                                             gctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTActtcgcttcg                                                                                                                                               >  2:154908‑M2/137‑251 (MQ=255)                                                                                actgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGcc                                                                                                                                             <  1:144058‑M2/117‑1 (MQ=255)                                                                                actgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacgtgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGcc                                                                                                                                             <  2:64453‑M2/117‑1 (MQ=255)                                                                                     cggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGttttt                                                                                                                                       >  1:366484‑M2/129‑251 (MQ=255)                                                                                         acaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAANCGGGAATTACTTCGCTTCGCCGTTTTTCAtt                                                                                                                                   >  2:231979‑M2/125‑251 (MQ=255)                                                                                           aaccagcaaacgggcattgtgttctgaaatcnaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTc                                                                                                                                 <  2:238212‑M2/129‑1 (MQ=255)                                                                                           aaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTc                                                                                                                                 >  2:108518‑M2/123‑251 (MQ=255)                                                                                               agcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatga                                                                                                                                                                                                       <  2:157291‑M2/59‑1 (MQ=255)                                                                                                                     aaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCACATAGCTTTCAGCTTTTTgg                                                                                                          >  1:307473‑M2/97‑248 (MQ=255)                                                                                                                       atccaggctattgattcaacttcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATa                                                                                                                                                                                                                          <  2:213688‑M2/40‑1 (MQ=255)                                                                                                                       atccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATa                                                                                                                                                                                                                          >  1:213688‑M2/95‑134 (MQ=255)                                                                                                                                      tcaacgtcagcgacagatttgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTAct                                                                                                                                                       >  2:416825‑M2/80‑186 (MQ=255)                                                                                                                                      tcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTAct                                                                                                                                                       <  1:416825‑M2/107‑1 (MQ=255)                                                                                                                                              agcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTc                                                                              >  1:347107‑M2/72‑251 (MQ=255)                                                                                                                                               gcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCg                                                                             <  2:281889‑M2/181‑1 (MQ=255)                                                                                                                                                                  acagaacttgtaaaacaacgtaatagaattatccggaatggtggcgcctatAATTTCGCTTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGCGAGCTCGGCATTTTTGATGATGTAAt                                                          >  2:470832‑M2/52‑251 (MQ=255)                                                                                                                                                                     gaacttgtaaaacaacgtaatagaattatacggaatgatggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCACATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGTGATCGGCATTTTTGATGATGTACTCGt                                                       >  1:448224‑M2/49‑251 (MQ=255)                                                                                                                                                                              aaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGAt                                              <  1:12210‑M2/212‑1 (MQ=255)                                                                                                                                                                                          agaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCAc                                  >  2:487611‑M2/28‑251 (MQ=255)                                                                                                                                                                                                  ccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCAt                           <  1:415986‑M2/231‑1 (MQ=255)                                                                                                                                                                                                     gaatggtggcgactatAATTTCGATTGGTTCTCGGCTGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTg                       >  2:169618‑M2/17‑251 (MQ=255)                                                                                                                                                                                                      aatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCCCCTTTCATGCTgg                      >  2:121764‑M2/16‑251 (MQ=255)                                                                                                                                                                                                          gtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCgat                  >  2:96407‑M2/12‑251 (MQ=255)                                                                                                                                                                                                             gcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTc                                                                              >  2:45876‑M2/9‑188 (MQ=255)                                                                                                                                                                                                             gcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTc                                                                              <  1:45876‑M2/180‑1 (MQ=255)                                                                                                                                                                                                                          cGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCt  <  1:96407/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        GGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGCCCTGGATTGCT  >  NC_000913/2078718‑2079185

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