Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,799,390	(T)8→7	coding (807/852 nt)	waaZ ←	lipopolysaccharide KdoIII transferase, lipopolysaccharide core biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,799,383	0	T	.	100.0%	148.9 / NA	36	coding (814/852 nt)	waaZ	lipopolysaccharide KdoIII transferase, lipopolysaccharide core biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/20);  total (16/20)

AAAGGAAATTAATCTTGATTTTGTTTATATAAAATCTTTTCTTTTTGTGAAAAGTTCCTAAGCGCATTTTTATACCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCATTTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGC  >  NC_000913/3799246‑3799513                                                                                                                                          |                                                                                                                                   aaagGAAATTAATCTTGATTTTGTTTATATAAAATCTTTTCTTTTTGTGAAAAGTTCCTAAGCGCATTTTTATACCATATTATTTTAGTTTCAGGATGATTAATAAGCATAAGAATCGGTTGAAGTTCTATATTGCAtt                                                                                                                                   >  1:190381/1‑139 (MQ=255)                       gTTTATATAAAATCTTTTCTTTTTGTGAAAAGTTCCTAAGCGCATTTTTATACCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGa                                                                                                            <  2:204645/139‑1 (MQ=255)                       gTTTATATAAAATCTTTTCTTTTTGTGAAAAGTTCCTAAGCGCATTTTTATACCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGa                                                                                                            <  1:71529/139‑1 (MQ=255)                                 aaTCTTTTCTTTTTGTGAAAAGTTCCTAAGCGCATTTTTATACCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCa                                                                                                  >  2:159288/1‑139 (MQ=255)                                                                  aTTTTTATACCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACa                                                                 <  1:302690/139‑1 (MQ=255)                                                                     tttATACCATATTATTTTAGTTTCAGGATGAGTAATAAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCg                                                              <  2:190381/139‑1 (MQ=255)                                                                       tataCCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           >  2:750137/1‑124 (MQ=255)                                                                       tataCCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           <  1:750137/124‑1 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTg                                                                                               >  1:552789/1‑97 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           >  2:163154/1‑117 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           <  1:163154/117‑1 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACa                                                                 >  2:418139/1‑127 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACa                                                                 <  1:418139/127‑1 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATa                                                     >  2:331785/1‑139 (MQ=255)                                                                              tattatTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATa                                                     >  1:609056/1‑139 (MQ=255)                                                                              tattatTTTAGTTTCAGGATCATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTg                                                                                               <  2:552789/97‑1 (MQ=255)                                                                               attattTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAg                                                    <  1:159288/139‑1 (MQ=255)                                                                                 tattTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCt                                                  <  1:331785/139‑1 (MQ=255)                                                                                        gTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGa                                                                                                            <  1:477542/74‑1 (MQ=255)                                                                                        gTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGa                                                                                                            >  2:477542/1‑74 (MQ=255)                                                                                           tCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTg                                        >  1:225903/1‑139 (MQ=255)                                                                                             aGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGa                                                                                                            <  2:154991/69‑1 (MQ=255)                                                                                             aGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGa                                                                                                            >  1:154991/1‑69 (MQ=255)                                                                                                     taattaGCATAAGAATTGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCAGTATCATCTGATAAATcata                              <  1:523936/139‑1 (MQ=255)                                                                                                     taattaGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATtata                              <  1:402317/139‑1 (MQ=255)                                                                                                         taGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTg                                                                                               <  1:269208/70‑1 (MQ=255)                                                                                                         taGCATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTg                                                                                               >  2:269208/1‑70 (MQ=255)                                                                                                            cATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           <  1:704465/87‑1 (MQ=255)                                                                                                            cATAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           >  2:704465/1‑87 (MQ=255)                                                                                                              tAAGAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAat                     <  2:118889/139‑1 (MQ=255)                                                                                                                 gAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGaaa                                                       <  1:12510/102‑1 (MQ=255)                                                                                                                 gAATCGGTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGaaa                                                       >  2:12510/1‑102 (MQ=255)                                                                                                                      ggTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGAt                                      >  2:34887/1‑114 (MQ=255)                                                                                                                      ggTTGAAGTTCTATATTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGAt                                      <  1:34887/114‑1 (MQ=255)                                                                                                                                 tataTTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTACTAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATAATGc  >  2:692861/1‑139 (MQ=255)                                                                                                                                 tataTTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           >  2:479912/1‑66 (MQ=255)                                                                                                                                 tataTTGCA‑TTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATaaa                                                                           <  1:479912/66‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   AAAGGAAATTAATCTTGATTTTGTTTATATAAAATCTTTTCTTTTTGTGAAAAGTTCCTAAGCGCATTTTTATACCATATTATTTTAGTTTCAGGATGATTAATTAGCATAAGAATCGGTTGAAGTTCTATATTGCATTTTTTTTCTTGTCATATCTTTTGATTCTTCTGCACTGATGTCTTGAAGTTTAATAAATGGAATTACATCGTAAGAAATAGCTGTATCATCTGATAAATTATAAATATTAATATCTTTTACATTATCATGC  >  NC_000913/3799246‑3799513

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