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

GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601                                              |                                               gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:4432073/48‑1 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:3017752/48‑1 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:443179/1‑49 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             <  1:3847317/49‑1 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:3132465/1‑50 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:3367539/1‑50 (MQ=255)  cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCa                                            <  1:1065693/49‑1 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:3925023/1‑49 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:575860/50‑1 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:5228327/50‑1 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:5886596/50‑1 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       >  1:4467221/1‑50 (MQ=255)       gAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      >  1:422587/1‑50 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACgg                                     >  1:7033547/1‑50 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:1528709/49‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:2660/49‑1 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:6529094/1‑49 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTg                                   <  1:5029148/50‑1 (MQ=255)            attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:6037551/1‑50 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 <  1:1230470/49‑1 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:1106822/49‑1 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:43973/49‑1 (MQ=255)                tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              <  1:2630794/49‑1 (MQ=255)                 ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:51158/1‑50 (MQ=255)                  tcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:4318073/1‑49 (MQ=255)                   cGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          >  1:680026/1‑50 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:6453631/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:5580569/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:2506078/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:6634402/50‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:7030583/1‑50 (MQ=255)                     aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        >  1:4582797/1‑50 (MQ=255)                      gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:4787933/50‑1 (MQ=255)                       cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCa                      <  1:498057/50‑1 (MQ=255)                        gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:796583/1‑48 (MQ=255)                         cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:4434034/1‑47 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:3512137/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:2702146/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:1102292/1‑49 (MQ=255)                           gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:3054331/50‑1 (MQ=255)                           gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:1845971/50‑1 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:6728315/50‑1 (MQ=255)                            cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:906961/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:264300/1‑50 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:6144626/1‑50 (MQ=255)                              tCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               >  1:4708160/1‑50 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:6083745/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:6116936/48‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:4863196/1‑50 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:7078487/1‑50 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             >  1:5560916/1‑50 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:5636275/49‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:6604679/49‑1 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            >  1:2704986/1‑50 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:3663113/50‑1 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           <  1:4859384/50‑1 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           <  1:685709/50‑1 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           >  1:4117414/1‑50 (MQ=255)                                   atGGACGACAGTCGCAATAACGGTGAGCTTTCCACGCAAGAAATTCgcgc          >  1:4487644/1‑50 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:5603163/1‑50 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:4567665/1‑50 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:5803903/50‑1 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     <  1:1188551/50‑1 (MQ=255)                                        cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     >  1:110310/1‑50 (MQ=255)                                         gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    <  1:6464701/50‑1 (MQ=255)                                          acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   >  1:6711141/1‑50 (MQ=255)                                          acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:6349211/50‑1 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:6657437/50‑1 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  >  1:3466312/1‑50 (MQ=255)                                              |                                               GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601

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