Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,895,977 A→C 6.4% D404A (GAC→GCC)  PP_1700 → acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,895,9770AC 6.4% 171.0 / 2.4 58D404A (GAC→GCC) PP_1700acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base C (2/2);  total (29/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 = 8.31e-01

GCTGCAGCAGCGCCGCACGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAATATCGCG  >  NC_002947/1895832‑1896122
                                                                                                                                                 |                                                                                                                                                 
gCTGCAGCAGCGCCGCACGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGa                                                                                                                                                <  2:195718/149‑1 (MQ=255)
      gcagcGCCGCACGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCAAGGCCCAGGACgagcgagc                                                                                                                                          >  1:185077/1‑149 (MQ=255)
            ccgcACGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCTGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCCCGGCCGAGGCCGAGCGCGCCAgggg                                                                                                                                    >  6:296387/1‑148 (MQ=255)
              gcACGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCGCGCTCGTCCTCGGCCGTGCGGGCCTggatca                                                                                                                                  <  2:156762/149‑5 (MQ=38)
                aCGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATc                                                                                                                                <  6:28217/149‑1 (MQ=255)
                    ctggctggCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCgg                                                                                                                            <  3:114249/149‑1 (MQ=255)
                         tggcACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCGCGCTCGTCCTCGGCCGTGCGGGCCTggatcaac                                                                                                                                >  3:122278/1‑134 (MQ=37)
                         tggcACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCGCGCTCGTCCTCGGCCGTGCGGGCCTggatcaac                                                                                                                                <  4:122278/140‑7 (MQ=37)
                         tggcACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCcaa                                                                                                                        >  8:172957/1‑148 (MQ=255)
                            cACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCaacaac                                                                                                                    <  2:78964/149‑1 (MQ=255)
                               tgctgGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATc                                                                                                                 >  5:217970/1‑149 (MQ=255)
                                       aTGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGc                                                                                                         >  5:269302/1‑149 (MQ=255)
                                       aTGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGc                                                                                                         >  2:74595/1‑149 (MQ=255)
                                         gAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCt                                                                                                       >  7:48477/1‑149 (MQ=255)
                                          aGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCtt                                                                                                      >  3:164635/1‑149 (MQ=255)
                                           gCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTg                                                                                                     >  3:76963/1‑149 (MQ=255)
                                                     gtggTTCTTGCTGTGGTTCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTTCCACTGTATGCGTTGATCCAGGCCCGCACGTCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCAtggtg                                                                                            <  2:185077/148‑1 (MQ=255)
                                                     gtggATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCAtggtgg                                                                                           >  8:113514/1‑149 (MQ=255)
                                                     gtggATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCAtggtgg                                                                                           <  4:76963/149‑1 (MQ=255)
                                                         aTCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAgcgc                                                                                                                                            <  2:246658/96‑1 (MQ=255)
                                                         aTCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTAATCCAGGCCCGCACGGCCGAGGACGAgcgc                                                                                                                                            >  1:246658/1‑96 (MQ=255)
                                                          tCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTcg                                                                                      >  8:229395/1‑149 (MQ=255)
                                                                ggTCGCTCATTGGCCGGGGGGTGGTCGGCGGCTTCTACCTCGTACCCCTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGt                                                                                 <  5:296387/147‑1 (MQ=255)
                                                                        aTTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACca                                                                         <  8:48477/148‑1 (MQ=255)
                                                                            gCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCatcatc                                                                    <  2:150477/149‑1 (MQ=255)
                                                                             ccTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCTCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGCTGTCGGCGGTGCTCACCATCATcc                                                                   >  1:129513/1‑149 (MQ=255)
                                                                                 ggTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGc                                                                >  4:9092/1‑148 (MQ=255)
                                                                                  gtgtgtTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGc                                                                >  3:136511/1‑147 (MQ=255)
                                                                                     tgtTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGc                                                           <  7:113514/149‑1 (MQ=255)
                                                                                     tgtTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGc                                                           <  6:217970/149‑1 (MQ=255)
                                                                                          ggcggcTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGc                                                      <  3:9092/149‑1 (MQ=255)
                                                                                               cTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGc                                                 >  7:176445/1‑149 (MQ=255)
                                                                                                ttCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCt                                                >  7:7777/1‑149 (MQ=255)
                                                                                                  cTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCt                                                >  6:224031/1‑147 (MQ=255)
                                                                                                   tACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAg                                             >  1:21904/1‑149 (MQ=255)
                                                                                                    aCATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCAtg                                                                                               <  4:149303/98‑1 (MQ=255)
                                                                                                     cATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGc                                            <  4:136511/148‑1 (MQ=255)
                                                                                                     cATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCa                                           >  7:279044/1‑149 (MQ=255)
                                                                                                       tCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATc                                         <  7:172957/149‑1 (MQ=255)
                                                                                                          tACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCg                                      >  3:130304/1‑149 (MQ=255)
                                                                                                           aCCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCg                                      <  4:164635/148‑1 (MQ=255)
                                                                                                            ccACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCa                                    <  1:69444/149‑1 (MQ=255)
                                                                                                                   aTGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAg                                                             <  6:81247/117‑1 (MQ=255)
                                                                                                                   aTGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAg                                                             >  5:81247/1‑117 (MQ=255)
                                                                                                                     gCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGATAGGCGTTTACAAGCTGAGCATCCCCAAGCTTTTACt                           >  8:280021/1‑149 (MQ=255)
                                                                                                                      cGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCtg                          <  7:229395/149‑1 (MQ=255)
                                                                                                                       gTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCt                                                <  6:187644/126‑1 (MQ=255)
                                                                                                                       gTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCt                                                >  5:187644/1‑126 (MQ=255)
                                                                                                                              cAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTc                  >  3:12183/1‑149 (MQ=255)
                                                                                                                              cAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTc                  <  4:130304/149‑1 (MQ=255)
                                                                                                                              cAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTc                  >  2:225415/1‑149 (MQ=255)
                                                                                                                                    cGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCATTTCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCgctgct            <  5:4158/149‑1 (MQ=255)
                                                                                                                                      cACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTc           >  6:150291/1‑148 (MQ=255)
                                                                                                                                         ggCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAata       <  1:225415/149‑1 (MQ=255)
                                                                                                                                           ccGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAATATc     <  4:12183/149‑1 (MQ=255)
                                                                                                                                            cGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAATATcg    <  1:74595/149‑1 (MQ=255)
                                                                                                                                              aGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAATATcgcg  <  8:176445/149‑1 (MQ=255)
                                                                                                                                              aGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAATATcgcg  >  8:39380/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCTGCAGCAGCGCCGCACGACTGGCTGGCACTGCTGGGCATGAGCCAGGCCTGGTGGATCTTGCTGTCGATCATTGGCCTGGGTGTGTTCGGCGGCTTCTACATCGTACCACTGTATGCGTTGATCCAGGCCCGCACGGCCGAGGACGAGCGCGCCAGGGTGATCGCGGCCAACAACATCCTCAATGCCTTGTTCATGGTGGTGTCGGCGGTGCTCACCATCATCCTGCTAGGCGTGGCCAAGCTGAGCATCCCGCAGCTTTTCCTGGTGGTGTCGCTGCTCAATATCGCG  >  NC_002947/1895832‑1896122

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