Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,548,130	C→A	P130Q (CCA→CAA)	fimF →	type 1 fimbriae minor subunit FimF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,548,130	0	C	A	100.0%	83.0 / NA	25	P130Q (CCA→CAA)	fimF	type 1 fimbriae minor subunit FimF
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (11/14);  total (11/14)

TTGGGTTTACTGGCGTTGCAGATAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACT  >  NC_000913/4548001‑4548257                                                                                                                                  |                                                                                                                                ttGGGTTTACTGGCGTTGCAGATAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCcgcg                                                                                                                       <  1:405177/140‑1 (MQ=255)       ttACTGGCGTTGCAGATAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCg                                                                                                                 >  2:471408/1‑140 (MQ=255)       ttACTGGCGTTGCAGATAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCg                                                                                                                 <  1:474805/140‑1 (MQ=255)                      aTAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACg                                                                                                  <  1:424590/140‑1 (MQ=255)                      aTAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACg                                                                                                  <  2:115240/140‑1 (MQ=255)                          ccACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCgg                                                                                              >  1:76542/1‑140 (MQ=255)                             caATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTa                                                                                           >  1:276661/1‑140 (MQ=255)                             caATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTa                                                                                           >  2:197833/1‑140 (MQ=255)                                  ccAACCTGCTTGCACTTGAAAATACGGGGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCa                                                                                      >  2:97055/1‑140 (MQ=255)                                          cTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCt                                                                              >  1:443222/1‑140 (MQ=255)                                                           ggTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGc                                                             >  1:432286/1‑140 (MQ=255)                                                             tgtCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTa                                                           >  1:57670/1‑140 (MQ=255)                                                                cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAACCCAAATCCGCTGAATTTTTACGCCCGGCTAATg                                                        <  2:276661/140‑1 (MQ=255)                                                                       ttCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGacac                                                 <  1:97055/140‑1 (MQ=255)                                                                         cGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAg                                               <  2:122312/140‑1 (MQ=255)                                                                           ggACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGt                                             <  2:57670/140‑1 (MQ=255)                                                                                      aCAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTg                                  <  1:190831/140‑1 (MQ=255)                                                                                        aGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCg                                <  1:89953/140‑1 (MQ=255)                                                                                           ttCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCgggg                             >  1:453186/1‑140 (MQ=255)                                                                                               gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCata                         >  2:220820/1‑140 (MQ=255)                                                                                                      cagcaAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGc                  <  2:443222/140‑1 (MQ=255)                                                                                                        gcaAAATCAAATACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCa                <  1:81464/140‑1 (MQ=255)                                                                                                                  aTACCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCtt      <  2:453186/140‑1 (MQ=255)                                                                                                                    aCCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCa    <  1:42501/140‑1 (MQ=255)                                                                                                                    aCCCCTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCa    <  1:332164/140‑1 (MQ=255)                                                                                                                      cccTTAATGCTCAATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCAct  >  1:285805/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                TTGGGTTTACTGGCGTTGCAGATAGCCACAATGCCAACCTGCTTGCACTTGAAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGCTACCTTCACT  >  NC_000913/4548001‑4548257

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