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%	104.0 / NA	33	*329* (TAG→TAA)	lpxK	lipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/20);  total (13/20)

CGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGT  >  NC_000913/969231‑969476                                                                                                                          |                                                                                                                             cGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGt                                                                                                             <  1:944842/139‑1 (MQ=255)     gAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGa                                                                                                         <  2:1222137/139‑1 (MQ=255)        aTGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCg                                                                                                          <  1:1496439/135‑1 (MQ=255)        aTGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCg                                                                                                          >  2:1496439/1‑135 (MQ=255)           ccGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                   <  1:288053/139‑1 (MQ=255)                  tttGTAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCAt                                                                                            >  1:1307639/1‑139 (MQ=255)                    tGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGa                                                                                          >  1:1608745/1‑139 (MQ=255)                           gaaAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGc                                                                                   <  2:42488/139‑1 (MQ=255)                              aaTTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGc                                                                                <  2:1594754/139‑1 (MQ=255)                                   gtggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAc                                                                           <  2:1608745/139‑1 (MQ=255)                                   gtggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAc                                                                           >  2:43811/1‑139 (MQ=255)                                   gtggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACACAACTAACCTTGCTGGCATCTGGCAACTAATATCGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCCGCCGCAc                                                                           <  1:1406090/139‑1 (MQ=255)                                     ggtATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACct                                                                         >  2:1288310/1‑139 (MQ=255)                                        aTTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCCTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCtctc                                                                      <  1:1488801/139‑1 (MQ=255)                                                 tAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGa                                                             >  1:734445/1‑139 (MQ=255)                                                 tAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGa                                                             >  2:279445/1‑139 (MQ=255)                                                 tAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGa                                                             >  1:1404083/1‑139 (MQ=255)                                                       cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                       <  1:43811/139‑1 (MQ=255)                                                       cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                       <  2:734445/139‑1 (MQ=255)                                                                      atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAc                                                                           >  2:1351686/1‑104 (MQ=255)                                                                      atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAc                                                                           <  1:1351686/104‑1 (MQ=255)                                                                      atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGc                                        <  2:1307639/139‑1 (MQ=255)                                                                          aCCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGca                                    <  1:27922/139‑1 (MQ=255)                                                                                     tGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGtt                         <  2:1404083/139‑1 (MQ=255)                                                                                                            gCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAg                                >  2:1532072/1‑109 (MQ=255)                                                                                                            gCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAg                                <  1:1532072/109‑1 (MQ=255)                                                                                                            gCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGcgt  <  2:1211962/139‑1 (MQ=255)                                                                                                              ttCTGGCAACTAATTACGCCGCGGCAGCGTTCGAt                                                                                                        <  1:641475/35‑1 (MQ=255)                                                                                                              ttCTGGCAACTAATTACGCCGCGGCAGCGTTCGAt                                                                                                        >  2:641475/1‑35 (MQ=255)                                                                                                                  ggCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCtctc                                                                      >  2:374283/1‑65 (MQ=255)                                                                                                                  ggCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCtctc                                                                      <  1:374283/65‑1 (MQ=255)                                                                                                                      aCTAATTACGCCGCGGCAGCGTTCGATAGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                       <  1:872364/76‑1 (MQ=255)                                                                                                                      aCTAATTACGCCGCGGCAGCGTTCGATAGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                       >  2:872364/1‑76 (MQ=255)                                                                                                                          |                                                                                                                             CGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGCGT  >  NC_000913/969231‑969476

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