Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	150,562	A→G	I131T (ATA→ACA)	yadC ←	fimbrial tip‑adhesin YadC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	150,562	0	A	G	100.0%	111.0 / NA	33	I131T (ATA→ACA)	yadC	fimbrial tip‑adhesin YadC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

GGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATAGA  >  NC_000913/150485‑150637                                                                              |                                                                            ggttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGtt                                                                           <  1:36325/80‑1 (MQ=255)  gttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTaa                                                                         >  1:36351/1‑81 (MQ=255)   ttggAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAAc                                                                        <  1:36326/81‑1 (MQ=255)               ccGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATa                                                             >  2:36362/1‑80 (MQ=255)               ccGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATa                                                             >  1:36373/1‑80 (MQ=255)                 gATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAg                                                          <  2:36327/81‑1 (MQ=255)                  atatAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGa                                                         >  1:36354/1‑81 (MQ=255)                    ataGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGAc                                                        >  2:36350/1‑80 (MQ=255)                       gATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAgg                                                    <  2:36328/81‑1 (MQ=255)                        aTTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGa                                                   <  1:36329/81‑1 (MQ=255)                          tCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGaa                                                  >  2:36349/1‑80 (MQ=255)                            cTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGaa                                                  >  1:36356/1‑78 (MQ=255)                            cTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAg                                               <  2:36330/81‑1 (MQ=255)                              ggCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGAt                                             >  2:36368/1‑81 (MQ=255)                               gCGTTTTAATGTCAGTTTTTGTTTCGTAGGCAGACCAGACCCTTGATGTAAACAGGGTGAAAAACAGACAAGAAACAGATg                                            <  2:36334/81‑1 (MQ=255)                                cGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGt                                           <  1:36331/81‑1 (MQ=255)                                 gATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTa                                          >  2:36374/1‑81 (MQ=255)                                 gATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTa                                          >  1:36379/1‑81 (MQ=255)                                     gAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTa                                       >  2:36372/1‑80 (MQ=255)                                      aaTGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTaa                                      <  2:36332/80‑1 (MQ=255)                                       aTGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAAt                                    <  1:36333/81‑1 (MQ=255)                                           cAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGtt                                >  1:36360/1‑81 (MQ=255)                                           cAGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGtt                                >  2:36365/1‑81 (MQ=255)                                            aGTTATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGttt                               <  2:36335/81‑1 (MQ=255)                                               tATTGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATg                            >  2:36367/1‑81 (MQ=255)                                                 ttGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTc                          <  2:36336/81‑1 (MQ=255)                                                  tGTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTcc                         >  2:36364/1‑81 (MQ=255)                                                   gTATCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCg                        >  1:36371/1‑81 (MQ=255)                                                      tCGTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGcc                      <  1:36337/80‑1 (MQ=255)                                                        gTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCAt                    <  2:36339/80‑1 (MQ=255)                                                        gTATGCAGACCAGACCCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCAt                    <  2:36338/80‑1 (MQ=255)                                                               gaccagaccCTTGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTa            >  2:36369/1‑81 (MQ=255)                                                                         ttGATGTTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATAGa  >  1:36377/1‑81 (MQ=255)                                                                              |                                                                            GGTTGGAAGGATCTCCGATATAGATTCCTGGCGATTGAATGTCAGTTATTGTATCGTATGCAGACCAGACCCTTGATATTAACATGGTGTAATACAGACCAGGAACAGATGTATTAAATAGTTTATGTCCGCCATAATCTTTACCTGAATAGA  >  NC_000913/150485‑150637

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