Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,567,602 A→C 9.3% L109R (CTG→CGG)  PP_2250 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,567,6020AC 9.3% 153.1 / 4.8 53L109R (CTG→CGG) PP_2250LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base C (2/3);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GCGGGCATCGACCTGCAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCTGTGCGG  >  NC_002947/2567457‑2567749
                                                                                                                                                 |                                                                                                                                                   
gcggGCATCGACCTGCAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCcagca                                                                                                                                                  <  5:75665/149‑1 (MQ=255)
       tCGACCTGCAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAgcgggcg                                                                                                                                           <  1:169384/149‑1 (MQ=255)
          aCCTGCAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCA‑GGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGa                                                                                                                                       <  7:66345/149‑1 (MQ=255)
           ccTGCAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGa                                                                                                                                       >  6:304764/1‑149 (MQ=255)
               cAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGATCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATc                                                                                                                                     <  8:67586/147‑1 (MQ=255)
                  gCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATCACCGCCCGCTGCTGGCGGCCATCgccaa                                                                                                                                <  5:19525/149‑5 (MQ=255)
                    cGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCg                                                                                                                              >  6:22741/1‑149 (MQ=255)
                      gCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAgcg                                                                                                                                               <  2:215235/130‑1 (MQ=255)
                      gCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAgcg                                                                                                                                               >  1:215235/1‑130 (MQ=255)
                       cAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGt                                                                                                                           >  3:308682/1‑149 (MQ=255)
                        aaaGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTg                                                                                                                          >  5:105872/1‑149 (MQ=255)
                          aGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATCACCGCCCGCTgctggcggccatc                                                                                                                                     <  8:340577/136‑6 (MQ=255)
                          aGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATCACCGCCCGCTgctggcggccatc                                                                                                                                     >  7:117324/1‑131 (MQ=255)
                          aGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATCACCGCCCGCTgctggcggccatc                                                                                                                                     >  7:340577/1‑131 (MQ=255)
                          aGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATCACCGCCCGCTgctggcggccatc                                                                                                                                     <  8:117324/136‑6 (MQ=255)
                           ggCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCa                                                                                                                       >  4:85279/1‑149 (MQ=255)
                              tACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTg                                                                                                                    >  3:142334/1‑149 (MQ=255)
                                   agagGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAgg                                                                                                               >  8:199382/1‑149 (MQ=255)
                                    gagGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAggg                                                                                                              <  3:56914/149‑1 (MQ=255)
                                        tCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAg                                                                                                          >  2:346974/1‑149 (MQ=255)
                                          ggCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAgg                                                                                                         <  4:28482/148‑1 (MQ=255)
                                               ggCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGctc                                                                                                   >  6:130219/1‑149 (MQ=255)
                                                gCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGctct                                                                                                  <  1:132043/149‑1 (MQ=255)
                                                  cGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCtctc                                                                                                 >  7:49998/1‑148 (MQ=255)
                                                   gCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCa                                                                                               <  7:270973/149‑1 (MQ=255)
                                                   gCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCa                                                                                               <  5:86953/149‑1 (MQ=255)
                                                          aGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTg                                                                                        <  4:76590/149‑1 (MQ=255)
                                                          aGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTg                                                                                        <  4:81651/149‑1 (MQ=255)
                                                           gCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCATGGTCAGGCTGCTCTCCAGCCCTTgg                                                                                       <  4:308682/149‑1 (MQ=255)
                                                            cAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTgg                                                                                       >  3:243965/1‑148 (MQ=255)
                                                                  cAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCt                                                                                <  1:243014/149‑1 (MQ=255)
                                                                   aGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTc                                                                               >  4:236127/1‑149 (MQ=255)
                                                                   aGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTc                                                                               <  7:40327/149‑1 (MQ=255)
                                                                       cttcttCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCaga                                                                           <  4:19542/149‑1 (MQ=255)
                                                                       cttcttCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCaga                                                                           <  5:22741/149‑1 (MQ=255)
                                                                       cttcttCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCaga                                                                           <  6:105872/149‑1 (MQ=255)
                                                                       cttcttCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCa                                                                             >  3:19542/1‑147 (MQ=255)
                                                                             ccTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATg                                                                      <  1:346974/148‑1 (MQ=255)
                                                                              cTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGaa                                                                    <  7:199382/149‑1 (MQ=255)
                                                                                            gCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGccc                                                       >  2:19355/1‑148 (MQ=255)
                                                                                                ggaaaTATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGACAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCAGCAGCAGCCCCAgta                                                  >  6:112883/4‑148 (MQ=255)
                                                                                                    atatCGAGCATTGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGGAGCGGGAGGTGATCGATGTCTGGTGCCACTGCCAGGGTCAGGCTGCTCTCCTGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTc                                              <  5:112883/149‑1 (MQ=255)
                                                                                                          aGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCaa                                         >  5:212332/1‑148 (MQ=255)
                                                                                                          aGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAAt                                        >  5:141843/1‑149 (MQ=255)
                                                                                                                 ggTGGCGCTTACCTAGGTTGGCGATGGCCACCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTg                                 <  7:249181/149‑1 (MQ=255)
                                                                                                                      ctcTTACCTAGGTTGGCGATTGCCGACATCAGCGGGAGGTGATCGATGTCGGTTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTACAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCgtg                            <  4:243965/147‑1 (MQ=255)
                                                                                                                      cgcTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCgtg                            <  6:307963/149‑1 (MQ=255)
                                                                                                                        cTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGtgtg                          <  1:19355/149‑1 (MQ=255)
                                                                                                                           aCCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTgg                         <  4:142334/147‑1 (MQ=255)
                                                                                                                               aGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAg                   >  3:260449/1‑149 (MQ=255)
                                                                                                                                           gcagccAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCt        <  6:212332/145‑1 (MQ=255)
                                                                                                                                            tcgccAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCtgt      <  2:153767/148‑1 (MQ=255)
                                                                                                                                            ccgccAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCtgt      >  8:333927/1‑149 (MQ=255)
                                                                                                                                             cgccAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCtgtg     >  3:92282/1‑149 (MQ=255)
                                                                                                                                                cagcagCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCTGTGCgg  >  1:285200/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGGGCATCGACCTGCAGGCCGGCAAAGGCTACGCAGAGGTCGGCACGGCCGCTGTCCAGCAAGTGCAGGGCTTCTTCCTGCGGGGCGCTGAGCAGGGCAATATCCAGCAATGGGTGGCGCTTACCTAGGTTGGCGATGGCCGCCAGCAGCGGGCGGTGATCGATGTCGGGTGCCACTGCCAGGGTCAGGCTGCTCTCCAGCCCTTGGGACAGCTCCAGAGCATGAACCTGCAGCAGCCCCAGTTGTTCGGCAATCAGCCTGGCGTGTGGCTCCAGGGCCAAGGCCTGTGCGG  >  NC_002947/2567457‑2567749

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