Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,289 G→T intergenic (‑241/‑98) araB ← / → araC L‑ribulokinase/transcriptional regulator for ara operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,2890GT100.0% 343.4 / NA 114intergenic (‑241/‑98)araB/araCL‑ribulokinase/transcriptional regulator for ara operon
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (65/49);  total (65/49)

Only 100 of 114 total aligned reads displayed.

GTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTT  >  NC_000913/70218‑70360
                                                                       |                                                                       
 tttcgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTa                                                                    <  2:3804016/73‑1 (MQ=255)
 ttACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTACTAAGCGGGGTTACCGGTTTGGTTa                                                                    >  2:110422/1‑76 (MQ=255)
 ttACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTa                                                                    <  2:1864258/76‑1 (MQ=255)
   tcgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGc                                                                  <  2:2320674/75‑1 (MQ=255)
   tcgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGc                                                                  <  1:2229122/75‑1 (MQ=255)
    cgcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCg                                                                 <  1:4371788/76‑1 (MQ=255)
     gcgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCga                                                                <  1:4367375/76‑1 (MQ=255)
      cgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgag                                                               >  2:4516014/1‑76 (MQ=255)
      cgTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgag                                                               >  2:1556979/1‑76 (MQ=255)
       gTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCgaga                                                              >  2:175028/1‑76 (MQ=255)
       gTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGAGGGGTTACCGGTTTCGTTAGCgaga                                                              >  2:857321/1‑76 (MQ=255)
        tttttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaa                                                             <  1:4554321/76‑1 (MQ=255)
          tttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga                                                           <  2:4329640/76‑1 (MQ=255)
          tttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGagaaga                                                           <  2:1418126/76‑1 (MQ=255)
           ttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAatag                                                          >  2:208312/1‑76 (MQ=255)
           ttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag                                                          >  2:649289/1‑76 (MQ=255)
           ttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag                                                          >  2:144657/1‑76 (MQ=255)
           ttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag                                                          >  1:3080062/1‑76 (MQ=255)
           ttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAagag                                                          >  2:115831/1‑76 (MQ=255)
            tGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGc                                                         <  2:3767865/76‑1 (MQ=255)
             gTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGcc                                                        >  1:551589/1‑76 (MQ=255)
              gcATGGCTTTGGTCGCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCa                                                       <  2:951414/75‑1 (MQ=255)
                aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGCGCCAgc                                                     >  2:746572/1‑75 (MQ=255)
                aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGt                                                     >  2:4416279/1‑76 (MQ=255)
                aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGt                                                     <  1:4728712/76‑1 (MQ=255)
                aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGt                                                     <  1:4832868/76‑1 (MQ=255)
                  tggtttGGTGCGGTTTTGGTACAGAGTGGTTTTAATAAGGGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaa                                                   <  2:46258/73‑1 (MQ=255)
                  ggCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaa                                                   >  1:3781/1‑76 (MQ=255)
                  ggCTTGGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaa                                                   <  2:4302768/76‑1 (MQ=255)
                   gCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGCGAAGAGCCAGTaaa                                                  >  2:950668/1‑76 (MQ=255)
                   gCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaaa                                                  >  1:3213125/1‑76 (MQ=255)
                   gCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTaaa                                                  >  2:1277621/1‑76 (MQ=255)
                    gtttGGTCCCGGTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACGGGTTTGGTTAGCGAGAAGAGCCAGTaaaa                                                 <  2:1232866/75‑1 (MQ=255)
                     tttGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAg                                                >  2:4176495/1‑76 (MQ=255)
                     tttGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAg                                                >  1:2731717/1‑76 (MQ=255)
                       tGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGAc                                              <  1:2966650/76‑1 (MQ=255)
                        ggTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACg                                             >  1:4293734/1‑76 (MQ=255)
                         gTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGc                                            >  1:1366626/1‑76 (MQ=255)
                          tCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCa                                           <  1:1004668/76‑1 (MQ=255)
                          tCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCa                                           >  2:863687/1‑76 (MQ=255)
                           cccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAg                                          >  1:4614528/1‑76 (MQ=255)
                           cccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAg                                          >  2:348774/1‑76 (MQ=255)
                           cccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAg                                          <  1:1445401/76‑1 (MQ=255)
                           cccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAg                                          <  2:4100258/76‑1 (MQ=255)
                             cGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTg                                        >  1:4470064/1‑76 (MQ=255)
                             cGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTg                                        >  1:1880611/1‑76 (MQ=255)
                             cGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTg                                        <  1:4618656/76‑1 (MQ=255)
                              gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGa                                       >  1:2435027/1‑76 (MQ=255)
                              gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGa                                       >  2:4693889/1‑76 (MQ=255)
                                ttttttACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGAca                                     <  2:3863385/76‑2 (MQ=255)
                                tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACg                                     >  1:2198741/1‑76 (MQ=255)
                                tttGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACg                                     >  2:4162553/1‑76 (MQ=255)
                                  tGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGc                                   >  1:2058732/1‑76 (MQ=255)
                                  tGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGc                                   >  2:2917258/1‑76 (MQ=255)
                                   gTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCa                                  <  1:3183204/76‑1 (MQ=255)
                                    ttACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCaa                                 <  2:1957246/76‑1 (MQ=255)
                                     tACAGAATGCTTTTAATAAGGGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAAt                                <  2:1370302/76‑1 (MQ=255)
                                     tACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAAt                                <  2:3218581/76‑1 (MQ=255)
                                     tACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAAt                                >  1:4551814/1‑76 (MQ=255)
                                      aCAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATg                               <  1:4663470/76‑1 (MQ=255)
                                      aCAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATg                               >  1:1378788/1‑76 (MQ=255)
                                       aaGAATGCTTTTAATAGGCGGGGTTACCGGTTTGGTTAGCGGGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGt                              <  2:2036128/75‑1 (MQ=255)
                                          aaTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTg                           >  2:3040602/1‑76 (MQ=255)
                                          aaTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTg                           >  1:4291274/1‑76 (MQ=255)
                                           aTGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGa                          <  2:2105089/76‑1 (MQ=255)
                                            tGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGAt                         <  1:1708577/76‑1 (MQ=255)
                                             gCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATg                        <  2:3275054/76‑1 (MQ=255)
                                               ttttAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCa                      >  2:3822921/1‑76 (MQ=255)
                                               ttttAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCa                      <  1:3230041/76‑1 (MQ=255)
                                                 ttAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGCTGCAat                    >  2:4168190/1‑76 (MQ=255)
                                                 ttAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAat                    >  1:2588624/1‑76 (MQ=255)
                                                   aatagGCGGGGTTCCCGGTTTGGTTAGCGATAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAatat                  <  1:1534900/76‑1 (MQ=255)
                                                   aataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAGAAGACGCAGTGACGGCAATGTCTGATGCAatat                  <  2:1740161/76‑1 (MQ=255)
                                                    ataaGCGGGGTTACCGGTTTGGTTGGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATg                 <  2:3766631/76‑1 (MQ=255)
                                                    ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATg                 >  1:1514031/1‑76 (MQ=255)
                                                    ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATg                 <  1:4218752/76‑1 (MQ=255)
                                                     taaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATgg                >  1:908028/1‑76 (MQ=255)
                                                     taaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATgg                >  2:1619403/1‑76 (MQ=255)
                                                     taaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATgg                <  2:4434650/76‑1 (MQ=255)
                                                     aaaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATgg                <  1:4516014/75‑1 (MQ=255)
                                                      aaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGa               >  2:3853101/1‑76 (MQ=255)
                                                        ggggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACa             <  2:1056342/74‑1 (MQ=255)
                                                        gCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACa             >  1:3501506/1‑76 (MQ=255)
                                                          ggggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAt           >  2:3926993/1‑76 (MQ=255)
                                                            ggTTGCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTg         <  1:4420507/76‑1 (MQ=255)
                                                            ggTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTg         <  1:283985/76‑1 (MQ=255)
                                                             gTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACCATTgg        >  1:296409/1‑76 (MQ=255)
                                                             gTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTgg        >  2:4215134/1‑76 (MQ=255)
                                                             gTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTgg        <  1:3049395/76‑1 (MQ=255)
                                                              ttACCCGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGt       >  1:2552772/1‑76 (MQ=255)
                                                                aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGttt     <  2:2836546/76‑1 (MQ=255)
                                                                  cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttct   >  1:2751724/1‑76 (MQ=255)
                                                                  cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttct   >  2:1285314/1‑76 (MQ=255)
                                                                  cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttct   <  1:70488/76‑1 (MQ=255)
                                                                  cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttct   >  1:3675740/1‑76 (MQ=255)
                                                                  cGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAACATGGCCAATTGGGttc  >  2:2448324/1‑75 (MQ=255)
                                                                   ggTTTGGTTAGCGAGAAGAGCCGGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttctt  <  2:3777301/76‑1 (MQ=255)
                                                                   ggTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttctt  >  1:2130704/1‑76 (MQ=255)
                                                                   ggTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGAAATATGGACAATTGGTttctt  >  1:347256/1‑76 (MQ=255)
                                                                   ggTTTGGTTAGCGAGAAGAGCCAGCAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTttctt  >  2:4488124/1‑76 (MQ=255)
                                                                       |                                                                       
GTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTT  >  NC_000913/70218‑70360

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