Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,909,167 C→G 11.3% A42A (GCG→GCC PP_1709 ← fumarylacetoacetate hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,909,1670CG11.3% 151.5 / 11.3 53A42A (GCG→GCCPP_1709fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base G (4/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

GATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAG  >  NC_002947/1909023‑1909314
                                                                                                                                                 |                                                                                                                                                    
gctGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGgcgccgcgcc                                                                                                                                                 >  2:1200653/3‑150 (MQ=255)
gATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGgcgccgcgcc                                                                                                                                                 >  2:455395/1‑150 (MQ=255)
      ggCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                           <  1:859182/150‑1 (MQ=255)
           cgccTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCg                                                                                                                                      <  1:800982/150‑1 (MQ=255)
             ccTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGc                                                                                                                                    >  1:655141/1‑150 (MQ=255)
               tttGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                           >  1:48792/1‑141 (MQ=255)
               tttGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                           <  2:48792/141‑1 (MQ=255)
                   gCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACggg                                                                                                                              <  2:57364/150‑1 (MQ=255)
                       tGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGt                                                                                                                          >  1:777892/1‑150 (MQ=255)
                           cgcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCg                                                                                                                      >  2:499183/1‑150 (MQ=255)
                           cgcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCg                                                                                                                      <  2:655141/150‑1 (MQ=255)
                            gcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCgg                                                                                                                     >  1:849552/1‑150 (MQ=255)
                                  aGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAgg                                                                                                               <  1:676469/150‑1 (MQ=255)
                                    cTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                             <  1:1098806/150‑1 (MQ=255)
                                     tCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCCCGTTGCTTGCCGGTTTCATGAAAAACCCAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGccg                                                                                                            >  2:202294/1‑148 (MQ=255)
                                                gCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGcc                                                                                                >  1:1364193/1‑150 (MQ=255)
                                                      cTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCgg                                                                                          <  2:945922/150‑1 (MQ=255)
                                                       tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAg                                                                                                                <  1:872122/128‑1 (MQ=255)
                                                       tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCggg                                                                                         >  1:1038015/1‑150 (MQ=255)
                                                       tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTACGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAg                                                                                                                >  2:872122/1‑128 (MQ=255)
                                                          tcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAg                                                                                      <  2:494965/150‑1 (MQ=255)
                                                           cgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGTATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAgc                                                                                     <  1:175642/150‑1 (MQ=255)
                                                                            ggCAGTTCGCCCAGTGCCTCGACCACTTTGCTTGCCGGTTTCATGAAAAACACCAAGCCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACGCCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCgg                                                                    >  1:886386/1‑150 (MQ=255)
                                                                             gCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGa                                                                   <  2:777892/150‑1 (MQ=255)
                                                                             gCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGCTTAATGAAACCAACGCGCTCCTGCGGCGCCGCGCCGTAGGATTACGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGa                                                                   >  1:1096857/1‑150 (MQ=255)
                                                                                    gCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAg                                                            <  2:29989/150‑1 (MQ=255)
                                                                                      ccAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                             <  1:1094457/100‑1 (MQ=255)
                                                                                      ccAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                             >  2:1094457/1‑100 (MQ=255)
                                                                                             cTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGAt                                                   >  2:1271716/1‑150 (MQ=255)
                                                                                             cTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGAt                                                   >  2:347640/1‑150 (MQ=255)
                                                                                              tCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                  <  2:849552/150‑1 (MQ=255)
                                                                                              tCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                  <  2:1096857/150‑1 (MQ=255)
                                                                                                     cGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCcggcg                                           <  1:1164827/150‑1 (MQ=255)
                                                                                                       ttgcttgcCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCggg                                         <  1:1271716/150‑1 (MQ=255)
                                                                                                       ttgcttgcCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCggg                                         >  1:1215953/1‑150 (MQ=255)
                                                                                                           ttgcCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGGGGCCCATCCGGGATGCCCTCCTGGGCAGCTCAGCGGTGCCGCGGATCGGCAGGCTGACGATAGCCGGGGGGGTGa                                     >  1:525498/1‑150 (MQ=255)
                                                                                                           ttgcCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCg                                                                <  2:1018280/123‑1 (MQ=255)
                                                                                                           ttgcCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCg                                                                >  1:1018280/1‑123 (MQ=255)
                                                                                                            tgcCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGaa                                    <  1:1200653/150‑1 (MQ=255)
                                                                                                               cGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACAcc                                                                                                     <  1:167251/82‑1 (MQ=255)
                                                                                                               cGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACAcc                                                                                                     >  2:167251/1‑82 (MQ=255)
                                                                                                                   ttCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTaa                             <  2:188869/150‑1 (MQ=255)
                                                                                                                       tGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                  <  2:1413540/125‑1 (MQ=255)
                                                                                                                       tGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                  >  1:1413540/1‑125 (MQ=255)
                                                                                                                       tGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGtga                         >  1:1172274/1‑148 (MQ=255)
                                                                                                                       tGAAAAACACAGGCTCCTGCGGCGCCACGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCa                         >  1:831764/1‑150 (MQ=255)
                                                                                                                         aaaaaCACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATg                       >  2:232724/1‑150 (MQ=255)
                                                                                                                             acacAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTg                   <  1:499183/150‑1 (MQ=255)
                                                                                                                               gccggaatCCTACGGCGCGGCGCCGCAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTgg                  >  2:1184896/8‑150 (MQ=255)
                                                                                                                                    gCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGAtcct             >  1:297358/1‑150 (MQ=255)
                                                                                                                                    aatCCTACGGCGCGGCGCCGCAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGAtcct             >  1:1293549/3‑150 (MQ=255)
                                                                                                                                       ccTGCGGCGCCGCGCCGTAGGATTCCGGCACCGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATcctcct          <  2:337599/150‑1 (MQ=255)
                                                                                                                                       ccTACGGCGCGGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATcctcct          >  1:360228/1‑150 (MQ=255)
                                                                                                                                         tGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACa                                  >  2:278747/1‑124 (MQ=255)
                                                                                                                                         tGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACa                                  <  1:278747/124‑1 (MQ=255)
                                                                                                                                               gccgcgccgTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAg  <  2:297358/150‑1 (MQ=255)
                                                                                                                                               gccgcgccgTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCCCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAg  <  1:455395/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATGTTGGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACACAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGAACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAG  >  NC_002947/1909023‑1909314

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