Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,547,328	G→A	A745T (GCC→ACC)	fimD →	type I fimbriae usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,547,328	0	G	A	100.0%	33.7 / NA	14	A745T (GCC→ACC)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/8);  total (6/8)

ACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCAC  >  NC_000913/4547103‑4547507                                                                                                                                                                                                                                  |                                                                                                                                                                                    aCGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGCTGATATTAAGCAGCTCTATTACGGAGTACGTGGTGGGGTCCTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGAAAAAGACGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTCTACCGtgctgcgatcggaagagcaccc                                                                                                                                                            >  1:174853/1‑235 (MQ=255)                   ngtTACGGCACTGACAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGgc                                                                                                                                          >  1:203810/2‑251 (MQ=255)                            aaTGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATg                                                                                                                                                                      <  1:9735/214‑1 (MQ=255)                                             cAGCCATAGCGATGATATTAAGCAGCGCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTTTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATGc                                                                                                                                                                     <  1:38416/198‑1 (MQ=255)                                                  aTAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGGACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCtn                                                                                                                       <  1:75824/239‑2 (MQ=255)                                                                  gcagcTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAAGGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGGACCGACTGGCGTGGTTATACCGTGCTGCCTTATGCCACGGAATATCGGGAAAATAGAGTTGCGCTGGATACCAATGCCCTGGCTGATAACGTCGATTTAGATAACGCgg                                                                                           <  1:225068/251‑1 (MQ=255)                                                                                gAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGGGCAAAAGATGCGAAAGTCGAAAACCAGACGGGGGTGGGTACCGACTGGCGGGGTGATACCGTGCTGCCTTAGGCCAGTGAATAGCGGTAAAATAGAGTGGCGCTTGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACgttg                                                                               <  1:251215/249‑1 (MQ=255)                                                                                            gggTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCAGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATGACACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTggg                                                                 >  1:97605/1‑251 (MQ=255)                                                                                                                                 gCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATGCCACTGAATATTGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTANc                                                                                   <  1:199438/196‑1 (MQ=255)                                                                                                                                   cGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCTCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGc                                                               <  1:191126/214‑1 (MQ=255)                                                                                                                                       aaaCGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATAACGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATAACAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGCTCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAAc                      >  1:90571/1‑251 (MQ=255)                                                                                                                                                            tAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTCCCCACTGGCGTCCTTCTACCGTGCTGCCTTATGCCACTGAATCTCGAGAAAACAGAGTGGAGCTCGATACCCATCCCCTGGCTGCTACCGTAGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCATGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCcac  >  1:141353/1‑250 (MQ=255)                                                                                                                                                                  gcCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATACCGTGCTGCCTTACGCCACTGAATATCGGGAACATAGAGTGGCGCTGGGTACCAATACCCTGGCTGATAACGTCGCTTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTgg                                                                  >  1:129926/1‑180 (MQ=255)                                                                                                                                                                                                                       ggCGTGGTTATACCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGTCGCTGTATGCCAATACCCTGGCTGATAACGTCGATTGAGATAACTCGGGTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAgcgcgc                                   <  1:181807/158‑1 (MQ=255)                                                                                                                                                                                                                       ggCGTGGTNANNNCNNGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGNTACCACTACCCTGGCNGANAACGNCGANTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGCGTTTAAAgcgcgc                                   >  2:181807/1‑158 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                    ACGCTGAATTATCGCGGTGGTTACGGCAATGCCAATATCGGTTACAGCCATAGCGATGATATTAAGCAGCTCTATTACGGAGTCAGCGGTGGGGTACTGGCTCATGCCAATGGCGTAACGCTGGGGCAGCCGTTAAACGATACGGTGGTGCTTGTTAAAGCGCCTGGCGCAAAAGATGCAAAAGTCGAAAACCAGACGGGGGTGCGTACCGACTGGCGTGGTTATGCCGTGCTGCCTTATGCCACTGAATATCGGGAAAATAGAGTGGCGCTGGATACCAATACCCTGGCTGATAACGTCGATTTAGATAACGCGGTTGCTAACGTTGTTCCCACTCGTGGGGCGATCGTGCGAGCAGAGTTTAAAGCGCGCGTTGGGATAAAACTGCTCATGACGCTGACCCAC  >  NC_000913/4547103‑4547507

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