Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,266,921 G→C 8.9% P356R (CCG→CGG)  quiP ← acyl‑homoserine lactone acylase QuiP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,266,9210GC 8.9% 165.5 / 7.5 56P356R (CCG→CGG) quiPacyl‑homoserine lactone acylase QuiP
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base C (3/2);  total (28/28)
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.53e-01

CTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGCACAATTGGCAGGCCCG  >  NC_002947/1266778‑1267068
                                                                                                                                               |                                                                                                                                                     
cTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGc                                                                                                                                                  >  2:187699/1‑149 (MQ=255)
cTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCCGGCCCGCGCCCGCAGc                                                                                                                                                  <  4:5347/149‑1 (MQ=255)
       cTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTAc                                                                                                                                           >  6:216545/1‑149 (MQ=255)
          gggCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTAcc                                                                                                                                          >  5:204161/1‑147 (MQ=255)
          gggCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTAcc                                                                                                                                          <  6:204161/147‑1 (MQ=255)
          gggCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCAt                                                                                                                                        >  2:283638/1‑149 (MQ=255)
            gCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCg                                                                                                                                      >  2:267738/1‑149 (MQ=255)
            gCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGGGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCg                                                                                                                                      <  2:175568/149‑1 (MQ=255)
            gCAGGTGCCGCGCCAGGCCCAGCCTGGGGTTTTCGGGCGGGCCCAGCAGGGGGCCGGGGCGGGGGTCGTACCTCACTTCGCGCACGGGGCGTTGCCCACGTACAAAGTTGGTTTCGCTGCGGGCGCGGGCGGGCAGCCATTTACCATCg                                                                                                                                      <  6:293738/149‑1 (MQ=255)
             gaGGTTCAGCGCCAGGCACAGGCTTGCGTTTTCTGGCTGGGCCAGCAGGGTGGCGTGGCTGGTGTAGTTCATCACTTAGCGCACGGGGCGTTGGCCACGTACATAGTAGGTTTCGGTTCTGGCGAGGGCCGTCAGCCATTTACGATCgg                                                                                                                                     <  6:229804/148‑1 (MQ=255)
               ggTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGcc                                                                                                                                   <  1:283638/149‑1 (MQ=255)
                gTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGTGTGCCGTGGCGGGAGTAGTACATCACTTCGCGCACGGGGCGTTGGCGACGTGCAAAGGAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATATACCCTGGGc                                                                                                                                    >  6:207905/1‑147 (MQ=255)
                     gcgcCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAg                                                                                                                            <  1:267738/149‑1 (MQ=255)
                       gcCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCt                                                                                                                          <  3:168373/149‑1 (MQ=255)
                               caggcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATGTACCATCGGCCAGGTAGCTCAGCTgcc                                                                                                                  >  6:303376/1‑147 (MQ=255)
                               caggcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCCGGCCCGCGCCCGCa                                                                                                                                                    <  5:38004/116‑1 (MQ=255)
                               caggcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCCGGCCCGCGCCCGCa                                                                                                                                                    >  6:38004/1‑116 (MQ=255)
                                  gcTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCTCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTTACTGc                                                                                                               <  4:244608/149‑1 (MQ=255)
                                        gTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTggc                                                                                                         >  2:244322/1‑149 (MQ=255)
                                        gTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTggc                                                                                                         <  3:250593/149‑1 (MQ=255)
                                         ttttCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTggcg                                                                                                        >  5:140113/1‑149 (MQ=255)
                                                    ggCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAgctg                                                                                               <  1:244322/147‑1 (MQ=255)
                                                     gCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTc                                                                                            >  3:202947/1‑149 (MQ=255)
                                                              gTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAg                                                                                  <  3:196675/149‑1 (MQ=255)
                                                              gTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAg                                                                                  <  2:114031/149‑1 (MQ=255)
                                                                 ccGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTc                                                                               <  3:179288/149‑1 (MQ=255)
                                                                 ccGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTc                                                                               >  3:340236/1‑149 (MQ=255)
                                                                 ccGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTc                                                                               <  4:340236/149‑1 (MQ=255)
                                                                     ggCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGa                                                                           <  4:202947/149‑1 (MQ=255)
                                                                      gCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGAt                                                                          >  2:273046/1‑149 (MQ=255)
                                                                           tgtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCg                                                                     <  6:140113/149‑1 (MQ=255)
                                                                           tgtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCg                                                                     <  5:216545/149‑1 (MQ=255)
                                                                            gtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGa                                                                           >  1:155180/1‑142 (MQ=255)
                                                                            gtCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGa                                                                           <  2:155180/142‑1 (MQ=255)
                                                                             tCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGc                                                                   >  8:138657/1‑149 (MQ=255)
                                                                                   aTCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGAc                                                             <  1:187699/149‑1 (MQ=255)
                                                                                      aCTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACcgc                                                          >  5:197664/1‑149 (MQ=255)
                                                                                                   ggCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAg                                             <  7:245655/149‑1 (MQ=255)
                                                                                                           ccACGTACAAAGTAGGTTTCGCTGCGTGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTc                                                                      <  5:368412/116‑1 (MQ=255)
                                                                                                           ccACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTc                                                                      >  6:368412/1‑116 (MQ=255)
                                                                                                           ccACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCtt                                     <  7:138657/149‑1 (MQ=255)
                                                                                                                tACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTATACTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCa                                 >  1:98363/1‑149 (MQ=255)
                                                                                                                 aCAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCAtt                               >  3:221494/1‑149 (MQ=255)
                                                                                                                    aaGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACACGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGt                             <  1:273046/149‑1 (MQ=255)
                                                                                                                     aGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAg                           >  1:217198/1‑149 (MQ=255)
                                                                                                                       tAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGcc                         >  1:106434/1‑149 (MQ=255)
                                                                                                                         ggTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCgg                       <  4:221494/149‑1 (MQ=255)
                                                                                                                           tttCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGc                      >  8:134645/1‑148 (MQ=255)
                                                                                                                            ttCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGc                      >  1:330833/1‑147 (MQ=255)
                                                                                                                               gCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGca                 <  6:197664/149‑1 (MQ=255)
                                                                                                                               gCTGCCGGCCCGCGCCCGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGca                 >  6:13023/1‑149 (MQ=255)
                                                                                                                                cTGCCGGCCCGCGCCCGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGca                 >  3:350549/1‑148 (MQ=255)
                                                                                                                                         cgggcCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGa                                                              <  6:261508/94‑1 (MQ=255)
                                                                                                                                         cgggcCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGa                                                              >  5:261508/1‑94 (MQ=255)
                                                                                                                                            gcCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGCACAATTGGCAGGcc    >  8:272909/1‑149 (MQ=255)
                                                                                                                                              cGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGCACAATTGGCAGGCCCg  >  7:74822/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
CTTTGAACTGGGGCAGGTTCAGCGCCAGGCCCAGGCTGGCGTTTTCTGGCTGGGCCAGCAGGGTGCCGTGGCGGGTGTCGTACATCACTTCGCGCACGGGGCGTTGGCCACGTACAAAGTAGGTTTCGCTGCGGGCGCGGGCCGGCAGCCATTTACCATCGGCCAGGTAGCTCAGCTGACTGCCTTGGCGGCGCAGCTGCTCGAGGAACAGGTCCTGATTGTCGGCCATGACCGCGCTGCTGCTCCAGGCCAGCTTGCCATTGTAGCCGGCCAGCACAATTGGCAGGCCCG  >  NC_002947/1266778‑1267068

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: