Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,145,707 C→T intergenic (+38/‑157) dnaG → / → rpoD DNA primase/RNA polymerase sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,145,7070CT100.0% 688.2 / NA 175intergenic (+38/‑157)dnaG/rpoDDNA primase/RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (138/37);  total (138/37)

Only 100 of 234 total aligned reads displayed.

CTCTGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACC  >  REL606/3145634‑3145780
                                                                          |                                                                          
ctctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGcctcc                                                                         >  2:241830/1‑76 (MQ=255)
ctctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGcctcc                                                                         >  2:3737506/1‑76 (MQ=255)
 tctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCg                                                                        <  2:2671889/76‑1 (MQ=255)
  ctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGGCTCCGa                                                                       >  1:4551905/1‑76 (MQ=255)
  ctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGa                                                                       >  1:7817029/1‑76 (MQ=255)
  ctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGa                                                                       >  2:3591786/1‑76 (MQ=255)
  ctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGa                                                                       >  1:2790342/1‑76 (MQ=255)
   tGGACCTTAAACCAGCAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAACCCTCCGAc                                                                      <  2:5725597/76‑1 (MQ=255)
    ggACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCAACa                                                                     <  2:5620088/76‑1 (MQ=255)
      aCATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGc                                                                   >  1:7223904/1‑76 (MQ=255)
        ataaaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCGGAAGAGCGATCGGAAAGCCTCCGACAGCCg                                                                 <  2:2869628/73‑1 (MQ=255)
         ttAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGc                                                                >  2:1910931/1‑76 (MQ=255)
         ttAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGc                                                                >  2:6693449/1‑76 (MQ=255)
         ttAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGc                                                                >  1:1090502/1‑76 (MQ=255)
         atAAACCAAGGAGCGGGCAAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCAATCGGGAAGCCTCCGACAGCCg                                                                 <  2:7304089/75‑1 (MQ=255)
          tAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCa                                                               >  2:3142983/1‑76 (MQ=255)
           aaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc                                                              >  1:2657103/1‑76 (MQ=255)
           aaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc                                                              >  1:7459627/1‑76 (MQ=255)
           aaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc                                                              >  2:7254947/1‑76 (MQ=255)
           aaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc                                                              >  2:3499244/1‑76 (MQ=255)
           aaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc                                                              <  1:4137979/76‑1 (MQ=255)
            aaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACt                                                             <  1:1827609/76‑1 (MQ=255)
             aCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTg                                                            >  2:4909338/1‑76 (MQ=255)
             aCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTg                                                            <  1:1257850/76‑1 (MQ=255)
               ccgggAGCGGGCGAAAAATTGATTTAGCGGCTTAGGTGCCGAAGAGCGATGGGGAAGCCTCCGACAGCCGCACTga                                                           <  2:6143249/73‑1 (MQ=255)
               cAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTgag                                                          >  2:3249195/1‑76 (MQ=255)
                 aGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTgaga                                                         >  1:3083498/1‑76 (MQ=255)
                 aGGAGCTGGCGAAAAAGTGATTAAACGTCTTAAGTGCCGAAGAGCGATGGGGAAGCCTCCGACAGCCGCACTgaga                                                         <  2:1033010/76‑1 (MQ=255)
                    ggCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGc                                                      >  1:6522613/2‑76 (MQ=255)
                    aGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGc                                                      >  1:6049654/1‑76 (MQ=255)
                    aGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGc                                                      <  1:3238201/76‑1 (MQ=255)
                     gcgggCAAAAAAGTGTTTTAACGGCTAAAGTGCCAAAGAGCGATCGGAAAGCCTCCGACAGCCGCACTGAGAGGCa                                                     <  2:5334927/73‑1 (MQ=255)
                     gCTGGCGAAAAAGTGATTTAACGGCTTAAGTCCCGAAGGGCGAGCGGGAAGCCTCCGACAGCCGCACTGAGAGGCa                                                     <  1:4519564/76‑1 (MQ=255)
                      cTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGGGAGGCAg                                                    >  2:2996448/1‑76 (MQ=255)
                      cTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAg                                                    >  1:6661731/1‑76 (MQ=255)
                        ggCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCg                                                  >  1:1910623/1‑76 (MQ=255)
                        ggCAAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCg                                                  <  1:728741/76‑1 (MQ=255)
                        ggCAAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCg                                                  >  1:2053572/1‑76 (MQ=255)
                          cGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGc                                                >  2:4516233/1‑76 (MQ=255)
                          cGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGc                                                >  1:893775/1‑76 (MQ=255)
                          cGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGc                                                >  1:6507339/1‑76 (MQ=255)
                           gAAAAATTATTTTACCGGTTAAAGTCCTGAAGAGGACACGGAAACCCCCCGACAGCCGCACTGAGAGGCAGCGGCa                                               <  2:3542966/76‑1 (MQ=255)
                           aaaaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCa                                               <  1:2862279/75‑1 (MQ=255)
                            aaaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATGGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCaa                                              <  1:2908427/76‑1 (MQ=255)
                            aaaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCaa                                              >  2:5866253/1‑76 (MQ=255)
                            aaaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCaa                                              >  2:4744111/1‑76 (MQ=255)
                            aaaaaGTAATTAAACGTTTAAGGTGCCAAAGAGCACCGGGAAACCCTCCGACAGCCGCACTGAGAGGCAGCGGCaa                                              <  2:4373269/76‑1 (MQ=255)
                              aaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAat                                            >  1:3319902/1‑76 (MQ=255)
                               aaGTGATTTAACGGCTTAAGTGCAGAAGAGCGATGGGAAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAata                                           <  2:5475181/76‑1 (MQ=255)
                                 gTGATTCAACGGCTTAAGTGCGGAAGAGCGATGGGGAACCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatata                                         <  2:530464/76‑1 (MQ=255)
                                  tattttAACGGTTAAAGGCCCAAAGGGCAACGGGAAACCCCCCGCCAGCCGCACTGAGAGGCAGCGGCAAatatat                                        <  2:6977933/73‑1 (MQ=255)
                                  tGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatatat                                        >  1:1233532/1‑76 (MQ=255)
                                  tGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatatat                                        >  1:3887105/1‑76 (MQ=255)
                                  tGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatatat                                        >  2:6849900/1‑76 (MQ=255)
                                   gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAAtatata                                       >  2:7358199/1‑76 (MQ=255)
                                    aTTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATaa                                      >  1:1488895/1‑76 (MQ=255)
                                     tttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAg                                     >  2:5056322/1‑76 (MQ=255)
                                     tttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAg                                     >  1:4398693/1‑76 (MQ=255)
                                      ttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGt                                    >  2:7542780/1‑76 (MQ=255)
                                      ttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGt                                    >  1:2559484/1‑76 (MQ=255)
                                       tAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTa                                   >  1:4259171/1‑76 (MQ=255)
                                        aaCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTAc                                  >  1:1672450/1‑76 (MQ=255)
                                        aaCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTAc                                  >  2:4423793/1‑76 (MQ=255)
                                         aCGGCTTAAGTTCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACg                                 >  1:7857746/1‑76 (MQ=255)
                                         aCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACg                                 >  2:3307137/1‑76 (MQ=255)
                                          cGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTAAGc                                >  2:7076302/1‑76 (MQ=255)
                                            gCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAACTATATGAGTACGccc                              >  2:2971721/1‑76 (MQ=255)
                                              ttAAGTGCCGAAGAGCGATGGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc                            <  2:1847049/76‑1 (MQ=255)
                                              ttAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATACGTACGCCCtg                            >  1:5408290/1‑75 (MQ=255)
                                              ttAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc                            >  1:4273150/1‑76 (MQ=255)
                                              ttAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc                            >  2:5824128/1‑76 (MQ=255)
                                              ttAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc                            >  2:6271426/1‑76 (MQ=255)
                                               tAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCCAATATATAAGTACGCCCTCg                           >  1:6775389/1‑76 (MQ=255)
                                               tAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCg                           >  2:246353/1‑76 (MQ=255)
                                                 agggCCGAAGAGCAATGGGAAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTa                         <  2:4662316/73‑1 (MQ=255)
                                                 aGTGCCGAAGAGCGAGCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTa                         <  1:5797282/76‑1 (MQ=255)
                                                   tGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAAt                       <  1:3736389/76‑1 (MQ=255)
                                                    tccGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt                      >  2:7739594/2‑76 (MQ=255)
                                                    gCCAAAGAGCGATGGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt                      <  1:5502187/76‑1 (MQ=255)
                                                     ccGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa                     >  1:5368317/1‑76 (MQ=255)
                                                      cGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt                    <  1:1716358/76‑1 (MQ=255)
                                                      cGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt                    >  2:2720823/1‑76 (MQ=255)
                                                      cGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt                    >  2:2281315/1‑76 (MQ=255)
                                                      cGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt                    >  2:5553312/1‑76 (MQ=255)
                                                       gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATc                   >  2:2157872/1‑76 (MQ=255)
                                                        aaGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCg                  >  1:3711815/1‑76 (MQ=255)
                                                        aaGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCg                  >  1:6036195/1‑76 (MQ=255)
                                                         agagCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGt                 >  2:2600173/1‑76 (MQ=255)
                                                         agagCGAGGGGAAACCCCCCGCCAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGt                 <  2:1220883/76‑1 (MQ=255)
                                                              gATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg            >  2:1505432/1‑76 (MQ=255)
                                                              gATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg            >  1:484949/1‑76 (MQ=255)
                                                                tCGGGAAGCCTCCGCCAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGt          <  1:6814009/76‑1 (MQ=255)
                                                                ccGGGACGCCCCGCCCAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGt          <  2:4041063/75‑1 (MQ=255)
                                                                 cGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTa         >  2:3836562/1‑76 (MQ=255)
                                                                  gggAACCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTaa        <  1:4625924/76‑1 (MQ=255)
                                                                    gAACCCCCCGCCAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaac      <  1:2437469/76‑1 (MQ=255)
                                                                     aaGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaaca     >  1:5846428/1‑76 (MQ=255)
                                                                      acccccccACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaa    <  2:1320970/74‑1 (MQ=255)
                                                                      aGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaa    <  1:980152/76‑1 (MQ=255)
                                                                       ccccccGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaacaac   <  1:4867973/75‑1 (MQ=255)
                                                                          |                                                                          
CTCTGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACC  >  REL606/3145634‑3145780

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