Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,691,360 G→C Y416* (TAC→TAG bcsB ← cellulose synthase periplasmic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,691,3600GC100.0% 245.3 / NA 63Y416* (TAC→TAGbcsBcellulose synthase periplasmic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (32/31);  total (32/31)

CTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGC  >  NC_000913/3691227‑3691490
                                                                                                                                     |                                                                                                                                   
cTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCtaatta                                                                                                                               <  2:249939/140‑1 (MQ=255)
cTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATCGTGTAGCGCtaatta                                                                                                                               <  2:1117279/140‑1 (MQ=255)
  tGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCtaattaat                                                                                                                             >  1:335817/1‑140 (MQ=255)
    tAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAatat                                                                                                                           >  1:374833/1‑140 (MQ=255)
      aTACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATcc                                                                                                                         <  2:1066288/140‑1 (MQ=255)
          cGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCtaatta                                                                                                                               <  2:275278/130‑1 (MQ=255)
          cGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCtaatta                                                                                                                               >  1:275278/1‑130 (MQ=255)
          cGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCatat                                                                                                                     >  2:975531/1‑140 (MQ=255)
          cGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCatat                                                                                                                     >  2:1145034/1‑140 (MQ=255)
           ggAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATc                                                                                                                    <  2:886835/140‑1 (MQ=255)
           ggAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATc                                                                                                                    <  2:374833/140‑1 (MQ=255)
           ggAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATc                                                                                                                    <  1:1145034/140‑1 (MQ=255)
                 cGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATc                                                                                                                          >  1:320200/1‑128 (MQ=255)
                 cGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATc                                                                                                                          <  2:320200/128‑1 (MQ=255)
                   cagcagcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCgg                                                                                                            >  2:201224/1‑140 (MQ=255)
                      cagcagGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGACATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTAc                                                                                                         >  1:113797/1‑140 (MQ=255)
                              ggTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATc                                                                                                <  1:946671/140‑1 (MQ=255)
                              ggTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATc                                                                                                <  1:975531/140‑1 (MQ=255)
                                ttCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAg                                                                                              >  2:376576/1‑140 (MQ=255)
                                          gTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCa                                                                                               >  1:434751/1‑129 (MQ=255)
                                          gTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCa                                                                                               <  2:434751/129‑1 (MQ=255)
                                               ctgctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcggcgg                                                                               <  1:376576/140‑1 (MQ=255)
                                               ctgctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGACATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcggcgg                                                                               <  2:113797/140‑1 (MQ=255)
                                                  ctCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAg                                                                            >  1:591331/1‑140 (MQ=255)
                                                     aGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGtt                                                                         >  2:464621/1‑140 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCCATCAGGTAGAGATCCGGCGGCAGGTTTAGCg                                                                    <  2:1150055/134‑1 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCCATCAGGTAGAGATCCGGCGGCAGGTTTAGCg                                                                    >  1:1150055/1‑134 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcggcgg                                                                               <  2:826534/122‑1 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCcggcgg                                                                               >  1:826534/1‑122 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGtt                                                                         <  1:964524/128‑1 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGtt                                                                         >  2:964524/1‑128 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGtt                                                             >  1:559354/1‑140 (MQ=255)
                                                                 tGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGtt                                                             >  1:502638/1‑140 (MQ=255)
                                                                    aGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAAt                                                          <  2:591331/140‑1 (MQ=255)
                                                                      gAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCg                                                        >  2:42160/1‑140 (MQ=255)
                                                                         cTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgctg                                                     >  1:213383/1‑140 (MQ=255)
                                                                                ttCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTc                                              >  1:931515/1‑140 (MQ=255)
                                                                                  cAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCaa                                            <  1:464621/140‑1 (MQ=255)
                                                                                    ggCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGa                                          <  2:559354/140‑1 (MQ=255)
                                                                                     gCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGAc                                         >  1:670813/1‑140 (MQ=255)
                                                                                     gCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGAc                                         >  1:986648/1‑140 (MQ=255)
                                                                                     gCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGATTAGCGAAACGTTAATCGCTGCTGGCTCAAGAc                                         >  1:852234/1‑140 (MQ=255)
                                                                                      cTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGAcc                                        >  2:395227/1‑140 (MQ=255)
                                                                                       tGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCg                                       <  2:986648/140‑1 (MQ=255)
                                                                                       tGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCg                                       <  2:852234/140‑1 (MQ=255)
                                                                                        gATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGc                                      <  1:86904/140‑1 (MQ=255)
                                                                                        gATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGc                                      <  2:618159/140‑1 (MQ=255)
                                                                                         atatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgctgc                                                    <  2:228316/125‑1 (MQ=255)
                                                                                         atatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCgctgc                                                    >  1:228316/1‑125 (MQ=255)
                                                                                         atatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCt                                     >  2:147125/1‑140 (MQ=255)
                                                                                          tatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTg                                    <  1:42160/140‑1 (MQ=255)
                                                                                          tatCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTg                                    <  2:931515/140‑1 (MQ=255)
                                                                                               atccGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTg                                    >  2:548049/1‑135 (MQ=255)
                                                                                               atccGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTg                                    <  1:548049/135‑1 (MQ=255)
                                                                                               atccGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTg                               >  1:284059/1‑140 (MQ=255)
                                                                                               atccGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTg                               >  1:1096283/1‑140 (MQ=255)
                                                                                                  cgcgAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTaa                            <  2:670813/140‑1 (MQ=255)
                                                                                                   gcgAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAAc                           >  2:859977/1‑140 (MQ=255)
                                                                                                     gAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTg                         <  2:502638/140‑1 (MQ=255)
                                                                                                     gAACTGTCTTTCACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTg                         <  1:147125/140‑1 (MQ=255)
                                                                                                                cACCGGCGGCATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAgg              >  2:429713/1‑140 (MQ=255)
                                                                                                                     gcggcATGGTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTc         >  1:10000/1‑140 (MQ=255)
                                                                                                                            ggTGTAGCGCTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGc  <  1:201224/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CTTGTAATACCGGAATCCGCAGCAGCAGGCGGTTCGCCTCCTGTTTGCTGCTCAGGTTGAAGGATTGCAGGAACTGGTTATTCAGGCTGATATCCATCCGCGAACTGTCTTTCACCGGCGGCATGGTGTAGCGGTAATTAATATCCATATCAATGCCGGTACTGCGCATCAGGTAGAGATCCGGCGGCAGGTTTAGCGAAACGTTAATCGCTGCTGGCTCAAGACCGCTGGATTGTAACTGTTCTTCATAGGTTTTCAGTTCGC  >  NC_000913/3691227‑3691490

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