Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	305,596	G→A	L408F (CTC→TTC)	ecpD ←	fimbrial adhesin EcpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	305,596	0	G	A	100.0%	50.4 / NA	20	L408F (CTC→TTC)	ecpD	fimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/14);  total (6/14)

GGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAG‑GTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGG  >  NC_000913/305351‑305798                                                                                                                                                                                                                                                       |                                                                                                                                                                                                           ngTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAG‑GTTTTTGTAGCTCCGTTGCGGGTAATAAAGGCAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGACATCGGAGCTGAAGACACAATAGAAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTCCGATATAGCCGAAATCGAATGCCg                                                                                                                                                                                                       >  1:53329/2‑251 (MQ=255) ggTCAACCACAGCGCATCGGTTATATCCCGCCAGCTATCATCGCACCCGGCATCGTAG‑GTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCTTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCg                                                                                                                                                                                                       <  1:56864/251‑1 (MQ=255)  gTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAG‑GTTTTTGTAGCTCCGTTGCGGGTAATAACGGAAAGAGTCGCCGGGAACGGTACTTTCGCCTTACCGTCCTCGGCGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCg                                                                                                                                                                                                       >  1:200827/1‑250 (MQ=255)             ccctccGGTCATACCCCGCCAGCTATCATCGCCCCCGGCTTCGTGG‑GGTTTTGTAGCCCCGTTGCGGGTATAAAAGGAAAGTGTCGCCGGAAACGTTCCTTTCGCCTTCCCGTCATCGGAGCTGAAGACACAATAGGAGCCCCCGCCAATCTCCGGCGCGGGGCCGTTCGCCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTTTTTCCGATATAGCCGAAATCGAATGCCGACTCGCCGCTGc                                                                                                                                                                                           <  1:299747/246‑1 (MQ=255)                             acgtgtgctctgccaatctcGGCATCGTAGTGTTTTT‑TAGCTCCGTTGCGGGTAAGAAAGGAAAGTGTCGCCGGGAACGGTACTGTCGCCTTACCGTCATCGGAGCGGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTGTTACCTCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCcg                                                                                                                                                                           <  1:78011/232‑1 (MQ=255)                                           gCACCCGGCATCGTAG‑GTTTTTGTAGCTCCGGTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACGTGCGCCTTACCGTCATCGGATCTGAAGACACAAAAGGAGCGCCCGCCAATCACCTTCGCGGGTCCGGTCACCTTGATCAGCACTGCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACGCGCCGCTGCTAACATAGCCCTCCCGCGACGGTGTCTGAg                                                                                                                                                             <  1:431208/251‑1 (MQ=255)                                             aCCCGGCATCGTAG‑GTTTTTGTAGCTCCGTTGCGGGTCATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGCCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGAGTCTGAGt                                                                                                                                                            >  1:228929/1‑250 (MQ=255)                                                       gTAG‑GTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATAtccctcccn                                                                                                                                                                           <  1:44975/225‑2 (MQ=255)                                                                                           aaaGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATAtccctcccn                                                                                                                                                                           <  1:84100/190‑2 (MQ=255)                                                                                                                     gCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATTAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTc                                                                                    <  1:397858/251‑1 (MQ=255)                                                                                                                                  cGGAGCGGAAGACACAATATGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTGTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATAGCCCTCCCGCGACGTGGTCTGAGTATATGCATCGGAGACAATACTGATGCTGAAGCCACGGGGGTTGATATGCAGCGTGTTGGAGGTCGAAAATTCATACCAGCCAGAc                                                                       <  1:226270/251‑1 (MQ=255)                                                                                                                                            gACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGAtg                                                             <  1:17027/251‑1 (MQ=255)                                                                                                                                                             cgccAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGTTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATAtccctcc                                                                                                                                                                             <  1:301993/122‑1 (MQ=255)                                                                                                                                                                  aTCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGGTGGAGGTTAGGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTACCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGtt                                       >  1:173819/1‑251 (MQ=255)                                                                                                                                                                                               cAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATAGAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTCAATAGAACTTTGGTGTTGATATCGCTGa          <  1:139018/251‑1 (MQ=255)                                                                                                                                                                                                 gCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGTTAGAAAATTCATACCAGCCAGACTCCGGTGTTGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCgn             <  1:50153/246‑2 (MQ=255)                                                                                                                                                                                                            gCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAAGCACGGGGTTTGATAAGCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGAGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAgg  <  1:391195/246‑1 (MQ=255)                                                                                                                                                                                                                  gTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAATGCCGACTCGACGCTGACAACATATCACTCACGCGCCGGGGTCTGAGTATATTCATCGGAGATAATAatgatg                                                                                                                                   >  1:221326/1‑111 (MQ=255)                                                                                                                                                                                                                                      gATATAGCCGAAATCGAATGCCGACTCGCCGCTGCCAACATATCCCGCCCGCGACGGTGTATGAGTATAGTCAACGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAAGTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTgg                         <  1:19539/197‑1 (MQ=255)                                                                                                                                                                                                                                      gATATAGCCGAAATCGAATGCCGACTCGCAGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCCTCGGAGATAATACTGNNGCTGAANNc                                                                                                                           >  2:19539/1‑99 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                           GGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAG‑GTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGTACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGGTTTGATAATCAGCGTGTTGGAGGTAGAAAATTCATACCAGCCAGACTCCGGTGATGTGGTGTTCTGAAAGCGAAAGTTGAATAGATCTTTGGTGTTGATATCGCTGATCCCGAGG  >  NC_000913/305351‑305798

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