Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	372,465	Δ1 bp	coding (351/810 nt)	mhpD →	2‑hydroxypentadienoate hydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	372,463	0	C	.	100.0%	107.9 / NA	26	coding (349/810 nt)	mhpD	2‑hydroxypentadienoate hydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/14);  total (12/14)

ATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAG  >  NC_000913/372383‑372539                                                                                 |                                                                             aTAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGc                                                                              <  1:76639/81‑1 (MQ=255)  tAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGcc                                                                             >  2:76676/1‑81 (MQ=255)    aCGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGccg                                                                            >  1:76675/1‑79 (MQ=255)       aaaTCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAAc                                                                       >  1:76685/1‑81 (MQ=255)         aTCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAAcc                                                                      >  2:76673/1‑80 (MQ=255)           cATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGa                                                                    >  2:76684/1‑80 (MQ=255)             ttCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGatat                                                                 <  2:76640/81‑1 (MQ=255)              tCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGatat                                                                 >  1:76686/1‑80 (MQ=255)                   tttCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCt                                                            <  1:76641/80‑1 (MQ=255)                    ttCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCtt                                                           <  1:76643/80‑1 (MQ=255)                    ttCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCtt                                                           <  2:76642/80‑1 (MQ=255)                    ttCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTc                                                          >  1:76687/1‑81 (MQ=255)                      cccGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTcga                                                        >  1:76678/1‑81 (MQ=255)                            ttCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAAtt                                                  <  2:76644/81‑1 (MQ=255)                                 cAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGt                                                <  1:76645/78‑1 (MQ=255)                                 cAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGt                                                <  1:76646/78‑1 (MQ=255)                                    ccccGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACg                                          <  1:76647/81‑1 (MQ=255)                                     cccGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGc                                         <  1:76648/81‑1 (MQ=255)                                        gCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCAt                                      <  2:76649/81‑1 (MQ=255)                                         cATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCAtt                                     >  1:76689/1‑81 (MQ=255)                                           ttGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTg                                    >  1:76694/1‑80 (MQ=255)                                             gAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAAt                                 >  1:76696/1‑81 (MQ=255)                                                       tCGCACTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTcc                       >  1:76691/1‑81 (MQ=255)                                                           aCTGGTGTTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGgcg                   >  1:76693/1‑81 (MQ=255)                                                               gtgtTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTg                <  2:76650/80‑1 (MQ=255)                                                                tgtTGAACCGCGATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGa              <  1:76651/81‑1 (MQ=255)                                                                        cgcgATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGtggtgg      <  2:76652/81‑1 (MQ=255)                                                                          cgATTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGTGGTgggg    <  2:76653/81‑1 (MQ=255)                                                                            aTTTG‑CCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAg  >  1:76692/1‑81 (MQ=255)                                                                                 |                                                                             ATAACGAAATCATTCCTTTTTCCCGTGTTCTGCAACCCCGCATTGAAGCGGAGATCGCACTGGTGTTGAACCGCGATTTGCCCGCAACCGATATCACCTTCGACGAATTGTATAACGCCATTGAATGGGTACTTCCGGCGCTGGAAGTGGTGGGGAG  >  NC_000913/372383‑372539

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