Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,695,029 C→G 7.3% A332G (GCC→GGC)  wspC → biofilm formation methyltransferase WspC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,695,0290CG 7.3% 169.6 / 3.6 55A332G (GCC→GGC) wspCbiofilm formation methyltransferase WspC
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

GCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCG  >  NC_002947/1694893‑1695170
                                                                                                                                          |                                                                                                                                              
gCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGGCACCGGCATTGGCCAGCCGGGCAATGCCggc                                                                                                                                      <  4:189552/149‑2 (MQ=38)
gCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGc                                                                                                                                    >  6:174137/1‑149 (MQ=255)
 cACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGTCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑ACGGCATTGGCCAGCCGGGCAATGCCggccag                                                                                                                                   <  1:273260/149‑5 (MQ=255)
                    cAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTa                                                                                                                  <  4:111179/146‑1 (MQ=255)
                    cAGTGTGACGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTa                                                                                                                  >  3:111179/1‑146 (MQ=255)
                       tgtgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAAt                                                                                                            <  2:226521/149‑1 (MQ=255)
                         tgCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGc                                                                                                          <  5:134447/149‑1 (MQ=255)
                           ccGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGcc                                                                                                         >  1:209400/1‑148 (MQ=255)
                                gCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGc                                                                                                    >  4:199568/1‑148 (MQ=255)
                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGc                                                                                                                <  1:218542/134‑1 (MQ=255)
                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGc                                                                                                                >  2:218542/1‑134 (MQ=255)
                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCt                                                                                                   >  6:278410/1‑147 (MQ=255)
                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCt                                                                                                   <  5:278410/147‑1 (MQ=255)
                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 <  3:70095/149‑1 (MQ=255)
                                  cTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 >  5:25147/1‑149 (MQ=255)
                                     gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAGCCAGGCCCGTAGCGAATGCCAACGCTAc                                                                                                 <  1:14636/147‑1 (MQ=255)
                                     gTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCt                                                                                               >  8:310934/1‑148 (MQ=255)
                                      ttCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGa                                                                                             >  5:210946/1‑149 (MQ=255)
                                          cggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTCAGCCa                                                                                         <  7:310934/149‑1 (MQ=255)
                                           ggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAAt                                                                                                            >  7:58145/1‑129 (MQ=255)
                                           ggcgTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAAt                                                                                                            <  8:58145/129‑1 (MQ=255)
                                                   ctgcCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACCCTACCTGAGCCAGTACccgcc                                                                                >  8:149299/1‑149 (MQ=255)
                                                              ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑Ca                                                                                                                          >  6:38834/1‑96 (MQ=255)
                                                              ggccaggccaGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑Ca                                                                                                                          <  5:38834/96‑1 (MQ=255)
                                                               gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                    <  6:25147/149‑1 (MQ=255)
                                                               gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                    <  6:164757/149‑1 (MQ=255)
                                                               gccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgt                                                                    >  2:2979/1‑149 (MQ=255)
                                                                ccaggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtg                                                                   <  3:199568/149‑1 (MQ=255)
                                                                 caggccagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGgtgt                                                                  <  2:209400/149‑1 (MQ=255)
                                                                      cagCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                             >  1:337222/1‑149 (MQ=255)
                                                                         cgcgAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActgg                                                          <  6:210946/149‑1 (MQ=255)
                                                                            gAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTActggctg                                                       >  7:123598/1‑149 (MQ=255)
                                                                              gCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTggg                                                     >  8:297403/1‑149 (MQ=255)
                                                                                ccGCGAGGACAGCCTCGTGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCATCGAG‑CAGGGCCGTAGCTAATGCCAACGCTACCTGAGCCAGTACCGGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                                   <  1:146461/149‑1 (MQ=255)
                                                                                ccGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGtt                                                   >  3:149924/1‑149 (MQ=255)
                                                                                    gAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCt                                               >  8:134527/1‑149 (MQ=255)
                                                                                      gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                             <  5:75482/149‑1 (MQ=255)
                                                                                      gCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCa                                             <  1:28510/149‑1 (MQ=255)
                                                                                              aCGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATAcc                                     >  3:136798/1‑149 (MQ=255)
                                                                                                 tgggggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGaa                                  <  4:149924/148‑1 (MQ=255)
                                                                                                    cggCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGc                               >  2:38072/2‑149 (MQ=255)
                                                                                                        cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATg                           <  6:8783/149‑1 (MQ=255)
                                                                                                        cAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATg                           >  8:237912/1‑149 (MQ=255)
                                                                                                             aGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcag                      <  7:149299/149‑1 (MQ=255)
                                                                                                              gTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCgcag                      >  6:64279/1‑148 (MQ=255)
                                                                                                                  cTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAgg                  <  7:297403/148‑1 (MQ=255)
                                                                                                                     g‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCt              >  8:223263/1‑149 (MQ=255)
                                                                                                                         cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                            >  6:236827/1‑133 (MQ=255)
                                                                                                                         cGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGAt                            <  5:236827/133‑1 (MQ=255)
                                                                                                                           gCATTGGCCAGCCGGGCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc          >  2:224633/1‑149 (MQ=255)
                                                                                                                           gCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGc          <  7:223263/149‑1 (MQ=255)
                                                                                                                             aTTGGCCAGCCGGGCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCa        >  1:137765/1‑149 (MQ=255)
                                                                                                                                  ccGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                             >  1:58624/1‑91 (MQ=255)
                                                                                                                                  ccGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGtactac                                                             <  2:58624/91‑1 (MQ=255)
                                                                                                                                    ggCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCg  >  7:67119/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                              
GCACCGGTTTATCCGCAGCCCAGTGTGCCGTTGCCTCGTTCGCGGCGTGTGCTGCCGGTGGCGGCCAGGCCAGCGCGAGCCCGCGAGCACAGCCACGAGGGGGCCAGCGAGTTGCTGG‑‑CCGGCATTGCCCGGCTGGCCAATGCCGGTGCCAGCGAG‑CAGGCCCGTAGCGAATGCCAACGCTACCTGAGCCAGTACCCGCCGTCGGCACAGGTGTACTACTGGCTGGGTTTGCTCAGTGATACCGAAGGCGATGCGCAGCAGGCGCTGAGCCACTATCG  >  NC_002947/1694893‑1695170

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: ‑