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_0009131,382,5550CT100.0% 261.4 / NA 63S203S (AGC→AGTompGouter membrane porin G
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/28);  total (35/28)

CGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382511‑1382601
                                            |                                              
cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:2289868/1‑50 (MQ=255)
cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:3586264/1‑50 (MQ=255)
 gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:570782/1‑49 (MQ=255)
 gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:4356761/50‑1 (MQ=255)
  aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          >  1:2707305/1‑49 (MQ=255)
  aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:4093479/50‑1 (MQ=255)
   gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        >  1:161309/1‑50 (MQ=255)
   gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:2245003/50‑1 (MQ=255)
   gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        >  1:2496978/1‑50 (MQ=255)
   gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:2133889/50‑1 (MQ=255)
   gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:3746915/50‑1 (MQ=255)
    tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       >  1:3863236/1‑50 (MQ=255)
    tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:4346888/50‑1 (MQ=255)
      aaCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACgg                                     >  1:940149/1‑50 (MQ=255)
        cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:3059794/1‑49 (MQ=255)
        cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTg                                   <  1:1000767/50‑1 (MQ=255)
          attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:1442652/1‑50 (MQ=255)
           ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                <  1:130471/50‑1 (MQ=255)
           ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                <  1:1131335/50‑1 (MQ=255)
            tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:1816585/1‑50 (MQ=255)
            tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:282042/1‑50 (MQ=255)
            tatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               >  1:2367372/1‑50 (MQ=255)
             atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              >  1:2191236/1‑50 (MQ=255)
             atCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtt                               <  1:1397525/49‑1 (MQ=255)
               ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:2643149/1‑50 (MQ=255)
               ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:2505512/1‑50 (MQ=255)
                  gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTcc                           <  1:2381950/48‑1 (MQ=255)
                  gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:3903862/1‑50 (MQ=255)
                  gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:3892637/1‑50 (MQ=255)
                  gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         >  1:3605796/1‑50 (MQ=255)
                  gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:3600972/50‑1 (MQ=255)
                  gAGCGCGGCTTCAATATGGACGACAATCGCAATAACGGTGAGTTTTcc                           <  1:1419925/48‑1 (MQ=38)
                   aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:3094849/50‑1 (MQ=255)
                     cgcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:2814973/1‑49 (MQ=255)
                       cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:3347565/50‑1 (MQ=255)
                        ttcTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:3931743/3‑49 (MQ=255)
                        ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    >  1:4224693/1‑49 (MQ=255)
                         gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:4450059/50‑1 (MQ=255)
                          cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:2375986/50‑1 (MQ=255)
                          cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:3322769/1‑50 (MQ=255)
                           ttCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                >  1:2997689/1‑50 (MQ=255)
                             cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:96835/49‑1 (MQ=255)
                             cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               <  1:1000572/49‑1 (MQ=255)
                             cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:2164783/48‑1 (MQ=255)
                             cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:3258527/1‑50 (MQ=255)
                              aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:2285468/49‑1 (MQ=255)
                              aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:1592376/49‑1 (MQ=255)
                              aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              <  1:1357049/49‑1 (MQ=255)
                               atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:3545225/50‑1 (MQ=255)
                               atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:62930/50‑1 (MQ=255)
                                 atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:949799/1‑50 (MQ=255)
                                 atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:1948608/1‑50 (MQ=255)
                                 atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:2030045/1‑50 (MQ=255)
                                 atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:1405117/1‑50 (MQ=255)
                                    gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  1:151836/1‑50 (MQ=255)
                                    gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       >  1:1028440/1‑50 (MQ=255)
                                     acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      >  1:2093679/1‑50 (MQ=255)
                                       gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:3111416/1‑50 (MQ=255)
                                       gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    <  1:119515/50‑1 (MQ=255)
                                        acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   >  1:2031785/1‑50 (MQ=255)
                                        acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   >  1:2991463/1‑50 (MQ=255)
                                         cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:3085091/50‑1 (MQ=255)
                                         cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:1780787/50‑1 (MQ=255)
                                            |                                              
CGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382511‑1382601

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