Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	371,588	(G)6→5	coding (350/867 nt)	mhpC →	2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	371,583	0	G	.	100.0%	114.4 / NA	28	coding (345/867 nt)	mhpC	2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/13);  total (15/13)

TCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATCA  >  NC_000913/371444‑371717                                                                                                                                            |                                                                                                                                       tCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATg                                                                                                                                        <  1:817160/140‑1 (MQ=255)   ctgctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                      >  1:773375/1‑140 (MQ=255)          aTTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATa                                                                                                                              <  1:832820/140‑1 (MQ=255)             gTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGtt                                                                                                                           <  1:359648/140‑1 (MQ=255)                cGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                      >  2:833438/1‑127 (MQ=255)                cGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                      <  1:833438/127‑1 (MQ=255)                     tGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                      >  2:373813/1‑122 (MQ=255)                     tGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATggg                                                                                                                                      <  1:373813/122‑1 (MQ=255)                       gggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGtt                                                                                                                 <  2:773375/140‑1 (MQ=255)                        ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGtt                                                                                                                 >  2:395269/1‑139 (MQ=255)                        ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGtt                                                                                                                 <  1:395269/139‑1 (MQ=255)                        ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTc                                                                                                                >  1:603195/1‑140 (MQ=255)                                  aTTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAAt                                                                                                      >  1:411659/1‑140 (MQ=255)                                        tCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCgg                                                                                                >  1:126818/1‑140 (MQ=255)                                                  tGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCgg                                                                                      >  2:340676/1‑140 (MQ=255)                                                    gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCa                                                                                    <  2:281277/140‑1 (MQ=255)                                                          gATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTgg                                                                              <  2:66218/140‑1 (MQ=255)                                                             cGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGc                                                                           <  1:340676/140‑1 (MQ=255)                                                                    tAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATggg                                                                    <  2:161913/140‑1 (MQ=255)                                                                      aTGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcg                                                                  >  1:719379/1‑140 (MQ=255)                                                                        gCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaaa                                                                                  >  1:48082/1‑122 (MQ=255)                                                                        gCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaaa                                                                                  <  2:48082/122‑1 (MQ=255)                                                                             aaTCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTAc                                                           >  2:729097/1‑140 (MQ=255)                                                                              aTCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACg                                                          >  2:575738/1‑140 (MQ=255)                                                                                ccTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACggg                                                        <  2:481743/140‑1 (MQ=255)                                                                                             gtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgtttg                                           >  1:53185/1‑140 (MQ=255)                                                                                                  tCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTAc                                      >  1:533288/1‑140 (MQ=255)                                                                                                     aCTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGcc                                   <  2:603195/140‑1 (MQ=255)                                                                                                        ggATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGAt                                >  1:40146/1‑140 (MQ=255)                                                                                                          atatCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGc                              >  1:696717/1‑140 (MQ=255)                                                                                                          atatCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGc                              >  2:50424/1‑140 (MQ=255)                                                                                                                ccAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                        >  1:854817/1‑140 (MQ=255)                                                                                                                  aaaaTCCACCTGCTGGGCAACTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGaa                      <  2:126818/140‑1 (MQ=255)                                                                                                                                     aCTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATc   <  1:575738/140‑1 (MQ=255)                                                                                                                                      cTCGAT‑GGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATCa  <  1:50424/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                       TCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATCA  >  NC_000913/371444‑371717

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