Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	705,972	(T)9→8	intergenic (+82/‑121)	nagE → / → glnS	N‑acetylglucosamine specific PTS enzyme IIABC component/glutamine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	705,964	0	T	.	100.0%	54.1 / NA	14	intergenic (+74/‑129)	nagE/glnS	N‑acetylglucosamine specific PTS enzyme IIABC component/glutamine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (7/7);  total (7/7)

AATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATACGTTGTTTACGCTTTGAGGAATCC  >  NC_000913/705878‑706089                                                                                       |                                                                                                                              aaTCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAAttt                                                                                >  1:56977/1‑133 (MQ=255)                         ttggATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCt                                                        >  1:117480/4‑133 (MQ=255)                            gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTc                                                     >  1:362090/1‑133 (MQ=255)                            gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTc                                                     <  1:426735/133‑1 (MQ=255)                            gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTc                                                     >  1:611478/1‑133 (MQ=255)                            gATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTc                                                     >  2:426735/1‑133 (MQ=255)                                        gAGCGTCGCGCCCAAGAATGACAATCCGAGGACATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGa                                         <  2:278917/133‑1 (MQ=255)                                              cGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATAc                                   <  2:67172/133‑1 (MQ=255)                                               gCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACg                                  <  2:611478/133‑1 (MQ=255)                                                        aaTGACAAGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATAc                         <  2:369881/133‑1 (MQ=255)                                                               aGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATAc                         <  1:374051/126‑1 (MQ=255)                                                               aGCGGTGGAGATCTTCTCTGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATAc                         >  2:374051/1‑126 (MQ=255)                                                                              ctctGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATACGTTGTTTACGCTTTGAGGAATc   >  2:618115/1‑133 (MQ=255)                                                                               tctGCCGC‑TTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATACGTTGTTTACGCTTTGAGGAATcc  <  2:56977/133‑1 (MQ=255)                                                                                       |                                                                                                                              AATCAAAAAGTAATCTGCTTTATGCCTGATGCGACGCTTGAGCGTCGCATCCAACAATGACAAGCGGTGGAGATCTTCTCTGCCGCTTTTTTTTTCATCAATCATCCCCATAATCCTTGTTAGATTATCAATTTTAAAAAACTAACAGTTGTCAGCCTGTCCCGCTTATAAGATCATACGCCGTTATACGTTGTTTACGCTTTGAGGAATCC  >  NC_000913/705878‑706089

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