Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 826,693 G→C 8.6% intergenic (‑95/‑117) PP_0710 ← / → ycaC‑I hypothetical protein/putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947826,6930GC 8.6% 169.6 / 7.6 58intergenic (‑95/‑117)PP_0710/ycaC‑Ihypothetical protein/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base C (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

CGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGGACAAGA  >  NC_002947/826548‑826840
                                                                                                                                                 |                                                                                                                                                   
cgccgcGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGaaa                                                                                                                                                  >  4:198957/1‑149 (MQ=255)
   cgcgACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAAttt                                                                                                                                               <  5:157206/149‑1 (MQ=255)
     cgacgaGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTtgt                                                                                                                                             <  5:239222/149‑1 (MQ=255)
          aGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTAGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACt                                                                                                                                        >  1:55193/1‑149 (MQ=255)
           ggTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCATGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTg                                                                                                                                       <  1:159792/149‑1 (MQ=255)
           ggTCGGCGGCGGCTTGCCGCACTTCAAGTGTCAGGACCACCGGAGTGCACTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCCAGGGCAAAAATCGCTGGCGGCAAATAATCAGTAAAGAAATTTGTCAACt                                                                                                                                        >  8:155528/1‑148 (MQ=255)
            gTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGa                                                                                                                                      >  8:382673/1‑149 (MQ=255)
            gTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCACAAATTTCTGAACTGa                                                                                                                                      <  4:177850/149‑1 (MQ=255)
                gcggcggcTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATggg                                                                                                                                  <  2:207182/149‑1 (MQ=255)
                gcggcggcTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATggg                                                                                                                                  <  6:109213/149‑1 (MQ=255)
                   gcggcTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCACAAATTTCTGAACTGAttatttg                                                                                                                               <  6:383906/149‑7 (MQ=255)
                    cggcTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGcac                                                                                                                              <  7:357506/149‑1 (MQ=255)
                        ttGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAggg                                                                                                                          >  7:212395/1‑149 (MQ=255)
                          gCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAggg                                                                                                                          <  7:382673/147‑1 (MQ=255)
                               aCTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTa                                                                                                                     >  8:78396/1‑147 (MQ=255)
                               aCTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGGGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATGTGTGAACTGATGGGGCACAGGGCAGTACa                                                                                                                   <  6:227520/149‑1 (MQ=255)
                                   gAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTtgt                                                                                                                                             <  8:376738/119‑1 (MQ=255)
                                   gAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTtgt                                                                                                                                             >  7:376738/1‑119 (MQ=255)
                                   gAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATg                                                                                                                                    >  6:368302/1‑128 (MQ=255)
                                   gAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATg                                                                                                                                    <  5:368302/128‑1 (MQ=255)
                                         tCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGcac                                                                                                                              >  5:113957/1‑128 (MQ=255)
                                         tCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGcac                                                                                                                              <  6:113957/128‑1 (MQ=255)
                                         tCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAgg                                                                                                                           >  8:63011/1‑131 (MQ=255)
                                             gtgcACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGAc                                                                                                        <  7:155528/143‑1 (MQ=255)
                                             gTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGATAAATCGCTGGCGGCAAATAATCAGTTCAGAAAATTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGa                                                                                                     >  2:82445/1‑149 (MQ=255)
                                             gTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGa                                                                                                     <  3:198957/149‑1 (MQ=255)
                                                cACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACgg                                                                                                  >  8:98703/1‑149 (MQ=255)
                                                 aCGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGAc                                                                                                        <  8:253826/142‑1 (MQ=255)
                                                 aCGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGAc                                                                                                        >  7:253826/1‑142 (MQ=255)
                                                   gggAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGtct                                                                                               <  8:212395/149‑1 (MQ=255)
                                                               aCGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTc                                                                                    <  2:55193/148‑1 (MQ=255)
                                                                  gTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATc                                                                                <  7:98703/149‑1 (MQ=255)
                                                                  gTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATc                                                                                >  7:73634/1‑149 (MQ=255)
                                                                        ttGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCACAAATTTCTGAACTGAt                                                                                                                                     >  4:98769/1‑90 (MQ=255)
                                                                        ttGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCACAAATTTCTGAACTGAt                                                                                                                                     <  3:98769/90‑1 (MQ=255)
                                                                              ttGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCt                                                                    >  5:324057/1‑149 (MQ=255)
                                                                                 tggtgAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCt                                                                    >  5:148733/1‑146 (MQ=255)
                                                                                 tggtgAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCt                                                                    <  6:148733/146‑1 (MQ=255)
                                                                                           aTTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAAcc                                                       <  7:78396/149‑1 (MQ=255)
                                                                                                   ggCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAg                                               >  2:288819/1‑149 (MQ=255)
                                                                                                       gTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCTACCCGCTGAAGCAAg                                           <  6:324057/149‑1 (MQ=255)
                                                                                                       gTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAg                                           <  1:288819/149‑1 (MQ=255)
                                                                                                          ggCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAg                                           >  3:204875/1‑146 (MQ=255)
                                                                                                                   aaaaaTCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACtgt                               >  2:294016/1‑149 (MQ=255)
                                                                                                                    aaaaTCGCTGGCGGCAAATAATCAGTTCAGAAAATTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACtgtg                              <  1:82445/149‑1 (MQ=255)
                                                                                                                          gCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTc                        <  1:294016/149‑1 (MQ=255)
                                                                                                                           cTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTcc                       <  4:204875/149‑1 (MQ=255)
                                                                                                                             ggcggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCaa                     >  7:299672/1‑149 (MQ=255)
                                                                                                                              gcggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCaa                     >  3:195403/1‑148 (MQ=255)
                                                                                                                                ggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGc                                                                                                            <  2:27506/59‑1 (MQ=255)
                                                                                                                                ggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGc                                                                                                            >  1:27506/1‑59 (MQ=255)
                                                                                                                                ggcAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATa                   >  8:149786/1‑148 (MQ=255)
                                                                                                                                    aataatCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGc              <  3:374163/149‑1 (MQ=255)
                                                                                                                                     ataatCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTc                                                             >  6:26412/1‑101 (MQ=255)
                                                                                                                                     ataatCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTc                                                             <  5:26412/101‑1 (MQ=255)
                                                                                                                                      ccatCAGTTCACAAATTTCTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGt            >  1:135286/3‑149 (MQ=255)
                                                                                                                                                aGAAATTTGTTAACTGATGGTGCACAGGGCAGTACAGGCAAGCCTACTGAAGTGCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCACCCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGGACAAGa  <  8:369750/149‑1 (MQ=255)
                                                                                                                                                aGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGGACAAg   >  5:357269/1‑148 (MQ=255)
                                                                                                                                                 gAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGGACAAGa  >  8:377460/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGCCGCGACGAGGTCGGCGGCGGCTTGCCGCACTTGAAGTGTCTGGTCCACGGGAGTGGCCTCACGGTTTTTTTGTGATTGTGGTGAAGCGATTCTGGTGGCGGTGGGCGAGGGGAAAAATCGCTGGCGGCAAATAATCAGTTCAGAAATTTGTGAACTGATGGGGCACAGGGCAGTACAGGCAAGCCGACTGACGGTCTCTCCCACCCTCCATCCTTGCCTACGCTGATACTCCCAACCCGCTGAAGCAAGGAGATCCACTGTGGCCTTCCAATACAAGCGTCTGGACAAGA  >  NC_002947/826548‑826840

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: