Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,685	Δ1 bp	coding (2226/2481 nt)	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,682	0	A	.	100.0%	109.7 / NA	26	coding (2229/2481 nt)	yehB	putative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/10);  total (16/10)

TTGACTTCATAACCAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATCAAAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCTTCGCTATC  >  NC_000913/2188602‑2188756                                                                                 |                                                                           ttGACTTCATAACCAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATCa                                                                            <  1:382201/81‑1 (MQ=255)   gACTTCATAACCAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATCaa                                                                           <  2:382202/80‑1 (MQ=255)     cTTCATAACCAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATCaaatt                                                                        <  1:382203/81‑3 (MQ=255)        cATAACCAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAAc                                                                    >  1:382241/1‑81 (MQ=255)             ccAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTa                                                               >  2:382246/1‑81 (MQ=255)               aaaCGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTAc                                                              >  1:382243/1‑80 (MQ=255)               aaaCGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTAc                                                              >  1:382233/1‑80 (MQ=255)               aaaCGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACa                                                             <  1:382204/81‑1 (MQ=255)                           ttGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGa                                                 >  1:382238/1‑81 (MQ=255)                             gCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGAt                                                >  2:382232/1‑80 (MQ=255)                             gCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGAt                                                <  1:382205/80‑1 (MQ=255)                                 aTCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAggg                                           >  1:382239/1‑81 (MQ=255)                                     gCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGggcgg                                       >  1:382253/1‑81 (MQ=255)                                           aaCGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAAt                                   <  1:382206/79‑1 (MQ=255)                                             cGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTc                               <  2:382207/81‑1 (MQ=255)                                             cGCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTc                               <  2:382208/81‑1 (MQ=255)                                              gCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTcc                              >  1:382242/1‑81 (MQ=255)                                              gCTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTcc                              >  1:382248/1‑81 (MQ=255)                                               cTTTTATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCg                             <  1:382209/81‑1 (MQ=255)                                                  ttATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGtt                          >  2:382237/1‑81 (MQ=255)                                                  ttATAAACCATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGt                           <  2:382210/80‑1 (MQ=255)                                                         ccATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGt                   <  1:382211/81‑1 (MQ=255)                                                          cATGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGt                   >  1:382249/1‑80 (MQ=255)                                                            tGGCTTGCGCTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAAt                >  1:382251/1‑81 (MQ=255)                                                                  gcgcTGATCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCt          >  1:382261/1‑81 (MQ=255)                                                                         tCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCTTCGCTAt   <  2:382213/81‑1 (MQ=255)                                                                         tCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCTTCGCTAt   <  2:382212/81‑1 (MQ=255)                                                                         tCGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCTTCGCTAt   >  2:382268/1‑81 (MQ=255)                                                                          cGGTATC‑AAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCTTCGCTATc  >  2:382255/1‑81 (MQ=255)                                                                                 |                                                                           TTGACTTCATAACCAAACGTTAATGATTGCCCATCTGCTCTTAACGCTTTTATAAACCATGGCTTGCGCTGATCGGTATCAAAATTAACCAGTACAACCGCGCCGCGATAAGGGGCGGCAATTTTCCGGTTGCCACGTAATTCTGCTTCGCTATC  >  NC_000913/2188602‑2188756

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