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

TTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382508‑1382601                                                |                                               ttGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              >  1:3626677/1‑50 (MQ=255)  tGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              >  1:1096944/1‑49 (MQ=255)   gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:1380210/48‑1 (MQ=255)   gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:5295573/48‑1 (MQ=255)   gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:1398917/1‑49 (MQ=255)   gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:5568807/1‑49 (MQ=255)    cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:4065664/1‑50 (MQ=255)     gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:911407/1‑49 (MQ=255)     gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:5101395/1‑49 (MQ=255)     gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:2928728/50‑1 (MQ=255)     gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:515890/50‑1 (MQ=255)      aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          >  1:1298702/1‑49 (MQ=255)      aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:2492927/50‑1 (MQ=255)       gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:1355367/50‑1 (MQ=255)        tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:5926642/50‑1 (MQ=255)        tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:2772329/50‑1 (MQ=255)          aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACg                                      <  1:6073010/49‑1 (MQ=255)            cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:4164529/1‑49 (MQ=255)            cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:2981578/1‑49 (MQ=255)             tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:1072579/1‑48 (MQ=255)             tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:537545/1‑48 (MQ=255)             tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:5236801/1‑48 (MQ=255)              attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:6002867/1‑50 (MQ=255)               ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGa                                  <  1:1624282/48‑1 (MQ=255)               ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGa                                  <  1:5697311/48‑1 (MQ=255)               ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 <  1:4126577/49‑1 (MQ=255)                 atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              >  1:136513/1‑50 (MQ=255)                 atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:2709595/49‑1 (MQ=255)                 atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:2216254/49‑1 (MQ=255)                  tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              <  1:2878730/49‑1 (MQ=255)                  tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              <  1:331996/49‑1 (MQ=255)                  tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              <  1:3586445/49‑1 (MQ=255)                   ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:5454733/1‑50 (MQ=255)                    tcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:6080067/1‑49 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCTCAATAACGGTGAGTTTTCCa                          <  1:5223901/49‑1 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:5435806/46‑1 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:6078788/46‑1 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          <  1:3351794/49‑1 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          <  1:449807/49‑1 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:5691697/50‑1 (MQ=255)                      gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:957043/50‑1 (MQ=255)                       aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        >  1:2317854/1‑50 (MQ=255)                       aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:2909277/50‑1 (MQ=255)                       aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:4385201/50‑1 (MQ=255)                        gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:3603560/50‑1 (MQ=255)                        gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:4651321/50‑1 (MQ=255)                         cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:6493093/1‑49 (MQ=255)                         cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:5762086/1‑49 (MQ=255)                          gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:3649822/1‑48 (MQ=255)                          gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCaa                     <  1:6083929/50‑1 (MQ=255)                           cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:3545302/50‑1 (MQ=255)                            ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:2552361/1‑49 (MQ=255)                            ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   <  1:2225985/50‑1 (MQ=255)                            ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   <  1:5689301/50‑1 (MQ=255)                             gCTTCAATATTGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:3588013/50‑1 (MQ=255)                             gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:4816593/50‑1 (MQ=255)                             gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:5218918/50‑1 (MQ=255)                              cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:1665456/50‑1 (MQ=255)                              cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:4712696/1‑50 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:1499900/49‑1 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:1563386/49‑1 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:6480497/49‑1 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:5482192/49‑1 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:6725557/48‑1 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:1203287/48‑1 (MQ=255)                                 cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:5614224/1‑50 (MQ=255)                                  aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:6060764/49‑1 (MQ=255)                                  aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:137600/49‑1 (MQ=255)                                   atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:4199594/50‑1 (MQ=255)                                   atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            >  1:5724677/1‑50 (MQ=255)                                   atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:4232835/50‑1 (MQ=255)                                    tatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgcg           >  1:5196060/1‑50 (MQ=255)                                      tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          <  1:6490638/49‑1 (MQ=255)                                      tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          <  1:2905481/49‑1 (MQ=255)                                      tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:1276535/1‑50 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:855931/49‑1 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:6062549/49‑1 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:2991446/49‑1 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:6285740/49‑1 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:3565979/50‑1 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  1:3058194/1‑50 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  1:1430766/1‑50 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  1:4374818/1‑50 (MQ=255)                                        gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:1954652/50‑1 (MQ=255)                                         acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      >  1:6255448/1‑50 (MQ=255)                                         acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      <  1:2599883/50‑1 (MQ=255)                                           gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:5825509/1‑50 (MQ=255)                                           gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    <  1:5214917/50‑1 (MQ=255)                                           gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:5545139/1‑50 (MQ=255)                                           gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:2815316/1‑50 (MQ=255)                                           gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:2620890/1‑50 (MQ=255)                                            acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:54333/50‑1 (MQ=255)                                             cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  >  1:438788/1‑50 (MQ=255)                                             cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:2546345/50‑1 (MQ=255)                                                |                                               TTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382508‑1382601

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