Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,522,146 T→A D650V (GAT→GTT)  yjhG ← putative dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,522,1460TA100.0% 228.5 / NA 59D650V (GAT→GTT) yjhGputative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (33/26);  total (33/26)

CGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATT  >  NC_000913/4522005‑4522291
                                                                                                                                             |                                                                                                                                                 
ctaTAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCaa                                                                                                                                            >  2:519105/3‑149 (MQ=255)
cGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCaa                                                                                                                                            >  1:902985/1‑149 (MQ=255)
  aTAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGcc                                                                                                                                        <  1:631316/151‑1 (MQ=255)
     aTTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCcgcg                                                                                                                                     >  2:660206/1‑151 (MQ=255)
      ttAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCgcgc                                                                                                                                    >  2:51986/1‑151 (MQ=255)
          tGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAAt                                                                                                                                >  1:854013/1‑151 (MQ=255)
              aTAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAAtttt                                                                                                                             >  1:536210/1‑150 (MQ=255)
                     ccTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCTGCGCCAATTTTGTTTAcat                                                                                                                     >  2:260994/1‑151 (MQ=255)
                     ccTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTAca                                                                                                                      <  1:742188/150‑1 (MQ=255)
                       tctcATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATGAAAACGCGCAAAGCCGCGCCAATTTTGTTTAcatc                                                                                                                    >  1:50483/1‑150 (MQ=255)
                         tcATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTAcatcat                                                                                                                  >  2:707932/1‑150 (MQ=255)
                          cATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTAcatca                                                                                                                   <  2:954116/148‑1 (MQ=255)
                             aaCAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAAt                                                                                                              <  1:946160/150‑1 (MQ=255)
                             aaCAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATa                                                                                                             <  1:248805/151‑1 (MQ=255)
                                     aGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGt                                                                                                     >  2:376863/1‑151 (MQ=255)
                                       cGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTcc                                                                                                   >  1:559392/1‑151 (MQ=255)
                                        gACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCa                                                                                                  <  2:546969/151‑1 (MQ=255)
                                         aCCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCAt                                                                                                 >  2:363309/1‑151 (MQ=255)
                                          cccATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATg                                                                                                <  2:503648/151‑1 (MQ=255)
                                            catcatTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTc                                                                                              >  2:278514/1‑151 (MQ=255)
                                                atTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCccac                                                                                          >  1:66452/1‑151 (MQ=255)
                                                          tGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGc                                                                                >  2:210461/1‑151 (MQ=255)
                                                            aGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCt                                                                              <  2:451661/151‑1 (MQ=255)
                                                             gCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTg                                                                             <  1:753963/151‑1 (MQ=255)
                                                                       aaTGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGcc                                                                   <  2:359715/151‑1 (MQ=255)
                                                                           atgCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCAt                                                                >  2:782569/1‑150 (MQ=255)
                                                                           atgCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATa                                                               <  2:258604/151‑1 (MQ=255)
                                                                             gCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAg                                                              <  1:812980/150‑1 (MQ=255)
                                                                             gCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAg                                                              <  1:97164/150‑1 (MQ=255)
                                                                               agagTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCg                                                           >  2:310661/1‑151 (MQ=255)
                                                                                 agTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCggg                                                         >  2:196451/1‑151 (MQ=255)
                                                                                    cGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTat                                                      <  2:74643/151‑1 (MQ=255)
                                                                                         aaTTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCt                                                 <  2:450980/151‑1 (MQ=255)
                                                                                          aTTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTg                                                <  1:850979/151‑1 (MQ=255)
                                                                                              aGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCaa                                            <  2:891749/151‑1 (MQ=255)
                                                                                                    cTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAgg                                      >  1:17336/1‑151 (MQ=255)
                                                                                                         tGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGAATCATCTGGCAATTCAGGATCg                                  >  1:162322/1‑150 (MQ=255)
                                                                                                            ttttATTTATGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGa                               >  2:569101/1‑150 (MQ=255)
                                                                                                             tttatttaTGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGa                               >  1:826861/1‑149 (MQ=255)
                                                                                                               tatttatGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGt                            >  2:317276/1‑150 (MQ=255)
                                                                                                                    atGCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATcc                      <  2:668348/151‑1 (MQ=255)
                                                                                                                      gCTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTg                    <  1:660168/151‑1 (MQ=255)
                                                                                                                       cTTTCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCACACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGTTCGGGAAGTAAATCATg                    >  1:314136/1‑150 (MQ=255)
                                                                                                                         ttCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAt                  <  2:311000/150‑1 (MQ=255)
                                                                                                                         ttCAGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATg                 >  2:896808/1‑151 (MQ=255)
                                                                                                                            aGTTTTTATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCt              <  2:386958/151‑1 (MQ=255)
                                                                                                                                 ttATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGt          <  1:740285/150‑1 (MQ=255)
                                                                                                                                 ttATTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTc         >  2:525204/1‑151 (MQ=255)
                                                                                                                                   aTTCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTg       >  2:510793/1‑151 (MQ=255)
                                                                                                                                    ttCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATAAAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCCTCTGGCAATTCAGGATTGGGCAGCAAATCCTGATGGCTGGATCTgg      >  1:174718/1‑151 (MQ=255)
                                                                                                                                     tCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGc     >  2:814165/1‑151 (MQ=255)
                                                                                                                                     tCATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGc     >  2:895921/1‑151 (MQ=255)
                                                                                                                                      cATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGc     <  2:738845/150‑1 (MQ=255)
                                                                                                                                      cATAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCa    <  1:196440/151‑1 (MQ=255)
                                                                                                                                       aTAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCAt   <  1:544413/151‑1 (MQ=255)
                                                                                                                                       aTAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCAt   >  2:884471/1‑151 (MQ=255)
                                                                                                                                        tAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCAtt  >  2:209582/1‑151 (MQ=255)
                                                                                                                                        tAAAAACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCAtt  >  1:730994/1‑151 (MQ=255)
                                                                                                                                         aaaaACGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCAtt  <  2:880819/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
CGATAATTAGTGGCATAACGTCCTCTCATAACAGGGTAGCGACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATT  >  NC_000913/4522005‑4522291

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