Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,313,600 T→C V38A (GTA→GCA)  rcsD → RcsD phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,313,6000TC100.0% 218.2 / NA 62V38A (GTA→GCA) rcsDRcsD phosphotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/33);  total (29/33)

AGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGTACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGCGACGACCTC  >  NC_000913/2313464‑2313735
                                                                                                                                        |                                                                                                                                        
aGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGcaca                                                                                                                                       <  2:460864/140‑1 (MQ=255)
     tAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAgcg                                                                                                                                  <  2:529148/140‑1 (MQ=255)
          ggTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGtt                                                                                                                             >  1:638856/1‑140 (MQ=255)
            tCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAgc                                                                                                                                   <  2:387989/132‑1 (MQ=255)
            tCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAgc                                                                                                                                   >  1:387989/1‑132 (MQ=255)
            tCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTaa                                                                                                                           >  1:617937/1‑140 (MQ=255)
            tCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTaa                                                                                                                           >  1:75452/1‑140 (MQ=255)
            tCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTaa                                                                                                                           >  2:533799/1‑140 (MQ=255)
               tAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGc                                                                                                                        <  2:242123/140‑1 (MQ=255)
                 aaCACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCt                                                                                                                      <  2:463036/140‑1 (MQ=255)
                      gAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAgcg                                                                                                                                  <  2:268977/123‑1 (MQ=255)
                      gAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAgcg                                                                                                                                  >  1:268977/1‑123 (MQ=255)
                      gAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTg                                                                                                                 >  1:303908/1‑140 (MQ=255)
                      gAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTg                                                                                                                 <  2:574783/140‑1 (MQ=255)
                        aTGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGaa                                                                                                               >  1:356927/1‑140 (MQ=255)
                             tCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATa                                                                                                          >  2:41460/1‑140 (MQ=255)
                                    gagaCAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCATAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTa                                                                                                   <  2:733449/140‑1 (MQ=255)
                                         aaCGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGa                                                                                              <  2:356927/140‑1 (MQ=255)
                                            ggCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTg                                                                                           <  1:553206/140‑1 (MQ=255)
                                            ggCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTg                                                                                           <  1:821763/140‑1 (MQ=255)
                                            ggCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTg                                                                                           <  2:608586/140‑1 (MQ=255)
                                             gCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGt                                                                                          <  2:617937/140‑1 (MQ=255)
                                               cACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCg                                                                                        >  2:646255/1‑140 (MQ=255)
                                                 cGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGAc                                                                                      <  1:164437/140‑1 (MQ=255)
                                                   aCCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACAt                                                                                    <  1:653393/140‑1 (MQ=255)
                                                      cGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTAc                                                                                 >  1:27631/1‑140 (MQ=255)
                                                       gCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTAcc                                                                                <  1:287261/140‑1 (MQ=255)
                                                           ttCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATg                                                                            <  1:560823/140‑1 (MQ=255)
                                                                     ccGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCaaaa                                                                 >  1:287292/1‑140 (MQ=255)
                                                                       ggggAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAgc                                                               <  1:41460/140‑1 (MQ=255)
                                                                        gggAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAgcg                                                              <  2:287292/140‑1 (MQ=255)
                                                                              aTTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCATATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGa                                                        >  2:644524/1‑140 (MQ=255)
                                                                               ttACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAgc                                                               <  2:79738/132‑1 (MQ=255)
                                                                               ttACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAgc                                                               >  1:79738/1‑132 (MQ=255)
                                                                                tACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATa                                                      >  2:376842/1‑140 (MQ=255)
                                                                                 aCCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATaa                                                     <  2:638856/140‑1 (MQ=255)
                                                                                     gCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTgg                                                 <  1:384054/140‑1 (MQ=255)
                                                                                       ttcttTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCCAAAAGCGCGTCGATAACTGGc                                                >  1:736408/1‑140 (MQ=255)
                                                                                       ttcttTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCg                                               >  1:146059/1‑140 (MQ=255)
                                                                                        tcttTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGt                                              <  2:443503/140‑1 (MQ=255)
                                                                                           ttttACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTAc                                           >  1:500311/1‑140 (MQ=255)
                                                                                               aCTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGa                                       <  1:460353/140‑1 (MQ=255)
                                                                                               aCTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGa                                       <  1:452054/140‑1 (MQ=255)
                                                                                                 tGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGAcc                                     >  1:703676/1‑140 (MQ=255)
                                                                                                   ttGATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTg                                   >  1:730592/1‑140 (MQ=255)
                                                                                                     gATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGc                                 <  2:730592/140‑1 (MQ=255)
                                                                                                     gATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGc                                 <  2:27631/140‑1 (MQ=255)
                                                                                                     gATCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGc                                 <  1:376842/140‑1 (MQ=255)
                                                                                                      aTCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGa                                       <  2:704878/133‑1 (MQ=255)
                                                                                                      aTCATTGTGTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGa                                       >  1:704878/1‑133 (MQ=255)
                                                                                                           tgtgTTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCt                           >  1:118962/1‑140 (MQ=255)
                                                                                                             tgtTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTAc                         >  2:630076/1‑140 (MQ=255)
                                                                                                             tgtTACTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTAc                         >  1:128677/1‑140 (MQ=255)
                                                                                                                  cTGGTGACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGacaa                    >  1:271514/1‑140 (MQ=255)
                                                                                                                       gACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGc                                 >  1:137062/1‑122 (MQ=255)
                                                                                                                       gACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGc                                 <  2:137062/122‑1 (MQ=255)
                                                                                                                       gACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTg               >  1:457901/1‑140 (MQ=255)
                                                                                                                       gACGATGGGTGTAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTg               <  2:703676/140‑1 (MQ=255)
                                                                                                                                tgtAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGcgacga      <  1:163738/140‑1 (MQ=255)
                                                                                                                                  tAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCATATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGCGACGAcc    <  1:644524/140‑1 (MQ=255)
                                                                                                                                  tAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGCGACGAcc    <  2:271514/140‑1 (MQ=255)
                                                                                                                                  tAATGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGCGACGAcc    >  1:795251/1‑140 (MQ=255)
                                                                                                                                    aTGGCACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGCGACGACCTc  >  2:794361/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
AGTTTTAACAGGTCATAAACACGAATGCGTCAGAAAGAGACAACGGCCACGACCCGCTTTTCACTCCTACCGGGGAGCATTACCCGCTTCTTTTTACTGTTGATCATTGTGTTACTGGTGACGATGGGTGTAATGGTACAAAGCGCCGTTAACGCCTGGCTGAAAGATAAAAGTTACCAGATTGTCGACATTACCCATGCTATCCAAAAGCGCGTCGATAACTGGCGTTACGTGACCTGGCAGATCTACGACAACATTGCCGCGACGACCTC  >  NC_000913/2313464‑2313735

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