Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	612,648	G→T	intergenic (‑154/‑89)	fepA ← / → fes	ferric enterobactin outer membrane transporter/enterochelin esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	612,648	0	G	T	100.0%	341.8 / NA	89	intergenic (‑154/‑89)	fepA/fes	ferric enterobactin outer membrane transporter/enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (32/57);  total (32/57)

GTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTT  >  NC_000913/612601‑612693                                                |                                              gtgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTActct                                             <  1:4791249/50‑1 (MQ=255)  tgcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTa                                            >  1:6032290/1‑50 (MQ=255)   gcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAt                                           >  1:425328/1‑50 (MQ=255)   gcagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAt                                           >  1:2443354/1‑50 (MQ=255)    cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAtt                                          <  1:926616/50‑1 (MQ=255)    cagtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTAt                                           <  1:2718583/49‑1 (MQ=255)     agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                         >  1:3714925/1‑50 (MQ=255)     agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                         <  1:1456452/50‑1 (MQ=255)     agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                         >  1:268859/1‑50 (MQ=255)     agtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtg                                         >  1:6484171/1‑50 (MQ=255)      gtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgc                                        <  1:99905/50‑1 (MQ=255)      gtTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgc                                        <  1:2182612/50‑1 (MQ=255)          gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                    <  1:3826960/50‑1 (MQ=255)          gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                    >  1:3020885/1‑50 (MQ=255)          gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                    >  1:1041874/1‑50 (MQ=255)          gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                    >  1:2067267/1‑50 (MQ=255)          gTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaat                                    <  1:906595/50‑1 (MQ=255)           ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                   <  1:5403117/50‑1 (MQ=255)           ttaattaaTGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                   <  1:1133810/50‑1 (MQ=255)            taattaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatg                                   >  1:6120569/1‑49 (MQ=255)              attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaa                                <  1:6520847/50‑1 (MQ=255)              attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaa                                <  1:4679103/50‑1 (MQ=255)              attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                 >  1:437949/1‑49 (MQ=255)              attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                 >  1:2758994/1‑49 (MQ=255)              attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                 <  1:1379565/49‑1 (MQ=255)              attaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgca                                 <  1:5098271/49‑1 (MQ=255)               ttaatGTCCGCGCTTCCCACGGCGCGCCATTACTCTATtgcaaatgcaaa                               <  1:2525950/50‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                             <  1:3244523/49‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                             <  1:4039588/49‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                             <  1:2972159/49‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                             <  1:2998052/49‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATa                             <  1:691738/49‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAg                            <  1:6593389/50‑1 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAg                            >  1:4040567/1‑50 (MQ=255)                  atGTCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAg                            <  1:2467081/50‑1 (MQ=255)                     tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGtt                          <  1:287191/49‑1 (MQ=255)                     tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGtt                          <  1:6480203/49‑1 (MQ=255)                     tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTa                         <  1:6403506/50‑1 (MQ=255)                     tCCGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTa                         <  1:3606320/50‑1 (MQ=255)                      ccGCGCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                        >  1:2465150/1‑50 (MQ=255)                       cgcgCTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTAt                        >  1:4596092/1‑49 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATc                       >  1:3230755/1‑49 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATc                       >  1:5367962/1‑49 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATc                       >  1:2058037/1‑49 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATc                       >  1:3738393/1‑49 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCa                      <  1:913179/50‑1 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCa                      <  1:1114320/50‑1 (MQ=255)                        gcgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCa                      <  1:2809459/50‑1 (MQ=255)                         cgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                     >  1:5450023/1‑50 (MQ=255)                         cgcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                     >  1:6563084/1‑50 (MQ=255)                          gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                    <  1:5759717/50‑1 (MQ=255)                          gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                    <  1:4622399/50‑1 (MQ=255)                          gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaat                    <  1:3365306/50‑1 (MQ=255)                          gcTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaa                     >  1:6149545/1‑49 (MQ=255)                           cTTCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaata                   <  1:5641191/50‑1 (MQ=255)                            ttCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataa                  >  1:5146376/1‑50 (MQ=255)                            ttCCCACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataa                  >  1:924196/1‑50 (MQ=255)                              cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCaataa                  <  1:5701904/48‑1 (MQ=255)                              cccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                >  1:237214/1‑50 (MQ=255)                               ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                <  1:2375544/49‑1 (MQ=255)                               ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAataata                <  1:3009978/49‑1 (MQ=255)                               ccACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAatat               <  1:2655324/50‑1 (MQ=255)                                cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt              <  1:6090062/50‑1 (MQ=255)                                cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt              <  1:2763307/50‑1 (MQ=255)                                cACGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatt              <  1:3236605/50‑1 (MQ=255)                                 aCGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtatta             >  1:5762008/1‑50 (MQ=255)                                  cGGCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat            <  1:418455/50‑1 (MQ=255)                                   ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAAtattat            >  1:1829604/1‑49 (MQ=255)                                   ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATc           <  1:4259319/50‑1 (MQ=255)                                   ggCGCGCCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATc           <  1:2392170/50‑1 (MQ=255)                                    gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa          <  1:323921/50‑1 (MQ=255)                                    gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa          <  1:6028028/50‑1 (MQ=255)                                    gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa          <  1:2581897/50‑1 (MQ=255)                                    gcgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCa          <  1:4063061/50‑1 (MQ=255)                                     cgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCaa         <  1:434640/50‑1 (MQ=255)                                     cgcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCaa         >  1:2662703/1‑50 (MQ=255)                                      gcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAat        <  1:3391739/50‑1 (MQ=255)                                      gcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAat        <  1:249031/50‑1 (MQ=255)                                      gcgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAat        <  1:6745328/50‑1 (MQ=255)                                       cgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAata       <  1:2785706/50‑1 (MQ=255)                                       cgcCATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAata       <  1:191650/50‑1 (MQ=255)                                         ccATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatata     >  1:3817801/1‑50 (MQ=255)                                         ccATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatata     >  1:3744122/1‑50 (MQ=255)                                          cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat    <  1:2453972/50‑1 (MQ=255)                                          cATTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAatatat    <  1:5997101/50‑1 (MQ=255)                                           aTTACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAtt   >  1:3982540/1‑50 (MQ=255)                                            ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt  >  1:3493267/1‑50 (MQ=255)                                            ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt  <  1:3439591/50‑1 (MQ=255)                                            ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt  <  1:1909017/50‑1 (MQ=255)                                            ttACTCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATAttt  <  1:4874824/50‑1 (MQ=255)                                                |                                              GTGCAGTTGGTTAATTAATGTCCGCGCTTCCCACGGCGCGCCATTACGCTATTGCAAATGCAAATAGTTATCAATAATATTATCAATATATTT  >  NC_000913/612601‑612693

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