Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,573,171 C→A 7.8% S170R (AGC→AGA pagL‑II → lipid A 3‑O‑deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,573,1710CA 7.8% 156.6 / 3.9 51S170R (AGC→AGApagL‑IIlipid A 3‑O‑deacylase
Reads supporting (aligned to +/- strand):  ref base C (23/24);  new base A (2/2);  total (25/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.08e-01

GGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTC  >  NC_002947/3573026‑3573317
                                                                                                                                                 |                                                                                                                                                  
ggTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTTCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTacagcaca                                                                                                                                                 <  1:257235/149‑1 (MQ=255)
      cggTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTc                                                                                                                                           <  1:197263/149‑1 (MQ=255)
        gTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAACGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCt                                                                                                                                          <  1:289570/148‑1 (MQ=255)
        gTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCt                                                                                                                                          >  2:289570/1‑148 (MQ=255)
                tGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTAc                                                                                                                                  <  6:336428/148‑1 (MQ=255)
                  aaGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACgg                                                                                                                                <  7:281668/148‑1 (MQ=255)
                  aaGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACgg                                                                                                                                <  1:351549/148‑1 (MQ=255)
                         tATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCAAAGAGTCGTACTCGCTGTTCTacagcaca                                                                                                                                                 <  6:16870/124‑1 (MQ=255)
                         tATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCAAAGAGTCGTACTCGCTGTTCTacagcaca                                                                                                                                                 >  5:16870/1‑124 (MQ=255)
                               cTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGc                                                                                                                  <  7:93323/149‑1 (MQ=255)
                                    gCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCgg                                                                                                             <  6:212131/149‑1 (MQ=255)
                                     ccTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCgg                                                                                                             >  1:224628/1‑148 (MQ=255)
                                           gTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAAtt                                                                                                        >  4:6326/1‑147 (MQ=255)
                                           gTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAAtt                                                                                                        <  3:6326/147‑1 (MQ=255)
                                           gTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCa                                                                                                      <  3:357626/149‑1 (MQ=255)
                                               ccGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGCGGCCACTGCGCGGCCCCGGCAATTCAGaa                                                                                                   <  5:240487/148‑1 (MQ=255)
                                                cGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAgg                                                                                                 <  8:46665/149‑1 (MQ=255)
                                                        cAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGa                                                                                          >  6:135652/1‑148 (MQ=255)
                                                          gAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAAc                                                                                       >  3:368526/1‑149 (MQ=255)
                                                              gTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCg                                                                                   <  7:179819/149‑1 (MQ=255)
                                                                 ggCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGc                                                                                 <  1:232666/148‑1 (MQ=255)
                                                                    aTCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTgcgc                                                                                                                     >  1:272836/1‑109 (MQ=255)
                                                                    aTCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTgcgc                                                                                                                     <  2:272836/109‑1 (MQ=255)
                                                                    aTCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTg                                                                              >  4:372106/1‑148 (MQ=255)
                                                                      ccGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCAc                                                                           >  2:91098/1‑149 (MQ=255)
                                                                            cATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTgcg                                                                     >  4:215550/1‑149 (MQ=255)
                                                                              tGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGcgcg                                                                   >  3:287213/1‑149 (MQ=255)
                                                                               gCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCgg                                                                  >  1:186463/1‑149 (MQ=255)
                                                                                 actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGc                                                                                                            <  3:30984/105‑1 (MQ=255)
                                                                                 actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGc                                                                                                            >  4:30984/1‑105 (MQ=255)
                                                                                 actactCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGt                                                                 <  2:186463/148‑1 (MQ=255)
                                                                                        cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCggg                                                         >  5:47674/1‑149 (MQ=255)
                                                                                        cAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGC‑‑GCACCGCGCACGGTgcaggccg                                                         <  1:330946/147‑7 (MQ=255)
                                                                                            gCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGtat                                                     <  7:347546/149‑1 (MQ=255)
                                                                                            gCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGta                                                      >  1:101239/1‑148 (MQ=255)
                                                                                                 gATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATa                                                <  5:135652/149‑1 (MQ=255)
                                                                                                     aaGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCt                                             <  3:372106/148‑1 (MQ=255)
                                                                                                       gcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGtt                                          <  2:101239/149‑1 (MQ=255)
                                                                                                       gcagcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGtt                                          >  1:308792/1‑149 (MQ=255)
                                                                                                          gcCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGaa                                       >  8:701/1‑149 (MQ=255)
                                                                                                           ccAAACGGCGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAAc                                      >  4:142318/1‑149 (MQ=255)
                                                                                                                     gCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGTGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                            <  3:215550/149‑1 (MQ=255)
                                                                                                                      cATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                            >  3:65501/1‑148 (MQ=255)
                                                                                                                      cATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGt                            >  5:49991/1‑148 (MQ=255)
                                                                                                                       aTAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCa                          >  6:157906/1‑149 (MQ=255)
                                                                                                                             tCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATAcccc                    <  2:224628/149‑1 (MQ=255)
                                                                                                                                agaacGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTc                              >  5:220699/5‑136 (MQ=255)
                                                                                                                                agaacGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTc                              <  6:220699/132‑1 (MQ=255)
                                                                                                                                  cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgtt               >  6:219476/1‑149 (MQ=255)
                                                                                                                                  cTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgt                <  4:65501/148‑1 (MQ=255)
                                                                                                                                    cGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCttgttg              >  1:201171/1‑148 (MQ=255)
                                                                                                                                    cGCTGTGCTGTAGAACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGa             <  2:201171/149‑1 (MQ=255)
                                                                                                                                           cTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGtt      >  6:262045/1‑149 (MQ=255)
                                                                                                                                           cTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGtt      >  7:257654/1‑149 (MQ=255)
                                                                                                                                             acagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGtttt    >  8:221697/1‑149 (MQ=255)
                                                                                                                                              cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCg                             <  2:10495/123‑1 (MQ=255)
                                                                                                                                              cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCg                             >  1:10495/1‑123 (MQ=255)
                                                                                                                                              cagcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGttt     <  4:368526/147‑1 (MQ=255)
                                                                                                                                               agcacagcGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTc  <  3:142318/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTCGGCGGTCAACTTTGAAGACCGTATCGGCTTTGGCCTGAAGTTGCCGGGAGAACAGAAGGTAGGCATCCGTGCCATGCACTACTCCAACGCCGGGATAAAGCAGCCAAACGACGGCATAGAGTCGTACTCGCTGTTCTACAGCACAGCGTTCTGACATTACGGGGCCACTGCGCGGCCCCGGCAATTCAGAAGGTATAGGAAACGCCGGCCTGCACCGTGCGCGGTGCGCCGGGGTATACATAGCTGTTGAACGCGCCCTCGTCATACCCCTTGTTGAACACGTTTTTC  >  NC_002947/3573026‑3573317

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: