Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 761,422 C→T A906V (GCC→GTC)  sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913761,4220CT98.2% 156.2 / NA 56A906V (GCC→GTC) sucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCAGCAGTTTGCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGCCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCCTGACCTGCCT  >  NC_000913/761283‑761557
                                                                                                                                           |                                                                                                                                       
ttGCAGCAGTTTGCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGtctc                                                                                                                                      >  2:807205/1‑143 (MQ=255)
ttGCAGCAGTTTGCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGtctc                                                                                                                                      >  2:265719/1‑143 (MQ=255)
   cagcagTTTGCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTgcgc                                                                                                                                   >  2:393833/1‑143 (MQ=255)
   cagcagTGTGCTCCGGTCAAGGATTTTGTCTGGTGCCAGGCAGGGCCGCTCAGCCAGGGCGCATGGTGCTGCAGCCGGCATCATGTCCGGGAGGGGATTCCGTTTGGGGCTTCTCTGCGTTTTGCAGGCCGCCCGGTCTCCGc                                                                                                                                   <  1:835509/143‑1 (MQ=255)
    agcagTTTGCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCACATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCGGcc                                                                                                                                  >  1:480301/1‑143 (MQ=255)
            gCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCcggc                                                                                                                          <  2:525629/143‑1 (MQ=255)
                 cGTCAAGGATGTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAg                                                                                                                     <  1:361959/143‑1 (MQ=255)
                  gTCAAGGATTTTGTCTGGTCCCCGGCCGCGCCCCCCAACCAGGTGGCGTGGGGCTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAgg                                                                                                                    <  2:353677/143‑1 (MQ=255)
                      aGGATTGTGTCTGGGGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGtat                                                                                                                <  2:613847/143‑1 (MQ=255)
                      aGGAGTTTGTCTGGTGGCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGtat                                                                                                                <  2:156767/143‑1 (MQ=255)
                       ggATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCGGCGACTCTCCAGCGCTAGGGtata                                                                                                               >  2:604951/1‑143 (MQ=255)
                            ttggccgggaccAGGAGGAGCCGCTCGACCAGGGCGCATGTTGCGGCAGCCAGCATCATCCCCGTGAAGTGAGTCTGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATAGGTcc                                                                                                          <  1:915573/133‑1 (MQ=255)
                              gTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAGGAGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGt                                                                                                        <  2:487410/143‑1 (MQ=255)
                                  ggTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGCAGGGTATATGTCCGTTCAc                                                                                                    >  2:245736/1‑143 (MQ=255)
                                   gTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCAcc                                                                                                   >  1:578949/1‑143 (MQ=255)
                                   gTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCAGGCGGTAGGGTATATGGCCGTTCAcc                                                                                                   >  2:253659/1‑143 (MQ=255)
                                    tGCCATGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCa                                                                                                  <  2:937061/143‑1 (MQ=255)
                                    tGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCa                                                                                                  <  2:262012/143‑1 (MQ=255)
                                    tGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGACAGCATCATTTACGTGAAGTGATTCCTTTTGGGGCTTCTCTGCGTTATGCAGCCCGCCCGGTCTCCGCCTATCCGGCGGTAGGGTAGATGACCGATCAccg                                                                                                  >  1:112420/1‑142 (MQ=255)
                                     gCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAg                                                                                                 >  1:485429/1‑143 (MQ=255)
                                      ccAGGAAGGGCCGCTCAACCAGGGCGCAGGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGa                                                                                                <  2:187466/143‑1 (MQ=255)
                                        aGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGaaa                                                                                              <  2:511013/143‑1 (MQ=255)
                                        aGGAAGAGCCGCTCAACCAGGGCGCATGGTACTCCAGCCAGCATCATTTCCGTGAAGTGATTCGGTATGGGGCGTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCCCTCAGACGGGCGGGTATATGTCCGGTGAccgggaa                                                                                              >  2:358779/1‑138 (MQ=255)
                                          gAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACa                                                                                            <  2:123443/143‑1 (MQ=255)
                                                        cctgggCGCATGGTCCTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTgg                                                                              <  2:8186/140‑1 (MQ=255)
                                                        ccggggCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTgg                                                                              <  1:648806/140‑1 (MQ=255)
                                                             gcgcCTGGTACTGCAGCCAGCATCATTCCCGTGAAGTGCTTCAGTTTGGCCCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAAt                                                                         <  1:457311/143‑1 (MQ=255)
                                                             gcgcATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTCTACGTCCGTTCACCAGCAACAGCAACCAGCTCCGGTTCAt                                                                         >  1:434649/1‑143 (MQ=255)
                                                              cggctGGTACTGCAGCCAGCATCATTTCCGTGGAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATg                                                                        <  1:92185/139‑1 (MQ=255)
                                                               gcATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATCTGTCCGTTCGCCAGAAACAGCACCACGaaacagataatag                                                                       >  2:256210/1‑131 (MQ=255)
                                                                        tGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCGGGTTAATGACGCGCTgga                                                              >  2:430747/1‑141 (MQ=255)
                                                                        ggCAGCCAGCATCTTTGCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGaa                                                              <  2:240507/142‑1 (MQ=255)
                                                                          cagccagcATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATCTGTCCGTTCCCCCGACACAGCAACAGGATCTGGTTAATGACGCGCTGCACg                                                            >  1:489936/1‑143 (MQ=255)
                                                                           agccagcATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCACCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGt                                                           <  1:720172/143‑1 (MQ=255)
                                                                            gccagcATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTc                                                          <  2:483737/143‑1 (MQ=255)
                                                                                 catcatTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCCGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGaata                                                     >  1:459933/1‑143 (MQ=255)
                                                                                 catcatTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGaata                                                     >  1:454707/1‑143 (MQ=255)
                                                                                 catcatTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGaata                                                     >  1:24541/1‑143 (MQ=255)
                                                                                    caattCCGTGAAGGGAGTCCGTCTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGaataaat                                                  <  1:953295/140‑1 (MQ=255)
                                                                                     ggttCGGTGGAGTGATTGCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGTGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGaataaata                                                 <  1:501058/141‑1 (MQ=255)
                                                                                      tttCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGGTAATGACGCGCTGCACGTCGAacaaataa                                                >  1:88166/1‑143 (MQ=255)
                                                                                         ccGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAgg                                             <  2:651760/143‑1 (MQ=255)
                                                                                              aaGTGGTTCCGTGTGGGGCTTCTGTGCGTTATGCGGGCCGCCCGGTCTCCGCCGCTCCGGCGGTAGGGTATATGTGCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATaca                                        <  1:951540/143‑1 (MQ=255)
                                                                                              aaGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGACCGGTCTCCGCCTCGCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATaca                                        >  2:491877/1‑143 (MQ=255)
                                                                                               agtcattttgTTTGGGGCTTCTCTGCGTTCTGCTGGGGCCCCGGTCTTCGCCTTTCCTGTGGTAGGGTATTTGTCCGTTCATCAGAAAGAGCAACGAGATCTGGTTAATGTCGCGCTTATCGTCGTATAAATATAGGATacac                                       <  1:614223/134‑1 (MQ=255)
                                                                                                    ttCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGa                                  >  1:552307/1‑143 (MQ=255)
                                                                                                         tttGGGGCTTCTCTGCGTTATGAAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACAAGAAACAGCAACAAGCTCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACAAAATCAGTAGc                             >  1:425557/1‑143 (MQ=255)
                                                                                                               gCTTCTCTGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAGCAGCACCAAGATCTGGTTAATGACGCGCTGAACGTCGAACAAATAAAGGATACACAATGCGTCGCGTAGat                       >  1:25340/1‑143 (MQ=255)
                                                                                                                    tctGCGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAAGATACACAATGAGAAGCGTAGATATACt                  >  1:476868/1‑143 (MQ=255)
                                                                                                                     cgtgtatgtggaGGACGCCCGGTCGCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTg                 <  1:505988/132‑1 (MQ=255)
                                                                                                                        cGTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTAGATGTCCGGTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAACGGATACACCATGAGCAGCGAAGATATTCAggac              >  1:953812/1‑141 (MQ=255)
                                                                                                                         gTTATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTcc             <  1:265719/143‑1 (MQ=255)
                                                                                                                          ttgtGCAGGCCGCCCGCTCTCCGCCTCGCCTGCGGTAGGGTATCTGTTCGTTCACCAGAAACAGCAACAAGTTCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTccc            <  1:124395/140‑1 (MQ=255)
                                                                                                                           tATGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCCt           <  2:24541/143‑1 (MQ=255)
                                                                                                                            aTGCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCCTg          <  2:449089/143‑1 (MQ=255)
                                                                                                                              gCAGGCCGCCCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCCTGAc        <  2:643480/143‑1 (MQ=255)
                                                                                                                                  gccgccCGGTCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCATGCcctgc    >  1:440236/1‑143 (MQ=255)
                                                                                                                                    cgccCGGCCTCCGCCGCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCCTGAcctgcct  <  1:450135/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TTGCAGCAGTTTGCTCACGTCAAGGATTTTGTCTGGTGCCAGGAAGAGCCGCTCAACCAGGGCGCATGGTACTGCAGCCAGCATCATTTCCGTGAAGTGATTCCGTTTGGGGCTTCTCTGCGTTATGCAGGCCGCCCGGCCTCCGCCTCTCCGGCGGTAGGGTATATGTCCGTTCACCAGAAACAGCAACAAGATCTGGTTAATGACGCGCTGAACGTCGAATAAATAAAGGATACACAATGAGTAGCGTAGATATTCTGGTCCCTGACCTGCCT  >  NC_000913/761283‑761557

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