Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,960,118:1 +C 66.0% coding (1676/2247 nt) ptsP ← fused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,960,1181.C66.0% 129.7 / 71.3 53coding (1676/2247 nt)ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base . (11/7);  new base C (18/17);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGG  >  CP009273/2959970‑2960260
                                                                                                                                                     |                                                                                                                                              
cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAc                                                                                                                                               >  2:711983/1‑150 (MQ=255)
  gATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAccc                                                                                                                                             <  2:631124/150‑2 (MQ=255)
      cGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAcccatcg                                                                                                                                         <  1:140468/150‑6 (MQ=255)
       ggCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCgg                                                                                                                                         >  2:523869/1‑150 (MQ=255)
        gcagcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAcc                                                                                                                                              >  2:916623/1‑143 (MQ=255)
         cagcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAccc                                                                                                                                             <  1:916623/143‑2 (MQ=255)
          agcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCaa                                                                                                                                      >  2:901948/1‑149 (MQ=255)
             aTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGACGCTTGTGACCCATCGGCaa                                                                                                                                      >  2:638359/1‑147 (MQ=255)
               aaaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGaa                                                                                                                                 >  1:1048622/1‑149 (MQ=255)
                   aCCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATAtt                                                                                                                             >  1:308890/1‑150 (MQ=255)
                    ccATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTc                                                                                                                            <  1:655459/150‑1 (MQ=255)
                             gCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGccc                                                                                                                   <  2:652214/150‑1 (MQ=255)
                                   ccAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAg                                                                                                               <  2:1063048/148‑1 (MQ=255)
                                    cAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCgg                                                                                                            >  1:636218/1‑149 (MQ=255)
                                     aGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCgg                                                                                                            >  2:583876/1‑149 (MQ=255)
                                      gCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCa                                                                                                          >  1:2624/1‑150 (MQ=255)
                                        gggatgCCGATACGTGTTTTGGGAATTTCGAACCCGATCATCTCCTCGATTTCCCGTCCGCCACTTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGGGGCTTGTGCCCATCGACAACAAAATATTCGGGTTGCCCGTAGGGGCATTa                                                                                                       <  2:297096/148‑1 (MQ=255)
                                                gATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACgcag                                                                                               >  1:706776/1‑150 (MQ=255)
                                                 aTACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACgcag                                                                                               <  2:899390/150‑1 (MQ=255)
                                                 aTACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACgcag                                                                                               <  1:563230/150‑1 (MQ=255)
                                                     gTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATc                                                                                           <  2:557056/150‑1 (MQ=255)
                                                        gTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATcaaa                                                                                        >  1:999731/1‑147 (MQ=255)
                                                        gTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGcc                                                                                        >  2:959047/1‑150 (MQ=255)
                                                             ggAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGAGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCAc                                                                                   >  2:1008788/1‑150 (MQ=255)
                                                                aTTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTg                                                                                <  1:437998/150‑1 (MQ=255)
                                                                 tttCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGa                                                                              >  2:182539/1‑150 (MQ=255)
                                                                  ttCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGc                                                                                     >  2:428227/1‑143 (MQ=255)
                                                                   tCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCa                                                                                    <  1:428227/143‑1 (MQ=255)
                                                                        aCCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaag                                                                         <  1:158822/149‑1 (MQ=255)
                                                                          cGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGa                                                                              >  2:1014074/1‑142 (MQ=255)
                                                                           gATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGAt                                                                             <  1:1014074/142‑1 (MQ=255)
                                                                             tcatcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtc                                                                   <  2:308890/150‑1 (MQ=255)
                                                                               atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtct                                                                  <  1:459670/149‑1 (MQ=255)
                                                                               atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTcc                                                                <  1:901948/150‑1 (MQ=255)
                                                                                 ctcctcGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCgg                                                              >  1:760339/1‑150 (MQ=255)
                                                                                   cctcGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGc                                                             >  1:360903/1‑150 (MQ=255)
                                                                                       gACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATc                                                        <  1:366940/150‑1 (MQ=255)
                                                                                            cACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGc                                                    >  1:207449/1‑150 (MQ=255)
                                                                                             aCGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGt                                                  <  1:711983/150‑1 (MQ=255)
                                                                                                tCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTaa                                                <  1:124201/150‑1 (MQ=255)
                                                                                                      aCGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGaa                                          >  2:440904/1‑150 (MQ=255)
                                                                                                       cGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATc                                        >  1:130065/1‑150 (MQ=255)
                                                                                                        gTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATcc                                       >  1:1099018/1‑150 (MQ=255)
                                                                                                         ttCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCTGCAACAGATTATTCAGGTTGCTCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATcc                                       >  2:584070/1‑149 (MQ=255)
                                                                                                            aaTCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCAc                                    <  1:836657/150‑1 (MQ=255)
                                                                                                                 ggtggcgTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCaa                               >  1:1016917/4‑150 (MQ=255)
                                                                                                                    ggcgTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAAcc                             >  1:459278/1‑149 (MQ=255)
                                                                                                                       gTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAgg                         <  2:65042/150‑1 (MQ=255)
                                                                                                                         gCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCa                       >  2:72808/1‑150 (MQ=255)
                                                                                                                         gCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGc                                                             >  2:253499/1‑111 (MQ=255)
                                                                                                                          cTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCAt                      <  1:959047/150‑1 (MQ=255)
                                                                                                                          cTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCt                                                            <  1:253499/111‑1 (MQ=255)
                                                                                                                             catcaACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGAt                  >  2:723425/1‑150 (MQ=255)
                                                                                                                                 aaCTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATtctc              <  2:1048622/150‑1 (MQ=255)
                                                                                                                                   cTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTctct             >  2:251880/1‑150 (MQ=255)
                                                                                                                                   cTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTctct             >  1:377075/1‑150 (MQ=255)
                                                                                                                                       aTCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCt        <  1:867796/150‑1 (MQ=255)
                                                                                                                                              cTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCgg  <  2:907495/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGG  >  CP009273/2959970‑2960260

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