Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,983,761 C→T 100% E437K (GAA→AAA)  araG ← L‑arabinose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,7610CT100.0% 280.3 / ‑6.7 82E437K (GAA→AAA) araGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/42);  minor base G (0/1);  total (39/43)

TCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983665‑1983856
                                                                                                |                                                                                               
tCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACATCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTtgtg                                                                                              <  2:2399997/100‑1 (MQ=255)
 ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             >  2:408908/1‑100 (MQ=255)
 ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                             <  2:1035744/100‑1 (MQ=255)
   ggTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCtt                                                                                           >  2:1281664/1‑100 (MQ=255)
    gTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTa                                                                                          <  1:742350/100‑1 (MQ=255)
      cGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAg                                                                                         <  2:280983/99‑1 (MQ=255)
       ggCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcg                                                                                       >  1:899048/1‑100 (MQ=255)
        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:850392/1‑100 (MQ=255)
        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:1413700/1‑100 (MQ=255)
        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                      >  1:1465628/1‑100 (MQ=255)
         cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                     <  2:1392896/100‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTTGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                  <  1:1057386/100‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTTGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                  <  1:1244675/100‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGTAAACACCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCaa                                                                                   <  2:1728075/99‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                  <  2:1717054/100‑1 (MQ=255)
             tccGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                 <  1:378622/99‑1 (MQ=255)
             gCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                 <  1:896365/100‑1 (MQ=255)
              ccGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACaa                                                                                >  2:378618/1‑99 (MQ=255)
                gAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCa                                                                              >  1:288642/1‑100 (MQ=255)
                 aGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCaa                                                                             <  1:1119933/100‑1 (MQ=255)
                    aCTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGc                                                                          >  1:305827/1‑100 (MQ=255)
                         aGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgc                                                                     <  1:602021/100‑1 (MQ=255)
                          ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                    >  1:309388/1‑100 (MQ=255)
                            tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                  <  1:448581/100‑1 (MQ=255)
                            tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                  <  2:303878/100‑1 (MQ=255)
                              aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                >  1:843409/1‑100 (MQ=255)
                              aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                <  2:1449361/100‑1 (MQ=255)
                               gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                               <  2:185136/100‑1 (MQ=255)
                                tCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATCTACTTTATGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                >  2:39257/1‑98 (MQ=255)
                                tCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                              <  2:954026/100‑1 (MQ=255)
                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            >  1:1830097/1‑100 (MQ=255)
                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                            <  1:1817337/100‑1 (MQ=255)
                                     ggcggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                         >  1:31967/4‑100 (MQ=255)
                                      gAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCa                                                        >  2:832552/1‑100 (MQ=255)
                                       aGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAg                                                       >  1:1593733/1‑100 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                      >  1:347190/1‑100 (MQ=255)
                                         gCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCa                                                     >  1:1100659/1‑100 (MQ=255)
                                         gCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCa                                                     >  2:718589/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  1:1983276/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  1:804781/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  2:537967/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                    >  2:1508653/1‑100 (MQ=255)
                                              cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                                <  2:2264477/100‑1 (MQ=255)
                                              cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                                >  1:11692/1‑100 (MQ=255)
                                               aGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGa                                               >  1:1206993/1‑100 (MQ=255)
                                                 cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             <  1:2496897/100‑1 (MQ=255)
                                                 cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                             <  1:648643/100‑1 (MQ=255)
                                                  ttccgatctcccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCt                                            <  1:383724/91‑1 (MQ=255)
                                                   ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtt                                           <  1:993785/100‑1 (MQ=255)
                                                     gCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCa                                         >  1:2274768/1‑100 (MQ=255)
                                                           ttatGCGCTGGCGGCGCATAAATTACGTTATATATTTTGTGCTTAGGGCTAACATCAATGTCACGAGAAGGTTCGTCGCGCACAATGACCTTCAt                                        >  2:25922/4‑95 (MQ=255)
                                                           cccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtagatcggaa                                   >  2:383720/1‑91 (MQ=255)
                                                           cccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCtt                                   >  1:779624/1‑100 (MQ=255)
                                                            ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                  >  2:1412807/1‑100 (MQ=255)
                                                            ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                  >  1:631725/1‑100 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  1:671727/100‑1 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  2:1792177/100‑1 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  2:850385/100‑1 (MQ=255)
                                                              tGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                 <  2:2373943/99‑1 (MQ=255)
                                                                   cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                           <  1:2379290/100‑1 (MQ=255)
                                                                   cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                           >  2:514042/1‑100 (MQ=255)
                                                                   cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                           <  1:144804/100‑1 (MQ=255)
                                                                    gccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAcc                          <  1:1452038/100‑1 (MQ=255)
                                                                      cAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAg                        <  1:1254207/100‑1 (MQ=255)
                                                                       aGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGc                       >  1:1043020/1‑100 (MQ=255)
                                                                        tctcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCTATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCg                      <  2:1634305/97‑1 (MQ=255)
                                                                          gcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGc                    <  1:353210/100‑1 (MQ=255)
                                                                          gcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGc                    <  1:1138417/100‑1 (MQ=255)
                                                                           cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgaga                   >  1:1634317/1‑97 (MQ=255)
                                                                            aTAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTATGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGccc                  >  2:2104663/1‑100 (MQ=255)
                                                                              aaaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAg                <  2:276664/100‑1 (MQ=255)
                                                                              aaaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAg                <  2:276679/100‑1 (MQ=255)
                                                                                aTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGaa              >  2:1668193/1‑100 (MQ=255)
                                                                                 ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt             <  2:265273/100‑1 (MQ=255)
                                                                                   aCGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATgn           <  2:2322106/100‑2 (MQ=255)
                                                                                   aCGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATgg           <  2:355637/100‑1 (MQ=255)
                                                                                      ttATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtt        <  2:597606/100‑1 (MQ=255)
                                                                                       tatataTTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttt       >  1:2100300/1‑100 (MQ=255)
                                                                                           tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  2:2072216/100‑1 (MQ=255)
                                                                                           tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  1:2208889/100‑1 (MQ=255)
                                                                                           tatTTGGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAACTGACCTTCATCTCCTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc   <  2:2426896/100‑1 (MQ=255)
                                                                                            atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  1:175212/1‑100 (MQ=255)
                                                                                            atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt  >  2:2187785/1‑100 (MQ=255)
                                                                                                |                                                                                               
TCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCT  >  NC_000913/1983665‑1983856

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