Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	371,542	C→T	Q102* (CAA→TAA)	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,542	0	C	T	100.0%	112.2 / NA	36	Q102* (CAA→TAA)	mhpC	2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/15);  total (21/15)

CTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGT  >  NC_000913/371419‑371673                                                                                                                            |                                                                                                                                     cTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGTGTGGTGGATTAACTGGATATCGCCa                                                                                                                       >  2:1169025/1‑139 (MQ=255)    gTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCaaaa                                                                                                                    >  1:898675/1‑139 (MQ=255)              gCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCTCctgct                                                                                                          >  1:405700/1‑139 (MQ=255)                  tCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGc                                                                                                      <  1:252543/139‑1 (MQ=255)                            ctgctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑gg                                                                                           >  1:736695/1‑139 (MQ=255)                              gctgGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑gggg                                                                                         <  2:898675/139‑1 (MQ=255)                                     tGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAAGAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATa                                                                                  <  2:736695/139‑1 (MQ=255)                                       tCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGt                                                                                >  2:304500/1‑139 (MQ=255)                                                 ggCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGt                                                                      >  1:850689/1‑139 (MQ=255)                                                  gCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGtt                                                                     <  2:850689/139‑1 (MQ=255)                                                  gCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGtt                                                                     <  2:767336/139‑1 (MQ=255)                                                   cAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGTGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTc                                                                    <  1:1169025/139‑1 (MQ=255)                                                     agagCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCAc                                                                  <  1:304500/139‑1 (MQ=255)                                                                           tGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGtt                                                                     <  2:224337/114‑1 (MQ=255)                                                                           tGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGtt                                                                     >  1:224337/1‑114 (MQ=255)                                                                             gTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATACACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCgg                                          >  2:760884/1‑139 (MQ=255)                                                                                 gcgATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTaaa                                                           <  2:1062368/118‑1 (MQ=255)                                                                                 gcgATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTaaa                                                           >  1:1062368/1‑118 (MQ=255)                                                                                     tCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGtt                                                                     <  1:437829/104‑1 (MQ=255)                                                                                     tCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGtt                                                                     >  2:437829/1‑104 (MQ=255)                                                                                             tAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTg                                   >  1:332837/1‑130 (MQ=255)                                                                                             tAATGCACGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTg                                   <  2:332837/130‑1 (MQ=255)                                                                                                   aCGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCa                                                                                                                >  2:1022158/1‑48 (MQ=255)                                                                                                   aCGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCa                                                                                                                <  1:1022158/48‑1 (MQ=255)                                                                                                   aCGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                    >  2:912147/1‑139 (MQ=255)                                                                                                   aCGAATCCTGAAAAGCGTGGTGGAGTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGggcggc                    >  1:303808/1‑139 (MQ=255)                                                                                                    cGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑g                                                                                            <  2:1137840/66‑1 (MQ=255)                                                                                                    cGAATCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑g                                                                                            >  1:1137840/1‑66 (MQ=255)                                                                                                       aTCCTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTa                <  1:912147/139‑1 (MQ=255)                                                                                                         ccTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACg              >  1:1097654/1‑139 (MQ=255)                                                                                                          cTGAAAAGCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACgg             >  2:844175/1‑139 (MQ=255)                                                                                                                 gCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACggg            >  1:905536/1‑133 (MQ=255)                                                                                                                 gCGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACggg            <  2:905536/133‑1 (MQ=255)                                                                                                                  cGTGGTGGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGATGGGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCAt      >  1:1063509/1‑139 (MQ=255)                                                                                                                   gtggtgGATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGa    >  2:1193197/1‑139 (MQ=255)                                                                                                                     ggtggATTAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAgt  >  2:198927/1‑139 (MQ=255)                                                                                                                            |                                                                                                                                     CTGGTAGAGGCGGGCTATCGGGTGATCCTGCTGGATTGTCCGGGTTGGGGCAAGAGCGATTCGGTCGTTAATAGTGGTTCGCGATCGGATCTTAATGCACGAATCCTGAAAAGCGTGGTGGATCAACTGGATATCGCCAAAATCCACCTGCTGGGCAACTCGAT‑GGGGGGCCATAGTTCTGTGGCGTTCACCCTTAAATGGCCGGAGCGCGTCGGCAAACTGGTGCTGATGGGCGGCGGTACGGGCGGCATGAGT  >  NC_000913/371419‑371673

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