Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,103,469 A→T 7.0% K219M (AAG→ATG)  PP_2713 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,103,4690AT 7.0% 171.2 / 3.2 57K219M (AAG→ATG) PP_2713two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base T (2/2);  total (30/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.08e-01

CGTGAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGT  >  NC_002947/3103324‑3103616
                                                                                                                                                 |                                                                                                                                                   
cGTGAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAgg                                                                                                                                                   >  1:225764/1‑148 (MQ=255)
   gAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCAt                                                                                                                                                >  2:289488/1‑148 (MQ=255)
    aGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGa                                                                                                                                              <  1:7352/149‑1 (MQ=255)
        ccTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTc                                                                                                                                          <  1:289488/149‑1 (MQ=255)
          tGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGc                                                                                                                                        >  3:266875/1‑149 (MQ=255)
            aCTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCATGCCTTGACCTCgac                                                                                                                                       <  7:174388/148‑3 (MQ=255)
            aCTCAACTGCTCTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCATGCCTTGACCTCga                                                                                                                                        <  6:307457/147‑2 (MQ=255)
                 aCTGCTGTATGGCTGGCTTGAAGAGCCGGAGAGCAAGACGCTGGAGGTGAAGATCTACCATTTGCCCATGATGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGGCGAGGTCAAGGGATGACCTCGCTTCGCCa                                                                                                                                 <  2:93644/149‑1 (MQ=255)
                   tgctgTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAgc                                                                                                                               >  7:284428/1‑149 (MQ=255)
                      tgTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCa                                                                                                                            >  8:121105/1‑149 (MQ=255)
                      tgTATGCCTGGGATGAAGAGCCGGAGAGCAACACGCTTGAGGTAAACATCTACCATTTGCTAAAGAAGCTTTTCAACGGACTGATTCGCACGGTGCGCGGAAATCGATATCTGTTAGAGGTCAAGGCATGACCTAACTTAGCCAgagaa                                                                                                                            >  8:246316/1‑145 (MQ=255)
                            gCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCtgt                                                                                                                      <  7:121105/149‑1 (MQ=255)
                             cTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCtgtg                                                                                                                     >  2:86375/1‑149 (MQ=255)
                              tGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTgg                                                                                                                    >  4:279609/1‑149 (MQ=255)
                                ggATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGgcg                                                                                                                  >  1:118458/1‑149 (MQ=255)
                                   tGAAGAGCCGGAGAGCAACAAGCTGCAGCTGACCAGCTACCATTTGCGCAAGAAGCTGTTCAACCGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGt                                                                                                               <  4:266875/149‑1 (MQ=255)
                                          ccGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGctc                                                                                                        >  4:340568/1‑149 (MQ=255)
                                           cGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTg                                                                                                              <  5:185438/142‑1 (MQ=255)
                                           cGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTg                                                                                                              >  6:185438/1‑142 (MQ=255)
                                               gagCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTtgct                                                                                                   <  2:225764/149‑1 (MQ=255)
                                                agCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTtgctg                                                                                                  <  4:330151/149‑1 (MQ=255)
                                                 gCAACACGCTGGAGGTGAACATCTACCATTTGCGAAAGAAGCTTTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTtgctgc                                                                                                 <  7:246316/149‑1 (MQ=255)
                                                      acGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTc                                                                                            <  3:340568/149‑1 (MQ=255)
                                                      acGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTc                                                                                            <  8:284428/149‑1 (MQ=255)
                                                        gCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGGGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCgg                                                                                          <  4:61070/149‑1 (MQ=255)
                                                               gTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGc                                                                                   <  3:279609/149‑1 (MQ=255)
                                                                        tACCATTTGCGCGAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                          >  2:258298/1‑149 (MQ=255)
                                                                        tACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                          >  8:222418/1‑149 (MQ=255)
                                                                         aCCATTTGCGCAAGAAGCTGTTCAACTGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGc                                                                                   <  4:151932/139‑1 (MQ=255)
                                                                         aCCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                          <  3:83912/148‑1 (MQ=255)
                                                                         aCCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTcgcg                                                                          >  4:83912/1‑148 (MQ=255)
                                                                         aCCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCgcgc                                                                         >  1:315214/1‑149 (MQ=255)
                                                                         aCCATTTGCGCAAGAAGCTGTTCAACCGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGc                                                                                   >  3:151932/1‑139 (MQ=255)
                                                                          ccATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCgcgc                                                                         <  1:47056/148‑1 (MQ=255)
                                                                                  gcAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAAc                                                                 >  8:221423/1‑148 (MQ=255)
                                                                                     agaagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTAcc                                                             <  6:59085/149‑1 (MQ=255)
                                                                                       cagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCAc                                                           <  4:34952/148‑1 (MQ=255)
                                                                                       aagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCa                                                            >  1:325715/1‑148 (MQ=255)
                                                                                       aagCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCa                                                            <  2:325715/148‑1 (MQ=255)
                                                                                          cTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGACGCTCtt                                                                                                      <  2:8098/103‑1 (MQ=255)
                                                                                          cTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGAAGCTCtt                                                                                                      >  1:8098/1‑103 (MQ=255)
                                                                                                            cGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTaa                                       <  7:221423/148‑1 (MQ=255)
                                                                                                                   tGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCGCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCCCGACAGCCGCCACGAAATTAACGAGGta                                >  6:46220/1‑148 (MQ=255)
                                                                                                                            ttGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGAcgcg                       >  5:125239/1‑148 (MQ=255)
                                                                                                                               gCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATaa                                                                  >  8:130673/1‑102 (MQ=255)
                                                                                                                               gCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATaa                                                                  <  7:130673/102‑1 (MQ=255)
                                                                                                                                 tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGa                                    >  6:150128/1‑130 (MQ=255)
                                                                                                                                 tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGa                                    <  5:150128/130‑1 (MQ=255)
                                                                                                                                 tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTg                 >  4:156089/1‑149 (MQ=255)
                                                                                                                                 tATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGAACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTg                 >  8:235265/1‑149 (MQ=255)
                                                                                                                                   tCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGc               >  3:275329/1‑149 (MQ=255)
                                                                                                                                        gcGAGGTCATGCCTTGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAGCGAGGTATATGACGCGCACCTGGCCCAg           >  2:115004/2‑148 (MQ=255)
                                                                                                                                         cGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGa          <  7:256463/148‑1 (MQ=255)
                                                                                                                                         cGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGa          >  8:256463/1‑148 (MQ=255)
                                                                                                                                          gAGGTCATGCCTTGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAAc        >  8:63849/1‑149 (MQ=255)
                                                                                                                                          gAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGaa         <  2:118458/148‑1 (MQ=255)
                                                                                                                                             gTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAAcgcg     >  5:274331/1‑149 (MQ=255)
                                                                                                                                                aaGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGt  >  7:138311/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGTGAGCGCCTGACTCAACTGCTGTATGGCTGGGATGAAGAGCCGGAGAGCAACACGCTGGAGGTGAACATCTACCATTTGCGCAAGAAGCTGTTCAACGGCCTGATTCGCACGGTGCGCGGAATTGGCTATCTGGTCGAGGTCAAGGCATGACCTCGCTTCGCCAGCGCACGCTGTGGCGCGTGATGCTCTTGCTGCTACTCGGTACCGGCTTGCTCGCGCTCTATAACTACCACGACAGCAGCCACGAAATTAACGAGGTATATGACGCGCACCTGGCCCAGAACGCGCGT  >  NC_002947/3103324‑3103616

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