Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,847,969 C→T G14R (GGA→AGA)  selD ← selenide, water dikinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,847,9690CT100.0% 119.1 / NA 36G14R (GGA→AGA) selDselenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/14);  total (22/14)

CGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGA  >  NC_000913/1847727‑1848204
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
cgcgAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGt                                                                                                                                                                                                                                     >  1:46189/1‑251 (MQ=255)
     aGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAAACGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTgg                                                                                                                                                                                                                                >  1:256549/1‑251 (MQ=255)
                                cAATGCGGCAAAAATCGAAAGTAGTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATg                                                                                                                                                                                                     <  1:307077/251‑1 (MQ=255)
                                 caTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGCAGACGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATgg                                                                                                                                                                                                    >  2:401342/2‑251 (MQ=255)
                                                      aTTATCAACGATCGGCATAAAGCCGTCGGTGGTACTGATAACGCTGGTGCCATTGCACAGATCGTACACCGCCTCATCGTCGCGGGGTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGTAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCCATGCAATc                                                                                                                                                                               >  2:210059/1‑251 (MQ=255)
                                                              cgagcggaaTAAAAAAGTCGGTGGTACTGATAACGCGGGTCCCTTTCCCAAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCGCATGGACATCTCCTGGCAATGCAATCCGGGTATg                                                                                                                                                                       <  2:221431/243‑1 (MQ=255)
                                                                 ccGGCATAAAGAAGTCGGTGGTACTGATAACCCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAGCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTa                                                                                                                                                                    <  1:301195/250‑1 (MQ=255)
                                                                   ggCATACAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAACCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGAACTCTCCTGTCACTGCAATCCGGGGTTGGTAAc                                                                                                                                                                  >  2:114997/1‑251 (MQ=255)
                                                                             agTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGACGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCg                                                                                                                                                        >  1:224379/1‑251 (MQ=255)
                                                                                      tACTGATAACGCTGGAGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAACCGAATCGGGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGGGAAGTgct                                                                                                                                               >  2:211530/1‑249 (MQ=255)
                                                                                                   ggTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTc                                                                                                                                  >  1:156246/1‑251 (MQ=255)
                                                                                                                    ggTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCttt                                                                                                                 <  1:299012/250‑1 (MQ=255)
                                                                                                                    cGTACACCGCCGCATCGTCGCTGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCttt                                                                                                                 <  1:30034/251‑1 (MQ=255)
                                                                                                                     gTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGCTGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTg                                                                                                                >  2:443994/1‑251 (MQ=255)
                                                                                                                      acccccgccgCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGGGTCTGAATTCCTGCGCCTTTGc                                                                                                               <  2:115731/247‑1 (MQ=255)
                                                                                                                         acCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCt                                                                                                                   >  1:192516/1‑244 (MQ=255)
                                                                                                                         acCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCt                                                                                                                   <  2:192516/244‑1 (MQ=255)
                                                                                                                                    gtcgCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGAc                                                                                                  >  1:380754/1‑250 (MQ=255)
                                                                                                                                                agtACCCACAAGCAAAATCGGATCAACAAACGCCGCCTGCCCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTGACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTTATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACaa                                                                                     <  2:550174/249‑1 (MQ=255)
                                                                                                                                                aTTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACa                                                                                      >  2:248480/1‑250 (MQ=255)
                                                                                                                                                               aTTCGGAACAACAAACTTCGCCTGCTCCCTCTTCAGGATGGTTTCCAACTCTTTTGGGGAAATTTTACATCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAAGCCAGACAGTTTCGCGACAATTATCAGAAATAAGt                                                                      <  2:35745/251‑1 (MQ=255)
                                                                                                                                                                           caaCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACCTCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGGGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACATTTATCAGCAATAAGTTACAAACGGCGt                                                          >  2:153845/2‑251 (MQ=255)
                                                                                                                                                                            aaCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTc                                                         <  1:206079/251‑1 (MQ=255)
                                                                                                                                                                                tCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCgg                                                       <  1:20937/249‑1 (MQ=255)
                                                                                                                                                                                         cACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGAc                                              <  2:156246/249‑1 (MQ=255)
                                                                                                                                                                                         cACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAAGCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGAc                                              <  2:493491/249‑1 (MQ=255)
                                                                                                                                                                                                 aGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGCTAAGTGAGAGTGTCTg                                                                                                                                >  2:165156/1‑159 (MQ=255)
                                                                                                                                                                                                 aGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTg                                                                                                                                <  1:165156/159‑1 (MQ=255)
                                                                                                                                                                                                   gATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAACCGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCCTTCCGTGAAGTGATAAGTGAGAGCGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGCCGTTAATCGACGt                                  >  2:10785/1‑251 (MQ=255)
                                                                                                                                                                                                      ggTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGt                                  >  2:251773/1‑248 (MQ=255)
                                                                                                                                                                                                      ggTTTCCAAAACTTTTGGGGCAATTTTACAGCCGCCACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACCTCTCCTGTCAATGCAATCAGGGTATGGTAACCCTCATTCAGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCCGACAGTTTCGCGACAATTATCAGAAATAAGTTACACACGGCGTCGGGTCCGGGACGTTAATCGACGt                                  >  2:427603/1‑248 (MQ=255)
                                                                                                                                                                                                                  ttttGGGGAAATTTTACAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGGGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGCCAGTTTCGCGACAATTATCAGAACTAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTg                   >  2:80606/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggAAATTTTCCAGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTgcggcgg             >  2:454255/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cAGCCGAAACCAGCTCTGTGGCTGTATTGGGTCACACGAATCGCGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTACCCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATACGTTACCAACGGCGTCGGGTCCGGGCCGTTAATCGACGTAGATGCTTTCAGCTGCGCCGTACCGAGGt     >  2:234861/1‑249 (MQ=255)
                                                                                                                                                                                                                                   aGCCGCAACCAGCTCTGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGt     >  1:42650/1‑248 (MQ=255)
                                                                                                                                                                                                                                   aGCCGCAACCAGCTCTGTGCCTGACTTGGGTAAAACGAATAGCGGTCTCGCTCATGGACATCTCCTGTCAATGCAATACGGGCATGGTAACACTCATTCCGTGAAGAGATAAGCGAGCGTGTCTGAATTCCTGCGCCATTGCTCACACTCCAGACAGTTTCGCGACAATTATCAGACATAAGTTACAAACGGCGTCGGGCCCGGGGCGATAATCGAAGTCGATCATTTCAGCTGCGGCGAACCGAgggagg  >  2:365106/1‑247 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
CGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCCGTGAAGTGATAAGTGAGAGTGTCTGAATTCCTGCGCCTTTGCTCACAATCCAGACAGTTTCGCGACAATTATCAGAAATAAGTTACAAACGGCGTCGGGTCCGGGACGTTAATCGACGTAGATGCTTTCAGCTGCGGCGTACCGAGGTAGA  >  NC_000913/1847727‑1848204

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: