Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	969,352	G→A	*329* (TAG→TAA)	lpxK →	lipid A 4'kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	969,352	0	G	A	100.0%	57.9 / NA	19	*329* (TAG→TAA)	lpxK	lipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/9);  total (10/9)

TGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTGG  >  NC_000913/969220‑969482                                                                                                                                     |                                                                                                                                   tGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTAcgc                                                                                                                              <  2:1009109/139‑1 (MQ=255)           gCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCg                                                                                                                    >  1:177550/1‑139 (MQ=255)                      ccGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                         <  2:177550/139‑1 (MQ=255)                         ggCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAg                                                                                                      >  1:543591/1‑139 (MQ=255)                                                ggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                               <  2:808698/139‑1 (MQ=255)                                                           gTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTg                                                                    <  1:220982/139‑1 (MQ=255)                                                            tAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGa                                                                   >  1:238804/1‑139 (MQ=255)                                                                      gCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCg                                                                                    >  1:555516/1‑112 (MQ=255)                                                                      gCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCg                                                                                    <  2:555516/112‑1 (MQ=255)                                                                        tttCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCg                                                                                                                >  1:379601/1‑82 (MQ=255)                                                                        tttCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCg                                                                                                                <  2:379601/82‑1 (MQ=255)                                                                                   ggaCCAGCGAAACTGCTTACGAAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATCTCGCTGCTGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCg                                            >  2:777064/3‑139 (MQ=255)                                                                                              aCTGCTTACGCAACTAACCTTGCTAGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTACTCTTCACCTTGCCGCACAAGGCCTg                                 >  2:876721/1‑139 (MQ=255)                                                                                                                   gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAg                                      <  1:1157000/113‑1 (MQ=255)                                                                                                                   gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAg                                      >  2:1157000/1‑113 (MQ=255)                                                                                                                    cTGGCTTCTGGCAACTAATTACGCCGCGTCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAg           <  1:777064/139‑1 (MQ=255)                                                                                                                            tGGCAACTAATTACGCCGCGGCAGCGTTCGATTGAt                                                                                                          >  2:657319/1‑36 (MQ=255)                                                                                                                            tGGCAACTAATTACGCCGCGGCAGCGTTCGATTGAt                                                                                                          <  1:657319/36‑1 (MQ=255)                                                                                                                             ggCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTgg  >  1:1083459/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   TGAAAAAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGTCGTTGG  >  NC_000913/969220‑969482

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