Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 478,533 C→T A29T (GCC→ACC)  ylaB ← putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,5330CT100.0% 180.1 / NA 56A29T (GCC→ACC) ylaBputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/27);  total (29/27)

AGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGT  >  NC_000913/478402‑478663
                                                                                                                                   |                                                                                                                                  
aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             >  2:47045/1‑139 (MQ=255)
aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             <  2:1474956/139‑1 (MQ=255)
aGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             >  1:1474956/1‑139 (MQ=255)
      ttGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCtt                                                                                                                       >  1:1345813/1‑139 (MQ=255)
            tCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGGGACGCGGGAGGAATAGGTATCCTGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGGCAGATGCTTAAGCCg                                                                                                                 >  2:1354790/1‑139 (MQ=255)
               ttCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACa                                                                                                              >  1:903914/1‑139 (MQ=255)
                tCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAgg                                                                                                            <  2:1683155/139‑1 (MQ=255)
                       ccTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                      >  1:190594/1‑139 (MQ=255)
                        cTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  1:1159503/114‑1 (MQ=255)
                        cTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  2:1159503/1‑114 (MQ=255)
                          gCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATAc                                                                                                   >  1:1666895/1‑139 (MQ=255)
                              cgcATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGaa                                                                                               <  2:190594/139‑1 (MQ=255)
                                      tCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAg                                                                                       <  2:903914/139‑1 (MQ=255)
                                        ccTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     <  2:1666895/139‑1 (MQ=255)
                                           tgtgTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  1:1416850/1‑95 (MQ=255)
                                           tgtgTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  2:1416850/95‑1 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             <  1:953204/94‑1 (MQ=255)
                                             tgtCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             >  2:953204/1‑94 (MQ=255)
                                               tCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGa                                                                              >  2:8092/1‑139 (MQ=255)
                                                  ccACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  1:666279/1‑88 (MQ=255)
                                                  ccACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  2:666279/88‑1 (MQ=255)
                                                         aTTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  2:1354176/1‑81 (MQ=255)
                                                         aTTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  1:1354176/81‑1 (MQ=255)
                                                              cTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                               <  1:8092/139‑1 (MQ=255)
                                                               tCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGa                                                              >  1:997929/1‑139 (MQ=255)
                                                                  aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         <  1:1464491/93‑1 (MQ=255)
                                                                  aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         >  2:1464491/1‑93 (MQ=255)
                                                                  aaTAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGGAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGt                                                           >  2:855374/1‑139 (MQ=255)
                                                                      gCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCgtg                                                       <  1:47045/139‑1 (MQ=255)
                                                                       cGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCgtgt                                                      >  1:371568/1‑139 (MQ=255)
                                                                           gcgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                               >  1:1565402/1‑126 (MQ=255)
                                                                           gcgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAg                                                               <  2:1565402/126‑1 (MQ=255)
                                                                            cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTa                                                                                                                      <  1:1293848/70‑1 (MQ=255)
                                                                            cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTa                                                                                                                      >  2:1293848/1‑70 (MQ=255)
                                                                            cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                      <  2:1663311/86‑1 (MQ=255)
                                                                            cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCaa                                                                                                      >  1:1663311/1‑86 (MQ=255)
                                                                            cgGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCa                                                 <  2:997929/139‑1 (MQ=255)
                                                                             gggAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCAc                                                >  2:1371034/1‑139 (MQ=255)
                                                                                 ggaATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         >  1:987589/1‑78 (MQ=255)
                                                                                 ggaATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGgcag                                                                                                         <  2:987589/78‑1 (MQ=255)
                                                                                      aGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATc                                       <  1:1308025/139‑1 (MQ=255)
                                                                                           tCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCttat                                  >  1:1365323/1‑139 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             <  2:1292090/46‑1 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGtcagccag                                                                                                                             >  1:1292090/1‑46 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTg                                                                                                 <  1:463141/74‑1 (MQ=255)
                                                                                             cAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTg                                                                                                 >  2:463141/1‑74 (MQ=255)
                                                                                               gCTCTTCAATAAACGATGTTTCTACCTGCTAATGGGTCAGCCATATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGAACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAac                              <  1:1354790/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:1345813/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:612376/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  1:855374/139‑1 (MQ=255)
                                                                                                ctctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  2:1365323/139‑1 (MQ=255)
                                                                                                 tctTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  1:417284/41‑1 (MQ=255)
                                                                                                 tctTCAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  2:417284/1‑41 (MQ=255)
                                                                                                 tctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     <  2:603611/82‑1 (MQ=255)
                                                                                                 tctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTa                                                                                     >  1:603611/1‑82 (MQ=255)
                                                                                                  ctTCAATAAACGATGTTTCTACCTGCTGATGGGTCAGGCAGATACTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGTCCAGATGTCGTGTTCTCACCAGCATATCCTTATTAAcaca                           >  1:1130952/1‑139 (MQ=255)
                                                                                                   ttCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCa                                                 <  1:118608/116‑1 (MQ=255)
                                                                                                   ttCAATAAACGATGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCa                                                 >  2:118608/1‑116 (MQ=255)
                                                                                                     cAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              >  2:722927/1‑37 (MQ=255)
                                                                                                     cAATAAACGATGTTTCTACCTGCTGATGGgtcagcca                                                                                                                              <  1:722927/37‑1 (MQ=255)
                                                                                                               tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAc                                                                  <  1:1665262/87‑1 (MQ=255)
                                                                                                               tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGAc                                                                  >  2:1665262/1‑87 (MQ=255)
                                                                                                               tGTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGt              >  2:100816/1‑139 (MQ=255)
                                                                                                                gTTTCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGta             >  2:709647/1‑139 (MQ=255)
                                                                                                                  ttCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCg                                                                                                                 <  1:822563/37‑1 (MQ=255)
                                                                                                                  ttCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCg                                                                                                                 >  2:822563/1‑37 (MQ=255)
                                                                                                                  ttCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             >  2:387076/1‑121 (MQ=255)
                                                                                                                  ttCTACCTGCTGATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAaca                             <  1:387076/121‑1 (MQ=255)
                                                                                                                          gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAg                                                                        <  2:1256175/70‑1 (MQ=255)
                                                                                                                          gctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAg                                                                        >  1:1256175/1‑70 (MQ=255)
                                                                                                                           ctgATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATtgt  >  2:477052/1‑139 (MQ=255)
                                                                                                                             gATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCAc                                                <  1:1539282/91‑1 (MQ=255)
                                                                                                                             gATGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCAc                                                >  2:1539282/1‑91 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATtg   <  1:1309343/135‑1 (MQ=255)
                                                                                                                              aTGGGTCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATtg   >  2:1309343/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                  
AGCGCCTTGCCATCTTTCCAGCTCCTGCAGCGCATCTTTCCCTTGTGTCGCCACCTTATTGGCTCGAATAGCGACGCGGGAGGAATAGGTATCCAGCTCTTCAATAAACGATGTTTCTACCTGCTGATGGGCCAGCCAGATGCTTAAGCCGACAGGCAGCAATACTGAAAGAATCAGTACTCCCGAAATAAGGCCGACCAGATGTCGTGTTCTCACCAGCATATCCTTATTAACACAGTGTAGGGAAAGTATATCTGATTGT  >  NC_000913/478402‑478663

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