Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 960,601 T→C V113A (GTT→GCT)  ycaL → periplasmic protease YcaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913960,6010TC100.0% 137.3 / NA 39V113A (GTT→GCT) ycaLperiplasmic protease YcaL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/18);  total (21/18)

AGTAAATCCAGCAAGTACGGTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGTTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGCTTCTTATGCGATCGTTGCCGCACG  >  NC_000913/960465‑960736
                                                                                                                                        |                                                                                                                                       
aGTAAATCCAGCAAGTACGGTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCg                                                                                                                                      >  2:927395/1‑140 (MQ=255)
           cAAGTACGGTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTg                                                                                                                           <  1:23803/140‑1 (MQ=255)
              gTACGGTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGcc                                                                                                                        <  2:535554/140‑1 (MQ=255)
              gTACGGTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGCAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGcc                                                                                                                        >  1:348178/1‑140 (MQ=255)
                   gTAAACGTCTGGCAAAACTCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATTGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATg                                                                                                                   <  2:856946/140‑1 (MQ=255)
                   gTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATg                                                                                                                   <  1:266961/140‑1 (MQ=255)
                    tAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATgg                                                                                                                  >  1:160970/1‑140 (MQ=255)
                       aCGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTGCAGTGGCCTGATGGACa                                                                                                               <  1:651710/140‑1 (MQ=255)
                            tGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACatgatg                                                                                                          >  1:271656/1‑140 (MQ=255)
                                 aaaaTCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGatgaatga                                                                                                     <  2:251509/140‑1 (MQ=255)
                                   aaTCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACa                                                                                                   >  1:272485/1‑140 (MQ=255)
                                   aaTCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACa                                                                                                   >  2:40119/1‑140 (MQ=255)
                                                ttGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAgg                                                                                      <  2:664449/140‑1 (MQ=255)
                                                   ggTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGt                                                                                   >  1:844999/1‑140 (MQ=255)
                                                             ttAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCAt                                                                         >  2:690788/1‑140 (MQ=255)
                                                               aaTGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTGCAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGa                                                                       >  1:23870/1‑140 (MQ=255)
                                                                aTGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGaa                                                                      >  1:858730/1‑140 (MQ=255)
                                                                   gCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTg                                                                   >  1:951405/1‑140 (MQ=255)
                                                                     aCGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGcc                                                                           >  1:40331/1‑130 (MQ=255)
                                                                     aCGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGcc                                                                           <  2:40331/130‑1 (MQ=255)
                                                                     aCGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGaa                                                                      <  1:126929/135‑1 (MQ=255)
                                                                     aCGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGaa                                                                      >  2:126929/1‑135 (MQ=255)
                                                                      cGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGc                                                                <  1:73820/140‑1 (MQ=255)
                                                                       gCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGaa                                                                      >  1:618659/1‑133 (MQ=255)
                                                                       gCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGaa                                                                      <  2:618659/133‑1 (MQ=255)
                                                                          ggTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACg                                                            <  2:348178/140‑1 (MQ=255)
                                                                           gTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGt                                                           <  2:845844/140‑1 (MQ=255)
                                                                               aCTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTcgcg                                                       <  1:673241/140‑1 (MQ=255)
                                                                                           atatGACCAGCGACGTCAACGCATGGGCGATGGCGAAAGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCg                                           >  2:849180/1‑140 (MQ=255)
                                                                                                 ccAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCgctggct                                     >  2:105874/1‑140 (MQ=255)
                                                                                                  cAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCgctggctg                                    <  1:105874/140‑1 (MQ=255)
                                                                                                    gCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCtgaa                                  >  1:325900/1‑140 (MQ=255)
                                                                                                     cGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCtgaat                                 >  2:32905/1‑139 (MQ=255)
                                                                                                     cGACGTCAACGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCtgaaa                                 <  1:927395/140‑1 (MQ=255)
                                                                                                            aaCGCATGGGCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATCGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGc                          <  2:844999/140‑1 (MQ=255)
                                                                                                                   gggCGATGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGCTTCTTAt                   >  1:152971/1‑140 (MQ=255)
                                                                                                                         tGGCGAACGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGCTTCTTATGCGATc             >  1:333875/1‑140 (MQ=255)
                                                                                                                              aaCGGCTGTGCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGCTTCTTATGCGATCGTTGc        >  1:817444/1‑140 (MQ=255)
                                                                                                                                    tgtgCTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGCTTCTTATGCGATCGTTGCCGCAcg  <  2:749033/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AGTAAATCCAGCAAGTACGGTAAACGTCTGGCAAAAATCGCCAAAGCATTGGGTAACAATATTAATGGCACGCCGGTCAACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGTTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACTCGCTGGCTGAAATGAAAGCTTCTTATGCGATCGTTGCCGCACG  >  NC_000913/960465‑960736

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