Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,700,758 T→C 6.9% intergenic (‑71/+79) kdpD ← / ← kdpC sensor protein/potassium‑transporting ATPase subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,700,7580TC 6.9% 171.2 / 2.3 58intergenic (‑71/+79)kdpD/kdpCsensor protein/potassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base C (2/2);  total (32/26)
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.82e-01

GAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGA  >  NC_002947/4700613‑4700904
                                                                                                                                                  |                                                                                                                                                  
gAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTAcc                                                                                                                                                 <  4:215043/149‑1 (MQ=255)
   cacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                              <  6:165202/149‑1 (MQ=255)
      cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           >  5:65268/1‑149 (MQ=255)
      cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGaa                                                                                                                                           >  3:94981/1‑149 (MQ=255)
      cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATACTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCAGCGaa                                                                                                                                           >  6:140453/1‑149 (MQ=255)
       ttCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           >  2:363653/1‑148 (MQ=255)
       acgcaccggccggggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCACCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                           >  5:105033/6‑148 (MQ=255)
             cggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCa                                                                                                                                    <  8:334573/149‑1 (MQ=255)
                ccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGATAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGt                                                                                                                                  <  1:82845/149‑1 (MQ=255)
                 cccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTgg                                                                                                                                >  5:214542/1‑149 (MQ=255)
                 cccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTgg                                                                                                                                >  7:330640/1‑149 (MQ=255)
                      ccggTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAACCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcacca                                                                                                                           <  1:359426/149‑1 (MQ=255)
                               gggAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcacc                                                                                                                            >  1:164828/1‑139 (MQ=255)
                               gggAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcacc                                                                                                                            <  2:164828/139‑1 (MQ=255)
                                     ttGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                            <  1:12085/149‑1 (MQ=255)
                                        gCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTg                                                                                                         <  7:162198/149‑1 (MQ=255)
                                            aCAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCACCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAAcc                                                                                                     >  8:295636/1‑149 (MQ=255)
                                                 gcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTc                                                                                                <  6:214542/149‑1 (MQ=255)
                                                 gcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAACCCGCTc                                                                                                <  4:94981/149‑1 (MQ=21)
                                                    tCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTc                                                                                                <  7:119772/146‑1 (MQ=255)
                                                    tCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTc                                                                                                >  8:119772/1‑146 (MQ=255)
                                                       gcgcGGGCGGAGTCACTCATGTTTCATCCATTACTTTGTCAAACAGTCAGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTAACCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAg                                                                                          >  1:88224/1‑149 (MQ=255)
                                                         gcgggcggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGc                                                                                        >  4:180585/1‑149 (MQ=255)
                                                            ggcggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCACCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTc                                                                                     <  7:295636/149‑1 (MQ=255)
                                                                gAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGc                                                                                 <  3:276485/149‑1 (MQ=255)
                                                                      cTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                            <  7:44996/116‑1 (MQ=18)
                                                                      cTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                            >  8:44996/1‑116 (MQ=18)
                                                                                 tCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCt                                                                >  8:332105/1‑149 (MQ=255)
                                                                                  ccTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcac                                                                                                                             <  4:75530/87‑1 (MQ=255)
                                                                                  ccTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcac                                                                                                                             >  3:75530/1‑87 (MQ=255)
                                                                                    tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATa                                                                                                                                         >  5:343040/1‑73 (MQ=255)
                                                                                    tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATa                                                                                                                                         <  6:343040/73‑1 (MQ=255)
                                                                                        cTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                  <  1:374802/140‑1 (MQ=255)
                                                                                        cTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                  >  2:374802/1‑140 (MQ=255)
                                                                                         tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATgg                                                        >  7:53304/1‑149 (MQ=255)
                                                                                         tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATgg                                                        >  1:84644/1‑149 (MQ=255)
                                                                                             tCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAg                                                    >  6:251930/1‑149 (MQ=255)
                                                                                              ccAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGc                                                   >  6:72549/1‑149 (MQ=255)
                                                                                               cAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTg                                                                                  >  1:80616/1‑117 (MQ=255)
                                                                                               cAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTg                                                                                  <  2:80616/117‑1 (MQ=255)
                                                                                               cAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCGTCGATTACCTGACAATGGCGAGCt                                                  >  6:189440/1‑149 (MQ=255)
                                                                                                 aCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcaccac                                                                                                                          >  1:300933/1‑75 (MQ=255)
                                                                                                 aCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcaccac                                                                                                                          <  2:300933/75‑1 (MQ=255)
                                                                                                   aGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTc                                              <  1:363653/149‑1 (MQ=255)
                                                                                                     tCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGAGTCGATTACCGGACAATGGCGAGCTTATcc                                             >  8:291473/1‑148 (MQ=255)
                                                                                                       cgccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAg                                          <  6:105033/149‑1 (MQ=255)
                                                                                                       cgccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAg                                          >  5:228169/1‑149 (MQ=255)
                                                                                                           gcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCt                                      <  7:332105/149‑1 (MQ=255)
                                                                                                                 gggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTc                                 >  5:182301/1‑148 (MQ=255)
                                                                                                                   gAAAGCCCCCTTAACGTGTAGGAGCGGTTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCTGGTTAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTTACAATGGCTAGCTTTTCCAAGGCCTGGTTCAAc                              <  5:189440/149‑1 (MQ=255)
                                                                                                                        gCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCa                          <  5:251930/149‑1 (MQ=255)
                                                                                                                                aaCGTGTAGGAGCGGGTTTTCCCGCGAATACGTGAGTGTCACCACCACCGCATTCGCGTGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGtt                   <  2:88224/148‑1 (MQ=255)
                                                                                                                                 aCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGa                 >  5:196160/1‑149 (MQ=255)
                                                                                                                                        ggAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcgg          <  6:196160/149‑1 (MQ=255)
                                                                                                                                         gAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggc         >  8:32747/1‑149 (MQ=255)
                                                                                                                                          aGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggcg        >  7:140778/1‑149 (MQ=255)
                                                                                                                                           gCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggcgg       >  5:308511/1‑149 (MQ=255)
                                                                                                                                            cGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCggg      >  4:181825/1‑149 (MQ=255)
                                                                                                                                                 ttACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGa  <  6:228169/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGA  >  NC_002947/4700613‑4700904

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: