Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,837	G→A	D158N (GAT→AAT)	lpxK →	tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,837	0	G	A	96.9%	102.2 / NA	32	D158N (GAT→AAT)	lpxK	tetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/16);  minor base C (1/0);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGAT  >  NC_000913/968593‑969078                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  tATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGGTTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGCTCGTTCTGATGCGGGAAAAAGCC‑TTCTGGCGCAACACCCCGATGGGCAGATCCTCGTAACCGACGACGGGTTACAGCATTAACGTCTGGCGCGGAATGTgg                                                                                                                                                                                                                                             >  1:5983/1‑251 (MQ=255)                 gTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGa                                                                                                                                                                                                                             <  2:5983/251‑1 (MQ=255)                 gTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAAACTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTCTCCCGGTCGTTCTGATGCGG‑TAAAAGCCATTCTGGGGCAACACCCTGATGTGCAGATCATCGTAACCGACGAAGGGTTACAGCCTTAACGGCTGGCGCGGAATGTGGAAATTGTCGTTATTg                                                                                                                                                                                                                              >  1:89618/1‑250 (MQ=255)                                ggtggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCtt                                                                                                                                                                                                              <  1:325174/251‑1 (MQ=255)                                    gtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATTAACCTGTGGTGATTTTTCAACGCAATGATGCGCCTGTTGCGGTTTTTCCCCTTAGTTCTGATGCGG‑TAAAAGACCTTCTGGCGGAACACCCTTATGTGCAGATCATCGTAACCGACAACGGTTTAAAGCCATACCGTCTGGGGCGTAATGTGGAAATTGTCGTTTTTGTTGGGGTGCGTCGCttttgc                                                                                                                                                                                                          >  2:229705/1‑248 (MQ=255)                                                  tATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCtggtg                                                                                                                                                                                               <  1:330850/248‑1 (MQ=255)                                                            tATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCggg                                                                                                                                                                                   <  1:257766/250‑1 (MQ=255)                                                               tGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTTTCAAAGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACCACGGGTTACAGCATTTACGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCgggggcc                                                                                                                                                                               >  1:38350/1‑248 (MQ=255)                                                                              caacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGagcgagcg                                                                                                                                                                <  2:80480/251‑1 (MQ=255)                                                                              caacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGagcgagcg                                                                                                                                                                <  1:163892/251‑1 (MQ=255)                                                                                     acaGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCt                                                                                                                                                         <  2:38350/251‑1 (MQ=255)                                                                                                 tgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTAGTTCTGATGCGG‑TAAAAGACATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGGGGGGGGCTTTAAGGCgggtg                                                                                                                                             >  1:51871/1‑248 (MQ=255)                                                                                                              aTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCGGGGGCAAAGGGTGAGCGGGCGGGGGGGTTTAAGTCGGTTGATGCggggatcgt                                                                                                                                >  2:95216/1‑244 (MQ=255)                                                                                                                  aTCAACGCACTGATGCGCCTGTTGGGGTTTCTCCCGTTCGTTCGGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAAc                                                                                                                            <  2:228659/251‑1 (MQ=255)                                                                                                                            tGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGGGGGGGGGGGTTAAAGGCGGTTGAGGCGGGAAACGTCACcggcggggctc                                                                                                                  >  2:337087/1‑246 (MQ=255)                                                                                                                                  gactggtGCGGTGTCGCCCGTTAGTGCTGATGCGG‑TAAAAGCCAGTCTGGCGCAACACCCTGATGTGCAGATCATCGGAACCGACGACGGTGTACAGCAGTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATGGATGGTGGGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCa                                                                                                            <  2:414095/245‑1 (MQ=255)                                                                                                                                           ggTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGGCCCTCGCCGCggggga                                                                                                    >  1:192921/1‑246 (MQ=255)                                                                                                                                               tctcCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGGTTCGGTAATCGTCAACGGCGGGGGCCCTCGCAGCGGGGGAATccc                                                                                               >  1:303286/1‑251 (MQ=255)                                                                                                                                               tctcCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGTAATCGTCAACGGCGGTGTCCCTCCCAGCGGGGGAATccc                                                                                               >  1:43463/1‑251 (MQ=255)                                                                                                                                                     gttcgttcTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTCATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGCGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGCGGGGTCCCCCGCAGCGGTGAAATcccccatc                                                                                          >  1:135926/1‑246 (MQ=255)                                                                                                                                                                          aGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGAGCGGGGCGCTTAAAGTCGGGTGGTGGGGGAATCGTCAACGGCGGTGTCCCCCCCAGGGGGGGAATCCCCAAGCATCTGGTTGCGGGTCAggggg                                                                     >  1:122144/1‑248 (MQ=255)                                                                                                                                                                                      gcAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCCCCGCGGGGGAATCCCCCTGCATCTGCTGCCGGGGCAGGGGGGGAATTTACGt                                                          >  1:135969/1‑250 (MQ=255)                                                                                                                                                                                      gcAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTa                                                         <  2:92066/251‑1 (MQ=255)                                                                                                                                                                                           acCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCACTGCGGGAGCGAGGGGGGGGCTTAAAGTCGGTTGATGCGGGAATCGGCAACGGCGGGGGCCCCCCCAGCGGGGGAATACCCCCGCATCTTCTGCCGGGTCAAGCGGGGGATTTAcctatccgg                                                    >  2:362085/1‑243 (MQ=255)                                                                                                                                                                                              cTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCTGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGGCCGGCGGGGAATTTACGTACCGGGAcg                                                 >  1:348002/1‑251 (MQ=255)                                                                                                                                                                                                gATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGGGCCTAAAGTCGGGTGCTGCGGCAATCGGCAACGGCGGGGGCCCCCGCAGCGGGGGAATCCCCATGAATCTGCTGGCGGGGCAGGCGGGGGATTTAccgaccgggaccac                                               >  1:246611/1‑238 (MQ=255)                                                                                                                                                                                                   gtgCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGAGGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGtt                                             <  2:335875/250‑1 (MQ=255)                                                                                                                                                                                                           catcGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCTGGTCAGGCTGTGAATTTACGTACCTGTACGCGTTGTGActtt                                     <  2:51871/250‑4 (MQ=255)                                                                                                                                                                                                               gTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGGGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCAGTACGCGTTGTGACGTTGCTCa                                <  2:267458/251‑1 (MQ=255)                                                                                                                                                                                                                       cgacgGTTTACAGCATTACCGTATGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGGCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAAc                        <  2:135926/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCg             <  1:229705/251‑1 (MQ=255)                                                                                                                                                                                                                                               ggCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCACTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGAt  <  1:159890/249‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGG‑TAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGAT  >  NC_000913/968593‑969078

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