Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,505,423	G→A	C46Y (TGT→TAT)	yjhC →	KpLE2 phage‑like element; putative oxidoreductase YjhC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,505,423	0	G	A	100.0%	63.1 / NA	18	C46Y (TGT→TAT)	yjhC	KpLE2 phage‑like element; putative oxidoreductase YjhC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/11);  total (7/11)

TGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCC  >  NC_000913/4505345‑4505492                                                                               |                                                                      tGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGtata                                                                     <  1:905487/81‑1 (MQ=255)      atGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCaa                                                                 >  2:905516/1‑80 (MQ=255)         cATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAatat                                                              <  2:905488/80‑1 (MQ=255)                 tGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCtt                                                     <  1:905489/81‑1 (MQ=255)                     aaaaTTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGAt                                                 <  1:905490/81‑1 (MQ=255)                      aaaTTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATg                                                <  2:905491/81‑1 (MQ=255)                                tgtATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGtct                                       <  1:905492/80‑1 (MQ=255)                                   ataCGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAg                                   >  2:905524/1‑81 (MQ=255)                                      cGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTaa                                <  1:905493/81‑1 (MQ=255)                                            tGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGt                          >  1:905521/1‑81 (MQ=255)                                              aaaaTGGAGAAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTc                         <  1:905494/80‑1 (MQ=255)                                                     agaAAATATTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGt                 <  2:905495/81‑1 (MQ=255)                                                          atatTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCg            >  1:905518/1‑81 (MQ=255)                                                          atatTGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCg            >  2:905523/1‑81 (MQ=255)                                                             ttGCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAg         >  1:905530/1‑81 (MQ=255)                                                               gCCCGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGc        <  2:905496/80‑1 (MQ=255)                                                                  cGTGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCAcc    >  1:905519/1‑81 (MQ=255)                                                                    tGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCAcccc  <  1:905498/81‑1 (MQ=255)                                                                    tGAATTGCAGTATATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCAcccc  <  1:905497/81‑1 (MQ=255)                                                                               |                                                                      TGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATGGAGAAAATATTGCCCGTGAATTGCAGTGTATCAATATGTCAAGCTTGGATGCTTTAGTCTCAAGTAAATTAGTCGATTGCGTGATCGTAGCCACCCC  >  NC_000913/4505345‑4505492

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: ‑