Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,755,639 G→A 100% intergenic (‑498/‑59) ydhZ ← / → pykF uncharacterized protein/pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,755,6390GA100.0% 281.6 / NA 82intergenic (‑498/‑59)ydhZ/pykFuncharacterized protein/pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/48);  total (34/48)

GAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755544‑1755735
                                                                                               |                                                                                                
gagaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGt                                                                                              >  1:508659/1‑100 (MQ=255)
    cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATCCCGGATTAGCTTTCCGACAGTgcgc                                                                                          <  1:438173/100‑1 (MQ=255)
     acaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          <  2:68550/99‑1 (MQ=255)
     acaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGcc                                                                                         <  2:1526641/100‑1 (MQ=255)
      caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                        <  1:1431394/100‑1 (MQ=255)
      caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                        >  2:1955338/1‑100 (MQ=255)
       aaGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCa                                                                                       <  2:1196137/100‑1 (MQ=255)
       aaGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCa                                                                                       <  1:1939818/100‑1 (MQ=255)
           acacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaaa                                                                                   <  2:2042420/100‑1 (MQ=255)
            cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  >  1:1281609/1‑100 (MQ=255)
              caTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCa                                                                                <  2:2216659/100‑1 (MQ=255)
                 tCCTCTGCACGCTTTTTTTTTTTTTTTTTTCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAg                                                                                      >  2:999348/1‑91 (MQ=255)
                  ccTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                            >  1:966753/1‑100 (MQ=255)
                  ccTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTGTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                            <  1:1145422/100‑1 (MQ=255)
                  ccTCTGCACGCTTGTTAGATGTCACCTATCCTTAGAGCTAGGCACCACCACTTTCGTATTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                            <  2:1283730/100‑1 (MQ=255)
                   ctctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGttt                                                                           >  1:1519251/1‑100 (MQ=255)
                       gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCAGCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  2:1654059/100‑1 (MQ=255)
                       gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  1:1638142/100‑1 (MQ=255)
                       gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  2:1905437/100‑1 (MQ=255)
                       gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  1:725561/100‑1 (MQ=255)
                       gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  1:682814/100‑1 (MQ=255)
                       gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       >  1:553360/1‑100 (MQ=255)
                        cACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTccc                                                                      <  2:2290286/100‑1 (MQ=255)
                          cGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCAt                                                                    >  1:1622444/1‑100 (MQ=255)
                           gCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATc                                                                   <  1:999353/100‑1 (MQ=255)
                           gCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATc                                                                   >  1:1170093/1‑100 (MQ=255)
                             tttttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCt                                                                 <  2:1929263/100‑1 (MQ=255)
                               tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               >  1:284706/1‑100 (MQ=255)
                               tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               <  1:1955352/100‑1 (MQ=255)
                               tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               >  2:2172247/1‑100 (MQ=255)
                                   gATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCaa                                                           >  1:836747/1‑100 (MQ=255)
                                    aTGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAAc                                                          <  2:1493096/100‑1 (MQ=255)
                                    aTGTCACCGATCCGTAGAGCGAGGCACCACCACGTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAAc                                                          <  2:2290656/100‑1 (MQ=255)
                                     tGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACt                                                         <  1:159479/100‑1 (MQ=255)
                                      gatctcTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACtt                                                        <  1:1537168/95‑1 (MQ=255)
                                        cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      >  2:227068/1‑100 (MQ=255)
                                        cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      >  2:1671900/1‑100 (MQ=255)
                                        cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      >  2:870419/1‑100 (MQ=255)
                                          ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    <  2:377649/100‑1 (MQ=255)
                                          ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    <  1:382709/100‑1 (MQ=255)
                                           cTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTagatc                                                   >  2:1537159/1‑96 (MQ=255)
                                            tATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagac                                                  <  1:272716/100‑1 (MQ=255)
                                              tCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagacta                                                <  1:1473805/100‑1 (MQ=255)
                                               ccTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaa                                               >  2:1020638/1‑100 (MQ=255)
                                               ccTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaa                                               >  1:793239/1‑100 (MQ=255)
                                                cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGTAAGCAATTTTCTCCCATCCTTCTCAACTTAaagactaag                                              >  2:1654368/1‑100 (MQ=255)
                                                cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                              >  2:115078/1‑100 (MQ=255)
                                                  tAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagac                                            >  1:987247/1‑100 (MQ=255)
                                                     agcttgtccccaccacTTTCGTAATACCGACGTCGTTTTCCGACATTGCGCCCAGAAACCAAGTTTATCCCATCTTTCTCAACTTAAAGACTAAAACTGt                                         <  2:2233458/91‑1 (MQ=255)
                                                      gCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTc                                        >  1:1913591/1‑100 (MQ=255)
                                                       cGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCa                                       <  1:944040/100‑1 (MQ=255)
                                                         aGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATg                                     >  2:1290437/1‑100 (MQ=255)
                                                         aGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATg                                     <  2:404674/100‑1 (MQ=255)
                                                          ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                    <  2:276292/100‑1 (MQ=255)
                                                              ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                                >  1:975566/1‑100 (MQ=255)
                                                               caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAg                               <  1:697667/100‑1 (MQ=255)
                                                                 ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAc                             <  1:561114/100‑1 (MQ=255)
                                                                 ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAc                             <  2:793236/100‑1 (MQ=255)
                                                                 ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAc                             <  1:348141/100‑1 (MQ=255)
                                                                  cacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAcc                            <  1:1460531/100‑1 (MQ=255)
                                                                     tttCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCaaa                         >  2:2030159/1‑100 (MQ=255)
                                                                         gTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTg                     >  1:1458942/1‑100 (MQ=255)
                                                                          tAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGt                    >  1:1956674/1‑100 (MQ=255)
                                                                            aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                  <  1:648234/100‑1 (MQ=255)
                                                                            aTACCGGATTCGCTTTCCGACAGTGCGCCAAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                  <  2:350329/100‑1 (MQ=255)
                                                                             tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 <  2:1993958/100‑1 (MQ=255)
                                                                             tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 <  2:1191656/100‑1 (MQ=255)
                                                                              aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc                >  1:2224731/1‑100 (MQ=255)
                                                                               ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCa               >  1:1673255/1‑100 (MQ=255)
                                                                                cGGATTCTCTTTCCGACATTTCGCCCAGAAAGCAAGGTTCGCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGATAAAGACCAAAATTGTTTGCa               <  2:374228/99‑1 (MQ=255)
                                                                                cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAc              <  1:1636053/100‑1 (MQ=255)
                                                                                    ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAGAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATc          <  2:455679/100‑1 (MQ=255)
                                                                                     tCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAATACCAAAATTGTTTGCACCATCg         <  2:806332/100‑1 (MQ=255)
                                                                                     tCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCg         >  1:590207/1‑100 (MQ=255)
                                                                                     tCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCg         >  1:590195/1‑100 (MQ=255)
                                                                                        cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      >  1:1137091/1‑100 (MQ=255)
                                                                                        cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      <  2:1519242/100‑1 (MQ=255)
                                                                                        cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      >  1:1520315/1‑100 (MQ=255)
                                                                                         tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAATACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc     <  2:1281603/100‑1 (MQ=255)
                                                                                          ttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGTCCAAAATTGTTTGCACCATCGGACCg    <  1:1438592/100‑1 (MQ=255)
                                                                                          ttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCg    <  2:584052/100‑1 (MQ=255)
                                                                                            ccGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGaa  >  2:381561/1‑100 (MQ=255)
                                                                                               |                                                                                                
GAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755544‑1755735

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