Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,665,120 G→A *205* (TAG→TAA dmsD → twin‑argninine leader‑binding protein for DmsA and TorA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,665,1200GA100.0% 249.7 / NA 76*205* (TAG→TAAdmsDtwin‑argninine leader‑binding protein for DmsA and TorA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/38);  total (38/38)

CACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGC  >  NC_000913/1664984‑1665241
                                                                                                                                        |                                                                                                                         
cacGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGa                                                                                                                         <  2:1780323/139‑1 (MQ=255)
 acGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                        <  2:1845523/139‑1 (MQ=255)
 acGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                        <  1:1691746/139‑1 (MQ=255)
        cTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTa                                                                                                                 <  2:696866/139‑1 (MQ=255)
        cTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTa                                                                                                                 >  1:1542373/1‑139 (MQ=255)
                         aaaaGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                >  2:91098/1‑139 (MQ=255)
                                        tttCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCaa                                                                                   <  2:755605/137‑1 (MQ=255)
                                        tttCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCaa                                                                                   >  1:755605/1‑137 (MQ=255)
                                             aCCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCa                                                                            >  1:502932/1‑139 (MQ=255)
                                              ccGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                <  1:1029565/118‑1 (MQ=255)
                                              ccGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                >  2:1029565/1‑118 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGa                                                                                                                         >  2:974928/1‑91 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGa                                                                                                                         <  1:974928/91‑1 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGa                                                                                                      <  2:1844490/110‑1 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGa                                                                                                      >  1:1844490/1‑110 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAAt                                                                                  >  2:449088/1‑130 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAAt                                                                                  <  1:449088/130‑1 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                         >  2:1672866/1‑139 (MQ=255)
                                                gTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaata                                                                         <  1:1031792/139‑1 (MQ=255)
                                                  gCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                        >  1:370478/1‑90 (MQ=255)
                                                  gCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                        <  2:370478/90‑1 (MQ=255)
                                                  gCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCagtaat                                                                       <  2:1542373/139‑1 (MQ=255)
                                                  gCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaataat                                                                       <  2:322886/139‑1 (MQ=255)
                                                        ggTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGTGCAATCCGGCAATAATAGGTTa                                                                 >  1:544372/1‑139 (MQ=255)
                                                        ggTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTa                                                                 >  1:737443/1‑139 (MQ=255)
                                                         gTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTACTTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                >  2:1651434/1‑107 (MQ=255)
                                                         gTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                <  1:1651434/107‑1 (MQ=255)
                                                           gAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAata                                                                                                              <  2:1689974/91‑1 (MQ=255)
                                                           gAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAata                                                                                                              >  1:1689974/1‑91 (MQ=255)
                                                             aCTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaataa                                                                        >  1:239695/1‑127 (MQ=255)
                                                             aCTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaataa                                                                        <  2:239695/127‑1 (MQ=255)
                                                             aCTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAgtg                                                            >  1:1069400/1‑139 (MQ=255)
                                                             aCTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAgtg                                                            >  2:233734/1‑139 (MQ=255)
                                                                ggCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGc                                                                             <  1:247719/119‑1 (MQ=255)
                                                                ggCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGc                                                                             >  2:247719/1‑119 (MQ=255)
                                                                  cTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACg                                                       >  2:1380232/1‑139 (MQ=255)
                                                                   tCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                        >  1:29934/1‑73 (MQ=255)
                                                                   tCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGAt                                                                                                                        <  2:29934/73‑1 (MQ=255)
                                                                         aaCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                <  2:502932/139‑1 (MQ=255)
                                                                         aaCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                <  1:233734/139‑1 (MQ=255)
                                                                         aaCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                <  1:1380232/139‑1 (MQ=255)
                                                                           cGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTAt                                              <  2:1069400/139‑1 (MQ=255)
                                                                           cGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTAt                                              <  1:1672866/139‑1 (MQ=255)
                                                                                gcgcAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                >  2:97766/1‑84 (MQ=255)
                                                                                gcgcAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                <  1:97766/84‑1 (MQ=255)
                                                                                gcgcAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAc                                                        <  1:531141/124‑1 (MQ=255)
                                                                                gcgcAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCAc                                                        >  2:531141/1‑124 (MQ=255)
                                                                                 cgcAATGGCAATCACACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                              <  1:467495/100‑1 (MQ=255)
                                                                                 cgcAATGGCAATCACACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                              >  2:467495/1‑100 (MQ=255)
                                                                                              acaACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAataata                                                                      >  1:1990251/1‑96 (MQ=255)
                                                                                              acaACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAataata                                                                      <  2:1990251/96‑1 (MQ=255)
                                                                                                     ttAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTa                                               >  2:1609516/1‑112 (MQ=255)
                                                                                                     ttAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTa                                               <  1:1609516/112‑1 (MQ=255)
                                                                                                            cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                           <  1:1603228/77‑1 (MQ=255)
                                                                                                            cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCaa                                                                           >  2:1603228/1‑77 (MQ=255)
                                                                                                            cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAgcacg                                 >  1:1542961/1‑119 (MQ=255)
                                                                                                            cTGTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAgcacg                                 <  2:1542961/119‑1 (MQ=255)
                                                                                                              gTCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc           >  2:353828/1‑139 (MQ=255)
                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTAcg                                                                                        <  1:832059/61‑1 (MQ=255)
                                                                                                               tCGCGGTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTAcg                                                                                        >  2:832059/1‑61 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTAcgcg                                                                                      <  1:8378/58‑1 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTAcgcg                                                                                      >  2:8378/1‑58 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAAt                                                                                  <  1:1672483/62‑1 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAAt                                                                                  >  2:1672483/1‑62 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                              >  2:1528199/1‑66 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCgg                                                                              <  1:1528199/66‑1 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTa                                               >  2:1864713/1‑97 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTa                                               <  1:1864713/97‑1 (MQ=255)
                                                                                                                    gTTAAACCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc     >  1:348494/1‑139 (MQ=255)
                                                                                                                       aaaCCGCTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc  >  1:1170091/1‑139 (MQ=255)
                                                                                                                             cTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                >  2:741074/1‑87 (MQ=255)
                                                                                                                             cTGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                <  1:741074/87‑1 (MQ=255)
                                                                                                                              tGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                >  1:18085/1‑38 (MQ=255)
                                                                                                                              tGTTTCGATAAGATCACCGTAATATTGCCGGATGGTGa                                                                                                <  2:18085/38‑1 (MQ=255)
                                                                                                                               gTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGc                                                                             <  2:689361/56‑1 (MQ=255)
                                                                                                                               gTTTCGATAAGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGc                                                                             >  1:689361/1‑56 (MQ=255)
                                                                                                                                        |                                                                                                                         
CACGTTTTCTTGATGTTTTTATCGAAAAAGCAGAACACCCTTTCTACCGTGCACTGGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGC  >  NC_000913/1664984‑1665241

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