Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,574,228 C→G 6.7% R300R (CGC→CGG) ‡ glgX → glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,574,2280CG 6.7% 190.3 / 2.9 59R300R (CGC→CGG) ‡glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base C (32/23);  new base G (2/2);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAA  >  NC_002947/4574090‑4574376
                                                                                                                                          |                                                                                                                                                    
cTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGGGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACa                                                                                                                                            <  6:268100/149‑1 (MQ=255)
     tGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGcc                                                                                                                                       >  6:194392/1‑149 (MQ=255)
                 aGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGc                                                                                                                           >  4:358290/1‑149 (MQ=255)
                      gCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGa                                                                                                                       <  7:117757/148‑1 (MQ=255)
                       ccGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATg                                                                                                                     >  4:342592/1‑149 (MQ=255)
                            ttCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGga                                                                                                                <  7:324502/149‑1 (MQ=255)
                               aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACATCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacga                                                                                                             <  1:240908/149‑1 (MQ=255)
                               aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacga                                                                                                             >  3:164295/1‑149 (MQ=255)
                                 ggAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCg                                                                                                                  >  5:367019/1‑142 (MQ=255)
                                 ggAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCg                                                                                                                  <  6:367019/142‑1 (MQ=255)
                                  gagaTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                            >  6:130826/1‑147 (MQ=255)
                                      tGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGc                                                                                                                    <  4:203425/135‑1 (MQ=255)
                                      tGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGc                                                                                                                    >  3:203425/1‑135 (MQ=255)
                                      tGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgcgc                                                                                                      <  4:168582/149‑1 (MQ=255)
                                       cgTCGTGCAGGTGCGCGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgcgc                                                                                                      >  3:3158/2‑148 (MQ=255)
                                            cgcACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCtactac                                                                                                <  5:194392/149‑1 (MQ=255)
                                             gcACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACa                                                                                               <  7:50761/149‑1 (MQ=255)
                                               aCCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACATCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTAcc                                                                                                                             >  1:265238/1‑117 (MQ=255)
                                               aCCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTAcc                                                                                                                             <  2:265238/117‑1 (MQ=255)
                                               aCCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATc                                                                                             <  3:264716/149‑1 (MQ=255)
                                                    cACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                          <  3:342592/147‑1 (MQ=255)
                                                    cACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTAATGCCGGACGACAAGCGCTACTACCTCACCGa                                                                                        >  6:262999/1‑149 (MQ=255)
                                                     acgacgCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGAc                                                                                       >  7:37747/1‑149 (MQ=255)
                                                        acgCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTcc                                                                                    <  5:164474/149‑1 (MQ=255)
                                                                  tGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcaccggc                                                                           >  4:202906/1‑148 (MQ=255)
                                                                    aGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAac                                                                        <  5:130826/149‑1 (MQ=255)
                                                                      gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGAt                                                                                                                      >  2:72658/1‑101 (MQ=255)
                                                                      gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGAt                                                                                                                      <  1:72658/101‑1 (MQ=255)
                                                                       tGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACAcc                                                                     <  5:319640/149‑1 (MQ=255)
                                                                               gACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCt                                                             >  3:34747/1‑149 (MQ=255)
                                                                                aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcaccg                                                                             <  6:311106/132‑1 (MQ=255)
                                                                                aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccggcaccg                                                                             >  5:311106/1‑132 (MQ=255)
                                                                                 cGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                           <  3:202906/149‑1 (MQ=255)
                                                                                    ggTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                        >  2:228760/1‑101 (MQ=255)
                                                                                    ggTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                        <  1:228760/101‑1 (MQ=255)
                                                                                     gTCTACAACCACACCGCCGAAGGCAACGCGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGCCTCCGGCACCGGCAACACCCTGGACCTGCGCCa                                                       >  7:122187/1‑149 (MQ=255)
                                                                                       cTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                          >  2:116477/1‑112 (MQ=255)
                                                                                       cTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                          <  1:116477/112‑1 (MQ=255)
                                                                                        tACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                    >  3:109223/1‑149 (MQ=255)
                                                                                         aCAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACCCCCTGGACCTGAGCCAcccc                                                   >  6:289041/1‑149 (MQ=255)
                                                                                          cAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCt                                                  >  1:228680/1‑149 (MQ=255)
                                                                                             ccACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCCCCCCTGCg                                               >  4:93053/1‑149 (MQ=255)
                                                                                                  ccgccgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                          <  4:164295/149‑1 (MQ=255)
                                                                                                   cgcagAAGGCAACGAGAGCGGCCCGACCCTGTCGAAGCGCCGCATCGACAACGCCTCGTACTACAGCCTGATGCCGGACGACAAGCGCTACTACAACAACGACTCCGGCACCGGCAACACCCTGGACCAGACCCACCCCTCCGtgatg                                          >  2:82758/1‑148 (MQ=255)
                                                                                                        aaGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCGTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGTCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                    <  1:82758/149‑1 (MQ=255)
                                                                                                           gCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTc                                 <  8:37747/149‑1 (MQ=255)
                                                                                                                     cgGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgc                       <  2:228680/149‑1 (MQ=255)
                                                                                                                        cccGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttc                    >  6:373889/1‑149 (MQ=255)
                                                                                                                          cGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcgc                  >  2:118600/1‑149 (MQ=255)
                                                                                                                              ccTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACt              >  3:294569/1‑149 (MQ=255)
                                                                                                                                tGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAATACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg            >  7:35693/1‑149 (MQ=255)
                                                                                                                                tGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg            >  8:41891/1‑149 (MQ=255)
                                                                                                                                 gTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCCCCGACTCGCTTCGCTACTggg           >  8:146930/1‑149 (MQ=255)
                                                                                                                                 gTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTggg           >  3:12181/1‑149 (MQ=255)
                                                                                                                                  tCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGc          >  2:306146/1‑149 (MQ=255)
                                                                                                                                  tCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg            >  6:372007/1‑147 (MQ=255)
                                                                                                                                  tCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg            <  5:372007/147‑1 (MQ=255)
                                                                                                                                  tCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTCGGc          >  4:143900/1‑149 (MQ=255)
                                                                                                                                    gATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCg        >  5:138079/1‑149 (MQ=255)
                                                                                                                                          cgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGaaa  <  4:109223/149‑1 (MQ=255)
                                                                                                                                          cgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGaaa  <  4:34747/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAA  >  NC_002947/4574090‑4574376

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: