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

GCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382510‑1382601                                              |                                               gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:1048561/48‑1 (MQ=255) gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:3228057/48‑1 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:3432554/1‑49 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:2689171/1‑49 (MQ=255)   gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:628384/50‑1 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          >  1:4949017/1‑49 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          >  1:833571/1‑49 (MQ=255)    aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:5426157/50‑1 (MQ=255)     gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:1033143/50‑1 (MQ=255)     gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:3174739/50‑1 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       >  1:3810924/1‑50 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:3360174/50‑1 (MQ=255)      tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:2817052/50‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACgg                                     <  1:6777316/50‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACgg                                     >  1:1906731/1‑50 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:6213536/49‑1 (MQ=255)        aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:3578466/49‑1 (MQ=255)         aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:1218534/1‑50 (MQ=255)         aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:7026924/1‑50 (MQ=255)          cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:4324320/1‑49 (MQ=255)             ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 <  1:3103118/49‑1 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:2969567/49‑1 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:708433/49‑1 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:4447814/49‑1 (MQ=255)               atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:4534019/49‑1 (MQ=255)                tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             >  1:4679461/1‑50 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:2982217/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:5503762/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:1936900/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:3833202/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:6572676/46‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                           <  1:6425226/48‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                           <  1:6624509/48‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          <  1:5150868/49‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          <  1:5846426/49‑1 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:7456409/1‑50 (MQ=255)                    gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:6500678/1‑50 (MQ=255)                     aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:7125650/50‑1 (MQ=255)                     aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        >  1:2025581/1‑50 (MQ=255)                       cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:4216343/1‑49 (MQ=255)                       cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:1645533/1‑49 (MQ=255)                       cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:6342077/1‑49 (MQ=255)                       cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCa                      <  1:1422365/50‑1 (MQ=255)                        gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:669233/1‑48 (MQ=255)                        gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:5364818/1‑48 (MQ=255)                         cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:1411209/1‑47 (MQ=255)                         cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:3750530/50‑1 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:6025089/1‑49 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   <  1:617265/50‑1 (MQ=255)                          ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   <  1:5532051/50‑1 (MQ=255)                           gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:5360778/50‑1 (MQ=255)                             ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:4197968/1‑50 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:2181850/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:4066066/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:6298468/49‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:3090092/48‑1 (MQ=255)                               cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:242725/48‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             <  1:4512558/50‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:5517658/49‑1 (MQ=255)                                aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:3956580/49‑1 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:4161587/50‑1 (MQ=255)                                 atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:2297533/50‑1 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           <  1:6850780/50‑1 (MQ=255)                                  tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           <  1:1973519/50‑1 (MQ=255)                                    tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         <  1:916368/50‑1 (MQ=255)                                     ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:5135502/1‑49 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:6693460/50‑1 (MQ=255)                                      gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  1:5912814/1‑50 (MQ=255)                                       acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      <  1:1817605/50‑1 (MQ=255)                                         gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    <  1:7147058/50‑1 (MQ=255)                                         gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:2678369/1‑50 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:436508/50‑1 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:3814861/50‑1 (MQ=255)                                           cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  >  1:129412/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: ‑