Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,383,741	T→C	intergenic (+121/‑98)	epmA → / → yjeM	EF‑P‑lysine lysyltransferase/putative transporter YjeM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,383,741	0	T	C	100.0%	110.7 / NA	33	intergenic (+121/‑98)	epmA/yjeM	EF‑P‑lysine lysyltransferase/putative transporter YjeM
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/20);  total (13/20)

ATGCCTGACTGTTCCGGCGTGGCATTAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTT  >  NC_000913/4383513‑4383981                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 aTGCCTGACTTGTCCGGCGTGGCATTAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATAGTTGATATCATCCAGGCATCAATTCTTTATTGGttttt                                                                                                                                                                                                                             <  2:338157/250‑1 (MQ=255)   gCCTGACTGTTCCGGCGTGGCATTAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACg                                                                                                                                                                                                                          >  1:49607/1‑251 (MQ=255)           gTTCCGGCGTGGCATTAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGCTAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGctgcctgc                                                                                                                                                                                                                  >  2:186188/1‑251 (MQ=255)                      gCATTAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTgcg                                                                                                                                                                                                       <  1:186188/251‑1 (MQ=255)                          tAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcagt                                                                                                                                                                                                   >  1:263596/1‑251 (MQ=255)                                     cGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTg                                                                                                                                                                                        <  1:170749/251‑1 (MQ=255)                                             gATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACgtg                                                                                                                                                                                >  1:62487/1‑251 (MQ=255)                                                  tGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTtat                                                                                                                                                                            <  1:265495/250‑1 (MQ=255)                                                   ggCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTtata                                                                                                                                                                           <  2:201542/250‑1 (MQ=255)                                                                        ggCTGAAGTCAGAGCCTTTAGCGTTGAACGGGGATAATTCTGAAAAATGACTGAATGAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTTGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGAt                                                                                                                                                       <  1:31175/249‑1 (MQ=255)                                                                                 cATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGcc                                                                                                                                            <  2:271624/251‑1 (MQ=255)                                                                                             cGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATaaa                                                                                                                                <  1:251923/251‑1 (MQ=255)                                                                                                   ccGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGAt                                                                                                                          <  1:272715/251‑1 (MQ=255)                                                                                                      ggCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATg                                                                                                                         >  1:210987/1‑249 (MQ=255)                                                                                                        cATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTc                                                                                                                     <  2:213082/251‑1 (MQ=255)                                                                                                           aaTTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACAcg                                                                                                                                                                       >  1:70069/1‑198 (MQ=255)                                                                                                           aaTTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACAcg                                                                                                                                                                       <  2:70069/198‑1 (MQ=255)                                                                                                           aaTTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTg                                                                                                                   >  1:2702/1‑250 (MQ=255)                                                                                                            aTTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGAAATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGCAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGAt                                                                                                                 >  2:328382/1‑251 (MQ=255)                                                                                                               cTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGTCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAgg                                                                                                              >  2:34135/1‑251 (MQ=255)                                                                                                                         aCTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGAGAGGACTCATAtt                                                                                                     <  1:123384/250‑1 (MQ=255)                                                                                                                            gAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATtgatg                                                                                                 <  2:328180/251‑1 (MQ=255)                                                                                                                            gAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATtgatg                                                                                                 >  1:14444/1‑251 (MQ=255)                                                                                                                             aaTTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATtgatga                                                                                                <  1:111683/251‑1 (MQ=255)                                                                                                                                                                             tttCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTAc                                                 <  2:263596/250‑1 (MQ=255)                                                                                                                                                                             tttCAGTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTAGTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGGGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGAGAAAAAAGATGAGTCTGATAGGACTCATATTGATTATCTTTACTTCCGTATTTGGATTTTCGGATAGCCCATTGGCTTATTAc                                                 <  2:355559/250‑1 (MQ=255)                                                                                                                                                                                    tttCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGATTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCAAGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATAACTTAATggg                                         >  2:332637/1‑251 (MQ=255)                                                                                                                                                                                            aaaGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTg                                 <  2:14444/251‑1 (MQ=255)                                                                                                                                                                                                      aTTTAAGACTGGATATTGGATATCATCCAGGCATCAATTCTGTATTGGTTGTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTGTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTccatt                       <  2:75154/251‑4 (MQ=255)                                                                                                                                                                                                        ttAAGACTGGATATTTGATATCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCtttt                     >  2:105309/1‑251 (MQ=255)                                                                                                                                                                                                                         atatCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTg    <  1:357187/251‑1 (MQ=255)                                                                                                                                                                                                                         atatCATCCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTg    <  1:206292/251‑1 (MQ=255)                                                                                                                                                                                                                           atcatcCAGGCATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGtt  >  1:131687/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 ATGCCTGACTGTTCCGGCGTGGCATTAGGTGTTGATCGTCTGGTGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAAAAAAGATGAGTCTGATAGGACTCATATTGATGATCTTTACTTCCGTATTTGGATTTGCCAATAGCCCATCGGCTTATTACTTAATGGGTTATAGTGCGATTCCCTTTTATATATTTTCTGCATTGTT  >  NC_000913/4383513‑4383981

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