Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,164 G→T 9.7% intergenic (‑128/‑140) PP_5175 ← / → PP_5176 HlyD family secretion protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,903,1640GT 9.7% 144.6 / 4.9 56intergenic (‑128/‑140)PP_5175/PP_5176HlyD family secretion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/30);  new base T (2/4);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGCC  >  NC_002947/5903018‑5903308
                                                                                                                                                  |                                                                                                                                                
gCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCGCATAGCCGGCggg                                                                                                                                                >  7:307208/1‑147 (MQ=255)
gCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCATCTTCCCGGGCGac                                                                                                                                                >  1:288012/1‑148 (MQ=255)
  ggACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCACTTAGCGGGCGAGcc                                                                                                                                              >  7:168894/1‑149 (MQ=255)
     cGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcgc                                                                                                                                           <  8:18580/149‑1 (MQ=255)
      gCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcgc                                                                                                                                           >  7:296341/1‑148 (MQ=255)
      gCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCATAGCGGGCGAGCCcgc                                                                                                                                           >  8:80930/1‑148 (MQ=255)
        agcAACGCAATGCTGGATCCTTTCTTGTGCACATATGTTTTGTGGGATGCGACTGTGAGCGTTGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGaa                                                                                                                                        <  4:32662/147‑1 (MQ=255)
        aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGaa                                                                                                                                        <  5:329836/149‑1 (MQ=255)
        aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGAGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcgcg                                                                                                                                          >  5:272942/1‑147 (MQ=255)
         tCAACGCCATGCTGGCTCCTTGGTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGTTGCCTGCGGGGGCCTATTAGCGGGCGAGCCcgcgct                                                                                                                                        >  7:171557/1‑146 (MQ=255)
          cAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                      <  7:287243/149‑1 (MQ=255)
          cAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGAGGCGACTGTGAGCGTAGTTCAGGGCTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGCCCGCGAaga                                                                                                                                      <  3:223011/149‑1 (MQ=255)
           aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAaga                                                                                                                                      <  3:168715/148‑1 (MQ=255)
           aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAaga                                                                                                                                      <  2:128766/148‑1 (MQ=255)
               ccATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCt                                                                                                                                 >  6:165923/1‑149 (MQ=255)
                cATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCACTTCGCGGGCGAGCCCGCGAAGAGGCCTc                                                                                                                                >  7:7818/1‑149 (MQ=255)
                         tCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCCCCGGCGAGCCCGCGAAGAGGCCGGCACAGCCAt                                                                                                                       >  4:217430/1‑149 (MQ=255)
                           cttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                       <  3:217430/147‑1 (MQ=255)
                            gtgcttgTGCACAAATGATTTGTGGGTGTGGAAAGTGAGCGGTGTTGATGGTTGGATTAACAGCACAGAACGGGGCCAAGCGGGCAATTTTTGTTGCAGCCGGGGCCGCATAGCCCGCGCGCCCCCGAACAGGCCTCCCCAcccc                                                                                                                        >  2:52629/2‑144 (MQ=255)
                               cttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATc                                                                                                                  <  3:100321/148‑1 (MQ=255)
                               cttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCa                                                                                                                 <  6:94025/149‑1 (MQ=255)
                                  gtgCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAATCATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcg                                                                                                                                            <  6:127479/119‑1 (MQ=255)
                                  gtgCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAATCATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcg                                                                                                                                            >  5:127479/1‑119 (MQ=255)
                                      acaAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                       >  4:283124/1‑136 (MQ=255)
                                      acaAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                       <  3:283124/136‑1 (MQ=255)
                                         aaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCGCCACAGCCATGTATCAGTATCTGCCa                                                                                                       >  5:358388/1‑149 (MQ=255)
                                          aTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCa                                                                                                                 <  5:329691/138‑1 (MQ=255)
                                          aTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCAATTCGCCGGCGAGCCCGCGAAGAGGCCACCACAGCCAAGTATCa                                                                                                                 >  6:329691/1‑138 (MQ=255)
                                            gATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCata                                                                                                     <  8:7818/148‑1 (MQ=255)
                                               ttGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCGGGCATTTTTGGTGCCTGGGGGGGCCTCTTCGCGGGCTAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATcc                                                                                                  <  4:168586/148‑1 (MQ=255)
                                                  tgGGATGCTACTGTTACCGTAGTTCTGGGTTCTATAAACAACACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGTGGCCTCTTCGGGGGCGAGCCCGCTATGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGa                                                                                               <  5:245190/148‑1 (MQ=255)
                                                    ggATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCGGGCGCGCCCGCGAAGAGGCCTCCACAGCCATGTATAAGTATATTCCATATACAGAAAc                                                                                            >  3:277365/1‑149 (MQ=255)
                                                     gATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGCAGAGGCCTCCACAGCCATGTATCAGTATCTGGCATATGCAGAAAc                                                                                            >  4:105216/1‑148 (MQ=255)
                                                         cgtgtgttAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCTTTTTTGGTGGCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCa                                                                                       <  7:80930/145‑1 (MQ=255)
                                                         cGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGCGGCCTCCCCAGCCATGTATCCGTATCTGGCATATCCAGCCACCAt                                                                                         >  8:255128/1‑147 (MQ=255)
                                                         cGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCa                                                                                       <  7:296411/149‑1 (MQ=255)
                                                                  gCGTAGTTCAGGGTTCGATAAAGAGCACGGACGGGGGCCTTGCTGGCATTTTTGGGGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCgg                                                                               <  8:168894/148‑1 (MQ=255)
                                                                        ttCCGGGTTCGATAAACATCACTGACGGGGGCCATGCTGCCATTTTTGGGGCCTGGGGGGGCCGCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATcc                                                                         <  5:34385/148‑1 (MQ=255)
                                                                           tgggTCCCATAACCATCACTGCCGGGGGCCATGCGGGAATTTTGGGTGCCTGCGGGGGCCTCTTCGGGGGCGACCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCg                                                                      <  2:163187/147‑1 (MQ=255)
                                                                            gggTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGtt                                                                    <  7:88979/149‑1 (MQ=255)
                                                                                   aTAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGCGGGCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCCACCg                                                             >  5:293410/1‑149 (MQ=255)
                                                                                   aTAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCg                                                             >  6:391958/1‑149 (MQ=255)
                                                                                                   ggggCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGa                                             <  3:105216/149‑1 (MQ=255)
                                                                                                     ggCCATGCTGGAATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAgg                                           <  4:277365/149‑1 (MQ=255)
                                                                                                         aTGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTca                                       <  8:296341/149‑1 (MQ=255)
                                                                                                             tGGCATTTTTGGTGCCTGCGGGGGCCTATTAGCGGGCGAGGCCGCGAAGCGGCCTCCACAGCCATGTATCAGTATATGCCATTTACAGAAACCATCAGGCATCGGCCCTCCCCGGTCCAACCGCTACCTGTAACAGGGAGGCTcaccg                                    >  8:282198/1‑146 (MQ=255)
                                                                                                                cATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCTACCGCTACCTTTCACAGGGAGGCTCACATTCaa                                <  1:52629/149‑1 (MQ=255)
                                                                                                                 aTTTTTGGGGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAg                               <  1:378441/149‑1 (MQ=255)
                                                                                                                   ttttGTGGCCTGCGGGGGCCTCTGGGCGGTCTAGCCGGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAgg                              <  5:27367/148‑1 (MQ=255)
                                                                                                                    tttGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGa                             <  5:391958/148‑1 (MQ=255)
                                                                                                                     ttGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGaa                            <  7:57526/148‑1 (MQ=255)
                                                                                                                               cGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGc                 <  6:358388/149‑1 (MQ=255)
                                                                                                                                    gCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTa            <  3:220196/149‑1 (MQ=255)
                                                                                                                                       tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc         >  2:145812/1‑149 (MQ=255)
                                                                                                                                        ctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc         >  5:337479/1‑148 (MQ=255)
                                                                                                                                           cgcgGGCGAGCCCGCGCAGATGCCGCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATGAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGTTCAAAATGCCGTTAGCCGaa      <  2:288012/148‑1 (MQ=255)
                                                                                                                                              gggCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGcc  <  7:282198/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATGCC  >  NC_002947/5903018‑5903308

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