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%	325.2 / NA	80	S203S (AGC→AGT)	ompG	outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/45);  total (35/45)

GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601                                              |                                               gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:2435302/48‑1 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:1292323/1‑49 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:5488555/1‑49 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:1945886/1‑49 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:2532000/1‑50 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           <  1:4483575/50‑1 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCa                                            <  1:4292642/49‑1 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:3852897/50‑1 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          >  1:2615125/1‑49 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:5465657/50‑1 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:4821610/50‑1 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:4603151/50‑1 (MQ=255)     gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:5268582/50‑1 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       >  1:1077150/1‑50 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:2432659/50‑1 (MQ=255)       gAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      >  1:2405008/1‑50 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:3581355/49‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:4204453/49‑1 (MQ=255)         aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:2049860/1‑50 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTg                                   <  1:2233598/50‑1 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTg                                   <  1:2299328/50‑1 (MQ=255)            attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:2480946/1‑50 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGa                                  <  1:2440460/48‑1 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                >  1:3456257/1‑50 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                <  1:284770/50‑1 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                <  1:4174739/50‑1 (MQ=255)              tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:3270815/1‑50 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:3012594/49‑1 (MQ=255)                 ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:1382064/1‑50 (MQ=255)                  tcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:4691564/1‑49 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:4008794/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:197687/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          <  1:648990/49‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          <  1:3289167/49‑1 (MQ=255)                     aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:4187648/50‑1 (MQ=255)                     aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:4163146/50‑1 (MQ=255)                      gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:676132/50‑1 (MQ=255)                         cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:979981/1‑47 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:2562345/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:621572/1‑49 (MQ=255)                           gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   >  1:518026/1‑49 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:3166601/50‑1 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:3077908/50‑1 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:5034194/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:3015056/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:4255643/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:2631213/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:2799348/1‑50 (MQ=255)                              tCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               >  1:5538100/1‑50 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:1146059/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:4054633/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:4383235/48‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:4286683/48‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:4395420/48‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:1556049/1‑50 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             <  1:2093688/50‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             <  1:3817769/50‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             >  1:5220769/1‑50 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:130196/49‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:366536/49‑1 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            >  1:2623959/1‑50 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           >  1:3742537/1‑50 (MQ=255)                                   atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:1353916/1‑50 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:4411886/1‑50 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:5473564/1‑50 (MQ=255)                                     ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:4099228/1‑49 (MQ=255)                                     ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:3163148/1‑49 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:2766800/49‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:3527031/49‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:1729843/49‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:1848830/50‑1 (MQ=255)                                       acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      <  1:2760072/50‑1 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     <  1:1090333/50‑1 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     >  1:93168/1‑50 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     <  1:3177992/50‑1 (MQ=255)                                         gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:2836954/1‑50 (MQ=255)                                          acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:2348952/50‑1 (MQ=255)                                          acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:3436168/50‑1 (MQ=255)                                          acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   >  1:3148951/1‑50 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:2406698/50‑1 (MQ=255)                                              |                                               GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601

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