Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,316,912 G→C 7.6% intergenic (‑115/+58) ribC ← / ← treSA riboflavin synthase/trehalose synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,316,9120GC 7.6% 151.4 / 4.4 53intergenic (‑115/+58)ribC/treSAriboflavin synthase/trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base C (2/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AGGGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAA  >  NC_002947/3316765‑3317057
                                                                                                                                                   |                                                                                                                                                 
agcgCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCgg                                                                                                                                                  <  3:214708/146‑1 (MQ=255)
  ggCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCgggg                                                                                                                                                >  3:42868/1‑149 (MQ=255)
     ggACGGCCT‑GACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGc                                                                                                                                               >  8:177265/1‑146 (MQ=255)
     ggACGGCCT‑GACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGc                                                                                                                                               <  7:177265/146‑1 (MQ=255)
      gACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAgg                                                                                                                                            <  7:212680/149‑1 (MQ=255)
         ggCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCCCCGCGCAggcgc                                                                                                                                         <  4:24591/149‑2 (MQ=255)
           ccTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGgcggcg                                                                                                                                       >  1:342670/1‑149 (MQ=255)
               gACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCt                                                                                                                                   >  1:149082/1‑149 (MQ=255)
                    tGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTgg                                                                                                                              <  1:107833/149‑1 (MQ=255)
                       cGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCAATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGCCGGCGTGCTTGTGGGc                                                                                                                            >  8:274017/1‑148 (MQ=255)
                          tgtACATGGAGTTGCTCCTGGGCACGGAAAAGAACACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGaaa                                                                                                                        <  3:281109/149‑1 (MQ=255)
                               aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGc                                                                                                                            >  2:236793/1‑140 (MQ=255)
                               aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGc                                                                                                                            <  1:236793/140‑1 (MQ=255)
                               aTGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGt                                                                                                                   >  5:149005/1‑149 (MQ=255)
                                        cTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACgg                                                                                                          <  4:42868/149‑1 (MQ=255)
                                         tCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACgg                                                                                                          >  6:208696/1‑148 (MQ=255)
                                              ggCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTCCGACGGCATGGc                                                                                                    >  3:218322/1‑149 (MQ=255)
                                                          aaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCGCGCCAGAtt                                                                                         >  2:190603/1‑148 (MQ=255)
                                                          aaaaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGAtt                                                                                         >  6:156903/1‑148 (MQ=255)
                                                            aaCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATCCa                                                                                       >  5:6252/1‑148 (MQ=255)
                                                             aCTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                                     <  6:274427/149‑1 (MQ=255)
                                                               ttCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACa                                                                                   >  2:120511/1‑149 (MQ=255)
                                                               ttCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACa                                                                                   >  5:141726/1‑149 (MQ=255)
                                                                  cGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGc                                                                                <  6:149005/149‑1 (MQ=255)
                                                                             cccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGtt                                                                     >  2:278066/1‑149 (MQ=255)
                                                                              ccGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTg                                                                    <  2:149082/149‑1 (MQ=255)
                                                                               cGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGa                                                                   <  1:278066/149‑1 (MQ=255)
                                                                                   gcgcAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGAGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGa                                                                <  5:208696/148‑1 (MQ=255)
                                                                                      gggatgttGCAAGTGCTGGCTTCGTGGCTGGCAAGCCGGCTTGGGTGTGTGGCGGCTGCGGGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATc                                                              <  3:204659/140‑1 (MQ=255)
                                                                                        gTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCa                                                          <  6:141726/149‑1 (MQ=255)
                                                                                              gCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGca                                                      >  3:74636/1‑147 (MQ=255)
                                                                                               cAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                                     >  8:159475/1‑115 (MQ=255)
                                                                                               cAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCaaa                                                                                     <  7:159475/115‑1 (MQ=255)
                                                                                                   aCCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGcc                                               >  1:203386/1‑149 (MQ=255)
                                                                                                     cAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTc                                             >  6:51631/1‑149 (MQ=255)
                                                                                                            ccTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGGCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACgg                                       <  7:274017/148‑1 (MQ=255)
                                                                                                                  tcgcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGAGCTTGTGCGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAgg                                >  5:319365/3‑149 (MQ=255)
                                                                                                                    gcAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGtt                              <  2:342670/149‑1 (MQ=255)
                                                                                                                      aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCa                            <  5:156903/149‑1 (MQ=255)
                                                                                                                      aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCa                            >  4:95918/1‑149 (MQ=255)
                                                                                                                      aaGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTATTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCa                            >  7:329791/1‑149 (MQ=255)
                                                                                                                       aGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                           >  7:177224/1‑149 (MQ=255)
                                                                                                                        gCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCg                                            <  5:367669/131‑1 (MQ=255)
                                                                                                                        gCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCg                                            >  6:367669/1‑131 (MQ=255)
                                                                                                                        gCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCAc                           >  1:210674/1‑148 (MQ=255)
                                                                                                                           ggCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCAt                        >  3:133253/1‑148 (MQ=255)
                                                                                                                             cTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAg                     >  1:35208/1‑149 (MQ=255)
                                                                                                                                  gAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGcc                >  7:374646/1‑149 (MQ=255)
                                                                                                                                      acgccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTc            >  6:132089/5‑149 (MQ=255)
                                                                                                                                        gccgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCt          <  5:132089/147‑1 (MQ=255)
                                                                                                                                         gcgcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCt          <  5:51631/148‑1 (MQ=255)
                                                                                                                                          cgcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCg        <  2:203386/149‑1 (MQ=255)
                                                                                                                                          cgcCTGCCCCGCGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGCCGGCATCGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCg        >  7:108717/1‑149 (MQ=255)
                                                                                                                                           gcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAg                                                         <  6:271691/99‑1 (MQ=255)
                                                                                                                                           gcCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAg                                                         >  5:271691/1‑99 (MQ=255)
                                                                                                                                                cgcgGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATCCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCaaa  <  6:6252/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGGCGGACGGCCTGGACGATGCCGGTGTACATGGAGTTGCTCCTGGGCACGGAAAAGAAAACTTCCGAACGCTATGCCCGTAGCGCAGTAAATGCAAGACCAGGCTTCCTGGCTGGCAAGCCGGCTTGCGAGTCTGGCGCCTGCGCGGGGCAGGCGGCGTGCTTGTGGGCGAAACCGGTTTGCGACGGCATGGCCACGCCAGATTCAAACATGCCCGCTGCTGTTGACGATCCGCAGGCACAGCGCCTCGTACGGGTCCAGGTTCACCATCAGGCGCCCGTCCTCGGTCAAA  >  NC_002947/3316765‑3317057

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