Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,558,932 2 bp→GC coding (153‑154/831 nt) glpG ← GlpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,558,9320CG100.0% 336.8 / NA 127V52L (GTG→CTG) glpGGlpG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (63/64);  total (63/64)
*NC_0009133,558,9330GC100.0% 347.3 / NA 130D51E (GAC→GAGglpGGlpG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (65/65);  total (65/65)

Only 100 of 130 total aligned reads displayed.

CGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCACGTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAAT  >  NC_000913/3558863‑3559001
                                                                     ||                                                                    
 ggTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGcc                                                                >  2:1317866/1‑76 (MQ=255)
  gTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGccc                                                               <  1:3334724/76‑1 (MQ=255)
  gTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGccc                                                               >  2:3466391/1‑76 (MQ=255)
    tgggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGt                                                             <  2:3739715/74‑1 (MQ=255)
    aTGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGt                                                             >  1:1157528/1‑76 (MQ=255)
    aTGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGt                                                             <  1:3276468/76‑1 (MQ=255)
    aTGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGt                                                             >  1:3099547/1‑76 (MQ=255)
     tGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTa                                                            >  1:1830149/1‑76 (MQ=255)
      ggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAc                                                           <  1:4945198/76‑1 (MQ=255)
      ggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAc                                                           >  2:4228871/1‑76 (MQ=255)
      ggCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTAc                                                           <  1:270955/76‑1 (MQ=255)
          tgcctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACgcgc                                                       <  1:5578453/76‑1 (MQ=255)
          tgcctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGATTTCGAGAAAACGCGCCAGCTCCGCCCGTACgcgc                                                       >  1:1882057/1‑76 (MQ=255)
           gcctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCt                                                      <  2:3426422/76‑1 (MQ=255)
            cctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                     <  1:4842501/76‑1 (MQ=255)
            cctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                     <  1:4307776/76‑1 (MQ=255)
            cctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                     <  1:4307274/76‑1 (MQ=255)
            cctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                     <  1:5704783/76‑1 (MQ=255)
             ctgccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCg                                                    >  1:4464009/1‑76 (MQ=255)
             ctgccACCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCg                                                    <  1:55897/76‑1 (MQ=255)
               gccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGc                                                  >  1:4935587/1‑76 (MQ=255)
               gccAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGc                                                  >  2:2692730/1‑76 (MQ=255)
                 cAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                <  2:3608462/76‑1 (MQ=255)
                 cAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                >  1:4518686/1‑76 (MQ=255)
                 cAGCGCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                <  1:2396569/76‑1 (MQ=255)
                 cAGCGCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCt                                                <  1:1920165/76‑1 (MQ=255)
                  aGCTCGCCGCCAGCTAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTg                                               <  1:2193864/76‑1 (MQ=255)
                  aGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTg                                               >  2:5204841/1‑76 (MQ=255)
                  aGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGCACGCGCTCGGcacg                                               >  1:1216401/1‑73 (MQ=255)
                   gCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTgg                                              <  1:2425611/76‑1 (MQ=255)
                   gCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTgg                                              <  2:3912961/76‑1 (MQ=255)
                     tCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGa                                            >  1:5016997/1‑76 (MQ=255)
                     ccgccgcCAGATAACGCGGATCTGCCGGGTTTTCGAGAAACCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGa                                            <  1:2356749/75‑1 (MQ=255)
                      cgccgcCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGAc                                           <  1:5237680/76‑1 (MQ=255)
                      cgccgcCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGAc                                           >  1:1596285/1‑76 (MQ=255)
                      ccccgccAGATAACGCGGATCTCCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCCGCCCGGGAc                                           >  2:3906017/3‑76 (MQ=255)
                       gccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACt                                          >  1:5681514/1‑76 (MQ=255)
                        ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTc                                         >  1:1891513/1‑76 (MQ=255)
                        ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTc                                         <  2:1474299/76‑1 (MQ=255)
                        ccgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCGGGGACTc                                         >  2:5561185/1‑76 (MQ=255)
                        ccgccAGATAACGCGGATCTGCCGCGTTTGCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTAGGCCCGCGACTc                                         >  2:3109394/1‑76 (MQ=255)
                         cgcgAGATACAGCGGTTCTGCCGGGTCTTCGCGAAACCGCGTCAGCTGCGCCCGTCCGTGGTCGGCCTGGGACTCa                                        <  2:3583652/76‑1 (MQ=255)
                         cgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCa                                        <  2:3205501/76‑1 (MQ=255)
                         cgccAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCa                                        <  1:3804036/76‑1 (MQ=255)
                            cAGATAGCGGGTACCTGCCGGGTTGTCGAGGGAACGCGCCAGCTCCGCGCGTACGCGCTCGGCCTGGGACGCATcc                                     <  2:739776/76‑1 (MQ=255)
                            cAGATAAGGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATcc                                     <  2:5215298/76‑1 (MQ=255)
                            cAGAGAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATcc                                     <  2:5187874/76‑1 (MQ=255)
                             aGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                     >  2:4012462/1‑59 (MQ=255)
                             aGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTc                                                     <  1:4012463/59‑1 (MQ=255)
                              gATAACGCGGATCTGCCGGTTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgc                                   <  1:5764943/76‑1 (MQ=255)
                              gATAACGCGGATCTGCCGGGTTTTCTAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgc                                   <  1:4948561/76‑1 (MQ=255)
                              gATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgc                                   >  1:1202592/1‑76 (MQ=255)
                                tAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgcca                                 <  1:2416144/76‑1 (MQ=255)
                                tAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgcca                                 <  1:41455/76‑1 (MQ=255)
                                  aCGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagc                               <  2:3219980/76‑1 (MQ=255)
                                   cgcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcc                              >  2:2685631/1‑76 (MQ=255)
                                   cgcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTACGCCCGTACGCGCTCGGCCCGGGACTCCTCCcccagcc                              >  2:408700/1‑76 (MQ=255)
                                    ggggATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcca                             <  1:1295893/74‑1 (MQ=255)
                                    gcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGCCTCATCCgccagcca                             >  2:4156764/1‑76 (MQ=255)
                                    gcgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcca                             >  2:4580188/1‑76 (MQ=255)
                                    gcgGATACGCCGGGTGTCCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcca                             <  1:1644929/76‑1 (MQ=255)
                                    ccgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCgccagcca                             <  1:2887651/75‑1 (MQ=255)
                                     cgGATCTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGccagccag                            >  2:4831005/1‑76 (MQ=255)
                                      ggATCTGCCGGGTTGTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGa                           <  1:2180373/76‑1 (MQ=255)
                                          cTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATc                       <  2:4485044/76‑1 (MQ=255)
                                          cTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATc                       <  1:4358391/76‑1 (MQ=255)
                                          cTGCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATc                       >  1:998636/1‑76 (MQ=255)
                                            gCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGc                     <  1:3381825/76‑1 (MQ=255)
                                            gCCGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGc                     >  1:4658329/1‑76 (MQ=255)
                                             cggggTTTTGGGGAAACCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCt                    <  2:5681512/74‑1 (MQ=255)
                                             ccGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCt                    <  1:5193963/76‑1 (MQ=255)
                                              ngggTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCtt                   >  1:4118036/2‑76 (MQ=255)
                                              ggggTTTTCGATAAAACGCGCGAGCTCCTCCCGTACGGGCTCGTCCTGGGACTCATCCGCCAGCCAGACATCGCtt                   <  1:1983239/75‑1 (MQ=255)
                                              ggggTTTACGAGAAACCGCGCCAGCCCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCtt                   <  1:1342814/75‑1 (MQ=255)
                                              cGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTCCGCGCTCGGCCTGGGACTCATACGCCAGCCAGACATCGCtt                   >  1:4117997/1‑76 (MQ=255)
                                              cGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCCGACATCGCtt                   >  1:165772/1‑76 (MQ=255)
                                              cGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCtt                   >  1:4509638/1‑76 (MQ=255)
                                              cGGGTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCtt                   >  2:3046046/1‑76 (MQ=255)
                                                gtttttCGAGAAAGCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTg                 <  2:2654543/74‑1 (MQ=255)
                                                agTTTTCGCGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATGCGCAAGTGAGACATCGCTTTg                 >  1:4619083/2‑76 (MQ=255)
                                                 gTTTTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTgg                >  1:2907839/1‑76 (MQ=255)
                                                 gTTGTCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTgg                <  2:2819004/76‑1 (MQ=255)
                                                  ttttCGAGAAATCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGt               <  1:3961288/76‑1 (MQ=255)
                                                   tttCGAGAAAACGCGCCAGCTCCGCGCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGtt              <  1:1396003/76‑1 (MQ=255)
                                                   tttCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGtt              >  1:5027069/1‑76 (MQ=255)
                                                     tCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAt            <  1:3792475/76‑1 (MQ=255)
                                                     tCGAGAAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAt            >  1:2709659/1‑76 (MQ=255)
                                                       gagaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgt          >  2:2197243/1‑76 (MQ=255)
                                                       gagaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgt          >  2:456086/1‑76 (MQ=255)
                                                        agaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgtt         >  1:4551143/1‑76 (MQ=255)
                                                         gaAACCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttg        <  2:2130127/76‑1 (MQ=255)
                                                         gaAAACGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttg        >  2:2575246/1‑76 (MQ=255)
                                                          aaaaCGGGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttgt       <  1:5585526/76‑1 (MQ=255)
                                                          aaaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttgt       <  1:5063251/76‑1 (MQ=255)
                                                          aaaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttgt       >  1:2750646/1‑76 (MQ=255)
                                                          aaaaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTAtgttgt       >  1:5230828/1‑76 (MQ=255)
                                                           aaccgcgCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATgttgtt      <  1:2755337/73‑1 (MQ=255)
                                                            aaCGCGCCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGttgttg     >  2:2094861/1‑76 (MQ=255)
                                                              cgcgCCAGCTCCGCCGGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGaa   <  2:2346911/76‑1 (MQ=255)
                                                               gcgcCAGCTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAAt  <  1:5260439/76‑1 (MQ=255)
                                                                     ||                                                                    
CGGTATGGCCTGCCTGCCAGCTCGCCGCCAGATAACGCGGATCTGCCGGGTTTTCGAGAAAACGCGCCACGTCCGCCCGTACGCGCTCGGCCTGGGACTCATCCGCCAGCCAGACATCGCTTTGGTTATGTTGTTGAAT  >  NC_000913/3558863‑3559001

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