Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,462,961 T→C V13A (GTA→GCA)  gspK → Type II secretion system protein GspK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,462,9610TC100.0% 188.0 / NA 54V13A (GTA→GCA) gspKType II secretion system protein GspK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/21);  total (33/21)

CAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCG  >  NC_000913/3462824‑3463097
                                                                                                                                         |                                                                                                                                        
cAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCAc                                                                                                                                        >  1:382784/1‑140 (MQ=255)
       aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGGGGCGTCGCGCTGCTCATTGCACTCATGtt                                                                                                                                 >  1:228011/1‑140 (MQ=255)
       aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGtt                                                                                                                                 >  1:142371/1‑140 (MQ=255)
       aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGtt                                                                                                                                 >  1:416333/1‑140 (MQ=255)
       aGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGtt                                                                                                                                 >  1:378888/1‑140 (MQ=255)
        gCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTa                                                                                                                                <  1:112984/140‑1 (MQ=255)
             atGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGgc                                                                                                                           >  1:144490/1‑140 (MQ=255)
               ggggggTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGTGAATAATAACCATGAATAATGAACAACGAGGCGTCGGGCTGCTGATTGCAGTCATGTTACTGgcgc                                                                                                                         >  1:188627/1‑140 (MQ=255)
                    gTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATg                                                                                                                    <  2:142371/140‑1 (MQ=255)
                    gTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATg                                                                                                                    <  2:539227/140‑1 (MQ=255)
                      tACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGgc                                                                                                                  >  2:54050/1‑140 (MQ=255)
                      tACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGgc                                                                                                                  >  2:302436/1‑140 (MQ=255)
                       aCAGCGCCGTTTTGCATTACCTGAGCAGTTATCGACAGAGGAGTCGCCGGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGgct                                                                                                                 <  2:222685/140‑1 (MQ=255)
                              cGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTgg                                                                                                          >  2:108395/1‑140 (MQ=255)
                               gTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGgc                                                                                                         >  1:648037/1‑140 (MQ=255)
                                             gAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGc                                                                                           <  2:648037/140‑1 (MQ=255)
                                             gAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGc                                                                                           >  2:355200/1‑140 (MQ=255)
                                             gAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGc                                                                                           <  2:435385/140‑1 (MQ=255)
                                                     aGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                                   >  1:448025/1‑140 (MQ=255)
                                                     aGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                                   <  1:20188/140‑1 (MQ=255)
                                                     aGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                                   >  1:647145/1‑140 (MQ=255)
                                                     aGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGtttt                                                                                   >  1:571750/1‑140 (MQ=255)
                                                               gagTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCAcagccagc                                                                         >  2:565712/1‑140 (MQ=255)
                                                                    gCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAg                                                                    >  1:474100/1‑140 (MQ=255)
                                                                     ccAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCGGCAGACAGGACGCTCAGTTTTCACAGCCAGCTTCAGc                                                                   >  1:171212/1‑140 (MQ=255)
                                                                     ccAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGc                                                                   >  2:291292/1‑140 (MQ=255)
                                                                     ccAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGc                                                                   >  1:10744/1‑140 (MQ=255)
                                                                            aGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCg                                                            <  1:553451/140‑1 (MQ=255)
                                                                            aGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCg                                                            <  2:588769/140‑1 (MQ=255)
                                                                             gACGCAGGCCGGGAATAATAACGATGAATAATGAACATCGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGc                                                           <  1:298108/140‑1 (MQ=255)
                                                                             gACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGc                                                           <  1:565712/140‑1 (MQ=255)
                                                                                  aGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGt                                                      >  1:41534/1‑140 (MQ=255)
                                                                                        ggAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTGTTCACAGTCAGATTGAGCGAATCCGCCAGGTAAACca                                                >  1:317980/1‑140 (MQ=255)
                                                                                              ataaCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCt                                          >  1:523330/1‑140 (MQ=255)
                                                                                              ataaCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCt                                          >  2:372114/1‑140 (MQ=255)
                                                                                                aaCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGc                                        <  2:552082/140‑1 (MQ=255)
                                                                                                 aCCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCa                                       <  2:486390/140‑1 (MQ=255)
                                                                                                   cATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAAc                                     <  2:10744/140‑1 (MQ=255)
                                                                                                    aTGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTTCTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACg                                    >  2:430195/1‑140 (MQ=255)
                                                                                                    aTGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACg                                    >  1:167654/1‑140 (MQ=255)
                                                                                                     tGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACgg                                   <  2:378888/140‑1 (MQ=255)
                                                                                                       aataatGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCa                                 <  1:405781/140‑1 (MQ=255)
                                                                                                       aataatGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCa                                 >  2:464219/1‑140 (MQ=255)
                                                                                                          aatGAACAACGCGGCGTCGCGCTGCTCATTGCACTGATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAata                              <  2:382784/140‑1 (MQ=255)
                                                                                                          aatGAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAata                              <  2:144490/140‑1 (MQ=255)
                                                                                                             gAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAata                              >  2:100994/1‑137 (MQ=255)
                                                                                                             gAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAata                              <  1:100994/137‑1 (MQ=255)
                                                                                                             gAACAACGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGa                           <  1:54050/140‑1 (MQ=255)
                                                                                                                caacGCGGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCAGCCAGGTAAACGATAATCTGCAACGGCAATATGatat                        >  1:188269/1‑140 (MQ=255)
                                                                                                                   cgcgGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGa                     >  2:374132/1‑140 (MQ=255)
                                                                                                                   cgcgGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGa                     >  2:42553/1‑140 (MQ=255)
                                                                                                                   cgcgGCGTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGa                     >  2:662321/1‑140 (MQ=255)
                                                                                                                        cgTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTgg                >  1:341283/1‑140 (MQ=255)
                                                                                                                        cgTCGCGCTGCTCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTgg                >  1:46566/1‑140 (MQ=255)
                                                                                                                                 gctCATTGCACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAAc       <  2:341283/140‑1 (MQ=255)
                                                                                                                                  ctCATCGCACTCATGTTACTCGCGCTCACGGCTGCTCTGGCAATAGACATGACCCTCAGTTTTCACAGCCTGCTTCAGCCAACCCGCCAGGTAAACCATCATCTGCAACTGCAATATGATATTGAGCCGGCGgtgagaat      >  2:73375/1‑136 (MQ=255)
                                                                                                                                      ttGCACTCATGTTACTGGCGCTCATGGCTGCGCTGCCAGGAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTgcg  <  2:228011/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGTTACTGGCGCTCATGGCTGCTCTGGCAGCAGACATGACGCTCAGTTTTCACAGCCAGCTTCAGCGAACCCGCCAGGTAAACCATCATCTGCAACGGCAATATGATATTGAACTGGCGGAGAAACTTGCG  >  NC_000913/3462824‑3463097

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