Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,382,555	C→T	S203S (AGC→AGT)	ompG →	outer membrane porin G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,382,555	0	C	T	100.0%	344.5 / NA	88	S203S (AGC→AGT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (42/46);  total (42/46)

TGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTG  >  NC_000913/1382509‑1382600                                               |                                              tGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                             >  2:2394083/1‑49 (MQ=255)  gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                             <  2:645185/48‑1 (MQ=255)  gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                            <  1:4127329/49‑1 (MQ=255)  gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                            <  1:2489000/49‑1 (MQ=255)   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                          <  2:3989599/50‑1 (MQ=255)    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                          >  2:3227143/1‑49 (MQ=255)    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                          >  1:4514463/1‑49 (MQ=255)    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                          >  2:2291957/1‑49 (MQ=255)    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                         <  1:1611257/50‑1 (MQ=255)     aGTGAACTATTATCACGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                        <  2:3726682/50‑1 (MQ=39)      gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                       >  1:2635237/1‑50 (MQ=255)       tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                      <  1:1661733/50‑1 (MQ=255)       tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                      <  1:1145878/50‑1 (MQ=255)       tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                      >  2:872678/1‑50 (MQ=255)       tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                      <  1:825030/50‑1 (MQ=255)         aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                     <  1:3643978/49‑1 (MQ=255)          aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                   >  2:335862/1‑50 (MQ=255)          aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                   >  2:3976374/1‑50 (MQ=255)           cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                   >  2:2380126/1‑49 (MQ=255)              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                <  2:80927/49‑1 (MQ=255)              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                <  2:833101/49‑1 (MQ=255)              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                               >  2:3996490/1‑50 (MQ=255)              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                               >  2:1844691/1‑50 (MQ=255)                atCTCGAGCGCGGCTTCAATATGGACGACAGTCGTAATAACGGTGAGttt                             >  2:2489917/1‑50 (MQ=255)                atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                             >  1:3889696/1‑50 (MQ=255)                atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                              <  2:3435240/49‑1 (MQ=255)                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                            >  2:1303503/1‑50 (MQ=255)                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                            >  1:4472116/1‑50 (MQ=255)                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                             <  1:670118/49‑1 (MQ=255)                  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                           >  1:293686/1‑50 (MQ=255)                   gcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                          <  2:4073085/49‑1 (MQ=255)                    cGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         >  2:2614092/1‑50 (MQ=255)                    cGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         >  1:756811/1‑50 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                          <  1:579432/48‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         <  1:2234203/49‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         <  1:2128175/49‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         <  1:3387575/49‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                         <  1:652508/49‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                        >  1:3380875/1‑50 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                        <  2:1714424/50‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                        >  2:3535343/1‑50 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                        <  1:731955/50‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                        <  1:1506741/50‑1 (MQ=255)                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGCGAGTTTTCCAc                        >  1:1405604/1‑50 (MQ=255)                      aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                       <  2:945008/50‑1 (MQ=255)                      aGCGCGGCTGCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                       <  1:772271/50‑1 (MQ=255)                       gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                       >  2:3527255/1‑49 (MQ=255)                       gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                       >  2:3017857/1‑49 (MQ=255)                       gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                      <  2:3838494/50‑1 (MQ=255)                       gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                      <  1:2169330/50‑1 (MQ=255)                        cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCa                     <  1:3886304/50‑1 (MQ=255)                         gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                      >  1:2008735/1‑48 (MQ=255)                         gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                      >  2:4460890/1‑48 (MQ=255)                         gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCaa                    <  2:3680778/50‑1 (MQ=255)                          cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                      >  2:551213/1‑47 (MQ=255)                          cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                   <  1:4095956/50‑1 (MQ=255)                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                   >  2:1318592/1‑49 (MQ=255)                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                   >  1:4043319/1‑49 (MQ=255)                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                   >  2:2141821/1‑49 (MQ=255)                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                   >  2:1418942/1‑49 (MQ=255)                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                   >  2:409776/1‑49 (MQ=255)                            gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                 <  2:1995043/50‑1 (MQ=255)                             cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                >  2:4332648/1‑50 (MQ=255)                             cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                >  1:1570905/1‑50 (MQ=255)                              ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt               >  1:2568251/1‑50 (MQ=255)                                cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc             >  1:490943/1‑50 (MQ=255)                                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg            <  1:4176928/50‑1 (MQ=255)                                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg            >  1:1097363/1‑50 (MQ=255)                                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg            <  2:2699664/50‑1 (MQ=255)                                  atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc           <  1:1445752/50‑1 (MQ=255)                                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg          <  2:2827476/50‑1 (MQ=255)                                   tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg          <  2:3822896/50‑1 (MQ=255)                                    atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc         >  1:4071261/1‑50 (MQ=255)                                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc        >  2:2929749/1‑50 (MQ=255)                                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc        >  2:2154552/1‑49 (MQ=255)                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt       <  1:2035888/49‑1 (MQ=255)                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt       <  1:58568/49‑1 (MQ=255)                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa      >  1:2033616/1‑50 (MQ=255)                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa      >  2:3429679/1‑50 (MQ=255)                                        acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt     <  2:1519896/50‑1 (MQ=255)                                        acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt     >  2:2509570/1‑50 (MQ=255)                                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt    <  1:1576784/50‑1 (MQ=255)                                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt    <  2:324975/50‑1 (MQ=255)                                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt    <  2:2084446/50‑1 (MQ=255)                                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt   <  1:3343695/50‑1 (MQ=255)                                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt   <  2:1123339/50‑1 (MQ=255)                                           acAGTCGCAATAACGGTGAGTTTTCCGCGCAAGAAATTCGCGCCTTTTTg  <  1:1422521/50‑1 (MQ=255)                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg  <  1:1739919/50‑1 (MQ=255)                                               |                                              TGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTG  >  NC_000913/1382509‑1382600

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