Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	371,589	C→T	G117G (GGC→GGT)	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,589	0	C	T	100.0%	123.0 / NA	39	G117G (GGC→GGT)	mhpC	2‑hydroxy‑6‑ketonona‑2,4‑dienedioate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/17);  total (22/17)

CGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATCA  >  NC_000913/371459‑371717                                                                                                                                   |                                                                                                                                 cGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCa                                                                                                                                <  2:81835/133‑1 (MQ=255)    gTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAg                                                                                                                             >  1:397291/1‑133 (MQ=255)         ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGTTGGGCAACTCGATGGGGGGTCATAGTTCtg                                                                                                                        >  1:378339/1‑133 (MQ=255)         ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCtg                                                                                                                        >  2:405716/1‑133 (MQ=255)          gCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCtgt                                                                                                                       <  1:323632/133‑1 (MQ=255)                gCGATTCGGTCGTTAACAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGtt                                                                                                                 <  2:56173/133‑1 (MQ=255)                gCGATTCGGTAGTTAGTAGTGGTTAGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGtt                                                                                                                 <  2:204616/133‑1 (MQ=255)                 cGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTc                                                                                                                <  1:405716/133‑1 (MQ=255)                   aTTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCAc                                                                                                              >  2:56262/1‑133 (MQ=255)                        gTCGTTAAGAGTGGTCCGAGATAGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTa                                                                                                         <  1:215328/133‑1 (MQ=255)                                     gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcg                                                                                            >  1:139243/1‑133 (MQ=255)                                     gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcg                                                                                            >  2:578147/1‑133 (MQ=255)                                     gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcg                                                                                            >  2:575574/1‑133 (MQ=255)                                     gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcg                                                                                            <  2:139243/133‑1 (MQ=255)                                      ttCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcgc                                                                                           >  1:89975/1‑133 (MQ=255)                                      ttCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAgcgc                                                                                           >  1:93068/1‑133 (MQ=255)                                               ggATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCAACCTTCAATGGCCGGAGCTCGTCGGCaaa                                                                                  >  2:301804/1‑133 (MQ=255)                                               gcaTCTTAATTCACGGATCCTGACTAGCCTGGGGGATCAACGGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCaaa                                                                                  <  2:301627/131‑1 (MQ=255)                                                gATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAAc                                                                                 <  1:516980/133‑1 (MQ=255)                                                           aCGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTAGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGGAAActgatgctgatg                                                                      >  1:227832/1‑133 (MQ=255)                                                                          cGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAGACTAGTGCTGATAGGCGCCGGTACGggc                                                       >  2:189809/1‑133 (MQ=255)                                                                             ggtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGTGGCGGTACGggcggc                                                    <  2:89975/133‑1 (MQ=255)                                                                             ggtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGggcggc                                                    >  1:21164/1‑133 (MQ=255)                                                                              gtggATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCa                                                   <  1:56262/133‑1 (MQ=255)                                                                               tggATCAACTGGATATCGCCAAAATCCACCTGCTAGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACg                                                          >  1:436813/1‑125 (MQ=255)                                                                               tggATCAACTGGATATCGCCAAAATCCACCTGCTAGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACg                                                          <  2:436813/125‑1 (MQ=255)                                                                                            tatCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGACGGAACGGGCGGGATGAATTTGTTTACg                                     >  2:199408/1‑133 (MQ=255)                                                                                               cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCg                                  >  2:52924/1‑133 (MQ=255)                                                                                               cGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCg                                  >  1:124110/1‑133 (MQ=255)                                                                                                    aaaTCCACCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGcc                             <  1:52924/133‑1 (MQ=255)                                                                                                          aCCTGTTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGa                       <  2:378339/133‑1 (MQ=255)                                                                                                          aCCTGCTGGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGGAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGa                       <  2:93068/133‑1 (MQ=255)                                                                                                            ctgctgGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                        <  2:615656/130‑1 (MQ=255)                                                                                                            ctgctgGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCg                        >  1:615656/1‑130 (MQ=255)                                                                                                               ctgGGCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGGCGGTGCCAACCGAAGGTa                  >  1:251256/1‑133 (MQ=255)                                                                                                                 gggCAACTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATACCAACCCAAGGTAtt                >  1:260748/1‑133 (MQ=255)                                                                                                                      aCTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCg           >  2:609966/1‑133 (MQ=255)                                                                                                                       cTCGATGGGGGGTCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGa          >  2:555061/1‑133 (MQ=255)                                                                                                                               ggggTCATAGTTCTGTGGTGTTCACCCTTATATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATCa  <  1:199408/133‑1 (MQ=255)                                                                                                                                   |                                                                                                                                 CGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGTTTGTTTACGCCGATGCCAACCGAAGGTATTAAGCGACTGAATCA  >  NC_000913/371459‑371717

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