Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,588,191	C→A	A36S (GCT→TCT)	ydeR ←	putative fimbrial protein YdeR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,588,191	0	C	A	100.0%	195.2 / NA	56	A36S (GCT→TCT)	ydeR	putative fimbrial protein YdeR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (33/23);  total (33/23)

GTTTCAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGCGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAATTCAAGGGTA  >  NC_000913/1588055‑1588327                                                                                                                                         |                                                                                                                                         gTTTCAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGAc                                                                                                                                       >  2:473092/1‑140 (MQ=255)   ttCAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTa                                                                                                                                     <  2:614123/140‑1 (MQ=255)    tCAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTAc                                                                                                                                    >  1:1194738/1‑140 (MQ=255)     cAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTAcc                                                                                                                                   <  1:473092/140‑1 (MQ=255)     cAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTAcc                                                                                                                                   <  1:467916/140‑1 (MQ=255)       aCCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCg                                                                                                                                 <  2:319537/140‑1 (MQ=255)             cAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCAtt                                                                                                                           <  2:477608/140‑1 (MQ=255)                              aGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACAt                                                                                                          >  1:758121/1‑140 (MQ=255)                              aGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACAt                                                                                                          >  2:755279/1‑140 (MQ=255)                                tAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATcg                                                                                                        <  1:755279/140‑1 (MQ=255)                                tAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATcg                                                                                                        <  2:804873/140‑1 (MQ=255)                                 aaTATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATcgg                                                                                                       <  2:366141/140‑1 (MQ=255)                                  atatTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATcggc                                                                                                      >  1:653214/1‑140 (MQ=255)                                   tatTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATcggcg                                                                                                     >  1:1191288/1‑140 (MQ=255)                                     ttGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCggcggc                                                                                                   <  1:44918/140‑1 (MQ=255)                                        tgCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTg                                                                                                >  2:112091/1‑140 (MQ=255)                                        tgCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTg                                                                                                >  1:765956/1‑140 (MQ=255)                                        tgCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTg                                                                                                >  1:1227400/1‑140 (MQ=255)                                        tgCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTg                                                                                                >  2:338402/1‑140 (MQ=255)                                        tgCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTg                                                                                                >  1:1205240/1‑140 (MQ=255)                                         gccagccaGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTgg                                                                                               >  2:709096/1‑140 (MQ=255)                                         gccagccaGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTgg                                                                                               >  2:111042/1‑140 (MQ=255)                                          ccagccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGga                                                                                              >  2:371426/1‑140 (MQ=255)                                          ccagccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGga                                                                                              >  1:1188411/1‑140 (MQ=255)                                          ccagccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGga                                                                                              >  2:939743/1‑140 (MQ=255)                                          ccagccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGga                                                                                              >  1:643333/1‑140 (MQ=255)                                          ccagccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGga                                                                                              >  1:1137838/1‑140 (MQ=255)                                           cagccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGgag                                                                                             >  2:405718/1‑140 (MQ=255)                                             gccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGagag                                                                                           <  2:620319/140‑1 (MQ=255)                                             gccagATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGagag                                                                                           <  2:379497/140‑1 (MQ=255)                                                              gttgttGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAAcgcg                                                                          <  2:758121/140‑1 (MQ=255)                                                                  ttgCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCaaa                                                                      <  2:996490/140‑1 (MQ=255)                                                                  ttgCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCaaa                                                                      >  1:21016/1‑140 (MQ=255)                                                                  ttgCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCaaa                                                                      >  2:827459/1‑140 (MQ=255)                                                                        aTTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAg                                                                <  2:928349/140‑1 (MQ=255)                                                                            cgcgTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTaa                                                            >  2:1005190/1‑140 (MQ=255)                                                                              cgTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACa                                                          <  1:586381/140‑1 (MQ=255)                                                                                tataAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAgg                                                        <  1:1005190/140‑1 (MQ=255)                                                                                 ataAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGa                                                       >  2:1059446/1‑140 (MQ=255)                                                                                      gATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACctc                                                  >  2:385478/1‑140 (MQ=255)                                                                                         ccccGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTtg                                               >  2:325583/1‑140 (MQ=255)                                                                                         ccccGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTtg                                               >  2:253043/1‑140 (MQ=255)                                                                                          cccGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTtgt                                              >  1:258824/1‑140 (MQ=255)                                                                                           ccGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTtgtg                                             <  1:709096/140‑1 (MQ=255)                                                                                               gaTTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGt                                         <  1:1223476/140‑1 (MQ=255)                                                                                                    aCGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCtttt                                    <  1:112091/140‑1 (MQ=255)                                                                                                     cGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTc                                   <  2:1194738/140‑1 (MQ=255)                                                                                                           cttcttTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGTTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCt                             <  2:1191288/140‑1 (MQ=255)                                                                                                           cttcttTGGTTTGAATAGTGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCt                             >  2:1079318/1‑140 (MQ=255)                                                                                                                             tGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAAt           >  1:683946/1‑140 (MQ=255)                                                                                                                             tGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAAt           >  1:1086113/1‑140 (MQ=255)                                                                                                                             tGCAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAAt           >  1:341184/1‑140 (MQ=255)                                                                                                                               cAGGGTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAATTc         <  1:939743/140‑1 (MQ=255)                                                                                                                                  ggTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAATTCAAg      >  1:408339/1‑140 (MQ=255)                                                                                                                                  ggTTTAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTAACTCTGATTATTGAAATTCAAg      >  1:204220/1‑140 (MQ=255)                                                                                                                                      tAGAGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAATTCAAGGGTa  <  1:1059446/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                         GTTTCAACCGGACAATCGGTTAATGACAAAGTAATATTGTGCCAGCCAGATGCAGAACCAGGTTGTTGCAGATTGCGCGTATAAAGATCCCCGAGATTAACGTTAGCTTCTTTGGTTTGAATAGTGCAGGGTTTAGCGACTACCCGACCATTAACAGTGATAGTGACATCGGCGGCCTGGAGAGTTGCTGTGAATAACGCGCAAAACGTAGCTAACAGGAACCTCTTGTGAAGTCTTTTCATATCTGACTCTGATTATTGAAATTCAAGGGTA  >  NC_000913/1588055‑1588327

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