Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,804,385 G→C 8.5% V66V (GTG→GTC rbsR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,804,3850GC 8.5% 171.6 / 3.1 55V66V (GTG→GTCrbsRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base C (2/3);  total (26/29)
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

CCATGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTGGAAAAGCGCA  >  NC_002947/2804244‑2804533
                                                                                                                                             |                                                                                                                                                    
ccATGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCTGTGACGCGCTCGCTGAAAGCGCTCAGCACGGCCACCATCGGCCTGTTGGTGCCCAAcc                                                                                                                                               >  5:80817/1‑148 (MQ=255)
ccATGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCaa                                                                                                                                                 >  5:164267/1‑147 (MQ=255)
ccATGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACa                                                                                                                                               <  4:329043/149‑1 (MQ=255)
  aTGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGGCACCAACAgg                                                                                                                                             <  6:407226/149‑2 (MQ=255)
   tGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGGCACCAACAggc                                                                                                                                            <  3:428368/149‑3 (MQ=255)
          aacaaGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAAccc                                                                                                                                     >  5:20213/1‑149 (MQ=255)
          aacaaGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAAccc                                                                                                                                     >  3:328328/1‑149 (MQ=255)
           acaaGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCg                                                                                                                                    >  4:29191/1‑149 (MQ=255)
             aaGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTa                                                                                                                                  <  7:141215/149‑1 (MQ=255)
              aGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTAc                                                                                                                                 <  3:265379/149‑1 (MQ=255)
                              cGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGgcgc                                                                                                                 >  7:243592/1‑149 (MQ=255)
                               gAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGcgcg                                                                                                                >  3:428177/1‑149 (MQ=255)
                                aGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCgg                                                                                                               <  4:223829/149‑1 (MQ=255)
                                 gCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCggg                                                                                                              <  3:29191/149‑1 (MQ=255)
                                 gCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCggg                                                                                                              <  4:37902/149‑1 (MQ=255)
                                 gCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCggg                                                                                                              <  4:328328/149‑1 (MQ=255)
                                        cGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAAGAGCGTCAACCCGTTCTTCGCCGAGCTGGCGCTGGGTATCg                                                                                                        <  4:102731/148‑1 (MQ=255)
                                        cGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGa                                                                                                       <  4:428177/149‑1 (MQ=255)
                                          gCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCg                                                                                                        >  4:157710/1‑146 (MQ=255)
                                           cTGAAAGTGGAAGCCGCCATCATCGAACTGGAGTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTCCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGa                                                                                                    >  4:361206/1‑149 (MQ=255)
                                              aaaGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGaa                                                                                                      >  7:380990/1‑144 (MQ=255)
                                              aaaGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGaa                                                                                                      <  8:380990/144‑1 (MQ=255)
                                              aaaGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGc                                                                                                 >  8:251502/1‑149 (MQ=255)
                                                    gAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGa                                                                                           <  4:276828/149‑1 (MQ=255)
                                                           ccATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAAc                                                                                    <  3:157710/149‑1 (MQ=255)
                                                           ccATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAAc                                                                                    <  5:145045/149‑1 (MQ=255)
                                                             atcatcGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACgg                                                                                  >  2:306282/1‑149 (MQ=255)
                                                                 tcGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTAt                                                                              >  1:361203/1‑149 (MQ=255)
                                                                  cGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTAtt                                                                             >  6:239818/1‑149 (MQ=255)
                                                                     aCTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCg                                                                          <  5:431507/149‑1 (MQ=255)
                                                                        ggATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCgg                                                                                                               >  2:168098/1‑109 (MQ=255)
                                                                        ggATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCgg                                                                                                               <  1:168098/109‑1 (MQ=255)
                                                                         gactACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTCCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTg                                                                        <  3:361206/144‑1 (MQ=255)
                                                                         gATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGa                                                                       <  3:162771/148‑1 (MQ=255)
                                                                               gTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCtgtg                                                                >  4:343963/1‑149 (MQ=255)
                                                                                  cccAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCaa                                                             <  7:251502/149‑1 (MQ=255)
                                                                                   ccAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAAc                                                            >  4:370069/1‑149 (MQ=255)
                                                                                    cAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACt                                                           >  8:148991/1‑149 (MQ=255)
                                                                                    cAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACt                                                           >  3:59740/1‑149 (MQ=255)
                                                                                     aGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTc                                                          >  7:402275/1‑149 (MQ=255)
                                                                                         cGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCcgac                                                      <  6:164267/149‑1 (MQ=255)
                                                                                            tGGCGCGCTCGCTGAAGGCGCGCAGCACGGCCACCATCCGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCcgacga                                                    >  1:416277/1‑148 (MQ=255)
                                                                                                    tcgctGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCa                                           >  8:177903/1‑149 (MQ=255)
                                                                                                      gctGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAg                                          <  5:239818/148‑1 (MQ=255)
                                                                                                             gcgcgcAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAagcagtg                                  <  8:243592/149‑3 (MQ=255)
                                                                                                             gcgcgcAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAgcagcg                                  <  3:343963/149‑1 (MQ=255)
                                                                                                             gcgcgcAGAACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAgcagcg                                  <  6:80817/149‑1 (MQ=255)
                                                                                                               gcgcAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGTAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAgcagcgca                                <  6:20213/149‑1 (MQ=255)
                                                                                                                    gcaCGGCCACCATCCGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTAc                           <  2:416277/149‑1 (MQ=255)
                                                                                                                                tCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGtgctg               <  1:306282/149‑1 (MQ=255)
                                                                                                                                    gcTGTTGGGCACCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGctgctg            >  8:382542/2‑148 (MQ=255)
                                                                                                                                     cTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTgg           >  2:291580/1‑148 (MQ=255)
                                                                                                                                     cTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTgg           >  2:117790/1‑148 (MQ=255)
                                                                                                                                           gTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTGGAAAAgc     <  8:402275/148‑1 (MQ=255)
                                                                                                                                            tGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTGGAAAAgcg    >  4:23394/1‑148 (MQ=255)
                                                                                                                                             gCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTGGAAAAGCGCa  >  1:110921/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CCATGTACTCAACAAGACCCGTCCGGTCAGCGAGCAGGTGCGGCTGAAAGTGGAAGCCGCCATCATCGAACTGGATTACGTGCCCAGTGCGGTGGCGCGCTCGCTGAAAGCGCGCAGCACGGCCACCATCGGCCTGTTGGTGCCCAACAGCGTCAACCCGTACTTCGCCGAGCTGGCGCGGGGTATCGAAGACGCCTGCGAGCGCAACGGCTATTGCGTGATCCTGTGCAACTCCGACGACAACCCGCAGAAGCAGCGCAGCTACCTGCGCGTGCTGCTGGAAAAGCGCA  >  NC_002947/2804244‑2804533

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