Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,731 C→A W357C (TGG→TGT mrdA ← transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606649,7310CA98.8% 691.6 / ‑4.0 171W357C (TGG→TGTmrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (87/82);  total (87/84)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

TTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGT  >  REL606/649658‑649803
                                                                         |                                                                        
ttGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGccaaca                                                                        <  1:4082076/76‑1 (MQ=255)
ttGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGccaaca                                                                        >  2:7110256/1‑76 (MQ=255)
 tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAg                                                                       <  1:5590563/76‑1 (MQ=255)
 tgtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAg                                                                       <  1:1806122/76‑1 (MQ=255)
  gtgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGc                                                                      <  1:6132551/76‑1 (MQ=255)
   tgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGcc                                                                     >  1:5845692/1‑76 (MQ=255)
   tgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGcc                                                                     <  1:2122850/76‑1 (MQ=255)
   tgACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGcc                                                                     <  2:486122/76‑1 (MQ=255)
     aCATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTg                                                                   <  2:6224486/76‑1 (MQ=255)
      cATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTgg                                                                  <  2:801707/76‑1 (MQ=255)
       aTTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTggg                                                                 <  2:6251772/76‑1 (MQ=255)
        ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTggat                                                                >  1:2458099/1‑74 (MQ=255)
        ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGt                                                                <  1:4478952/76‑1 (MQ=255)
        ttCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGt                                                                <  2:531971/76‑1 (MQ=255)
         tCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTc                                                               <  2:1296356/76‑1 (MQ=255)
         tCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTc                                                               <  1:2104767/76‑1 (MQ=255)
            gACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCaaa                                                            <  2:2901668/76‑1 (MQ=255)
             aCGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAAc                                                           <  1:6213969/76‑1 (MQ=255)
               gCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAg                                                         >  2:4571150/1‑76 (MQ=255)
               gCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAg                                                         <  2:903324/76‑1 (MQ=255)
                cccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGc                                                        >  1:5294841/1‑76 (MQ=255)
                cccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGc                                                        <  1:1092467/76‑1 (MQ=255)
                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGcg                                                       <  2:1958210/76‑1 (MQ=255)
                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGcg                                                       >  2:1483359/1‑76 (MQ=255)
                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGcg                                                       <  1:55389/76‑1 (MQ=255)
                 ccGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGcg                                                       <  2:7680640/76‑1 (MQ=255)
                   gTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGcgtc                                                     <  2:363012/76‑1 (MQ=255)
                     ggCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGcgtcgt                                                   >  1:7238811/1‑76 (MQ=255)
                        cccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTAtt                                                <  1:4360988/76‑1 (MQ=255)
                        cccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTAtt                                                <  1:4014129/76‑1 (MQ=255)
                         ccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTg                                               >  1:3976976/1‑76 (MQ=255)
                         ccATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTg                                               <  2:4366473/76‑1 (MQ=255)
                          cATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTgc                                              <  2:977711/76‑1 (MQ=255)
                           aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTgcg                                             <  1:4236192/76‑1 (MQ=255)
                           aTTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTgcg                                             >  2:1583680/1‑76 (MQ=255)
                             tttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGcgcg                                           <  1:11629/76‑1 (MQ=255)
                              ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGt                                          >  1:3968325/1‑76 (MQ=255)
                              ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGt                                          >  1:3014946/1‑76 (MQ=255)
                              ttttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGt                                          >  1:7891184/1‑76 (MQ=255)
                               tttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTg                                         >  1:7408021/1‑76 (MQ=255)
                                ttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGa                                        >  2:6887519/1‑76 (MQ=255)
                                ttCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGa                                        >  2:7479153/1‑76 (MQ=255)
                                 tCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGAt                                       >  1:2836425/1‑76 (MQ=255)
                                 tCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGAt                                       >  1:3364668/1‑76 (MQ=255)
                                  ccAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATc                                      >  1:664232/1‑76 (MQ=255)
                                  ccAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATc                                      >  1:3836082/1‑76 (MQ=255)
                                  ccAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATc                                      >  1:995802/1‑76 (MQ=255)
                                   cAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCa                                     <  2:8144483/76‑1 (MQ=255)
                                      tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCAccc                                  >  2:7196491/1‑76 (MQ=255)
                                      tCACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACACCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCAccc                                  <  1:6307402/76‑1 (MQ=255)
                                       cACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCAcccc                                 <  2:1372877/76‑1 (MQ=255)
                                       cACGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGACTGGGCCAAACAGCGTCGTATTGCGCGTGATCAcccc                                 >  1:6121588/1‑76 (MQ=255)
                                        aCGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCg                                <  1:812758/76‑1 (MQ=255)
                                         cGATAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTAATCACCCCgg                               >  2:8012484/1‑76 (MQ=255)
                                            tAACGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGgcgc                            >  1:3177874/1‑76 (MQ=255)
                                             aaCGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCt                           >  1:371003/1‑76 (MQ=255)
                                             aaCGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCt                           >  1:3618363/1‑76 (MQ=255)
                                             aaCGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCt                           >  1:4706571/1‑76 (MQ=255)
                                             aaCGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCt                           >  1:7501460/1‑76 (MQ=255)
                                             aaCGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCt                           >  1:4493836/1‑76 (MQ=255)
                                              aCGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTc                          <  1:3551079/76‑1 (MQ=255)
                                               cGTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCa                         <  2:3307265/76‑1 (MQ=255)
                                                gTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCaa                        >  1:7062253/1‑76 (MQ=255)
                                                gTTTTTCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCaa                        >  2:6775709/1‑76 (MQ=255)
                                                 tttttCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAAt                       <  2:7741221/76‑1 (MQ=255)
                                                 tttttCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCGCAAt                       >  2:2327868/1‑76 (MQ=255)
                                                  ttttCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATg                      <  2:5763732/76‑1 (MQ=255)
                                                   tttCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                     <  1:483544/76‑1 (MQ=255)
                                                   tttCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGc                     >  2:3245063/1‑76 (MQ=255)
                                                    ttCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGcc                    >  1:7958301/1‑76 (MQ=255)
                                                    ttCCGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGcc                    >  1:7060625/1‑76 (MQ=255)
                                                      ccGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGa                  >  2:6350661/1‑76 (MQ=255)
                                                      ccGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGa                  <  1:6918376/76‑1 (MQ=255)
                                                      ccGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGa                  >  1:1706591/1‑76 (MQ=255)
                                                       cGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaa                 >  2:6361484/1‑76 (MQ=255)
                                                       cGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaa                 <  2:3762837/76‑1 (MQ=255)
                                                       cGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaa                 <  2:7566304/76‑1 (MQ=255)
                                                       cGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaa                 <  2:7914132/76‑1 (MQ=255)
                                                       cGAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaa                 <  2:5210619/76‑1 (MQ=255)
                                                        gAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaaa                >  2:5884713/1‑76 (MQ=255)
                                                        gAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaaa                >  1:7209054/1‑76 (MQ=255)
                                                        gAACCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGaaa                >  2:4163023/1‑76 (MQ=255)
                                                         aaCCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAAc               >  1:2182088/1‑76 (MQ=255)
                                                          aCCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAAcc              <  1:5979820/76‑1 (MQ=255)
                                                          aCCTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAAcc              <  2:1619183/76‑1 (MQ=255)
                                                           ccTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCg             >  1:7632775/1‑76 (MQ=255)
                                                           ccTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCg             <  2:7836461/76‑1 (MQ=255)
                                                            cTGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGc            >  1:6621592/1‑76 (MQ=255)
                                                            cTGGCAGTTGCCAACAGCCTGGGGCAAACAGCGTCGTATTGCGCGTGATGACCCCGGCGCTCAATGCCGAAACCGc            >  2:5717870/1‑76 (MQ=255)
                                                             tGGCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGcc           >  2:8068741/1‑76 (MQ=255)
                                                              ggCAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCca          >  2:5679570/1‑76 (MQ=255)
                                                                cAGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCcaca        >  1:464233/1‑76 (MQ=255)
                                                                 aGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt       >  2:3588568/1‑76 (MQ=255)
                                                                 aGTTGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt       >  2:439327/1‑76 (MQ=255)
                                                                 aGTGGCCAACAGCCGGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACAt       <  1:7028906/76‑1 (MQ=255)
                                                                   ttGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAg     >  1:890263/1‑76 (MQ=255)
                                                                    tGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACTGCCACATAgg    >  2:2069983/1‑76 (MQ=255)
                                                                    tGCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAgg    >  1:7673162/1‑76 (MQ=255)
                                                                     gCCAACAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAggg   >  1:173291/1‑76 (MQ=255)
                                                                      ccaacaGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGt  <  1:941715/76‑1 (MQ=255)
                                                                         |                                                                        
TTGTGACATTCAGACGCCCGTGGCCCCATTTTTTCCAGTCACGATAACGTTTTTCCGAACCTGGCAGTTGCCACCAGCCTGGGTCAAACAGCGTCGTATTGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGT  >  REL606/649658‑649803

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