Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,545,466	T→C	V124A (GTC→GCC)	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,466	0	T	C	100.0%	118.6 / NA	34	V124A (GTC→GCC)	fimD	type I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/14);  total (20/14)

TCGCGTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCG  >  NC_000913/4545329‑4545600                                                                                                                                          |                                                                                                                                       tCGCGTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCg                                                                                                                                      >  2:178823/1‑140 (MQ=255)  cgcgTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGc                                                                                                                                     <  1:178823/140‑1 (MQ=255)  cgcgTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGc                                                                                                                                     >  2:575445/1‑140 (MQ=255)  cgcgTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGc                                                                                                                                     <  1:29777/140‑1 (MQ=255)  cgcgTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGc                                                                                                                                     <  2:163229/140‑1 (MQ=255)            tCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAAt                                                                                                                           <  2:543362/140‑1 (MQ=255)                      taatGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCgg                                                                                                                 >  1:183307/1‑140 (MQ=255)                        atGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGgat                                                                                                               >  2:662948/1‑140 (MQ=255)                              atatGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGcc                                                                                                         >  1:434294/1‑140 (MQ=255)                                  ggCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCtgtg                                                                                                     <  2:230960/140‑1 (MQ=255)                                    cAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCtgtgtg                                                                                                   <  2:183307/140‑1 (MQ=255)                                     aaCGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCCGATGATGCCTGTGTGc                                                                                                  >  1:290319/1‑140 (MQ=255)                                      aCGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGcc                                                                                                 <  2:250400/140‑1 (MQ=255)                                       cgcgTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCtgt                                                                                                      <  2:232200/134‑1 (MQ=255)                                       cgcgTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCtgt                                                                                                      >  1:232200/1‑134 (MQ=255)                                         cgTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCAtt                                                                                              >  1:331719/1‑140 (MQ=255)                                             aTGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAAcc                                                                                          >  1:667602/1‑140 (MQ=255)                                                 cacaTTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACaa                                                                                      <  1:575445/140‑1 (MQ=255)                                                      ttAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTc                                                                                 <  2:434294/140‑1 (MQ=255)                                                                             gggATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGa                                                          >  1:400876/1‑140 (MQ=255)                                                                             gggATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGa                                                          >  1:465546/1‑140 (MQ=255)                                                                                          gcctgACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTcagcag                                             >  1:73630/1‑140 (MQ=255)                                                                                              gACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGAc                                         <  2:465546/140‑1 (MQ=255)                                                                                                     gcgcAACTCGCCAGTATGGGGCTGAATATGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCt                                  >  2:559347/1‑140 (MQ=255)                                                                                                           cTCGCCAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGAt                            >  2:595130/1‑140 (MQ=255)                                                                                                               ccAGTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCt                        >  2:153685/1‑140 (MQ=255)                                                                                                                  gTATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAg                     <  1:302613/140‑1 (MQ=255)                                                                                                                   tATGGGGCTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAgg                    <  1:662948/140‑1 (MQ=255)                                                                                                                          cTGAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTAt             >  2:248059/1‑140 (MQ=255)                                                                                                                            gAATACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGa           >  2:132164/1‑140 (MQ=255)                                                                                                                             aaTACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAg          >  1:198695/1‑140 (MQ=255)                                                                                                                             aaTACGACTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAg          >  2:13335/1‑140 (MQ=255)                                                                                                                               tACGGCTTCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTa        >  1:344925/1‑140 (MQ=255)                                                                                                                                     ttCTGCCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATcgcg  <  1:559347/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       TCGCGTCGATATCTATTTGAATAATGGTTATATGGCAACGCGTGATGTCACATTTAATACGGGCGACAGTGAACAAGGGATTGTTCCCTGCCTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCG  >  NC_000913/4545329‑4545600

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