Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,201,960 A→C D200E (GAT→GAG iclR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,201,9600AC98.7% 653.9 / ‑3.0 154D200E (GAT→GAGiclRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (83/69);  total (83/71)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 159 total aligned reads displayed.

CTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAG  >  REL606/4201887‑4202031
                                                                         |                                                                       
 tGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGctct                                                                      <  1:3282471/76‑1 (MQ=255)
 tGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtca                                                                      <  2:2928309/76‑2 (MQ=255)
  gCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAGTAACCCCGTTTGCGCGTTTGGGCGAGCtctt                                                                     <  2:3544363/76‑1 (MQ=255)
  gCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtctt                                                                     <  2:1433344/76‑1 (MQ=255)
   cAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtcttc                                                                    >  2:3493491/1‑76 (MQ=255)
    aaGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtcttct                                                                   <  1:1005882/76‑1 (MQ=255)
    aaGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtcttct                                                                   >  1:2470269/1‑76 (MQ=255)
    aaGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtcttct                                                                   >  1:484814/1‑76 (MQ=255)
    aaGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCtctcct                                                                   >  2:4999663/1‑76 (MQ=255)
     aGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTcttctt                                                                  <  2:2048356/76‑1 (MQ=255)
      ggCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCttt                                                                 <  2:1056942/76‑1 (MQ=255)
       gCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCtttt                                                                <  2:723374/76‑1 (MQ=255)
        cAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTa                                                               <  1:685044/76‑1 (MQ=255)
         aaCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTaa                                                              <  2:4349314/76‑1 (MQ=255)
         aaCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTaa                                                              <  1:3772948/76‑1 (MQ=255)
          aCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTaaa                                                             >  1:2231895/1‑76 (MQ=255)
           cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAAt                                                            >  1:1067545/1‑76 (MQ=255)
            gTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATg                                                           >  2:624147/1‑76 (MQ=255)
             tAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGc                                                          <  1:2683573/76‑1 (MQ=255)
              aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGca                                                         >  1:1863041/1‑76 (MQ=255)
              aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGca                                                         >  1:617953/1‑76 (MQ=255)
              aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGca                                                         >  1:493170/1‑76 (MQ=255)
              aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGca                                                         >  1:1278386/1‑76 (MQ=255)
              aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGca                                                         >  2:1913227/1‑76 (MQ=255)
              aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGca                                                         >  2:419286/1‑76 (MQ=255)
                ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGcaca                                                       <  2:674283/76‑1 (MQ=255)
                ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGcaca                                                       >  2:2749079/1‑76 (MQ=255)
                 cccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAg                                                      <  1:4172669/76‑1 (MQ=255)
                  ccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAACTCTTCTTTTAAATGCACAgg                                                     <  1:4660783/76‑1 (MQ=255)
                   cAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGga                                                    >  2:1186693/1‑76 (MQ=255)
                    aGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGgag                                                   >  2:747054/1‑76 (MQ=255)
                    aGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGgag                                                   <  2:2319680/76‑1 (MQ=255)
                    aGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGgag                                                   >  1:1283386/1‑76 (MQ=255)
                      cgcATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGac                                                 <  2:4480668/76‑1 (MQ=255)
                      cgcATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGac                                                 >  2:2441711/1‑76 (MQ=255)
                       gcATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGaca                                                <  2:1731454/76‑1 (MQ=255)
                        cATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGacac                                               >  1:397175/1‑76 (MQ=255)
                        cATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGacac                                               >  1:641653/1‑76 (MQ=255)
                        cATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGacac                                               >  2:3362567/1‑76 (MQ=255)
                          tGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCa                                             >  1:746656/1‑76 (MQ=255)
                           gTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAg                                            >  1:4528958/1‑76 (MQ=255)
                           gTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAg                                            <  2:2467550/76‑1 (MQ=255)
                           gTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAg                                            >  1:1179115/1‑76 (MQ=255)
                           gTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAg                                            >  1:1467354/1‑76 (MQ=255)
                             tCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCg                                          >  1:1412189/1‑76 (MQ=255)
                             tCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCg                                          >  2:568946/1‑76 (MQ=255)
                               cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGtt                                        <  2:2751335/76‑1 (MQ=255)
                                tCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTg                                       <  2:2204200/76‑1 (MQ=255)
                                tCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTg                                       <  1:3805966/76‑1 (MQ=255)
                                   tCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGt                                    >  1:3522426/1‑76 (MQ=255)
                                    cGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTg                                   <  2:822588/76‑1 (MQ=255)
                                     gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGGGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg                                  >  1:2349215/1‑76 (MQ=255)
                                     gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg                                  >  1:4050219/1‑76 (MQ=255)
                                     gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTgg                                  >  2:2427313/1‑76 (MQ=255)
                                      tCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTggg                                 <  2:3877688/76‑1 (MQ=255)
                                      tCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTggg                                 >  2:247233/1‑76 (MQ=255)
                                       cAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGt                                <  2:1230270/76‑1 (MQ=255)
                                       cAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGt                                <  1:4166119/76‑1 (MQ=255)
                                        aaaTGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGta                               >  1:2416908/1‑76 (MQ=255)
                                        aaaTGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGta                               <  1:2299279/76‑1 (MQ=255)
                                          aTGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGtata                             <  1:464615/76‑1 (MQ=255)
                                           tGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAg                            >  2:3070142/1‑76 (MQ=255)
                                           tGAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAg                            <  1:3856869/76‑1 (MQ=255)
                                            gAATAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAgg                           >  2:1867983/1‑76 (MQ=255)
                                             aaTAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGc                          >  2:4839140/1‑76 (MQ=255)
                                             aaTAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGc                          <  2:1885017/76‑1 (MQ=255)
                                              aTAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCa                         >  1:4042782/1‑76 (MQ=255)
                                               tAACCCCGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCAt                        >  1:210439/1‑76 (MQ=255)
                                                   cccGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTaa                    >  2:4217907/1‑76 (MQ=255)
                                                   cccGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTaa                    <  2:1722703/76‑1 (MQ=255)
                                                   cccGTTTGCGCGTTTGGACGAGCTCTTCGTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTaa                    <  2:4342284/76‑1 (MQ=255)
                                                    ccGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAAc                   <  2:1528335/76‑1 (MQ=255)
                                                    ccGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAAc                   <  1:4547315/76‑1 (MQ=255)
                                                     cGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAAcc                  >  1:1668912/1‑76 (MQ=255)
                                                     cGTTTGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAAcc                  <  2:4656781/76‑1 (MQ=255)
                                                       tttGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCt                <  1:2671411/76‑1 (MQ=255)
                                                        ttGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCtt               >  2:4526916/1‑76 (MQ=255)
                                                        ttGCGCGTTTGGGCGAGCTCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCtt               <  2:1300165/76‑1 (MQ=255)
                                                         tGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCttt              <  2:300874/76‑1 (MQ=255)
                                                         tGCGCGTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCttt              >  1:2179880/1‑76 (MQ=255)
                                                           cgcgTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGc            >  1:4704828/1‑76 (MQ=255)
                                                           cgcgTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGCATAGGCATGGAACACTTTGc            >  1:3765794/1‑76 (MQ=255)
                                                            gcgTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCg           >  1:1880842/1‑76 (MQ=255)
                                                            gcgTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCg           <  2:1937819/76‑1 (MQ=255)
                                                             cgTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCgg          >  1:507295/1‑76 (MQ=255)
                                                              gTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGt         <  1:2994603/76‑1 (MQ=255)
                                                              gTTTGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGt         >  1:3195924/1‑76 (MQ=255)
                                                               tttGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTg        <  2:3946642/76‑1 (MQ=255)
                                                                ttGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgc       <  1:3413187/76‑1 (MQ=255)
                                                                ttGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgc       <  2:4410526/76‑1 (MQ=255)
                                                                ttGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgc       >  1:2487719/1‑76 (MQ=255)
                                                                ttGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgc       >  2:1314482/1‑76 (MQ=255)
                                                                ttGGGCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgc       >  1:1159122/1‑76 (MQ=255)
                                                                  gggCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGTTgcag     >  2:4982090/1‑76 (MQ=255)
                                                                  gggCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag     >  1:3602750/1‑76 (MQ=255)
                                                                  gggCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag     >  2:3616512/1‑76 (MQ=255)
                                                                  gggCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag     >  1:3386610/1‑76 (MQ=255)
                                                                  gggCGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag     >  2:4181853/1‑76 (MQ=255)
                                                                     cGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGcagcag  >  2:3424684/1‑76 (MQ=255)
                                                                     cGAGCTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGcagcag  >  1:3529525/1‑76 (MQ=255)
                                                                         |                                                                       
CTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAG  >  REL606/4201887‑4202031

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