Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,986,815 T→C 7.1% V47A (GTA→GCA)  flgM → flagellin biosynthesis negative regulator FlgM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,986,8150TC 7.1% 158.0 / 3.1 56V47A (GTA→GCA) flgMflagellin biosynthesis negative regulator FlgM
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base C (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

TTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGATGCCGGCCGGGTCGCCA  >  NC_002947/4986668‑4986958
                                                                                                                                                   |                                                                                                                                               
ttCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTa                                                                                                                                                >  6:368576/1‑149 (MQ=255)
      tCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCt                                                                                                                                           >  3:256966/1‑148 (MQ=255)
        aTGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAg                                                                                                                                        >  6:202263/1‑149 (MQ=255)
          ggTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGc                                                                                                                                       >  3:376260/1‑148 (MQ=255)
              aTCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGTGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCaa                                                                                                                                    >  3:301588/1‑147 (MQ=255)
                  aCTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGTGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGca                                                                                                                              <  4:301588/149‑1 (MQ=255)
                  aCTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGca                                                                                                                              >  3:73372/1‑149 (MQ=255)
                    ttCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGca                                                                                                                              >  2:411085/1‑147 (MQ=255)
                        gTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCAcagcag                                                                                                                        >  8:386065/1‑149 (MQ=255)
                           gTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTg                                                                                                                     <  7:82285/149‑1 (MQ=255)
                           gTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTg                                                                                                                     >  3:311036/1‑149 (MQ=255)
                            tttGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGc                                                                                                                    <  4:376260/149‑1 (MQ=255)
                               gAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGc                                                                                                                    <  1:299036/146‑1 (MQ=255)
                               gAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGc                                                                                                                    >  2:299036/1‑146 (MQ=255)
                               gAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCag                                                                                                                  <  4:265363/148‑1 (MQ=255)
                                    aCTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGAt                                                                                                             <  4:311036/148‑1 (MQ=255)
                                     ctctCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAaga                                                                                                                                  <  3:305812/126‑1 (MQ=255)
                                     ctctCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAaga                                                                                                                                  >  4:305812/1‑126 (MQ=255)
                                      tctcCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTg                                                                                                                     <  3:255331/138‑1 (MQ=255)
                                      tctcCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTg                                                                                                                     >  4:255331/1‑138 (MQ=255)
                                                  aCGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTgcg                                                                                              <  5:202263/149‑1 (MQ=255)
                                                  aCGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTgc                                                                                               <  7:386065/148‑1 (MQ=255)
                                                      gcggcgTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGcgac                                                                                          <  2:238033/149‑1 (MQ=255)
                                                      gcggcgTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGcgac                                                                                          <  7:50062/149‑1 (MQ=255)
                                                             gcgcGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCg                                                                                   <  4:314915/149‑1 (MQ=255)
                                                                    aaCACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCaa                                                                            >  7:125747/1‑149 (MQ=255)
                                                                      cacCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACa                                                                          <  4:256966/149‑1 (MQ=255)
                                                                       acCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGAGCGACGAGCCCGTGGTCAACAg                                                                         >  6:382922/1‑149 (MQ=255)
                                                                       acCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGCACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAg                                                                         >  1:142519/1‑149 (MQ=255)
                                                                       acCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGCACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAg                                                                         >  1:137835/1‑149 (MQ=255)
                                                                            gTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGccc                                                                    <  4:73372/149‑1 (MQ=255)
                                                                             tCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGAc                                                                                                      <  5:150104/114‑1 (MQ=255)
                                                                             tCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGAc                                                                                                      >  6:150104/1‑114 (MQ=255)
                                                                                ggCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCgtgt                                                                <  1:411085/149‑1 (MQ=255)
                                                                                  cAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGCGGTCAACAGCGCCCGTGTgg                                                              <  5:368576/149‑1 (MQ=255)
                                                                                    gcgcCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGcc                                                            >  4:359011/1‑149 (MQ=255)
                                                                                         gAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGtt                                                       <  4:49128/149‑1 (MQ=255)
                                                                                         gAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGCACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGtt                                                       <  2:137835/149‑1 (MQ=255)
                                                                                         gAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGCACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGtt                                                       <  2:142519/149‑1 (MQ=255)
                                                                                          aaaaaCCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCa                                                          <  8:299160/145‑1 (MQ=255)
                                                                                          aaaaaCCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCa                                                          >  7:299160/1‑145 (MQ=255)
                                                                                            aaaCCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGaa                                                    <  3:326268/149‑1 (MQ=255)
                                                                                                cGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGaa                                                    >  8:44435/1‑145 (MQ=255)
                                                                                                cGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGaa                                                    <  7:44435/145‑1 (MQ=255)
                                                                                                cGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAg                                                >  3:314126/1‑149 (MQ=255)
                                                                                                  gCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGc                                              >  4:228901/1‑149 (MQ=255)
                                                                                                          ccagGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCg                                      >  1:146180/1‑149 (MQ=255)
                                                                                                                 gCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAg                               <  3:359011/149‑1 (MQ=255)
                                                                                                                      gAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTAc                           >  6:368894/1‑148 (MQ=255)
                                                                                                                         aGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTgg                                                              <  7:327726/110‑1 (MQ=255)
                                                                                                                         aGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTgg                                                              >  8:327726/1‑110 (MQ=255)
                                                                                                                              ccagcgccagcgGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGAt                  <  5:410720/149‑1 (MQ=255)
                                                                                                                                     cagcgGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGAGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGATgccggcc           >  4:209727/1‑149 (MQ=255)
                                                                                                                                       gcgGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGATGccggccgg         <  5:382922/149‑1 (MQ=255)
                                                                                                                                           agagGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGATGCCGGCCGGGTCg     <  3:209727/149‑1 (MQ=255)
                                                                                                                                              ggCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGATGCCGGCCGGGTCGCCa  <  8:125747/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TTCTGATCATGGTCATCGACTTCAGTCGTTTGAATAACTCTCCGTCCGTCACGGGCGGCGTGCGCGGCAACACCGCGTCCGGCAGCGCCGAAAAACCGGCAGCCAGCCAGGAAGCGGTGAAAGACACCAGCGCCAGCGGAGAGGCGGTACACCTCAGCCAAGAGGCACAGCAGTTGCAGAAGATCAGCGACAAGCTGCGCGACGAGCCCGTGGTCAACAGCGCCCGTGTGGCCCAGTTGAAGCAGGCGATCGCCGACGGCAGCTACCAGGTCGATGCCGGCCGGGTCGCCA  >  NC_002947/4986668‑4986958

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: