Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	372,753	C→T	R213R (CGC→CGT)	mhpD →	2‑hydroxypentadienoate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	372,753	0	C	T	100.0%	199.7 / NA	56	R213R (CGC→CGT)	mhpD	2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/32);  total (24/32)

GCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCT  >  NC_000913/372624‑372882                                                                                                                                  |                                                                                                                                  gcAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCa                                                                                                                         <  2:698109/140‑1 (MQ=255)  cAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAg                                                                                                                        <  1:50829/140‑1 (MQ=255)  cAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAg                                                                                                                        <  1:62338/140‑1 (MQ=255)   aaCGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGt                                                                                                                       <  1:104808/140‑1 (MQ=255)            cggGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGaa                                                                                                              <  2:740916/140‑1 (MQ=255)               ggTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCg                                                                                                           <  2:1013040/140‑1 (MQ=255)               ggTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCg                                                                                                           >  1:331761/1‑140 (MQ=255)                 ttAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCt                                                                                                         >  1:383211/1‑140 (MQ=255)                 ttAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCt                                                                                                         <  1:444415/140‑1 (MQ=255)                      ccTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCa                                                                                                    <  1:844012/140‑1 (MQ=255)                          aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCt                                                                                                         <  2:706492/131‑1 (MQ=255)                          aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCt                                                                                                         >  1:706492/1‑131 (MQ=255)                          aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCa                                                                                                    <  2:1148583/136‑1 (MQ=255)                          aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCa                                                                                                    >  1:1148583/1‑136 (MQ=255)                          aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCgg                                                                                                <  1:431716/140‑1 (MQ=255)                          aaaaaCTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCgg                                                                                                <  1:146173/140‑1 (MQ=255)                               cTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGata                                                                                           >  2:622030/1‑140 (MQ=255)                                 gcgcCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATc                                                                                         >  1:825634/1‑140 (MQ=255)                                    ccATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCAtt                                                                                      >  2:725946/1‑140 (MQ=255)                                            aTGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCgg                                                                              <  1:315612/140‑1 (MQ=255)                                                      ataaCGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGt                                                                    >  1:165215/1‑140 (MQ=255)                                                           gAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGAt                                                               >  1:506935/1‑140 (MQ=255)                                                             agagGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGAtgg                                                             >  1:419888/1‑140 (MQ=255)                                                               agGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGAtggtg                                                           >  2:882591/1‑140 (MQ=255)                                                               agGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGAtggtg                                                           >  2:797210/1‑140 (MQ=255)                                                               agGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGAtggtg                                                           >  2:707248/1‑140 (MQ=255)                                                                      tAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGa                                                    <  1:622030/140‑1 (MQ=255)                                                                           ggCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCg                                               >  2:381097/1‑140 (MQ=255)                                                                           ggCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCg                                               >  2:1055005/1‑140 (MQ=255)                                                                            gcgcGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCgg                                              >  1:827453/1‑140 (MQ=255)                                                                                        aTGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCGCGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTg                                  <  2:825634/140‑1 (MQ=255)                                                                                         tGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGGATGCGGGCGATCGTTTTGa                                 <  2:165215/140‑1 (MQ=255)                                                                                         tGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGa                                 <  1:797210/140‑1 (MQ=255)                                                                                           ccTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAg                               <  1:725946/140‑1 (MQ=255)                                                                                            cTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAg                               <  1:327055/139‑1 (MQ=255)                                                                                            cTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAg                               >  2:327055/1‑139 (MQ=255)                                                                                            cTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGc                              <  2:419888/140‑1 (MQ=255)                                                                                            cTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGc                              <  1:237948/140‑1 (MQ=255)                                                                                            cTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGc                              <  1:771923/140‑1 (MQ=255)                                                                                             tGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGcc                             <  1:707248/140‑1 (MQ=255)                                                                                                   ctCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATAtt                       <  2:123404/139‑1 (MQ=255)                                                                                                    tCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGcc                             <  2:164937/133‑1 (MQ=255)                                                                                                    tCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGcc                             >  1:164937/1‑133 (MQ=255)                                                                                                    tCCGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTg                      <  2:506935/140‑1 (MQ=255)                                                                                                     ccGCTTAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGa                     >  1:660850/1‑140 (MQ=255)                                                                                                          tAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCa                >  2:1148413/1‑140 (MQ=255)                                                                                                          tAATGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCa                >  1:357446/1‑140 (MQ=255)                                                                                                            aTGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATa              <  2:331761/140‑1 (MQ=255)                                                                                                            aTGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATa              <  1:362689/140‑1 (MQ=255)                                                                                                            aTGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATa              <  2:827453/140‑1 (MQ=255)                                                                                                            aTGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATa              >  1:939563/1‑140 (MQ=255)                                                                                                            aTGCGGCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATa              <  1:735939/140‑1 (MQ=255)                                                                                                                 gCCGTCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTc         <  2:880777/140‑1 (MQ=255)                                                                                                                     tCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGtt     >  1:719657/1‑140 (MQ=255)                                                                                                                     tCTGGCTGGCACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGtt     >  2:431277/1‑140 (MQ=255)                                                                                                                        ggctggcACGTAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTtgct  <  1:431277/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                  GCAACGTCCGGCGGGGTTAGACCTGAAAAACTGCGCCATGAAGATGACGCGTAATAACGAAGAGGTTTCTAGCGGGCGCGGCAGCGAATGCCTGGGACATCCGCTTAATGCGGCCGTCTGGCTGGCACGCAAAATGGCCAGTCTGGGTGAACCGCTGCGCACCGGAGATATCATTCTTACCGGGGCATTAGGTCCGATGGTGGCGGTGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCT  >  NC_000913/372624‑372882

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