Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,762,213 T→G Y67S (TAT→TCT)  rpoS ← RNA polymerase sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,762,2130TG99.2% 1019.2 / ‑3.4 238Y67S (TAT→TCT) rpoSRNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (136/100);  total (138/100)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCT  >  REL606/2762140‑2762285
                                                                         |                                                                        
gcgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTgaagaa                                                                        <  1:3526046/76‑1 (MQ=255)
 cgAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAAc                                                                       >  2:6027470/1‑76 (MQ=255)
 cgAGAGGCGACATCTCCACGCAGGGCGTGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGCGGTGAAGAAc                                                                       >  1:7152923/1‑76 (MQ=255)
  gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAAcc                                                                      <  2:2782214/76‑1 (MQ=255)
  gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAAcc                                                                      <  1:6216609/76‑1 (MQ=255)
   agagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCa                                                                     >  2:5755905/1‑76 (MQ=255)
    gagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCaa                                                                    >  1:5664142/1‑76 (MQ=255)
    gagGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCaa                                                                    >  1:2558758/1‑76 (MQ=255)
     agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAAt                                                                   >  1:5225506/1‑76 (MQ=255)
     agGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAAt                                                                   <  2:4848616/76‑1 (MQ=255)
       gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAAtct                                                                 <  2:4379516/76‑1 (MQ=255)
       gCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAAtct                                                                 <  1:3868156/76‑1 (MQ=255)
        cGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAAtctc                                                                <  1:2226240/76‑1 (MQ=255)
        cGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAAtctc                                                                >  1:3371453/1‑76 (MQ=255)
        cGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAAtctc                                                                >  1:2313019/1‑76 (MQ=255)
         gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCa                                                               <  1:3578036/76‑1 (MQ=255)
         gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCa                                                               <  1:5007641/76‑1 (MQ=255)
         gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCa                                                               <  1:2287135/76‑1 (MQ=255)
           cATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCAca                                                             >  2:5448822/1‑75 (MQ=255)
            aTCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCa                                                            <  1:6636731/76‑1 (MQ=255)
               tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAgg                                                         <  2:2926062/76‑1 (MQ=255)
               tcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAgg                                                         >  2:7393813/1‑76 (MQ=255)
                ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGt                                                        >  2:6938358/1‑76 (MQ=255)
                 cACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTa                                                       <  2:2351880/76‑1 (MQ=255)
                  aCGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTaa                                                      <  2:1131120/76‑1 (MQ=255)
                   cGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTaaa                                                     <  1:6724648/76‑1 (MQ=255)
                   cGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTaaa                                                     >  1:5376181/1‑76 (MQ=255)
                    gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTCAAg                                                    >  2:7269783/1‑76 (MQ=255)
                    gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAg                                                    >  1:4910967/1‑76 (MQ=255)
                    gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAg                                                    >  2:822081/1‑76 (MQ=255)
                     cAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGc                                                   >  2:3315463/1‑76 (MQ=255)
                     cAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGc                                                   >  2:1190058/1‑76 (MQ=255)
                     cAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGCAAAGc                                                   >  2:7372310/1‑76 (MQ=255)
                      aGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCt                                                  <  1:3135283/76‑1 (MQ=255)
                       gTGCGCGACGCGCAAAATAAACTTCTTCTTCGGTCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTg                                                 <  2:3037763/76‑1 (MQ=255)
                       gTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTg                                                 <  1:1725711/76‑1 (MQ=255)
                       gTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTg                                                 <  2:5780837/76‑1 (MQ=255)
                       gTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTg                                                 <  2:3433203/76‑1 (MQ=255)
                          cgcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGt                                              >  2:5141553/1‑76 (MQ=255)
                          cgcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGt                                              <  1:2798545/76‑1 (MQ=255)
                           gcgACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTc                                             <  1:4173905/76‑1 (MQ=255)
                              aCGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTcgcg                                          <  1:1410243/76‑1 (MQ=255)
                              aCGCGCAAAATAAACTTCTTCTCCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTcgcg                                          <  1:4557648/76‑1 (MQ=255)
                                gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTc                                        >  1:1094695/1‑76 (MQ=255)
                                gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTc                                        >  1:6806813/1‑76 (MQ=255)
                                  gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCa                                      >  2:5800030/1‑76 (MQ=255)
                                    aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAac                                    >  2:366074/1‑76 (MQ=255)
                                    aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAac                                    >  2:7078107/1‑76 (MQ=255)
                                    aaaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAac                                    >  1:2887162/1‑76 (MQ=255)
                                     aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAaca                                   >  2:3852154/1‑76 (MQ=255)
                                     aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAaca                                   <  1:3919470/76‑1 (MQ=255)
                                      aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacac                                  >  2:3291613/1‑76 (MQ=255)
                                      aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacac                                  <  1:5145496/76‑1 (MQ=255)
                                      aaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAacac                                  <  1:7045278/76‑1 (MQ=255)
                                         aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                               <  2:4169055/76‑1 (MQ=255)
                                         aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                               >  2:3153051/1‑76 (MQ=255)
                                         aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACg                               >  2:6820030/1‑76 (MQ=255)
                                          aaCTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGc                              <  2:1653097/76‑1 (MQ=255)
                                           aCTTCTTCTTCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCt                             >  2:6174019/1‑76 (MQ=255)
                                            cttcttcttCGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtg                            >  1:4503186/1‑76 (MQ=255)
                                              tcttcttcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCtgtg                          <  2:1440025/76‑1 (MQ=255)
                                                  cttcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGc                      >  2:4318006/1‑76 (MQ=255)
                                                  cttcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGc                      >  1:3272872/1‑76 (MQ=255)
                                                  cttcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGc                      >  1:5813595/1‑76 (MQ=255)
                                                   ttcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCt                     >  2:3552839/1‑76 (MQ=255)
                                                   ttcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCt                     <  2:6358138/76‑1 (MQ=255)
                                                    tcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTc                    >  1:1973345/1‑76 (MQ=255)
                                                    tcGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTc                    >  1:1166802/1‑76 (MQ=255)
                                                     cGGCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTcc                   >  1:4317404/1‑76 (MQ=255)
                                                       gCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCt                 >  2:2837899/1‑76 (MQ=255)
                                                       gCCGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCt                 <  2:5900728/76‑1 (MQ=255)
                                                        ccGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTg                <  2:4535836/76‑1 (MQ=255)
                                                         cGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc               <  1:4759415/76‑1 (MQ=255)
                                                         cGTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGc               >  2:6959972/1‑76 (MQ=255)
                                                          gTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCttgg              >  2:3738862/1‑73 (MQ=255)
                                                          gTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg              >  1:3198749/1‑76 (MQ=255)
                                                          gTTAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCg              <  1:1015854/76‑1 (MQ=255)
                                                           ttAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa             >  1:1785966/1‑76 (MQ=255)
                                                           ttAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa             >  1:5550111/1‑76 (MQ=255)
                                                           ttAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGa             >  1:2133143/1‑76 (MQ=255)
                                                            tAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt            >  2:2230044/1‑76 (MQ=255)
                                                            tAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt            >  2:2622288/1‑76 (MQ=255)
                                                            tAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGAt            >  2:1281423/1‑76 (MQ=255)
                                                            tAACAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGACCAACACACGCTGTGTGGCTCCCTGCGAt            <  2:7505290/76‑1 (MQ=255)
                                                             aaCAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATa           >  2:1716306/1‑76 (MQ=255)
                                                             aaCAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATa           >  2:6723273/1‑76 (MQ=255)
                                                             aaCAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATa           >  2:4791299/1‑76 (MQ=255)
                                                             aaCAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATa           <  1:5949366/76‑1 (MQ=255)
                                                             aaCAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATa           <  1:2635876/76‑1 (MQ=255)
                                                              aCAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATaa          <  2:2522535/76‑1 (MQ=255)
                                                               cAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc         <  1:6819739/76‑1 (MQ=255)
                                                               cAGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAAc         <  1:2961920/76‑1 (MQ=255)
                                                                aGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACa        <  1:1031929/76‑1 (MQ=255)
                                                                aGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACa        <  1:193221/76‑1 (MQ=255)
                                                                aGTGGTGAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACa        <  2:3215832/76‑1 (MQ=255)
                                                                 gtggtgAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAg       <  1:871294/76‑1 (MQ=255)
                                                                  tggtgAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGt      >  2:1997548/1‑76 (MQ=255)
                                                                  tggtgAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGt      <  1:938621/76‑1 (MQ=255)
                                                                   ggtgAAGAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGtt     <  2:246654/76‑1 (MQ=255)
                                                                      gaagaaCCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCct  <  2:3919538/76‑1 (MQ=255)
                                                                         |                                                                        
GCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCT  >  REL606/2762140‑2762285

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