Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,961	A→G	V650V (GTT→GTC)	yehB ←	putative fimbrial usher protein YehB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,188,961	0	A	G	100.0%	87.4 / NA	26	V650V (GTT→GTC)	yehB	putative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/12);  total (14/12)

TCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGG  >  NC_000913/2188891‑2189037                                                                       |                                                                             tcaCAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTc                                                                    <  1:387368/81‑1 (MQ=255)   acaGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTcc                                                                   <  2:387369/80‑1 (MQ=255)     aGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGc                                                                <  2:387370/81‑1 (MQ=255)        aaaCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTg                                                             >  2:387402/1‑81 (MQ=255)        aaaCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTg                                                             >  2:387410/1‑81 (MQ=255)        aaaCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTg                                                             >  1:387399/1‑81 (MQ=255)          aCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTc                                                           <  1:387371/81‑1 (MQ=255)              ttCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGAt                                                        >  2:387398/1‑80 (MQ=255)               tCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCTCCCTGAGACACTGGCTCCAGCCTGTCGAt                                                        >  2:387416/1‑79 (MQ=255)                 ggAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAg                                                     >  1:387405/1‑80 (MQ=255)                      gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTc                                                >  1:387409/1‑80 (MQ=255)                       aCGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCga                                              <  1:387372/81‑1 (MQ=255)                           ttCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCgactga                                          >  1:387397/1‑81 (MQ=255)                             cGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCgactgac                                         <  2:387373/80‑1 (MQ=255)                             cGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCgactgac                                         >  2:387407/1‑80 (MQ=255)                              gCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACta                                       <  2:387374/81‑1 (MQ=255)                                       aCGCCACCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCa                              >  2:387408/1‑81 (MQ=255)                                            aCCCGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCAc                         <  1:387375/81‑1 (MQ=255)                                               cGACCAGGCGACAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGt                      <  2:387376/81‑1 (MQ=255)                                                gACCAGGCGACAATGCCCACTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTc                     >  2:387418/1‑81 (MQ=255)                                                   cAGGCGACAATGCCCGCTGAGACGCTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATGGCGGCCATTCACTGTCGCa                  <  2:387377/81‑1 (MQ=255)                                                         aCAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGgc            >  2:387419/1‑81 (MQ=255)                                                          cAATGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGgcg           <  2:387378/81‑1 (MQ=255)                                                           aaTGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGgcgc          >  2:387417/1‑81 (MQ=255)                                                            aTGCCCCCTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGcgcg         <  1:387379/81‑1 (MQ=255)                                                                   cTGAGACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAgg  >  2:387412/1‑81 (MQ=255)                                                                       |                                                                             TCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGG  >  NC_000913/2188891‑2189037

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