Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,590,212	G→A	intergenic (‑211/+642)	ydeT ← / ← hipA	fimbrial usher domain‑containing protein YdeT/serine/threonine‑protein kinase toxin HipA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,590,212	0	G	A	100.0%	33.2 / NA	14	intergenic (‑211/+642)	ydeT/hipA	fimbrial usher domain‑containing protein YdeT/serine/threonine‑protein kinase toxin HipA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/3);  total (11/3)

GCTTGACGGTGTAACCACTCATTCGACTGTAGGCACTGTCGGATAAGTTATAAAAACCTTGCGTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGGAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGG  >  NC_000913/1589980‑1590328                                                                                                                                                                                                                                         |                                                                                                                     gCTTGACGGTGTAACCACTCATTCGACTGTAGGCACTGTCGGATAAGTTATCAAAACCTTGCGTAGAATAGCAATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTc                                                                                                     >  1:228211/1‑250 (MQ=255)  cTTGACGGTGTAACCACTCATTCGAATGTAGGCACTGTCGGATAAGTTGTAAAAACCTTGCGTAGAATAGAGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAACTAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCa                                                                                                    >  1:241510/1‑250 (MQ=255)                             gTAGGCACTGTCGGATAAGTTATAAAAACCTTGCGTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACGTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTAcn                                                                                           <  1:156863/232‑2 (MQ=255)                                         ggATAAGTTATAAAAACCTTGCGTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCt                                                             >  1:174787/1‑250 (MQ=255)                                         ggATAAGTTATAAAAACCTTGCGTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAAACGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCt                                                             >  1:174777/1‑250 (MQ=255)                                                         ccTTGCGTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTCTTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTttgctt                                            >  1:85260/1‑251 (MQ=255)                                                               gTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGa                                       >  1:61595/1‑250 (MQ=255)                                                                                                   ncacaACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAgg  >  1:53586/2‑251 (MQ=255)                                                                                                                                  nagagTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCg                            >  1:64950/2‑194 (MQ=255)                                                                                                                                                  aCATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAAc                                    >  1:20489/1‑170 (MQ=255)                                                                                                                                                                                  caatATGCATGTTATTTATAAACCGTGACTACACAGAATCTTATTGCGAAAAAGGAAAGTGTCAGGAAAGGTCAAGAGATACCAATGGTTAATATg                                                                              >  1:89637/4‑96 (MQ=255)                                                                                                                                                                                         cATGTTATTTAGAAACCGTGACTAAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAagag                                                                                                <  1:166247/71‑1 (MQ=255)                                                                                                                                                                                         cATGTTATNNATNNANCNTGACTNAACAGAATATTATTGCGAACAAGGAAAGTGTCAGGAAAGGTCAagan                                                                                                >  2:166247/1‑70 (MQ=255)                                                                                                                                                                                                                         annattGCGNNCAAGGAAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAAc                                    <  2:20489/99‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                     GCTTGACGGTGTAACCACTCATTCGACTGTAGGCACTGTCGGATAAGTTATAAAAACCTTGCGTAGAATAGCGATATCCTGCGACCTGGATATTGGTTCCACAACAAGTTTATTGCTGCACTTGTTGAAGAGAGTATTTGCAAAAACATTTTATGTCGCTAAATAACATGACGATTTCATTATGCATGTTATTTATAAACCGTGACTAAACAGAATATTATTGCGAACAAGGGAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGG  >  NC_000913/1589980‑1590328

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