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% 336.4 / NA 84S203S (AGC→AGTompGouter membrane porin G
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/46);  total (38/46)

TTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382508‑1382601
                                               |                                              
ttGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              >  1:628310/1‑50 (MQ=255)
ttGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              >  1:2533798/1‑50 (MQ=255)
ttGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              >  1:204789/1‑50 (MQ=255)
  gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCg                                              <  1:462393/48‑1 (MQ=255)
  gCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGc                                             >  1:5143860/1‑49 (MQ=255)
   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           <  1:1250383/50‑1 (MQ=255)
   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCaa                                           >  1:5131760/1‑50 (MQ=255)
   cGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCa                                            <  1:924037/49‑1 (MQ=255)
    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:4452228/50‑1 (MQ=255)
    gAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          <  1:2548533/50‑1 (MQ=255)
     aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAAt                                          >  1:1541315/1‑49 (MQ=255)
     aGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATa                                         <  1:1108116/50‑1 (MQ=255)
      gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:4242905/50‑1 (MQ=255)
      gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        >  1:3283400/1‑50 (MQ=255)
      gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        >  1:5345770/1‑50 (MQ=255)
      gTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATaa                                        <  1:4219097/50‑1 (MQ=255)
       tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:3414575/50‑1 (MQ=255)
       tGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAAc                                       <  1:4865731/50‑1 (MQ=255)
          aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    <  1:4571283/50‑1 (MQ=255)
          aCTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:4742150/1‑50 (MQ=255)
           cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:3870088/1‑49 (MQ=255)
           cTATTATCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTg                                   <  1:1853419/50‑1 (MQ=255)
            tattatCTCGAGCGCGGCTTCAATATGGACTACAGTCGCAATAACGGt                                    >  1:136044/1‑48 (MQ=255)
            tattatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGt                                    >  1:5304841/1‑48 (MQ=255)
             attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:5983378/1‑50 (MQ=255)
             attatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 >  1:485028/1‑50 (MQ=255)
              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGa                                  <  1:4562716/48‑1 (MQ=255)
              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAg                                 <  1:2688848/49‑1 (MQ=255)
              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                >  1:1611321/1‑50 (MQ=255)
              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                <  1:2951189/50‑1 (MQ=255)
              ttatCTCGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGt                                <  1:171569/50‑1 (MQ=255)
                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             >  1:2958219/1‑50 (MQ=255)
                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             >  1:3719510/1‑50 (MQ=255)
                 tctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGttt                              <  1:3710271/49‑1 (MQ=255)
                  ctcGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            >  1:3470306/1‑50 (MQ=255)
                    cGAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCa                          >  1:2206576/1‑50 (MQ=255)
                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGtttt                             <  1:1841464/46‑1 (MQ=255)
                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTc                            <  1:4814658/47‑1 (MQ=255)
                     gAGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCAc                         <  1:4966310/50‑1 (MQ=255)
                      aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:2276041/50‑1 (MQ=255)
                      aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:5049764/50‑1 (MQ=255)
                      aGCGCGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACg                        <  1:2251753/50‑1 (MQ=255)
                       gcgcGGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       <  1:2675140/50‑1 (MQ=255)
                         gcgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:4433686/1‑48 (MQ=255)
                          cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGc                       >  1:4466446/1‑47 (MQ=255)
                          cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:2494890/50‑1 (MQ=255)
                          cgGCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAg                    <  1:1964276/50‑1 (MQ=255)
                           ggCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGa                   <  1:782078/50‑1 (MQ=255)
                            gCTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaa                  <  1:1082787/50‑1 (MQ=255)
                             cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:4348470/50‑1 (MQ=255)
                             cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:1131529/1‑50 (MQ=255)
                             cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 <  1:5009940/50‑1 (MQ=255)
                             cTTCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGaaa                 >  1:1598687/1‑50 (MQ=255)
                               tCAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAtt               >  1:3696761/1‑50 (MQ=255)
                                cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAAt                <  1:5831556/48‑1 (MQ=255)
                                cAATATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTc              >  1:2365509/1‑50 (MQ=255)
                                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             >  1:4671712/1‑50 (MQ=255)
                                 aaTATGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcg             >  1:4923973/1‑50 (MQ=255)
                                  atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            >  1:5534987/1‑50 (MQ=255)
                                  atatGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTcgc            <  1:2249928/50‑1 (MQ=255)
                                    atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:2128630/1‑50 (MQ=255)
                                    atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:590683/1‑50 (MQ=255)
                                    atGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          >  1:2235456/1‑50 (MQ=255)
                                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCgcgc          <  1:5978082/49‑1 (MQ=255)
                                     tGGACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         <  1:2134551/50‑1 (MQ=255)
                                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:1543671/1‑49 (MQ=255)
                                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:4308914/1‑49 (MQ=255)
                                      ggACGACAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGcc         >  1:4651546/1‑49 (MQ=255)
                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCt        <  1:4070852/49‑1 (MQ=255)
                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:2513213/50‑1 (MQ=255)
                                       gacgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTa       <  1:3494263/50‑1 (MQ=255)
                                        acgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAt      >  1:2337820/1‑50 (MQ=255)
                                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     >  1:5549488/1‑50 (MQ=255)
                                         cgacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAtt     >  1:5849327/1‑50 (MQ=255)
                                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    <  1:1689512/50‑1 (MQ=255)
                                          gacAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTAttt    >  1:3055530/1‑50 (MQ=255)
                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:5328317/50‑1 (MQ=255)
                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:4644152/50‑1 (MQ=255)
                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:4276570/50‑1 (MQ=255)
                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:4339775/50‑1 (MQ=255)
                                           acAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTg   <  1:673130/50‑1 (MQ=255)
                                            cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:2028118/50‑1 (MQ=255)
                                            cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:3455719/50‑1 (MQ=255)
                                            cAGTCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGc  <  1:2994647/50‑1 (MQ=255)
                                               |                                              
TTGCGAGTGAACTATTATCTCGAGCGCGGCTTCAATATGGACGACAGCCGCAATAACGGTGAGTTTTCCACGCAAGAAATTCGCGCCTATTTGC  >  NC_000913/1382508‑1382601

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