Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,799,622	C→T	G192D (GGT→GAT)	waaZ ←	lipopolysaccharide core biosynthesis protein WaaZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,799,622	0	C	T	100.0%	61.4 / NA	18	G192D (GGT→GAT)	waaZ	lipopolysaccharide core biosynthesis protein WaaZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/11);  total (7/11)

GATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTTATAAAACCACCTTTTTCCCGACGGTAA  >  NC_000913/3799406‑3799840                                                                                                                                                                                                                         |                                                                                                                                                                                                                           gATTCTTCTGCACTGATGTCTCGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAat                                                                                                                                                                                           >  1:119886/1‑250 (MQ=255)  aTTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAatat                                                                                                                                                                                         >  1:182382/1‑251 (MQ=255)          gCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGAtt                                                                                                                                                                                   <  2:118859/249‑1 (MQ=255)                          gTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCa                                                                                                                                                                 <  2:263944/251‑1 (MQ=255)                                         ttACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATg                                                                                                                                                    >  1:100380/1‑249 (MQ=255)                                                                         aaaTTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATc                                                                                                                   <  2:119886/250‑1 (MQ=255)                                                                                aaaTATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAg                                                                                                            <  1:209567/250‑1 (MQ=255)                                                                                 aaTATTAATATCTTTTACATGATCATGCATAAAACGAAAAAATGGCAATATTTTTAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAgac                                                                                                          <  1:375458/251‑1 (MQ=255)                                                                                          atCTTTTACATTATAATGCATAAAACGAAAAAATGGCAATATTTTTAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAGGGCAATAGTATGGCATGAGCAATAACCAAAACTAATGTCCTTGCAAAATCCAACCAGACGACCttttt                                                                                                 >  2:242584/1‑251 (MQ=255)                                                                                                                            ggCNATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTtgtg                                                               <  2:355177/251‑1 (MQ=255)                                                                                                                                         aTAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCaa                                                   >  1:210885/1‑250 (MQ=255)                                                                                                                                                   aCTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAAt                                         <  1:174982/250‑1 (MQ=255)                                                                                                                                                           ccGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATc                                                                                                                   <  1:2963/168‑1 (MQ=255)                                                                                                                                                           ccGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATc                                                                                                                   >  2:2963/1‑168 (MQ=255)                                                                                                                                                            cGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAAtcttct                               <  1:242584/251‑1 (MQ=255)                                                                                                                                                               gggCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATACGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATTAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTta                            <  2:67808/251‑1 (MQ=255)                                                                                                                                                                                  tcatcaTAGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACAAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTATTTATAAAACAACCTTTTTCCCg         >  2:167281/1‑251 (MQ=255)                                                                                                                                                                                         aGAAACGAGAACAGCTACCCGTTAAATCAAGATCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTTATAAAACCACCTTTTTCCCGACGGTaa  <  1:223112/251‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                           GATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGCATAAAACGAAAAAATGGCAATATTTTGAATAAATCCCTACTTAATTCCGAGGGCATGGGATTTTTATTCTCATCATAGAAACGAGAACAGCTACCCGTTAAATCAAGACCAGAACAAATAATCCGCGCATACTTAAGTGAATATGCGATTTGAATTGCAGCAAAGGCAATAGTATGGCATGAGCAATAACCAAGACTAATGTCCTTGCAAAATCCAACCAGACGACCTTTTTTAGACAACGGTACAGAGATCAGCAGTTCTTTGTGAATCTGTCTCAAAATATTAAATTTAATCTTCTTTATAAAACCACCTTTTTCCCGACGGTAA  >  NC_000913/3799406‑3799840

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