Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,851,762 C→T G404G (GGC→GGT sppA → protease IV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,851,7620CT100.0% 86.7 / NA 28G404G (GGC→GGTsppAprotease IV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/15);  total (13/15)

TCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAT  >  W3110S.gb/1851705‑1851812
                                                         |                                                  
tCGCGCTGAACTGGCAGCAGCCCGGGCATCGGGTAAGCCTGTGGTTGTATCGATGggtggca                                                <  1:2132563/62‑1 (MQ=255)
tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggca                                                <  1:4151759/62‑1 (MQ=255)
tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggca                                                <  1:3420603/62‑1 (MQ=255)
tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggca                                                <  1:1299507/62‑1 (MQ=255)
tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggca                                                <  1:2341275/62‑1 (MQ=255)
tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggc                                                 <  1:890497/61‑1 (MQ=255)
       gAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggcatggcggc                                         <  1:2891996/62‑1 (MQ=255)
        aaCTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtggcatggcggc                                         <  1:1567319/61‑1 (MQ=255)
          cTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATc                                      <  1:3641974/62‑1 (MQ=255)
           tGGTAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCt                                     <  1:1731022/62‑1 (MQ=255)
           tGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCt                                     <  1:178106/62‑1 (MQ=255)
           tGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCt                                     <  1:3368666/62‑1 (MQ=255)
            ggCAGCAGCCCGGGTAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:885955/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:4364347/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:3483818/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:3805902/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:3934073/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:3988960/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:3534434/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:462668/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:4654828/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:623015/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:1092507/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:16633/1‑48 (MQ=255)
            ggCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGggtgg                                                  >  1:1354130/1‑48 (MQ=255)
                     ccGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACt                           <  1:4428920/62‑1 (MQ=255)
                             cGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTcc                   >  1:4057028/1‑62 (MQ=255)
                             cGGGTAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTcc                   >  1:238988/1‑62 (MQ=255)
                                 tAAGCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCAc                 >  1:1239254/1‑60 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTTGTTACTGGATTTCCACGCCa             >  1:2500609/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:1289369/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:3588768/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:673881/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:1825904/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:4486221/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:2820637/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCa             >  1:3896525/1‑61 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCAACGcc              >  1:3512127/1‑60 (MQ=255)
                                    gCCTGTGGTTGTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTACACGCCa             >  1:19610/1‑61 (MQ=255)
                                              gTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt  <  1:2044203/62‑1 (MQ=255)
                                              gTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt  <  1:1497569/62‑1 (MQ=255)
                                              gTATCGATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt  <  1:1745422/62‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc           <  1:2833730/48‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc           <  1:4309117/48‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc           <  1:50024/48‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc           <  1:1798173/48‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc           <  1:747999/48‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGc           <  1:3535840/48‑1 (MQ=255)
                                                   gATGGGTGGCATGGCGGCATCTGGTGGGTACTGGATTTCCACGCCAGc           <  1:2760328/48‑1 (MQ=39)
                                                         |                                                  
TCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAT  >  W3110S.gb/1851705‑1851812

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: