Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 882,722 G→A 100% R5C (CGC→TGC)  deoR ← deoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913882,7220GA100.0% 257.9 / ‑6.7 81R5C (CGC→TGC) deoRdeoxyribose‑5‑phosphate‑inducible deoxyribose operon transcriptional repressor; repressor of nupG and tsx
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/44);  minor base T (0/1);  total (36/45)

GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGG  >  NC_000913/882625‑882816
                                                                                                 |                                                                                              
gTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAAc                                                                                              >  1:489195/1‑100 (MQ=255)
gTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAAc                                                                                              <  1:353319/100‑1 (MQ=255)
 tCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACg                                                                                             <  2:1417223/100‑1 (MQ=255)
   aTCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACgt                                                                                            >  1:598463/1‑99 (MQ=255)
     ctcCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGtt                                                                                         >  1:2233446/1‑100 (MQ=255)
       ccGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTc                                                                                       >  2:2042509/1‑100 (MQ=255)
        cGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTcc                                                                                      <  1:803759/100‑1 (MQ=255)
          aaaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCa                                                                                     <  1:2211253/99‑1 (MQ=255)
          aaaCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCAt                                                                                    <  1:1906530/100‑1 (MQ=255)
            aCCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATCCGCTCTTCGCAACGTGTTTCCATaa                                                                                  >  2:815412/1‑100 (MQ=255)
             ccccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAAt                                                                                 >  1:1557373/1‑100 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                <  2:656279/100‑1 (MQ=255)
              cccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                >  2:1548040/1‑100 (MQ=255)
              ccaaaGCAGGGCGGCGGCGCCTTTAAGATGTAACTTATCGCTACGTTATAACTCTTGCAGCAGCTGCCCGATACGCACTACGCAACGTGATTCCatcatc                                                                                >  2:81190/4‑100 (MQ=255)
               ccAAGCAGGGCGTCGGCGTCGTTAAGATGTAACTTATCGCTGCGGTTTGATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATc                                                                                <  1:1112837/99‑1 (MQ=255)
               ccAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATcc                                                                               >  2:1366740/1‑100 (MQ=255)
                  aGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctc                                                                            >  1:1434963/1‑100 (MQ=255)
                   gCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCctct                                                                           <  1:2079270/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                          <  1:815421/100‑1 (MQ=255)
                    cAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTg                                                                          >  2:636696/1‑100 (MQ=255)
                        gcggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAAt                                                                       <  2:1878306/99‑1 (MQ=255)
                        gcggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATa                                                                      >  1:1049604/1‑100 (MQ=255)
                          ggcggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTACATAATCCCTCTGAATAGt                                                                    >  2:528121/1‑100 (MQ=255)
                             ggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                 >  1:2291692/1‑100 (MQ=255)
                             ggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                 >  1:395834/1‑100 (MQ=255)
                             ggcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                 >  1:728151/1‑100 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAtt                                                                <  2:691950/100‑1 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTAt                                                                 <  1:374526/99‑1 (MQ=255)
                              gcgTCTTTAAGATGTAACTTATCGCTGCGGTTTAAGTCTTGCAGTAGCTGCCCGATACGCTCGTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTa                                                                  <  2:1618017/98‑1 (MQ=255)
                                gTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGa                                                              >  1:1955950/1‑100 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                            >  1:2145232/1‑100 (MQ=255)
                                  cTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAg                                                            <  1:6376/100‑1 (MQ=255)
                                       gatcttAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                       <  2:1158753/95‑1 (MQ=255)
                                       aGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGc                                                       <  2:2191763/100‑1 (MQ=255)
                                        gATGTAACTTATCGCTGCGAATTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTAAAGCGAGc                                                       >  2:562820/1‑99 (MQ=255)
                                         aTGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCTTCTTAATAGTTATTGAAGCGAGc                                                       <  2:1083409/98‑1 (MQ=255)
                                          tGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGc                                                    <  2:1147219/100‑1 (MQ=255)
                                           gTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCt                                                   >  1:1688850/1‑100 (MQ=255)
                                           gTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCt                                                   >  1:1524612/1‑100 (MQ=255)
                                            tAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGcagatc                                                  >  1:1158764/1‑95 (MQ=255)
                                             aaCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTc                                                  <  2:544317/99‑1 (MQ=255)
                                              aCTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCACCGTGTTTCCATACTCCCTCTGAATAGTTATTGAAACGAACCGCTCaa                                                >  1:1630391/1‑100 (MQ=255)
                                                  aTCGCTGCGTTTTAATTCTTGCAGCATCTTCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtact                                            <  2:225366/100‑1 (MQ=255)
                                                   tCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtacta                                           >  2:1666241/1‑100 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                          <  1:1399766/100‑1 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                          <  2:2291666/100‑1 (MQ=255)
                                                    cGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAAtactac                                          <  1:1841105/100‑1 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                         >  1:55509/1‑100 (MQ=255)
                                                     gCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTaca                                         <  2:100249/100‑1 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATCCGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                        >  1:2276440/1‑100 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                        <  2:1934490/100‑1 (MQ=255)
                                                      cTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTacac                                        >  2:719947/1‑100 (MQ=255)
                                                            tttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTa                                  <  1:119570/100‑1 (MQ=255)
                                                             ttAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAg                                 <  2:59805/100‑1 (MQ=255)
                                                               aaTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCACACAAAAATACACATTTTaaaa                               >  1:600517/1‑97 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                              <  2:1290112/100‑1 (MQ=255)
                                                                aTTCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCag                              <  2:2366151/100‑1 (MQ=255)
                                                                 ttCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCaga                             <  1:1366754/100‑1 (MQ=255)
                                                                 ttCTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCaga                             <  2:296452/100‑1 (MQ=255)
                                                                 ttCTTGCAGCAGCTGCACGATACGCTCTTCTCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCaga                             >  2:128853/1‑100 (MQ=255)
                                                                   cTTGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCagag                            >  1:1770188/1‑99 (MQ=255)
                                                                   cTTGCAGCAGCCGACCGATACGCTCTTCGCAACGTCTTTCCCAACTCTCTCCGAATAGTCATTGAAGCGCGCCGCTCACTACTACACGTGTTAGCatat                            <  2:2003344/99‑4 (MQ=255)
                                                                    ttGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGAt                          >  1:195814/1‑100 (MQ=255)
                                                                     tGCAGCAGCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATc                         <  1:1296146/100‑1 (MQ=255)
                                                                      gcagcaGCTGCCCGATACGCTCTTCGCAACGTGGGTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTCAGCAGAGATc                         >  2:2287488/1‑99 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                       <  1:2258510/100‑1 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                       <  2:653589/100‑1 (MQ=255)
                                                                       cagcagCTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAg                       >  1:463395/1‑100 (MQ=255)
                                                                             cTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGc                 >  2:235825/1‑100 (MQ=255)
                                                                             cTGCCCGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGc                 <  2:308036/100‑1 (MQ=255)
                                                                                cccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAcc              <  1:1824378/100‑1 (MQ=255)
                                                                                cccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAcc              <  2:2121148/100‑1 (MQ=255)
                                                                                 ccGATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCAccc             >  2:1889965/1‑100 (MQ=255)
                                                                                   gATACGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAg           <  2:2058526/100‑1 (MQ=255)
                                                                                     tACGCTCTTCGCTACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACGCAGcc         <  2:489192/100‑1 (MQ=255)
                                                                                       cGCTCTTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCtt       <  1:1484984/100‑1 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc    <  2:1065457/100‑1 (MQ=255)
                                                                                          tctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc    <  2:2069025/100‑1 (MQ=255)
                                                                                           ctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCATCCTTTGCg   >  2:518774/1‑100 (MQ=255)
                                                                                           ctTCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGc  <  1:2136153/100‑2 (MQ=255)
                                                                                            ttCGCAACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCgg  <  2:633692/100‑1 (MQ=255)
                                                                                                 |                                                                                              
GTCATCTCCGAAACCCCAAGCAGGGCGGCGGCGTCTTTAAGATGTAACTTATCGCTGCGTTTTAATTCTTGCAGCAGCTGCCCGATACGCTCTTCGCGACGTGTTTCCATAATCCCTCTGAATAGTTATTGAAGCGAGCCGCTCAATACTACACTTTTTAGCAGAGATCAGTCACGCACCCAGCCTTTGCGG  >  NC_000913/882625‑882816

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