Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,188,459	C→T	D818N (GAT→AAT)	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,459	0	C	T	100.0%	74.1 / NA	21	D818N (GAT→AAT)	yehB	putative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/5);  total (16/5)

CATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACT  >  NC_000913/2188382‑2188535                                                                              |                                                                             catTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAAt                                                                           >  2:382205/1‑81 (MQ=255)      atatAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCtt                                                                      >  1:382198/1‑81 (MQ=255)          aaaTAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTAcc                                                                  >  1:382196/1‑81 (MQ=255)            aTAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCg                                                                 >  1:382216/1‑80 (MQ=255)               gCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAgg                                                             >  2:382197/1‑81 (MQ=255)                  ttATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGa                                                          <  2:382182/81‑1 (MQ=255)                    aTGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGAt                                                         <  2:382183/80‑1 (MQ=255)                              ccATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGaa                                               <  2:382184/80‑1 (MQ=255)                                  aGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTcc                                          >  1:382218/1‑81 (MQ=255)                                     accaccTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCttg                                       >  2:382203/1‑81 (MQ=255)                                         ccTGACTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGc                                   >  1:382202/1‑81 (MQ=255)                                              cTATTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTAt                               >  1:382210/1‑80 (MQ=255)                                                aTTACTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCaa                            >  1:382200/1‑81 (MQ=255)                                                   gcTGGCAAATATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAAtt                          <  1:382185/79‑1 (MQ=255)                                                       gCAAATATAATTTCTACTTTCATTAATCTCTTTACTGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCca                      >  1:382215/1‑80 (MQ=255)                                                          aaTATAATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACAt                   >  2:382201/1‑80 (MQ=255)                                                              taATTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAAcc              <  2:382186/81‑1 (MQ=255)                                                               aaTTTCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAAcc              >  1:382207/1‑80 (MQ=255)                                                                   tCTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATgg         >  2:382217/1‑81 (MQ=255)                                                                    cTACTTTCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGc        >  2:382206/1‑81 (MQ=255)                                                                         ttCATTAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTAc   >  2:382209/1‑81 (MQ=255)                                                                          tcattaATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACt  >  2:382204/1‑81 (MQ=255)                                                                              |                                                                             CATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACT  >  NC_000913/2188382‑2188535

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