Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,650,490 C→T V600M (GTG→ATG)  yfhM ← bacterial alpha2‑macroglobulin colonization factor ECAM, anti‑host protease defense factor, periplasmic inner membrane‑anchored lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,650,4900CT100.0% 229.6 / NA 70V600M (GTG→ATG) yfhMbacterial alpha2‑macroglobulin colonization factor ECAM, anti‑host protease defense factor, periplasmic inner membrane‑anchored lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/37);  total (33/37)

CCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATC  >  NC_000913/2650360‑2650619
                                                                                                                                  |                                                                                                                                 
ccTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGAc                                                                                                                           <  2:69078/139‑1 (MQ=255)
  ttCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCg                                                                                                                         >  2:388642/1‑139 (MQ=255)
    catcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGcc                                                                                                                       >  1:1445124/1‑139 (MQ=255)
     atcaACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGccc                                                                                                                      >  2:19341/1‑139 (MQ=255)
          caaTGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTc                                                                                                                 <  1:928876/139‑1 (MQ=255)
             tGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTc                                                                                                                 <  2:154717/135‑1 (MQ=255)
             tGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTc                                                                                                                 >  1:154717/1‑135 (MQ=255)
               ggCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCa                                                                                                            >  1:140570/1‑139 (MQ=255)
               ggCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCa                                                                                                            <  1:370826/139‑1 (MQ=255)
               ggCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCa                                                                                                            >  2:569703/1‑139 (MQ=255)
                 cTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCaaa                                                                                                          <  1:819861/139‑1 (MQ=255)
                    tttCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTa                                                                                                       >  2:1241030/1‑139 (MQ=255)
                    tttCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTa                                                                                                       >  1:318256/1‑139 (MQ=255)
                        aCCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCt                                                                                                   <  2:588215/139‑1 (MQ=255)
                         ccGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTg                                                                                                  >  1:10664/1‑139 (MQ=255)
                            tGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaa                                                                                               >  2:58066/1‑139 (MQ=255)
                              cTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaaa                                                                                             >  2:435354/1‑139 (MQ=255)
                                gTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATa                                                                                           <  2:140570/139‑1 (MQ=255)
                                  cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGAGCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                         <  2:1445124/139‑1 (MQ=255)
                                  cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                         <  1:1028941/139‑1 (MQ=255)
                                  cAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                         <  2:318256/139‑1 (MQ=255)
                                           aaTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATgg                                                                                >  1:856182/1‑139 (MQ=255)
                                            aTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                               <  1:1474123/139‑1 (MQ=255)
                                            aTCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                               <  2:10664/139‑1 (MQ=255)
                                              cGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGa                                                                             <  2:1390554/139‑1 (MQ=255)
                                               gTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGaa                                                                            <  1:1411796/139‑1 (MQ=255)
                                                tAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGTTTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAAt                                                                           >  1:1298162/1‑139 (MQ=255)
                                                   aCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCg                                                                        <  2:103172/139‑1 (MQ=255)
                                                   aCCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCg                                                                        <  1:58066/139‑1 (MQ=255)
                                                       aTTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGt                                                                                     >  2:336618/1‑122 (MQ=255)
                                                       aTTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGt                                                                                     <  1:336618/122‑1 (MQ=255)
                                                       aTTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttc                                                                    >  1:301893/1‑139 (MQ=255)
                                                           cGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTAcc                                                                                                     >  1:1168786/1‑102 (MQ=255)
                                                           cGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTAcc                                                                                                     <  2:1168786/102‑1 (MQ=255)
                                                              ggCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCAc                                                             <  1:388642/139‑1 (MQ=255)
                                                                  aaCTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGc                                                         <  2:301893/139‑1 (MQ=255)
                                                                               tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCa                                                                                                               >  1:686524/1‑72 (MQ=255)
                                                                               tCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCa                                                                                                               <  2:686524/72‑1 (MQ=255)
                                                                                   cGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCATCAAACTACCCTGGaaa                                                                                              >  1:936158/1‑85 (MQ=255)
                                                                                   cGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGaaa                                                                                              <  2:936158/85‑1 (MQ=255)
                                                                                        aTGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCa                                                                                                            >  2:833958/1‑66 (MQ=255)
                                                                                        aTGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCa                                                                                                            <  1:833958/66‑1 (MQ=255)
                                                                                        aTGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttc                                                                    >  1:246340/1‑106 (MQ=255)
                                                                                        aTGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTttc                                                                    <  2:246340/106‑1 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCaaa                                                                                                          <  2:508066/67‑1 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCaaa                                                                                                          >  1:508066/1‑67 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTa                                                                                                       >  2:340349/1‑70 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTa                                                                                                       <  1:340349/70‑1 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                         <  1:412583/84‑1 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAAc                                                                                         >  2:412583/1‑84 (MQ=255)
                                                                                         tGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCgcgc                                  >  1:214125/1‑139 (MQ=255)
                                                                                          gCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGcc                                 <  2:1327913/139‑1 (MQ=255)
                                                                                                      cAGATAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCa                     >  1:818086/1‑139 (MQ=255)
                                                                                                          tAGCCTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGt                 >  1:880315/1‑139 (MQ=255)
                                                                                                              cTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGAt                                                                                                                   >  1:596379/1‑37 (MQ=255)
                                                                                                              cTGCTCAGCGCGGCGCGTCATCGGGCGACCGCCCGAt                                                                                                                   <  2:596379/37‑1 (MQ=255)
                                                                                                                 cTCAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGa          <  1:19341/139‑1 (MQ=255)
                                                                                                                   cAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                               >  1:143776/1‑68 (MQ=255)
                                                                                                                   cAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATggg                                                                               <  2:143776/68‑1 (MQ=255)
                                                                                                                   cAGCGCGGCGCGTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAAc        <  1:569703/139‑1 (MQ=255)
                                                                                                                      cgcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCtt                                                       <  1:620100/89‑1 (MQ=255)
                                                                                                                      cgcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCtt                                                       >  2:620100/1‑89 (MQ=255)
                                                                                                                      cgcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCgcgc                                  >  2:347560/1‑110 (MQ=255)
                                                                                                                      cgcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCgcgc                                  <  1:347560/110‑1 (MQ=255)
                                                                                                                      cgcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTc     >  2:581868/1‑139 (MQ=255)
                                                                                                                       gcggcgcgTCATCGGGCGACCGCCCGATTCCAGCaa                                                                                                           <  1:1252573/36‑1 (MQ=255)
                                                                                                                       gcggcgcgTCATCGGGCGACCGCCCGATTCCAGCaa                                                                                                           >  2:1252573/1‑36 (MQ=255)
                                                                                                                       gcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTatcatc                       <  1:192377/120‑1 (MQ=255)
                                                                                                                       gcggcgcgTCATCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTatcatc                       >  2:192377/1‑120 (MQ=255)
                                                                                                                         ggcgcgTCATCGGGCGACCGCCCGTTTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATc  <  2:1298162/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
CCTTCATCAACAATGGGCTGTTTCACCGTGCTGTCAGTACGATAATCGTAAACCGATTTCGAGGCGAACTGCGGACGGATCCCCGGCAATGCATCGGCAGGCCAGATAGCCTGCTCAGCGCGGCGCGTCACCGGGCGACCGCCCGATTCCAGCAAACTACCCTGGAAAATAACCTGTAATGGGGAATGCGTTTCCTTCCACTGGCTTTCTGTAGAAACTTCGCCGCGCCCTTTATCATCCAGCGTCAACTGAACTTCATC  >  NC_000913/2650360‑2650619

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