Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,789,047 C→T A316V (GCT→GTT) 
*345* (TAG→TAA
yibQ →
waaH ←
putative polysaccharide deacetylase
LPS(HepIII)‑glucuronic acid glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,789,0470CT100.0% 217.7 / NA 67A316V (GCT→GTT) 
*345* (TAG→TAA
yibQ
waaH
putative polysaccharide deacetylase
LPS(HepIII)‑glucuronic acid glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/33);  total (34/33)

CCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCG  >  NC_000913/3788914‑3789178
                                                                                                                                     |                                                                                                                                   
ccGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                                <  1:1023147/139‑1 (MQ=255)
ccGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAAAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                                <  1:683098/139‑1 (MQ=255)
                  gCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTAGTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcg                                                                                                              <  1:1010347/139‑1 (MQ=255)
                    tAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCATAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATGAAAGGCGCTAcgcg                                                                                                            <  2:1016133/139‑1 (MQ=255)
                       tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcg                                                                                                              <  2:1369361/134‑1 (MQ=255)
                       tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTAcg                                                                                                              >  1:1369361/1‑134 (MQ=255)
                       tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                         >  2:964387/1‑139 (MQ=255)
                       tCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTaa                                                                                                         >  1:431774/1‑139 (MQ=255)
                                       ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                         >  1:1371312/1‑139 (MQ=255)
                                       ttAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACcg                                                                                         >  2:171037/1‑139 (MQ=255)
                                        tAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTGTACTTCCAGCATCAGTGGCAAGGCTGGAGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCTACcgc                                                                                        >  2:1670690/1‑139 (MQ=255)
                                               aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAg                                                                                                                      >  1:1304152/1‑102 (MQ=255)
                                               aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAg                                                                                                                      <  2:1304152/102‑1 (MQ=255)
                                               aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                                 >  1:83761/1‑139 (MQ=255)
                                               aaGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                                 >  1:1261003/1‑139 (MQ=255)
                                                 gCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                               <  2:1371312/139‑1 (MQ=255)
                                                   aTCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGAGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                             <  1:1670690/139‑1 (MQ=255)
                                                     cAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGcc                                                                           >  2:338620/1‑139 (MQ=255)
                                                            agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                                <  1:557864/79‑1 (MQ=255)
                                                            agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAgcg                                                                                                                                >  2:557864/1‑79 (MQ=255)
                                                            agCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCa                                                                    >  1:81038/1‑139 (MQ=255)
                                                               aCGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATaa                                                                 >  1:402795/1‑139 (MQ=255)
                                                                   tGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTTATGCGATGGTGCGACTGCGTTTTGTCGCGCCATTGCCATAAGCgg                                                             <  1:964387/139‑1 (MQ=255)
                                                                      tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcg                                                                                                                   >  1:841387/1‑82 (MQ=255)
                                                                      tCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcg                                                                                                                   <  2:841387/82‑1 (MQ=255)
                                                                               tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcg                                                                                                                   <  2:782705/73‑1 (MQ=255)
                                                                               tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGgcg                                                                                                                   >  1:782705/1‑73 (MQ=255)
                                                                               tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            >  1:142186/1‑96 (MQ=255)
                                                                               tCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            <  2:142186/96‑1 (MQ=255)
                                                                                 cAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTaa                                               >  2:696602/1‑139 (MQ=255)
                                                                                  aGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAAc                                              >  1:954526/1‑139 (MQ=255)
                                                                                   gCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAg                                                                <  2:1245105/120‑1 (MQ=255)
                                                                                   gCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAg                                                                >  1:1245105/1‑120 (MQ=255)
                                                                                        aGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATaa                                                                 >  2:874504/1‑114 (MQ=255)
                                                                                        aGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATGAAAGGCGCTACGCGTAATGAGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATaa                                                                 <  1:874504/114‑1 (MQ=255)
                                                                                          tGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTa                                      >  1:1441653/1‑139 (MQ=255)
                                                                                             cAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   <  2:431774/139‑1 (MQ=255)
                                                                                                ggCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                                 >  2:169653/1‑90 (MQ=255)
                                                                                                ggCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGt                                                                                 <  1:169653/90‑1 (MQ=255)
                                                                                                 gCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAAc                                              <  1:856710/124‑1 (MQ=255)
                                                                                                 gCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAAc                                              >  2:856710/1‑124 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                               >  2:1469849/1‑81 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcg                                                                               <  1:1469849/81‑1 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                               >  2:174505/1‑97 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGc                                                               <  1:174505/97‑1 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgc                     >  2:53302/1‑139 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgc                     >  2:1398840/1‑139 (MQ=255)
                                                                                                           aCAGCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTgc                     >  1:208179/1‑139 (MQ=255)
                                                                                                              gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                  <  2:402795/139‑1 (MQ=255)
                                                                                                              gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                  <  1:338620/139‑1 (MQ=255)
                                                                                                              gCCCGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGca                  <  1:53302/139‑1 (MQ=255)
                                                                                                               cccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCt                       >  1:1623382/1‑133 (MQ=255)
                                                                                                               cccGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCt                       <  2:1623382/133‑1 (MQ=255)
                                                                                                                 cGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            >  2:1542094/1‑62 (MQ=255)
                                                                                                                 cGAAGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGa                                                                                            <  1:1542094/62‑1 (MQ=255)
                                                                                                                    aGCGGCGAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTgg            >  2:929821/1‑139 (MQ=255)
                                                                                                                        gcgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAAt        >  1:157121/1‑139 (MQ=255)
                                                                                                                         cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   <  2:1440317/111‑1 (MQ=255)
                                                                                                                         cgAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAAcc                                   >  1:1440317/1‑111 (MQ=255)
                                                                                                                          gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa      <  2:208179/139‑1 (MQ=255)
                                                                                                                          gAAGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCa      <  2:83761/139‑1 (MQ=255)
                                                                                                                            aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg    <  1:171037/139‑1 (MQ=255)
                                                                                                                            aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg    <  2:81038/139‑1 (MQ=255)
                                                                                                                            aGTTTAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAgg    <  2:1261003/139‑1 (MQ=255)
                                                                                                                              tttAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                             <  2:1365546/64‑1 (MQ=255)
                                                                                                                              tttAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTcgcg                                                                             >  1:1365546/1‑64 (MQ=255)
                                                                                                                              tttAACGTTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGgcg  <  1:1398840/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCAGCATCAGTGGCAAGGCTGGGGCAAACAGCCCGAAGCGGCGAAGTTTAACGCTAGCGCAAATTAAAGGCGCTACGCGTAATGCGATGGTGCGACCGCGTTTTGTCGCGCCATTGCCATAAGCGGAAAGACCACAGTAACGCCTGGTAACCGACCTTCACGCTGCGCACATTGGTAATCAGGCG  >  NC_000913/3788914‑3789178

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