Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,368,837 C→T R33H (CGT→CAT)  rpoA ← DNA‑directed RNA polymerase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,368,8370CT100.0% 830.8 / NA 206R33H (CGT→CAT) rpoADNA‑directed RNA polymerase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (87/119);  total (87/119)

Only 100 of 209 total aligned reads displayed.

GGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTA  >  REL606/3368764‑3368909
                                                                         |                                                                        
 gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGct                                                                       >  1:4917156/1‑76 (MQ=255)
   cACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGctct                                                                     >  1:2540551/1‑76 (MQ=255)
   cACCGAGCAACCCGGCATCGATGAGAGCATAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGctct                                                                     <  2:4087073/76‑1 (MQ=255)
    aCCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTa                                                                    <  2:787183/76‑1 (MQ=255)
    aCCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTa                                                                    <  2:6413016/76‑1 (MQ=255)
     ccGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTaa                                                                   <  1:4150205/76‑1 (MQ=255)
       gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAg                                                                 <  2:574874/76‑1 (MQ=255)
       gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAg                                                                 <  2:3379764/76‑1 (MQ=255)
        cgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAgg                                                                >  2:1440950/1‑76 (MQ=255)
         gcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCCAAGCCATGCTCTAAAGGt                                                               >  1:1024468/1‑76 (MQ=255)
         gcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGt                                                               <  2:1262944/76‑1 (MQ=255)
           aaCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTc                                                             <  2:134285/76‑1 (MQ=255)
           aaCCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTc                                                             >  1:3428966/1‑76 (MQ=255)
             cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCaa                                                           >  1:4592262/1‑76 (MQ=255)
               cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAgg                                                         <  2:5386795/76‑1 (MQ=255)
               cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAgg                                                         <  1:2698920/76‑1 (MQ=255)
                ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAggg                                                        <  2:3826181/76‑1 (MQ=255)
                ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAggg                                                        >  1:3044183/1‑76 (MQ=255)
                ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAggg                                                        <  2:2830758/76‑1 (MQ=255)
                 gCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGt                                                       >  1:2330383/1‑76 (MQ=255)
                   aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCa                                                     >  2:2207395/1‑76 (MQ=255)
                   aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCa                                                     >  2:2988539/1‑76 (MQ=255)
                   aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCa                                                     >  2:6640990/1‑76 (MQ=255)
                       aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCtt                                                 <  2:237031/76‑1 (MQ=255)
                       aTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCtt                                                 <  2:570816/76‑1 (MQ=255)
                        tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTg                                                <  2:6130117/76‑1 (MQ=255)
                        tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTg                                                <  2:4003546/76‑1 (MQ=255)
                        tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTg                                                <  2:4389342/76‑1 (MQ=255)
                        tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTg                                                <  2:874255/76‑1 (MQ=255)
                        tGAGAGCAGAACACGGCGCAGTGCATTGCCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTg                                                <  1:4298475/76‑1 (MQ=255)
                          agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgc                                              <  2:227338/76‑1 (MQ=255)
                          agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgc                                              <  2:4634178/76‑1 (MQ=255)
                          agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgc                                              <  1:6350698/76‑1 (MQ=255)
                          agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgc                                              >  2:6204855/1‑76 (MQ=255)
                           gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgcg                                             >  1:4042881/1‑76 (MQ=255)
                             gCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgcgtg                                           <  2:1486577/76‑1 (MQ=255)
                               aGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgcgtgcg                                         >  1:873905/1‑76 (MQ=255)
                                gAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGgcgtgcgt                                        >  1:3155530/1‑76 (MQ=255)
                                 aaTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTc                                       >  1:3592617/1‑76 (MQ=255)
                                 aaTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTc                                       >  1:2914263/1‑76 (MQ=255)
                                 aaTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTc                                       >  1:4643578/1‑76 (MQ=255)
                                  aTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCg                                      <  2:2452529/76‑1 (MQ=255)
                                  aTACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCg                                      <  2:5160284/76‑1 (MQ=255)
                                    aCGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGaa                                    >  1:4953867/1‑76 (MQ=255)
                                    aCGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGaa                                    >  2:3459463/1‑76 (MQ=255)
                                    aCGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGaa                                    <  1:2420175/76‑1 (MQ=255)
                                    aCGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGaa                                    >  1:1268038/1‑76 (MQ=255)
                                     cGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAAc                                   <  1:167768/76‑1 (MQ=255)
                                     cGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAAc                                   <  1:1345422/76‑1 (MQ=255)
                                      ggCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACt                                  >  2:2099982/1‑76 (MQ=255)
                                        cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                <  1:3645230/76‑1 (MQ=255)
                                          cAGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                              >  1:5237857/1‑76 (MQ=255)
                                           aGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                             >  1:1929982/1‑76 (MQ=255)
                                           aGTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                             <  1:3910142/76‑1 (MQ=255)
                                            gTGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTg                            <  2:5238517/76‑1 (MQ=255)
                                             tGCGTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGc                           >  2:191569/1‑76 (MQ=255)
                                                gTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                        <  1:2161556/76‑1 (MQ=255)
                                                gTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                        >  2:1004977/1‑76 (MQ=255)
                                                gTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                        <  1:2736789/76‑1 (MQ=255)
                                                gTTACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCg                        >  2:6396030/1‑76 (MQ=255)
                                                 ttACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGa                       >  1:1076445/1‑76 (MQ=255)
                                                 ttACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGa                       >  1:5781050/1‑76 (MQ=255)
                                                  tACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                      <  2:2634691/76‑1 (MQ=255)
                                                  tACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                      >  1:730704/1‑76 (MQ=255)
                                                  tACCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                      >  1:2718409/1‑76 (MQ=255)
                                                   aCCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                     <  2:3776594/76‑1 (MQ=255)
                                                   aCCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                     <  1:1790811/76‑1 (MQ=255)
                                                   aCCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                     >  2:3899693/1‑76 (MQ=255)
                                                   aCCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                     <  2:4883031/76‑1 (MQ=255)
                                                   aCCCAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                     >  1:5729690/1‑76 (MQ=255)
                                                     ccAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATc                   <  1:398888/76‑1 (MQ=255)
                                                     ccAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATc                   <  1:5371815/76‑1 (MQ=255)
                                                      cAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCa                  >  2:5659689/1‑76 (MQ=255)
                                                      cAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCa                  <  1:6382784/76‑1 (MQ=255)
                                                      cAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCa                  >  2:352827/1‑76 (MQ=255)
                                                      cAGAGTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCa                  >  2:6576573/1‑76 (MQ=255)
                                                       agagTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCaa                 <  2:2900407/76‑1 (MQ=255)
                                                       agagTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCaa                 >  2:1166759/1‑76 (MQ=255)
                                                       agagTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCaa                 <  1:3860468/76‑1 (MQ=255)
                                                         agTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc               <  1:162720/76‑1 (MQ=255)
                                                         agTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc               <  1:327896/76‑1 (MQ=255)
                                                         agTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc               <  1:5952389/76‑1 (MQ=255)
                                                          gTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa              <  1:5638821/76‑1 (MQ=255)
                                                          gTATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa              <  1:3197972/76‑1 (MQ=255)
                                                           tATGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg             <  2:5487798/76‑1 (MQ=255)
                                                            aTGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAgg            >  1:2120655/1‑76 (MQ=255)
                                                            aTGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAgg            <  1:1485368/76‑1 (MQ=255)
                                                            aTGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAgg            >  2:735054/1‑76 (MQ=255)
                                                             tGGCCAAAGCCATGCTCTAAAGGTTCAAGGGTCATCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgc           >  2:2515212/1‑76 (MQ=255)
                                                                ccAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg        <  2:6010601/76‑1 (MQ=255)
                                                                ccAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg        <  2:2026371/76‑1 (MQ=255)
                                                                 cAAAGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCgg       <  2:2706305/76‑1 (MQ=255)
                                                                  aaaGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGt      >  2:5540913/1‑76 (MQ=255)
                                                                  aaaGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGt      >  2:525210/1‑76 (MQ=255)
                                                                  aaaGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGt      >  2:4157631/1‑76 (MQ=255)
                                                                    aGCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt    <  2:3804453/76‑1 (MQ=255)
                                                                     gCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt   <  2:2561700/76‑1 (MQ=255)
                                                                     gCCATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt   <  2:6612585/76‑1 (MQ=255)
                                                                      ccATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTa  >  1:353405/1‑76 (MQ=255)
                                                                      ccATGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTa  >  2:1866043/1‑76 (MQ=255)
                                                                         |                                                                        
GGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGTTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTA  >  REL606/3368764‑3368909

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