Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,589,116	G→A	intergenic (‑37/+13)	yjiX ← / ← btsT	PF04328 family protein YjiX/pyruvate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,589,116	0	G	A	100.0%	94.1 / NA	30	intergenic (‑37/+13)	yjiX/btsT	PF04328 family protein YjiX/pyruvate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/13);  total (17/13)

AACAGCGCATACCGCCTTTACCATCGCCGCCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGCCGTCAGGTCAGTACCGTTAGC  >  NC_000913/4588881‑4589359                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aacaGCGCATACCGCCTTTACCATCGCCGCCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGtt                                                                                                                                                                                                                                       >  2:94694/1‑250 (MQ=255)                      cATCGCCGCCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCAc                                                                                                                                                                                                                 >  1:132869/1‑251 (MQ=255)                                                      ggaaagaaTTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCAACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCCCTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGAGCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGttt                                                                                                                                                                                 <  2:34238/248‑1 (MQ=255)                                                           gaaTTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTg                                                                                                                                                                            >  1:157988/1‑251 (MQ=255)                                                                 ttcATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTCCCAAACATAGTTCTCTCTGAGAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGt                                                                                                                                                                       <  1:40767/250‑1 (MQ=255)                                                                      aGCTCCTGTACGGCATGTCGGGACGGGTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGTTAAACACGGTTAGTGGTGCGAAGAGCTCTTCACGCCGCCTTCCGGGATTGTAACGTACGGTGTTTCTTTGTCAGTACGttt                                                                                                                                                                 >  2:182725/1‑251 (MQ=255)                                                                                 ggCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTAc                                                                                                                                                      >  2:204412/1‑251 (MQ=255)                                                                                        cGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGCGATCTTCACGACGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTACGTCAGTTTGTCGCTGATACGCAcctc                                                                                                                                               >  1:149687/1‑249 (MQ=255)                                                                                             gggttggtCTTCAGATGCTCGACGTAGTTTTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGATATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCTCTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGTATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGcc                                                                                                                                          <  1:181651/250‑1 (MQ=255)                                                                                                 tggtCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAGGCCAgg                                                                                                                                       <  1:85075/250‑1 (MQ=255)                                                                                                             gCTCGAAGGAGTTGTCGTAGGACGGAAAGCCAATCAACATCTTCGCCGACTAGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCAGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCAGGGATAGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGATTTGAAACCGt                                                                                                                           >  2:22704/1‑250 (MQ=255)                                                                                                                cGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTagaa                                                                                                                       <  1:206189/251‑1 (MQ=255)                                                                                                                          gTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATNCTg                                                                                                             >  2:253488/1‑251 (MQ=255)                                                                                                                                                 aCATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAAAGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAAt                                                                                      >  2:28139/1‑251 (MQ=255)                                                                                                                                                       tCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACt                                                                                   <  2:167545/248‑1 (MQ=255)                                                                                                                                                         gccgccTGGCCGAGATATTATTTTGCCCGTCCTAAGCTACTAAACATAGTTCACTATGATAAGAACAAAGCCCCGACGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGGACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGAc                                                                              >  2:169520/1‑251 (MQ=255)                                                                                                                                                                 gCCGAGATATGTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTg                                                                      <  1:28139/251‑1 (MQ=255)                                                                                                                                                                  ccGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTgg                                                                     <  1:253488/251‑1 (MQ=255)                                                                                                                                                                     agaTATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATTAGAACAAAGCCCCGCCGAAGAGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTGTCTTTGTCAGTACGTTTGTCGCAGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACACTCAGGAACAGAATACTCAGAccggcttaggagt                                                                  >  2:98152/1‑240 (MQ=255)                                                                                                                                                                       atatTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAg                                                                <  1:94694/251‑1 (MQ=255)                                                                                                                                                                                  ccTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAAc                                                     <  1:199930/251‑1 (MQ=255)                                                                                                                                                                                    tGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACg                                                    >  2:104794/1‑250 (MQ=255)                                                                                                                                                                                                 cAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATAtt                                      <  1:204412/251‑1 (MQ=255)                                                                                                                                                                                                          ttCACTCTGATAAGAACAAAGACCCGCAGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGAAAGCCAGGttt                                                                                                                                    >  2:130242/1‑148 (MQ=255)                                                                                                                                                                                                          tccACTCTGATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCCCGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGttt                                                                                                                                    <  1:130242/146‑1 (MQ=255)                                                                                                                                                                                                                  gATAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGCCg                     >  1:210538/1‑251 (MQ=255)                                                                                                                                                                                                                    tAAGAACAAAGCCCCGCCGAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGCAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGCCgtc                   >  2:237431/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccgccgAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGACGTCAGGTCAGTACCGt      >  2:127455/1‑251 (MQ=255)                                                                                                                                                                                                                                     cgAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGCCGTCAGGTCAGTACCGTTAGc  <  1:67974/251‑1 (MQ=255)                                                                                                                                                                                                                                     cgAAGCGAGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGACGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTAAGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGACTACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGCCGTCAGGTCAGTACCGTTAGc  >  1:111149/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    AACAGCGCATACCGCCTTTACCATCGCCGCCGTAGCGCGCGTTCTGGCGTTCGCGGAAGAATTCTTCATAGCTCATGTACGGCTTGTCGGGATGGTTGGTCTTCATATGCTCGACGTAGTTGTCGTAGTCTGGAATGCCAATCAACATCTTCGCCGCCTGGCCGAGATATTTTTTTGCCTGTCCTAAGTTACCAAACATAGTTCACTCTGATAAGAACAAAGCCCCGCCGAAGCGGGGCTAAACACGGTTAGTGGTGCGAAGAGATCTTCACGCCGCCTTCCGGGATTGGAACGTACGGTGTTTCTTTGTCAGTACGTTTGTCGCTGTTACGCACCGCAAGCCAGGTTTTGAAACCGTAGAAGATGATGCTGTACACCACAATCAGGAACAGAATACTCAGACCTGCGTTGGTGTAGTTGTTCACAACGATGTGGTTCATATTGGCAATCTGCTGCGCCGTCAGGTCAGTACCGTTAGC  >  NC_000913/4588881‑4589359

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