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

GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601                                              |                                               gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             <  1:5054113/49‑1 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             <  1:5057996/49‑1 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           <  1:72314/50‑1 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCa                                            <  1:485592/49‑1 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:6077501/1‑49 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:3543274/1‑49 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:4325150/1‑49 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:394677/50‑1 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:2334171/50‑1 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:4048378/50‑1 (MQ=255)     gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:2769931/50‑1 (MQ=255)     gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:2887340/50‑1 (MQ=255)       gAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      >  1:4571439/1‑50 (MQ=255)       gAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      >  1:4357223/1‑50 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACgg                                     >  1:4504277/1‑50 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:3873699/49‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:4315134/49‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:3118940/49‑1 (MQ=255)         aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    <  1:4946632/50‑1 (MQ=255)         aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    <  1:1428187/50‑1 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:2806762/1‑49 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:3644412/1‑49 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                >  1:4235393/1‑50 (MQ=255)              tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:2790276/1‑50 (MQ=255)              tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:4673971/1‑50 (MQ=255)              tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:3623588/1‑50 (MQ=255)                tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              <  1:6367801/49‑1 (MQ=255)                  tcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:1025854/1‑49 (MQ=255)                   cGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          >  1:766342/1‑50 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:5144375/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            <  1:3855983/47‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:1017234/1‑50 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:4347247/1‑50 (MQ=255)                      gcgcGGCTTCAATATGGAAGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:5006470/50‑1 (MQ=38)                       cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:2857251/1‑49 (MQ=255)                        gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:989973/1‑48 (MQ=255)                        gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:1192532/1‑48 (MQ=255)                         cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:8491/50‑1 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:4879702/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:561925/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:2476346/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   <  1:3603659/50‑1 (MQ=255)                           gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:389736/50‑1 (MQ=255)                           gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:2200750/50‑1 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:31988/1‑50 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:5599370/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:2355935/1‑50 (MQ=255)                              tCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               >  1:2280770/1‑50 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:4316417/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:4619065/1‑50 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:3637339/1‑50 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:1031147/49‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:3726939/49‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:6043566/49‑1 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            >  1:3797142/1‑50 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            >  1:3676674/1‑50 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:565506/50‑1 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           <  1:2594324/50‑1 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          <  1:5303528/49‑1 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         <  1:4736776/50‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:5518791/49‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:895150/49‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:822790/50‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:4564007/50‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:2941717/50‑1 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     <  1:3299066/50‑1 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     >  1:3880405/1‑50 (MQ=255)                                         gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:5617018/1‑50 (MQ=255)                                         gacAGTCGCAATAACGGCGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:3585462/1‑50 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:1958936/50‑1 (MQ=255)                                              |                                               GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601

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: ‑