Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,545,140	C→T	C15C (TGC→TGT)	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,545,140	0	C	T	100.0%	51.6 / NA	19	C15C (TGC→TGT)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/8);  total (11/8)

GGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGGATGATCCCCAGGCTGTGGCCGATTTATCGCGTTTTGAAAATGGGCAAGAATTACCGCCAGGGACGTATCGCGTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTC  >  NC_000913/4544909‑4545377                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              ggTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATCTTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAACATGACGGGCGTAATGGAATAACGCAGGGGGACTTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGGTTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATc                                                                                                                                                                                                                            >  1:169662/1‑251 (MQ=255)    gCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTACGCATCGtt                                                                                                                                                                                                                         >  1:166871/1‑251 (MQ=255)                 gAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCACCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGtttttt                                                                                                                                                                                                            >  1:112277/1‑251 (MQ=255)                   aaGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTg                                                                                                                                                                                                          <  1:77546/251‑1 (MQ=255)                     gCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTc                                                                                                                                                                                                        <  1:169115/251‑1 (MQ=255)                                                                                   ggCGCACTTACCCCCAAAATGACGGGCGTCATGGAATAACGCAGGGGGAATTTTTAGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCAcc                                                                                                                                                                  >  1:104363/1‑227 (MQ=255)                                                                                     cgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGAAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCCTCTGCCGACCTCTAtttt                                                                                                                                         >  1:288996/1‑250 (MQ=255)                                                                                                              gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTGTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTATGTCATCTGCCGACCTCGATTTTAATCCGCGCTTTTTAGCGGATGATcn                                                                                                               <  1:29481/251‑2 (MQ=255)                                                                                                                   ggAATAACGCAGGGGGAATTTTTCGCCTGAATAACAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGGATGATCCCCAgg                                                                                                          >  1:273880/1‑251 (MQ=255)                                                                                                                      aTAACGCAGGGGGAATTTTTCGCCTGAACAAAAAGACTTGACTGCCGGGCGATTTTAAGCCGGAGGAATACTGTCATACCAGAATTTAAGAATTTATCAGCGAAACACACCATGTTTGCATGATCGTAAGCATCGTTTGGCTGGCTTTTTTGTGCGCATCGTTGTCGCAtgtg                                                                                                                                                                                     >  1:278460/1‑173 (MQ=255)                                                                                                                          cGCAGGGGGAATTTGTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACATTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTGTGGCTTGTGGTTGTGTTCGACTCGTTGTCGCCTGTGCGTTTGCCTCACAGGCGCCTTGGTCATCTGCCTACCTCTATTTTAATCCGCGCTTCTTAGCGGATGATCCCCAGGCTGTGGc                                                                                                   <  1:304914/251‑1 (MQ=255)                                                                                                                                           gCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTcn                                                                                                                                                                                                 <  1:59487/140‑2 (MQ=255)                                                                                                                                                              aCTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGAcc                                                                                                                                                 <  1:183487/169‑1 (MQ=255)                                                                                                                                                                                                       gAATTTAAGACTTTACCAGCGAAACACACAATGTTTGCNTATGNGNNAGNANCGNTTGGc                                                                                                                                                                                                                     <  2:166871/60‑1 (MQ=255)                                                                                                                                                                                                                   ttACCAGCGAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCTCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCAGACCTCTATTTTAATc                                                                                                                                     >  1:279222/1‑128 (MQ=255)                                                                                                                                                                                                                           naaaCACACAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGt                                                                                                                                                                                                                 >  1:39343/2‑44 (MQ=255)                                                                                                                                                                                                                           gAAACACACAATGTTTGCATATTCGTAAGCATCGTTTGGATGGTTTTTTTGTCCGACTCATTGTCGACTGTGCTTTTGCCGCACAGGGCCTTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTCTTTGAGGCTGATACCCAGGCTGTGGCCGATTTATCGCGTTTTGAAAATGGGCAAGCATTACCGCCAGGGACGTAACGCGTCGATATCTATTTGAATAATGGTTCTATGGCAACGCGTGATGTc  >  1:226509/1‑251 (MQ=255)                                                                                                                                                                                                                              acacacAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGttttttt                                                                                                                                                                                                           >  1:19542/1‑47 (MQ=255)                                                                                                                                                                                                                              acacacAATGTTTGCATATTCGTAAGCATCGTTTGGCTGGttttttt                                                                                                                                                                                                           <  2:19542/47‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              GGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGGATGATCCCCAGGCTGTGGCCGATTTATCGCGTTTTGAAAATGGGCAAGAATTACCGCCAGGGACGTATCGCGTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTC  >  NC_000913/4544909‑4545377

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