Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	371,515	G→A	A93T (GCA→ACA)	mhpC →	2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	371,515	0	G	A	100.0%	172.3 / NA	52	A93T (GCA→ACA)	mhpC	2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/27);  total (25/27)

CGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGT  >  NC_000913/371370‑371658                                                                                                                                                  |                                                                                                                                                cGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACa                                                                                                                                               <  7:99083/148‑1 (MQ=255) cGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACa                                                                                                                                               <  1:355852/148‑1 (MQ=255)        ggTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCt                                                                                                                                       <  4:177334/149‑1 (MQ=255)               ctggctggGCGAACTTCAGCCGCAATATCGAGCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGc                                                                                                                                <  7:406615/149‑1 (MQ=255)                  gctggGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCgtg                                                                                                                             >  5:266611/1‑149 (MQ=255)                      tgCGAACTTCAGACGAAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGTGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCgttgtg                                                                                                                          <  6:266611/147‑1 (MQ=255)                      ggCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCgtggtg                                                                                                                          <  7:352513/148‑1 (MQ=255)                      ggCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCgtggt                                                                                                                           <  3:390006/147‑1 (MQ=255)                            cTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAAc                                                                                                                   >  2:412978/1‑149 (MQ=255)                            cTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAAc                                                                                                                   <  5:347772/149‑1 (MQ=255)                               cAGCCGCAATATCGATCCGCTGGTAGAGGCGGACTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACt                                                                                                                  <  2:424194/147‑1 (MQ=255)                                  ccGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCgtgg                                                                                                                            <  3:250037/134‑1 (MQ=255)                                  ccGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCgtgg                                                                                                                            >  4:250037/1‑134 (MQ=255)                                     cAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCg                                                                                                          >  7:144629/1‑149 (MQ=255)                                             aTCCGCGGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATc                                                                                                  <  7:79567/149‑1 (MQ=255)                                                  cTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTCGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACct                                                                                             >  4:299803/1‑149 (MQ=255)                                                     gTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAAc                                                                                                                   >  7:290085/1‑124 (MQ=255)                                                     gTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAAc                                                                                                                   <  8:290085/124‑1 (MQ=255)                                                     gTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACctgct                                                                                          >  5:406744/1‑149 (MQ=255)                                                     gTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACct                                                                                             >  1:239494/1‑146 (MQ=255)                                                       agagGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTgg                                                                                        >  3:9857/1‑149 (MQ=255)                                                        gagGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTgg                                                                                        <  1:22415/148‑1 (MQ=255)                                                            cGGGCTATCGGGTGATCGTGCTGGTTTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGAGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAAc                                                                                   <  7:132631/149‑1 (MQ=255)                                                                   tCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAt                                                                              <  1:412978/147‑1 (MQ=255)                                                                     gggTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGa                                                                                                               <  2:24471/112‑1 (MQ=255)                                                                     gggTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGa                                                                                                               >  1:24471/1‑112 (MQ=255)                                                                      ggTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggggg                                                                         >  4:75253/1‑149 (MQ=255)                                                                      ggTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCATAATACAAGAATCCTGAAAAGCGTGGTGGATCAACTGGAAATCGCCAAAATCCCCATGCTGGGCAACTCGATggggg                                                                         >  6:36450/1‑149 (MQ=255)                                                                        tGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTggg                                                                                       >  6:350224/1‑133 (MQ=255)                                                                        tGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTggg                                                                                       <  5:350224/133‑1 (MQ=255)                                                                        tGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGc                                                                       <  4:9857/149‑1 (MQ=255)                                                                             ctgctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAg                                                                  >  4:385404/1‑149 (MQ=255)                                                                                   gATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGTGCAACTCGATGGGGGGCCATAGTTCtgt                                                            <  3:75253/149‑1 (MQ=255)                                                                                         ccGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGt                                                       >  8:234459/1‑148 (MQ=255)                                                                                          cGGGTTGGGGCAAGAGCGAATCGGTCGTTAATAGTGGTTCGCGATCGGATCGTAATACACGAATCCGGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAGGGGGGGCCATAGTGCTGTGGCGTTc                                                     >  4:360854/1‑149 (MQ=255)                                                                                             ggtGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGAGCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCAcc                                                  <  6:195137/147‑1 (MQ=255)                                                                                             gTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCAcc                                                  >  2:56905/1‑149 (MQ=255)                                                                                                  ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCTTCTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTa                                              >  4:282052/1‑148 (MQ=255)                                                                                                      agagCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACTCGCATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATgg                                         >  6:394435/1‑149 (MQ=255)                                                                                                         gCGATTCGGTCGTTATTAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCg                                      <  3:299803/149‑1 (MQ=255)                                                                                                           gATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCa                                                    <  4:126547/133‑1 (MQ=255)                                                                                                             ttCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGCCCTTAGTTATGTGGCGTTCCCCGTTAAATGGCCGGAgc                                  >  2:125693/1‑149 (MQ=255)                                                                                                                ggTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGcgcg                               <  3:385404/149‑1 (MQ=255)                                                                                                                 gTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACt                                                                                  <  2:60423/97‑1 (MQ=255)                                                                                                                 gTCGTTAATAGTGGTTCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACt                                                                                  >  1:60423/1‑97 (MQ=255)                                                                                                                            tGGTTCGCGATCGGATCTTAATACACGATTCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCAc                                                   <  6:209853/117‑1 (MQ=255)                                                                                                                            tGGTTCGCGATCGGATCTTAATACACGATTCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCAc                                                   >  5:209853/1‑117 (MQ=255)                                                                                                                               ttCGCGATCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGc                <  8:144629/149‑1 (MQ=255)                                                                                                                                     aTCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATgg          >  1:329882/1‑149 (MQ=255)                                                                                                                                      tCGGATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATggg         >  2:359715/1‑149 (MQ=255)                                                                                                                                         gATCTTAATACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcgg      >  1:206864/1‑149 (MQ=255)                                                                                                                                               aaTACACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGGTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGt  <  5:36450/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                CGGCCCGGGTGCTACTGGCTGGGCGAACTTCAGCCGCAATATCGATCCGCTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGT  >  NC_000913/371370‑371658

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