Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 468,510 A→G E43G (GAG→GGG)  ybaO → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913468,5100AG100.0% 212.7 / NA 59E43G (GAG→GGG) ybaOputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/27);  total (32/27)

GAAGGGCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACGGT  >  NC_000913/468375‑468639
                                                                                                                                       |                                                                                                                                 
gAAGGGCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                <  2:1575851/139‑1 (MQ=255)
   gggCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacga                                                                                                                             <  1:755488/139‑1 (MQ=255)
      cTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACgg                                                                                                                          >  1:764286/1‑139 (MQ=255)
        aTGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTa                                                                                                                        >  1:1044670/1‑139 (MQ=255)
        aTGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTa                                                                                                                        <  2:293003/139‑1 (MQ=255)
           ttAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATcc                                                                                                                     >  2:1275517/1‑139 (MQ=255)
               aTAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTAt                                                                                                                 >  2:1413578/1‑139 (MQ=255)
               aTAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTAt                                                                                                                 >  2:1253361/1‑139 (MQ=255)
                tAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATc                                                                                                                >  2:656237/1‑139 (MQ=255)
                  aaaTTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCgg                                                                                                              <  2:1462666/139‑1 (MQ=255)
                         ccGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTc                                                                                                       >  2:1580104/1‑139 (MQ=255)
                          cGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCg                                                                                                      <  2:764286/139‑1 (MQ=255)
                           gTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGc                                                                                                     <  2:1255591/139‑1 (MQ=255)
                            tAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGcc                                                                                                    <  1:1413578/139‑1 (MQ=255)
                                 tgctgGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCctgct                                                                                               <  1:1580104/139‑1 (MQ=255)
                                    tgGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacg                                                                                                                              >  1:224237/1‑105 (MQ=255)
                                    tgGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacg                                                                                                                              <  2:224237/105‑1 (MQ=255)
                                          tACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGa                                                                                      >  2:531204/1‑139 (MQ=255)
                                           aCTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                <  2:486973/96‑1 (MQ=255)
                                           aCTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                >  1:486973/1‑96 (MQ=255)
                                               cagcagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGaaaaaa                                                                                 >  2:1218879/1‑139 (MQ=255)
                                                  cagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacga                                                                                                                             <  2:305629/92‑1 (MQ=255)
                                                  cagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacga                                                                                                                             >  1:305629/1‑92 (MQ=255)
                                                  cagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCgg                                                                                       <  1:550679/130‑1 (MQ=255)
                                                  cagGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCgg                                                                                       >  2:550679/1‑130 (MQ=255)
                                                          cACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacga                                                                                                                             <  2:153601/84‑1 (MQ=255)
                                                          cACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGgacga                                                                                                                             >  1:153601/1‑84 (MQ=255)
                                                                ctctTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAcc                                                                >  2:521928/1‑139 (MQ=255)
                                                                      gCAGGCACTGGTTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCtttt                                                          <  2:1044670/139‑1 (MQ=255)
                                                                         ggCACTGGCTGAAGCCGTTAATCTGACATCCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTtgtg                                                       <  2:964826/139‑1 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                >  2:536090/1‑64 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                <  1:536090/64‑1 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                >  2:468500/1‑64 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                <  1:468500/64‑1 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGcc                                                                           <  1:88631/117‑1 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGcc                                                                           >  2:88631/1‑117 (MQ=255)
                                                                           cACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCt                                                     >  2:863342/1‑139 (MQ=255)
                                                                             cTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGa                                                   >  2:479858/1‑139 (MQ=255)
                                                                               ggCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATa                                                 >  2:1444121/1‑139 (MQ=255)
                                                                                 cTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATaaa                                               <  2:400220/139‑1 (MQ=255)
                                                                                   gAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATaaaaa                                             <  1:479858/139‑1 (MQ=255)
                                                                                     aGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAAt                                                                                <  1:1388444/102‑1 (MQ=255)
                                                                                     aGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAAt                                                                                >  2:1388444/1‑102 (MQ=255)
                                                                                      gCCGTTAATCTGACAACCACCCCTTTCTTGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGTCCTGACTGCTTTTGTGCTGATAAAAACGc                                          <  1:119601/139‑1 (MQ=255)
                                                                                      gCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGc                                          >  2:9997/1‑139 (MQ=255)
                                                                                          ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACa                                      >  2:1254215/1‑139 (MQ=255)
                                                                                          ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACa                                      >  2:1015412/1‑139 (MQ=255)
                                                                                          ttAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACa                                      >  1:1260586/1‑139 (MQ=255)
                                                                                               cTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATcaca                                 <  1:531204/139‑1 (MQ=255)
                                                                                                tGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCAcag                                <  1:863342/139‑1 (MQ=255)
                                                                                                 gACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCAcagc                               <  1:521928/139‑1 (MQ=255)
                                                                                                  aCAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCAcagca                              <  1:1015412/139‑1 (MQ=255)
                                                                                                       caccCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                <  2:1373793/36‑1 (MQ=255)
                                                                                                       caccCCTTGCTGGAAGCGCCTGAAACGGCTggggga                                                                                                                                >  1:1373793/1‑36 (MQ=255)
                                                                                                       caccCCTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGaa                         <  2:1260586/139‑1 (MQ=255)
                                                                                                         ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGcc                                                                     <  2:1355034/93‑1 (MQ=255)
                                                                                                         ccccTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGcc                                                                     >  1:1355034/1‑93 (MQ=255)
                                                                                                            cTTGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTa                    >  1:1543391/1‑139 (MQ=255)
                                                                                                              tGCTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTAtt                  >  1:1254886/1‑139 (MQ=255)
                                                                                                                cTGGAAGCGCCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTGTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGc                >  1:1045552/1‑139 (MQ=255)
                                                                                                                       cgcCTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTg         >  2:1504405/1‑139 (MQ=255)
                                                                                                                         ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAc                                                                 >  1:1348486/1‑81 (MQ=255)
                                                                                                                         ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGAc                                                                 <  2:1348486/81‑1 (MQ=255)
                                                                                                                         ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGc                                          <  1:1458459/104‑1 (MQ=255)
                                                                                                                         ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGc                                          >  2:1458459/1‑104 (MQ=255)
                                                                                                                         ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTc       >  1:402977/1‑139 (MQ=255)
                                                                                                                         ccTGAAACGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTc       >  2:904107/1‑139 (MQ=255)
                                                                                                                              aaCGGCTGGGGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACggt  <  1:1504405/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
GAAGGGCTATGTTAGATAAAATTGACCGTAAGCTGCTGGCCTTACTGCAGCAGGATTGCACCCTCTCTTTGCAGGCACTGGCTGAAGCCGTTAATCTGACAACCACCCCTTGCTGGAAGCGCCTGAAACGGCTGGAGGACGACGGTATCCTTATCGGCAAAGTCGCCCTGCTGGATCCGGAAAAAATAGGCCTCGGCCTGACCGCTTTTGTGCTGATAAAAACGCAACATCACAGCAGCGAATGGTATTGCCGCTTTGTCACGGT  >  NC_000913/468375‑468639

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