Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,797,329	C→T	L128L (CTA→TTA)	waaL →	O‑antigen ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,797,329	0	C	T	100.0%	64.3 / NA	19	L128L (CTA→TTA)	waaL	O‑antigen ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/10);  total (9/10)

TATAATATCACTGCTATAGTTTGCTTATTGTCACTAATTTTACGTGGCAGACAAGAAAATTATAATATAAAAAACCTTATTCTTCCCCTTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGCTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTGGCGTTGCGATTCTTTATACTAAGAAAAATCATCCTTTTTTATTTTTATTA  >  NC_000913/3797086‑3797547                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                           tataATATCACTGCTATAGTTTGCTTATTGTCACTAATTTTACGTGGCAGACAAGAAAATTATAATATAAAAAACCTTATTCTTCCCCTTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATc                                                                                                                                                                                                                     >  2:196874/1‑251 (MQ=255)         cACTGCTATAGTTTGCTTATTGTCACTAATTTTACGTGGCAGACAATAAAATTATAAGAGAAAAAACCTTATTCTTCCCCGTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTAACATAGTTATTTAAATACTGCCAAAATATTTATATATGGTTATTTTATTGTTTTCTTGACACTAACTAGCCCGTTAAAAGCaaaaaaa                                                                                                                                                                                                              <  2:270325/250‑1 (MQ=255)                             tGTCACTAATTTTACGTGGCAGACAAGAAAAGTATAATATAAAAAACCTTATTCTTCCCCTTGCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTAGATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTAGACACttt                                                                                                                                                                                          <  2:153848/250‑1 (MQ=255)                                                             tataatataaAAAACCTTATTCTTCCCCTTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGCTAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGc                                                                                                                                                         >  1:328238/1‑251 (MQ=255)                                                                   ataaAAAACCTTATTCTTCCCCTTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGta                                                                                                                                                   <  1:313635/251‑1 (MQ=255)                                                                       aaaaCCTTATTCTTCCCCTTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGGGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATAtt                                                                                                                                               <  1:68179/251‑1 (MQ=255)                                                                                              atatTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCg                                                                                                                        >  1:252506/1‑251 (MQ=255)                                                                                                                         gtttGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAATTATTTATATTTGGTTCTTTTATTGTTTTCTTGCCACTAACTAGCCAGTTAAAATCATAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAAc                                                                                             >  2:69741/2‑251 (MQ=255)                                                                                                                          tttGGTATTCTGCGTTTAAAGTAGATACTTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATCCTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTATGTCATTTCTAATTGCTGGATATGCAATGTATATTCTTAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAAc                                                                                             >  2:201306/1‑250 (MQ=255)                                                                                                                             ggTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCACAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGCATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAACAgg                                                                                          >  1:150273/1‑250 (MQ=255)                                                                                                                                          tAAAGTAGATAATTCGCCATTTCGTTCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTc                                                                                                                                                                              <  1:105140/153‑1 (MQ=255)                                                                                                                                          tAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAACATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTc                                                                                                                                                                              >  2:105140/1‑153 (MQ=255)                                                                                                                                                                   gCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTg                                                    <  1:196874/250‑1 (MQ=255)                                                                                                                                                                           ccATAGTTATTTAAATACTGCCACAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCCGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTa                                                                                                                                              >  1:120129/1‑152 (MQ=255)                                                                                                                                                                           ccATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTAGTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGAGATGCAATGTATATTa                                                                                                                                              <  2:120129/152‑1 (MQ=255)                                                                                                                                                                                                                ggTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTGGCGTTGCGATTCTTTATACTAAGAAAAATCATCCTtttttattttt      <  2:6375/251‑1 (MQ=255)                                                                                                                                                                                                                      tttATTGTTTTCTTTACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATTCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGTTGTAGGAACGTCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTGGCGTTGCGATTCTTTATACTAAGAAAAATCATCCTTTtttatttttatt   <  1:323581/248‑1 (MQ=255)                                                                                                                                                                                                                                   tGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTGGCGTTGCGATTCTTTATACTAAGAAAAATCATCCTTTTTTATTTttatta  >  1:92824/1‑236 (MQ=255)                                                                                                                                                                                                                                   tGACACTAACTAGCCAGTTAAAATCAAAAAAAGAGAGTGGATTATACACTTTGTATTCTCTTTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGATCGCATTTCTTTTGGTGTAGGAACGGCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTGGCGTTGCGATTCTTTATACTAAGAAAAATCATCCTTTTTTATTTttatta  <  2:92824/236‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                           TATAATATCACTGCTATAGTTTGCTTATTGTCACTAATTTTACGTGGCAGACAAGAAAATTATAATATAAAAAACCTTATTCTTCCCCTTTCTATATTTTTAATAGGCTTGCTTGATTTAATTTGGTATTCTGCGTTTAAAGTAGATAATTCGCCATTTCGTGCTACTTACCATAGTTATTTAAATACTGCCAAAATATTTATATTTGGTTCTTTTATTGTTTTCTTGACACTAACTAGCCAGCTAAAATCAAAAAAAGAGAGTGTATTATACACTTTGTATTCTCTGTCATTTCTAATTGCTGGATATGCAATGTATATTAATAGCATTCATGAAAATGACCGCATTTCTTTTGGTGTAGGAACGGCAACAGGAGCAGCATATTCAACAATGCTAATAGGGATAGTTAGTGGCGTTGCGATTCTTTATACTAAGAAAAATCATCCTTTTTTATTTTTATTA  >  NC_000913/3797086‑3797547

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