Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,186,053 T→G 10.2% E136A (GAG→GCG)  aceK ← isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,186,0530TG10.2% 149.9 / 6.4 56E136A (GAG→GCG) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base G (4/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGAT  >  NC_002947/5185906‑5186186
                                                                                                                                                    |                                                                                                                                      
aggcaTTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGCCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTc                                                                                                                                       >  5:214631/1‑149 (MQ=255)
    aTTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTAAGGGTTTGCGCGGCGAGCTcgctg                                                                                                                                   >  8:125882/1‑147 (MQ=255)
        gcgcAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCcgc                                                                                                                              <  6:247006/149‑1 (MQ=255)
          gcAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCCTGTCGAACGGGTACTCGGCGAACACCGCGCGCTGCAAGTTCCTGAGACTGCCGGCAAGCCGATAGGTAGGGGTTTGCGCGGCGCGCGCGCGGCcgcgg                                                                                                                            >  4:266759/1‑148 (MQ=255)
                gttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATgg                                                                                                                        >  3:314622/1‑147 (MQ=255)
                gttCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATgg                                                                                                                        >  7:172791/1‑147 (MQ=255)
                   cttcCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGAGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGc                                                                                                                   >  5:215788/1‑149 (MQ=255)
                       cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                                <  3:241600/148‑1 (MQ=255)
                       cAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCgg                                                                                                                <  7:182007/148‑1 (MQ=255)
                        aGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGTGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCCGCGAGCGCGCCGCGCa                                                                                                                           >  3:295597/1‑136 (MQ=255)
                        aGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGTGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCCGCGAGCGCGCCGCGCa                                                                                                                           <  4:295597/136‑1 (MQ=255)
                        aGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCggg                                                                                                               >  6:158655/1‑148 (MQ=255)
                        aGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCggg                                                                                                               >  3:304534/1‑148 (MQ=255)
                           cgcgCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAgt                                                                                                           <  4:314622/149‑1 (MQ=255)
                              gcCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGtgtg                                                                                                        <  6:117740/149‑1 (MQ=255)
                               ccAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGtgtg                                                                                                        >  4:126542/1‑148 (MQ=255)
                               ccAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTgg                                                                                                       <  3:292791/149‑1 (MQ=255)
                                cAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGa                                                                                                      <  3:226079/149‑1 (MQ=255)
                                 aGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGAt                                                                                                     >  4:132965/1‑149 (MQ=255)
                                       ccTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAAc                                                                                                 >  5:142598/1‑147 (MQ=255)
                                       ccTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAAc                                                                                                 <  6:142598/147‑1 (MQ=255)
                                       ccTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACAt                                                                                               <  4:304534/149‑1 (MQ=255)
                                       ccTTCCAGGGCGCCATAGGGCATGTCGAAGGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACTGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACAt                                                                                               <  8:235112/149‑1 (MQ=255)
                                         ttCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCCGCGCGCGAGCGGCCGCGCATGGAAGGCCCGGGAGTGTGGATGAACATg                                                                                              >  2:310682/1‑148 (MQ=255)
                                           ccAGGTCGCCATAGCGCATGGCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTATGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCa                                                                                            <  8:220672/148‑1 (MQ=255)
                                              ggTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCg                                                                                        >  2:184950/1‑149 (MQ=255)
                                                  gCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCg                                                                                     <  3:126542/148‑1 (MQ=255)
                                                  gCCATAGGGCATGTCGAACGGGTAAGCGTCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGGCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGATGGCCGGGTAGTGTGGATGAACATGCAGCCGTCg                                                                                     <  7:125882/148‑1 (MQ=255)
                                                    cATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTgg                                                                                                       <  8:84662/128‑1 (MQ=255)
                                                    cATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTgg                                                                                                       >  7:84662/1‑128 (MQ=255)
                                                      tAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGAt                                                                                >  7:176232/1‑149 (MQ=255)
                                                      tAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGAt                                                                                <  6:215788/149‑1 (MQ=255)
                                                         ggCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAg                                                                             <  5:101944/149‑1 (MQ=255)
                                                         ggCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAg                                                                             <  8:176232/149‑1 (MQ=255)
                                                            aTGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTc                                                                          <  1:184950/149‑1 (MQ=255)
                                                                cGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATgg                                                                      <  8:172791/149‑1 (MQ=255)
                                                                  aaCGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCt                                                                    >  1:12119/1‑149 (MQ=255)
                                                                    cGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTg                                                                   >  5:108263/1‑148 (MQ=255)
                                                                          aTCTGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAgg                                                             <  1:310682/148‑1 (MQ=255)
                                                                                gAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaa                                                       >  8:97054/1‑148 (MQ=255)
                                                                                  acacCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCaaaa                                                     >  4:98566/1‑148 (MQ=255)
                                                                                  acacCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAAc                                                    <  3:179064/149‑1 (MQ=255)
                                                                                       gcgcgcAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTTCGCGTCGCGCTCGCGGCCGCTCATGGAAGTCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGa                                               <  5:158655/149‑1 (MQ=255)
                                                                                           gcagcaAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGtcg                                           <  7:97054/149‑3 (MQ=255)
                                                                                            cagcaAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTg                                           >  3:234781/1‑148 (MQ=255)
                                                                                            cagcaAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTg                                           >  3:291674/1‑148 (MQ=255)
                                                                                              gcaAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                        <  3:132965/149‑1 (MQ=255)
                                                                                              gcaAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTAc                                        <  6:214631/149‑1 (MQ=255)
                                                                                                       tGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTctct                               >  1:68936/1‑149 (MQ=255)
                                                                                                         agaCTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTctct                               <  2:12119/147‑1 (MQ=255)
                                                                                                                  tCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACtcgtc                    <  3:258227/149‑1 (MQ=255)
                                                                                                                         gATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAgg             <  4:291674/149‑1 (MQ=255)
                                                                                                                                  gggTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTc     >  1:48904/1‑148 (MQ=255)
                                                                                                                                     tttGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGAt  >  5:114721/1‑148 (MQ=255)
                                                                                                                                       tGCGCGGCCGCGAGCGCGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCg           <  6:142572/137‑1 (MQ=255)
                                                                                                                                       tGCGCGGCCGCGAGCGCGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCg           >  5:142572/1‑137 (MQ=255)
                                                                                                                                                    |                                                                                                                                      
AGGCATTCGCGCAGTTGTTCTTCCAGCCGCGCCAGGTCGCCTTCCAGGTCGCCATAGGGCATGTCGAACGGGTAATCGGCGAACACCGCGCGCAGCAAGTTCCTGAGACTGCCGTCAAGCCGATAGGTACGGGTTTGCGCGGCGCGCTCGCGGCCGCGCATGGAAGGCCGGGTAGTGTGGATGAACATGCAGCCGTCGCTGATCAGGTCATGGCTGAAGAGGCTGCAAAACAGCGAGTTGTACCAGGTCTCTGCCAACTCGTCGTCCAGGCGGGGGTCGAT  >  NC_002947/5185906‑5186186

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: