Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 283,759 T→G A445A (GCT→GCG yagF → putative dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913283,7590TG99.2% 1571.4 / ‑6.0 586A445A (GCT→GCGyagFputative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base G (295/285);  total (295/291)
Fisher's exact test for biased strand distribution p-value = 2.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 591 total aligned reads displayed.

TGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATACCACCAC  >  NC_000913/283686‑283831
                                                                         |                                                                        
    ccggccgccGGAGAAGGCAAAAGCGAAAGGGCTGCCCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGaa                                                                    <  2:1200335/73‑1 (MQ=255)
    ccTGCCGCCGGAGAAGGCAAAAGCGGAAGGGCTGACCGCGCCGGTCTGTTTCCCGACGGGCAACATCGCGCCGGaa                                                                    >  1:4343074/1‑76 (MQ=255)
     cTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAg                                                                   <  2:3051709/76‑1 (MQ=255)
      tGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCACCATCGCGCCGGAAgg                                                                  >  1:3540910/1‑76 (MQ=255)
      tGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAgg                                                                  >  1:3522834/1‑76 (MQ=255)
        ccgccgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGtt                                                                <  1:1227875/76‑1 (MQ=255)
        ccgccgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGtt                                                                <  1:4491163/76‑1 (MQ=255)
        ccgccgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGtt                                                                >  2:4354247/1‑76 (MQ=255)
         ctccgGAGAGTGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTc                                                               <  1:3421599/74‑1 (MQ=255)
         cgccgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTc                                                               >  2:3240470/1‑76 (MQ=255)
          gccgGCGACGGCAACAGCGACAGGTGTGACCGCGCCGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCg                                                              <  2:3402576/76‑1 (MQ=255)
           ccgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCgg                                                             <  2:5188393/76‑1 (MQ=255)
            cgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGt                                                            >  2:5065879/1‑76 (MQ=255)
            cgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGt                                                            <  2:3878926/76‑1 (MQ=255)
            cgGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGt                                                            <  1:142220/76‑1 (MQ=255)
              gagaAGGCAGAAGCGTAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGa                                                          <  1:3641705/76‑1 (MQ=255)
              gagaAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCg                                                              >  2:2894906/1‑72 (MQ=255)
               agaAGGCAGTGGCGAAAGTGCTGTCGTCGACGGTCTGCTTCCGGGCGGGCAACATCGCTGCGGAAGTTTCTGTGAt                                                         <  2:3537071/76‑1 (MQ=255)
               agaAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGAt                                                         >  1:1139550/1‑76 (MQ=255)
                gaAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATc                                                        <  2:4919738/76‑1 (MQ=255)
                 aaGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCa                                                       <  1:4621003/76‑1 (MQ=255)
                  aGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCaa                                                      <  1:1760253/76‑1 (MQ=255)
                  aGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCaa                                                      <  2:2521601/76‑1 (MQ=255)
                  aGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCaa                                                      >  1:3856195/1‑76 (MQ=255)
                   ggCAAAAGCGAAAGGGCTGCCCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAg                                                     <  1:2371973/76‑1 (MQ=255)
                   ggCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAg                                                     <  2:2010622/76‑1 (MQ=255)
                     cAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGc                                                   <  2:3762883/76‑1 (MQ=255)
                      aaaaGGGAAGGGGCTGGCCGCGACGGGCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGcc                                                  <  2:5731066/76‑1 (MQ=25)
                      aaaaGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGcc                                                  <  1:3034254/76‑1 (MQ=255)
                       aaaGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCa                                                 <  2:3353836/76‑1 (MQ=255)
                       aaaGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCa                                                 >  1:565031/1‑76 (MQ=255)
                        aaGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCAc                                                >  1:1870285/1‑76 (MQ=255)
                         aGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACg                                               <  2:4249036/76‑1 (MQ=255)
                         aGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACg                                               >  2:3636685/1‑76 (MQ=255)
                         aGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACg                                               <  1:3361653/76‑1 (MQ=255)
                         aGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACg                                               <  2:1717801/76‑1 (MQ=255)
                          gCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACgg                                              <  1:523710/76‑1 (MQ=255)
                          gCGAAAGCGCCGACCTAGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGCTCAAGGCCACgg                                              >  1:5416991/1‑76 (MQ=255)
                           cGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGc                                             >  1:2251999/1‑76 (MQ=255)
                            gAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                            >  1:1067811/1‑76 (MQ=255)
                            gAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                            <  1:3944653/76‑1 (MQ=255)
                            gAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                            >  1:3735343/1‑76 (MQ=255)
                            gAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCg                                            <  1:2013364/76‑1 (MQ=255)
                             aaaGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGa                                           <  1:4046803/76‑1 (MQ=255)
                             aaaGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGa                                           <  1:5717659/76‑1 (MQ=255)
                             aaaGGGCTGACCTCGACGGTCTGCTTCCCGACGCGCAACATCGCGCCGGGAGGCTCGGTGACGAAGGCCACAGCGa                                           >  2:58180/1‑76 (MQ=255)
                              gcgggCAGACGTCGACGGTCGCGTGCCGGACGGGCAACATCGCGCCGGAAGGTTCGGGGATCAAGGCCACGGCGAt                                          <  2:1925579/74‑1 (MQ=255)
                              aaGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGAt                                          >  1:745472/1‑76 (MQ=255)
                               aGGGCTGACCTCGAGTGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATc                                         <  2:3469128/76‑1 (MQ=38)
                                 ggCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGa                                       >  1:4895370/1‑76 (MQ=255)
                                  gCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAc                                      >  1:2579073/1‑76 (MQ=255)
                                  gCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAc                                      <  1:2865259/76‑1 (MQ=255)
                                  gCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATAAAGGCCACGGCGATCggc                                      >  1:1567881/1‑74 (MQ=255)
                                    tGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGAccc                                    >  1:5828152/1‑76 (MQ=255)
                                     gACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGCTCAAGGCCACGGCGATCGAcccc                                   >  2:2451198/1‑75 (MQ=255)
                                     gACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCCAGGCCACGGCGATCGACCCg                                   >  2:5014691/1‑76 (MQ=255)
                                      aCCTCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGt                                  <  2:1832554/76‑1 (MQ=255)
                                      aCCTCGACGGGCTGCTTCCCGACGGGCAACATCGTGCCGGAAGGTTCGGTGGTCGAGGACACGGCGATCGAGCCGt                                  >  1:209971/1‑76 (MQ=255)
                                       cgtcgtcgGTGTGCTACCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTc                                 <  1:2691122/74‑1 (MQ=38)
                                         tCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCgg                               >  2:4405471/1‑76 (MQ=255)
                                         tCGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGACACGGCGCCCGCCCCGTCgg                               >  2:1284019/1‑76 (MQ=37)
                                         ccGACGGTCTGCTTCCCGACGGTCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCgg                               <  1:1297567/75‑1 (MQ=39)
                                          cGACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggt                              <  2:4604848/76‑1 (MQ=255)
                                           gACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtg                             <  1:4809372/76‑1 (MQ=255)
                                           gACGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCgggg                             >  1:427264/1‑74 (MQ=255)
                                           gACGGTCTGCTGCCCTACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtt                             <  1:4986028/76‑2 (MQ=38)
                                           gACGGTCTGCTGCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtg                             <  2:3710377/76‑1 (MQ=255)
                                             cGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtggg                           >  1:3190584/1‑75 (MQ=255)
                                             cGGTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCggtggg                           >  2:3995853/1‑75 (MQ=255)
                                              ggTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGTCCACGGCTCTCGACCCCTCGgcggtg                          >  2:1536034/1‑76 (MQ=21)
                                              ggTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGgtggtg                          >  1:3157278/1‑76 (MQ=255)
                                              ggTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGgtggtg                          <  2:3111283/76‑1 (MQ=255)
                                               gTCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGtggtgg                         <  1:5044142/76‑1 (MQ=255)
                                                tCTGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGGGGCggg                        >  2:561510/1‑76 (MQ=38)
                                                  tGCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGGGGGCg                      >  2:2586474/1‑76 (MQ=255)
                                                   gCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCgn                     <  1:2420810/76‑2 (MQ=255)
                                                   gCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGa                     <  1:4161483/76‑1 (MQ=255)
                                                   gCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGGGGGCGa                     >  1:3280791/1‑76 (MQ=255)
                                                   gCTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGGGGGGGGCgg                     >  1:1017907/1‑75 (MQ=38)
                                                    cTTCCCGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCCATCGACCCGTCGGTGGTGGGCGaa                    <  1:809317/76‑1 (MQ=38)
                                                      tCCCGACGGGCAATATCGCGCCGGGAGGCCCGGTAATCAAGGCCACGGGGAGCCACTCGTCGAGGGTGCGCGgaga                  >  2:2980539/1‑76 (MQ=255)
                                                       cccGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGag                 >  1:4226290/1‑75 (MQ=255)
                                                         cGACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGGGGTGGGCGAAGATgg               >  1:5216328/1‑76 (MQ=255)
                                                          gACGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGc              <  2:4704676/76‑1 (MQ=255)
                                                           aCGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCg             <  1:290498/76‑1 (MQ=255)
                                                           aCGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGCGGGGGGCGCAGCGGGCg             >  2:5442552/1‑76 (MQ=25)
                                                            cGGGCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGt            >  2:5304724/1‑76 (MQ=255)
                                                             gggCAACATCGCGCCGGAAGGGTCGGGGATCAAGGCCACGCCGATCGACCCGCCGGTGGGGGGCGCAGATGGCgaa           >  2:1451038/1‑74 (MQ=255)
                                                              ggCAACATCGCGCCGGACGGTTCGGTGATCAATGCCACTGCGATCGACCCTTCGGTGGTGGGCGAAGATGGCGtat          <  1:855412/76‑1 (MQ=255)
                                                              ggCAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGCCGGGGGGGGGCGCGGGTGGCgcag          >  2:77849/1‑73 (MQ=25)
                                                                cAACATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATac        <  2:2354305/76‑1 (MQ=255)
                                                                cAACATCGCGCCGGAAGGTGCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATac        <  1:4118376/76‑1 (MQ=255)
                                                                 aaCATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGGAGATGGCCGATacc       >  2:4552341/1‑76 (MQ=255)
                                                                 aaCATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATacc       >  2:2742998/1‑76 (MQ=255)
                                                                 aaCATCGCGCCGGAAGGTGCGGTGATCAAGGCCACGGCGATCGCCCCGTCGGTGGTGGGCGAAGATGGCGTATacc       <  2:1236877/76‑1 (MQ=255)
                                                                  aCATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGGAGATGGCGTATacca      >  2:1075875/1‑76 (MQ=255)
                                                                  aCATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATaccc      >  2:4809922/1‑75 (MQ=255)
                                                                  aCATCGCGCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGGGGGCGAAGATGGCGTATacca      >  1:960997/1‑76 (MQ=255)
                                                                  aCAGCGCGCCGGAAGGTTCGGGGATCAAGGCCACGGCGATCGGCCCGTCGGTGGTGGGCGAAGATGGCGTATacca      <  2:228120/76‑1 (MQ=255)
                                                                     tCGCGCCGGAAGGTTCGGTGCTCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATAccacca   >  1:4682456/1‑76 (MQ=255)
                                                                         |                                                                        
TGATCCTGCCGCCGGAGAAGGCAAAAGCGAAAGGGCTGACCTCGACGGTCTGCTTCCCGACGGGCAACATCGCTCCGGAAGGTTCGGTGATCAAGGCCACGGCGATCGACCCGTCGGTGGTGGGCGAAGATGGCGTATACCACCAC  >  NC_000913/283686‑283831

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