Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,814,212	C→T	H116Y (CAC→TAC)	katE →	catalase HPII, heme d‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,814,212	0	C	T	100.0%	54.4 / NA	18	H116Y (CAC→TAC)	katE	catalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/10);  total (8/10)

CAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTT  >  NC_000913/1814088‑1814321                                                                                                                             |                                                                                                              cAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGTGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAg                                                                                                 >  2:1225487/1‑139 (MQ=255)   gTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAgcg                                                                                               <  1:1221075/139‑1 (MQ=255)   gTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAgcg                                                                                               <  1:724331/139‑1 (MQ=255)                              aGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCAcg                                                                    >  2:1021025/1‑139 (MQ=255)                              aGGGCGTGCACATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAGATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCAcg                                                                    >  2:373816/1‑139 (MQ=255)                                     gcgcATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTa                                                                                 >  2:59831/1‑119 (MQ=255)                                     gcgcATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTa                                                                                 <  1:59831/119‑1 (MQ=255)                                     gcgcATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCa                                                             >  2:663606/1‑139 (MQ=255)                                                acgaTCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACgg                                                  <  1:431535/139‑1 (MQ=255)                                                           tCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGc                                       <  1:1021025/139‑1 (MQ=255)                                                              cTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGt                                                 <  2:1132490/126‑1 (MQ=255)                                                              cTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGt                                                 >  1:1132490/1‑126 (MQ=255)                                                                     ccGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCt                                                       <  2:117423/113‑1 (MQ=255)                                                                     ccGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCt                                                       >  1:117423/1‑113 (MQ=255)                                                                              gTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCa                                                                      >  1:529060/1‑89 (MQ=255)                                                                              gTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCa                                                                      <  2:529060/89‑1 (MQ=255)                                                                                              ggAAGATTTTATTCTGCGCGAGAAGATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAt    <  1:373816/139‑1 (MQ=255)                                                                                                aaGATTTTATTCTGCGCGAGAAAATCACCTACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGAttt  <  1:663606/139‑1 (MQ=255)                                                                                                                             |                                                                                                              CAGTGAAAATTATGCGCTGACCACTAATCAGGGCGTGCGCATCGCCGACGATCAAAACTCACTGCGTGCCGGTAGCCGTGGTCCAACGCTGCTGGAAGATTTTATTCTGCGCGAGAAAATCACCCACTTTGACCATGAGCGCATTCCGGAACGTATTGTTCATGCACGCGGATCAGCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTT  >  NC_000913/1814088‑1814321

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