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	95.2%	60.3 / ‑6.3	21	*329* (TAG→TAA)	lpxK	lipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/10);  minor base T (0/1);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGC  >  NC_000913/969225‑969474                                                                                                                                |                                                                                                                           aaGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCgg                                                                                                                 <  2:607106/139‑1 (MQ=255)   gATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTATTTACGCCGCGGCa                                                                                                               <  2:518736/139‑1 (MQ=255)      gCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCg                                                                                                            <  2:174620/139‑1 (MQ=255)              aTGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTg                                                                                                    <  2:206508/139‑1 (MQ=255)                 ccGGGCCTTTGCGGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                 <  2:684276/139‑1 (MQ=255)                 ccGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATg                                                                                                 <  2:385667/139‑1 (MQ=255)                                                         gACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATg                                                         >  1:147655/1‑139 (MQ=255)                                                             cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                     >  2:744429/1‑139 (MQ=255)                                                             cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                     <  2:323164/139‑1 (MQ=255)                                                             cacaGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGcgcg                                                     >  1:715751/1‑139 (MQ=255)                                                                            atgaACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGc                                      >  1:107399/1‑139 (MQ=255)                                                                                   aGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACaa                               >  2:305343/1‑139 (MQ=255)                                                                                      gAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGc                            <  1:305343/139‑1 (MQ=255)                                                                                       aaaCTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAATAACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCAC‑‑CTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTg                         >  2:268630/1‑139 (MQ=255)                                                                                                  gCAACTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTaaacaaa                >  2:478481/1‑139 (MQ=255)                                                                                                    aaCTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCa                                            <  2:46994/109‑1 (MQ=255)                                                                                                    aaCTAACCTTGCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCa                                            >  1:46994/1‑109 (MQ=255)                                                                                                              gCTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGa    >  2:372621/1‑139 (MQ=255)                                                                                                               cTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCGTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAg   <  1:445690/139‑1 (MQ=255)                                                                                                               cTGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAg   <  1:679214/139‑1 (MQ=255)                                                                                                                tGGCTTCTGGCAACTAATTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGc  >  2:560215/1‑139 (MQ=255)                                                                                                                                |                                                                                                                           AAGATGCGGTGAAATGCCGGGCCTTTGCAGAAGAAAATTGGTGGTATTTGCCTGTAGACGCACAGCTTTCAGGTGATGAACCAGCGAAACTGCTTACGCAACTAACCTTGCTGGCTTCTGGCAACTAGTTACGCCGCGGCAGCGTTCGATTGATGGAGTCATGAATGTCGCTGCCGCACCTCTCCCTTGCTGATGCGCGTAATCTTCACCTTGCCGCACAAGGCCTGTTAAACAAACCCCGCCGTCGAGC  >  NC_000913/969225‑969474

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