Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 584,722 C→T 8.1% L178L (CTG→CTA ansA ← type 1 L‑asparaginase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947584,7220CT 8.1% 141.4 / 2.5 49L178L (CTG→CTAansAtype 1 L‑asparaginase
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base T (2/2);  total (24/25)
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.37e-01

CACACCGCTGTCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCAATGCGCCGCACAGGTTGTCCCAGGCATCG  >  NC_002947/584575‑584869
                                                                                                                                                     |                                                                                                                                                   
cacaCCGCTGTCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGgca                                                                                                                                                    <  7:119680/149‑1 (MQ=255)
   acCGCTGTCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGgcagcg                                                                                                                                                 >  1:27637/1‑149 (MQ=255)
   acCGCTGTCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGgcagcg                                                                                                                                                 <  2:37929/149‑1 (MQ=255)
    ccGCTGTCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGgcagcgc                                                                                                                                                >  7:97299/1‑149 (MQ=255)
          tCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCgcagca                                                                                                                                            >  6:68348/1‑147 (MQ=255)
             aTCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAgg                                                                                                                                        >  4:55084/1‑148 (MQ=255)
              tCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCa                                                                                                                                      <  2:12568/149‑1 (MQ=255)
                aaGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGACTACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATc                                                                                                                                    <  7:196383/149‑1 (MQ=255)
                aaGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATc                                                                                                                                    <  6:106521/149‑1 (MQ=255)
                      tGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGc                                                                                                                              >  5:41966/1‑149 (MQ=255)
                         aaaTGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTc                                                                                                                           <  6:251416/149‑1 (MQ=255)
                           aTGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGc                                                                                                                         >  3:60445/1‑149 (MQ=255)
                               ccGGCCTGAAGCCCCGGGAATACTGGCAGAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCg                                                                                                                     >  1:152818/1‑149 (MQ=255)
                                 ggCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGgcagc                                                                                                                                                  >  1:121025/1‑118 (MQ=255)
                                 ggCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGgcagc                                                                                                                                                  <  2:121025/118‑1 (MQ=255)
                                       aaGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTc                                                                                                              <  4:60445/148‑1 (MQ=255)
                                       aaGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTc                                                                                                              <  3:98969/148‑1 (MQ=255)
                                             cGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGccc                                                                                                       <  8:762/149‑1 (MQ=255)
                                                     cTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGcc                                                                                                <  5:79058/148‑1 (MQ=255)
                                                      tGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCAt                                                                                              >  5:94466/1‑149 (MQ=255)
                                                        gCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGcc                                                                                                >  6:184529/1‑145 (MQ=255)
                                                        gCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGcc                                                                                                <  5:184529/145‑1 (MQ=255)
                                                        gCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATg                                                                                             >  2:19618/1‑148 (MQ=255)
                                                          acaACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGCCGCCAGACAGCCATgcag                                                                                          >  1:252546/1‑149 (MQ=255)
                                                          acaACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATgcag                                                                                          <  6:104279/149‑1 (MQ=255)
                                                          acaACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATgcag                                                                                          <  3:55084/149‑1 (MQ=255)
                                                            aaCCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATgcag                                                                                          >  1:195186/1‑148 (MQ=255)
                                                            aaCCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATgcag                                                                                          >  6:161189/1‑147 (MQ=255)
                                                             aCCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGcagcag                                                                                       <  8:97299/149‑1 (MQ=255)
                                                              ccGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCagcagc                                                                                      <  2:152818/149‑1 (MQ=255)
                                                                   aGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGccc                                                                                 >  6:228085/1‑149 (MQ=255)
                                                                      ttGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTg                                                                              <  2:195186/149‑1 (MQ=255)
                                                                         aCCGGCTTGCGCGGCTGCTTGTAGCCCAGTTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATgc                                                                                            >  1:24869/1‑133 (MQ=255)
                                                                         aCCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATgc                                                                                            <  2:24869/132‑1 (MQ=255)
                                                                         aCCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGaa                                                                           <  1:19618/149‑1 (MQ=255)
                                                                         aCCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGaa                                                                           <  6:212451/149‑1 (MQ=255)
                                                                            ggCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTa                                                                        >  7:74957/1‑149 (MQ=255)
                                                                             gCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTAc                                                                       <  2:27637/149‑1 (MQ=255)
                                                                               ttGCGCTGCTTCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGTAGCGCAGCAAAGGCATCGAACGCTTCg                                                                                                                          <  2:74803/96‑1 (MQ=255)
                                                                               ttGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGTAGCGCAGCAAAGGCATCGAACGCTTCg                                                                                                                          >  3:141787/1‑96 (MQ=255)
                                                                               ttGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGTAGCGCAGCAAAGGCATCGAACGCTTCg                                                                                                                          <  4:141787/96‑1 (MQ=255)
                                                                               ttGAGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGTAGCGCAGCAAAGGCATCGAACGCTTCg                                                                                                                          >  1:74803/1‑96 (MQ=255)
                                                                                 gcgcGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAg                                                                     <  8:119422/147‑1 (MQ=255)
                                                                                              tAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTc                                                      >  2:36143/1‑149 (MQ=255)
                                                                                                   cAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTgg                                                                             >  3:261356/1‑121 (MQ=255)
                                                                                                   cAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTgg                                                                             <  4:261356/121‑1 (MQ=255)
                                                                                                    aGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGcc                                                >  3:37085/1‑149 (MQ=255)
                                                                                                             cAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGTGCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCg                                         <  1:81153/149‑1 (MQ=255)
                                                                                                              aGGGATCGTATCGGCCCGCTCAGCGTCATGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCACGCAGCAGCTGCCCGTgg                                                                             >  1:153115/1‑112 (MQ=255)
                                                                                                              aGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCACGCAGCAGCTGCCCGTgg                                                                             <  2:153115/112‑1 (MQ=255)
                                                                                                                ggATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGcc                                                                                  >  3:129025/1‑105 (MQ=255)
                                                                                                                ggATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGcc                                                                                  <  4:129025/105‑1 (MQ=255)
                                                                                                                                       tCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCAATGCGCCGCACAGGtt               >  1:79801/1‑149 (MQ=255)
                                                                                                                                          cGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCAATGCGCCGCACAGGTTGTc            >  7:176094/1‑149 (MQ=255)
                                                                                                                                          cGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCAATGCGCCGCACAGGTTGTc            >  2:11239/1‑149 (MQ=255)
                                                                                                                                                 gtggcAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCaa                              >  6:116050/1‑124 (MQ=255)
                                                                                                                                                 gtggcAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCaa                              <  5:116050/124‑1 (MQ=255)
                                                                                                                                                    gcagcgcagcAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCAATGCGCCGCACAGGTTGTCCCAGGCATCg  <  1:11239/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACACCGCTGTCTATCAAGGCCTGCAAATGCCCGGCCTGAAGCCCCGGGAATACTGGCACAACCGCCAGGTTGACCGGCTTGCGCGGCTGCTTGTAGCCCAGTTCGGCAGGGATCGTATCGGCCCGCTCAGCGTCACGGTGGCGTGGCAGCGCAGCAAAGGCATCGAACGCTTCGCTGCGTAGCTTCGACGCCCGACAGCCATGCAGCAGCTGCCCGTGGAAGTACAGCTGCACGCCATCTTCCAGGCCGTGTTCGAACTGGCGCAATGCGCCGCACAGGTTGTCCCAGGCATCG  >  NC_002947/584575‑584869

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: