Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	372,780	G→A	L222L (CTG→CTA)	mhpD →	2‑hydroxypentadienoate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	372,780	0	G	A	100.0%	113.6 / NA	33	L222L (CTG→CTA)	mhpD	2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

GGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATG  >  NC_000913/372536‑373022                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ggAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCAcc                                                                                                                                                                                                                                              <  1:239157/251‑1 (MQ=255)            cgcgACTGGTCGATTCAGTTTGTCGATACCGTGTCAGATAACGCCTCGTGTGGGTTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCAt                                                                                                                                                                                                                                   <  2:64295/251‑1 (MQ=255)            cgcgACTGGGCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGAGGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCa                                                                                                                                                                                                                                    <  2:227491/250‑1 (MQ=255)                      cGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCgg                                                                                                                                                                                                                          <  1:103888/250‑1 (MQ=255)                          tCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGGAACCTCCGGCGGGGTTAGACCTGCAAAACTGGGGCCTGAAGATGGCGCGTAATAACGAAGGGGTTTCTAGCGGGCGCGGCGGCGCAAGCCCGGGACAGCCGCTTAATGCGGCCGTCTGGCTGGGACGCAAAAAGGGCAGTCTGGGTGACCCCCGACGCACCGGCGATATCATTTTTACCGGGGCAt                                                                                                                                                                                                                     >  2:202823/1‑251 (MQ=255)                            ggttggTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGTCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTa                                                                                                                                                                                                                   <  2:214323/246‑1 (MQ=255)                             gTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAg                                                                                                                                                                                                                  <  1:24350/251‑1 (MQ=255)                                                                                     cGGCGCAACGCCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCCGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCAt                                                                                                                                                           >  2:113745/1‑250 (MQ=255)                                                                                               accGGCGGGGCTTGACCTGAAAAACTGCGCCATGAAGATGACGCGGAATAAAGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATACGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTAACGGGGCATTAGGTCCGATGGTGGCGGTGAAAGCGGGCGATCGTTTTGAAGCCCATATTGAAGGccacgtatcagttg                                                                                                                                                >  2:228045/2‑238 (MQ=255)                                                                                                       ggTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGAc                                                                                                                                         >  1:147367/1‑250 (MQ=255)                                                                                                         ttAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACAtt                                                                                                                                      >  1:51533/1‑251 (MQ=255)                                                                                                             aCCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATTGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCa                                                                                                                                  <  2:251443/251‑1 (MQ=255)                                                                                                              ccTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCGCCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCaa                                                                                                                                 <  1:198925/251‑1 (MQ=255)                                                                                                                 gAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCTAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCGGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGGCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGACCATATTGAAGGCATAGGTTCAGTTGCTGTGACATTCTCAAgc                                                                                                                               <  2:122649/250‑1 (MQ=255)                                                                                                                 gAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAgcg                                                                                                                              <  2:147367/251‑1 (MQ=255)                                                                                                                 gAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAgcg                                                                                                                              <  2:271524/251‑1 (MQ=255)                                                                                                                  aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGATCAGTTGCTGCGACATTTTCAAgcgc                                                                                                                             >  1:57819/1‑251 (MQ=255)                                                                                                                             cATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCaaa                                                                                                                    >  1:60593/1‑249 (MQ=255)                                                                                                                                            taataaCGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACAGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATg                                                                                                     >  1:275831/1‑249 (MQ=255)                                                                                                                                                         gTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTc                                                                                      >  1:178751/1‑251 (MQ=255)                                                                                                                                                                            gCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGc                                                                    >  1:198469/1‑250 (MQ=255)                                                                                                                                                                                  gCCTGGGACATCCGCTTAAGGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACAtt                                                              <  2:241543/250‑1 (MQ=255)                                                                                                                                                                                           aTCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGACCCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATc                                                    >  2:174173/1‑251 (MQ=255)                                                                                                                                                                                                  tAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTgatgat                                             >  1:27230/1‑251 (MQ=255)                                                                                                                                                                                                          acGTCAGGCTGGCACGCAAAAGGGCCAGTCTGGGTGAACCGCCACACACCGGAGATCTCATTCTTACAGGGGCATTAGGTCCGATGGTGGCGGAGAATACGGGAGCTCGTTTTGAAGCCCATATTGAAGGCATAGGTTAAGTTGCAACGACCTTTTCAAGCGACGCCCCAAAAGGAAATCTGGCAGGAGTCAGCATAAAGGCGCCATTATGCGTTCTGCCCACATTGGGACTGATCTGATGCTTAAAAttt                                     >  2:122228/2‑251 (MQ=255)                                                                                                                                                                                                             tCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGGTCTGATGCTTAAAATTTTgt                                  >  2:223641/1‑250 (MQ=255)                                                                                                                                                                                                                       aCGNAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGTCCCAAAAGGAAGTCTGGCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCCTCACGGTCa                        <  2:261981/251‑1 (MQ=255)                                                                                                                                                                                                                       aCGCAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTc                         >  2:23710/1‑250 (MQ=255)                                                                                                                                                                                                                          cAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGtt                                                                                                                                                      >  1:195815/1‑122 (MQ=255)                                                                                                                                                                                                                          cAAAATGGCCAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGtt                                                                                                                                                      <  2:195815/122‑1 (MQ=255)                                                                                                                                                                                                                                   cAGTCTGGGTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGCGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTgg                >  2:69689/1‑247 (MQ=255)                                                                                                                                                                                                                                         gggTGAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGc         <  1:69689/248‑1 (MQ=255)                                                                                                                                                                                                                                             gAACCGCTACGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATg  <  2:198469/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   GGAGCCGCATTCGCGACTGGTCGATTCAGTTTGTCGATACCGTGGCAGATAACGCCTCCTGTGGGGTGTATGTCATCGGCGGTCCGGCGCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATG  >  NC_000913/372536‑373022

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