BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,583,050	G→A	Q446* (CAA→TAA)	hsdR ←	endonuclease R

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,583,050	0	G	A	98.6%	430.9 / NA	147	Q446* (CAA→TAA)	hsdR	endonuclease R
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (63/82); total (63/84)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

CATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGG  >  NC_000913/4582978‑4583124

cgtCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                         <  2:802120/74‑1 (MQ=255)

cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                         <  2:4502693/76‑1 (MQ=255)

cATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACgg                                                                         >  2:4373107/1‑76 (MQ=255)

 aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                        <  1:1528843/76‑1 (MQ=255)

 aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                        >  1:2630657/1‑76 (MQ=255)

 aTCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGt                                                                        >  1:2652849/1‑76 (MQ=255)

  tCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGtt                                                                       <  1:751874/76‑1 (MQ=255)

    ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCg                                                                     >  2:3456659/1‑76 (MQ=255)

    ttCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCg                                                                     <  1:1932478/76‑1 (MQ=255)

     tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                    <  1:1187762/76‑1 (MQ=255)

     tCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGGTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGc                                                                    >  1:3125647/1‑76 (MQ=255)

      ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                   >  2:699107/1‑76 (MQ=255)

      ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                   <  1:1662770/76‑1 (MQ=255)

      ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                   <  1:4354627/76‑1 (MQ=255)

      ccTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCt                                                                   <  1:4304619/76‑1 (MQ=255)

       atGATACTAACGCAGGCCAAGCTCGCTCAGGCCGGGGTCAGCGGCAGACCACTGATTCTGTTTTTACGGTTCGccg                                                                  >  2:156292/2‑74 (MQ=255)

          ataataACGCAGGCCCAGTTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCg                                                               <  2:2097551/76‑1 (MQ=255)

          ataataACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGcca                                                               <  2:3709972/76‑2 (MQ=255)

           taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGa                                                              <  1:3029162/76‑1 (MQ=255)

           taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGa                                                              <  2:2309204/76‑1 (MQ=255)

           taataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACAGTTCGCTGCCGa                                                              >  1:44414/1‑76 (MQ=255)

             ataaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGc                                                            >  2:802744/1‑76 (MQ=255)

              taaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCa                                                           >  2:4117882/1‑76 (MQ=255)

               aaCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAt                                                          >  1:2570120/1‑76 (MQ=255)

                aCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAtt                                                         <  1:2958112/76‑1 (MQ=255)

                aCGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAtt                                                         >  1:2377437/1‑76 (MQ=255)

                 cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                        <  1:2767931/76‑1 (MQ=255)

                 cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                        <  2:1602953/76‑1 (MQ=255)

                 cGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCAttt                                                        <  2:2556124/76‑1 (MQ=255)

                   gaGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                      <  2:280650/75‑1 (MQ=255)

                   cAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTcc                                                      >  2:2226412/1‑76 (MQ=255)

                     ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcag                                                    >  2:3195888/1‑76 (MQ=255)

                     ggCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTGACGGTTCGCTGCCGAGCATTTCcag                                                    >  1:415576/1‑76 (MQ=255)

                      gCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagc                                                   <  2:3165582/76‑1 (MQ=255)

                       cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagca                                                  <  2:3428129/76‑1 (MQ=255)

                       cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagca                                                  <  1:4758734/76‑1 (MQ=255)

                       cccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagca                                                  <  2:883660/76‑1 (MQ=255)

                        ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                 >  1:2794067/1‑76 (MQ=255)

                        ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcag                                                 <  2:446192/76‑1 (MQ=255)

                        ccAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCcagcac                                                 >  2:396113/1‑75 (MQ=255)

                         cAGCTCGCTCATGCCAGGTTAATCGGCAAACCACTGATTCTGTTTTTACGATTCGCTGCCGAGCATTTCCagcagc                                                <  1:841415/76‑1 (MQ=255)

                         cAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCagcagc                                                >  2:4799555/1‑76 (MQ=255)

                          aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                               >  1:1237858/1‑76 (MQ=255)

                          aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                               >  1:2423422/1‑76 (MQ=255)

                          aGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGct                                               >  2:3823971/1‑76 (MQ=255)

                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                              >  1:1326786/1‑76 (MQ=255)

                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                              <  1:3795108/76‑1 (MQ=255)

                           gctcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGctc                                              <  1:844048/76‑1 (MQ=255)

                             tcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCtt                                            <  1:1474547/76‑1 (MQ=255)

                             tcgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCtt                                            >  2:1419915/1‑76 (MQ=255)

                              cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                           >  2:1770865/1‑76 (MQ=255)

                              cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                           <  1:1904085/76‑1 (MQ=255)

                              cgctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                           >  1:4473059/1‑76 (MQ=255)

                               gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTcc                                          >  1:2601910/1‑76 (MQ=255)

                               gctcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTcc                                          <  2:1399330/76‑1 (MQ=255)

                                 tcATGCCAGGTTTATCGTCAAACCACTGATTCTGGTGTGACGGTTCGCTGCCGAGCACTCCCAACAGCCCCTcccg                                        >  2:1732424/1‑74 (MQ=255)

                                 tcATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCgg                                        <  1:1445829/76‑1 (MQ=255)

                                  cATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCggg                                       <  2:4235571/76‑1 (MQ=255)

                                   ctGCCGGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                      <  1:1299689/75‑1 (MQ=255)

                                   aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                      >  2:4146035/1‑76 (MQ=255)

                                   aTGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGATCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGc                                      <  1:4247418/76‑1 (MQ=255)

                                     gCCGGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCgg                                    <  2:2140917/76‑1 (MQ=255)

                                     gCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCgg                                    >  2:3607573/1‑76 (MQ=255)

                                     gCCAGGGCTATAGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCCTTTACAGCACCTCTTTCGGCCgg                                    >  2:871871/1‑76 (MQ=255)

                                      ccAGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGt                                   <  1:3801018/76‑1 (MQ=255)

                                       cggggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTg                                  <  2:1652553/74‑1 (MQ=255)

                                        aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCCTTTCCCGCAGCTCTTCCGGGCGGTGc                                 >  2:3377668/1‑76 (MQ=255)

                                        aGGGTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGc                                 >  1:2287606/1‑76 (MQ=255)

                                         gggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGcc                                <  1:394299/76‑1 (MQ=255)

                                          ggTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCa                               >  2:3987625/1‑76 (MQ=255)

                                           gTTATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCAc                              >  1:4781709/1‑76 (MQ=255)

                                            ttATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACt                             <  2:3730624/76‑1 (MQ=255)

                                            ttATCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACt                             <  1:3271782/76‑1 (MQ=255)

                                              ggcGGCAAACCACTGATTCTGTTTGTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCg                           <  2:4521355/74‑1 (MQ=255)

                                              aTCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACCCg                           >  2:2012483/1‑76 (MQ=255)

                                               tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                           <  2:1095041/59‑1 (MQ=255)

                                               tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTc                                           >  1:1095042/1‑59 (MQ=255)

                                               tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCCCTCgg                          >  2:1896640/1‑76 (MQ=255)

                                               tCGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCgg                          >  2:4622534/1‑76 (MQ=255)

                                                cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                         >  1:3627253/1‑76 (MQ=255)

                                                cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                         <  1:1543006/76‑1 (MQ=255)

                                                cGGCAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCggg                         <  2:4717519/76‑1 (MQ=255)

                                                  ccAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTa                       <  1:4525207/75‑1 (MQ=255)

                                                   cAAACCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTaa                      >  2:3988522/1‑76 (MQ=255)

                                                    aaaCCACTGCTTCTGTCTTGACGGTGCGCTGCCTGGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                     <  1:485964/76‑1 (MQ=255)

                                                    aaaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                     >  1:1802218/1‑76 (MQ=255)

                                                    aaaCCACTGATTCTGTTTTCACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAg                     >  1:4841257/1‑76 (MQ=255)

                                                     aaCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAgg                    <  1:2699080/76‑1 (MQ=255)

                                                      aCCACTGGTTCTGGTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                   <  2:4607766/76‑1 (MQ=255)

                                                      aCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                   <  1:792129/76‑1 (MQ=255)

                                                      aCCACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGc                   <  1:699107/76‑1 (MQ=255)

                                                        cgcTGATTCTGTTTTTACGGTGCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCCCTCGGGTAAGGCtt                 <  2:2377133/74‑1 (MQ=255)

                                                        cacaGATCCTGTTTTTACGCTTCGCTGCCAAGCATTTCCCTCAGCTCTTCGGGGCGGTGCCACTCGGTTAAGGCtt                 <  2:2240493/76‑1 (MQ=255)

                                                        cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCtt                 >  2:629389/1‑76 (MQ=255)

                                                        cACTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGCGCCACTCGGGTAAGGCtt                 >  2:1077991/1‑76 (MQ=255)

                                                          cTGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTc               >  1:3536887/1‑76 (MQ=255)

                                                           tGATTCTGTTTTTACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCg              <  1:529611/76‑1 (MQ=255)

                                                             gtgcTGTTTTGACGGTTCGCTGCGGAGGATTTCCAGCGGCTCTTCCGGGCGGTGCCACTCGGGTAGGGCTTTCGAc            <  2:3911852/73‑1 (MQ=255)

                                                                   tttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTg      <  1:1162005/76‑1 (MQ=255)

                                                                    ttttACGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGc     <  2:1952303/76‑1 (MQ=255)

CATCTTCCTGATAATAACGCAGGCCCAGCTCGCTCATGCCAGGGTTATCGGCAAACCACTGATTCTGTTTTTGCGGTTCGCTGCCGAGCATTTCCAGCAGCTCTTCCGGGCGGTGCCACTCGGGTAAGGCTTTCGACATATTGCGGG  >  NC_000913/4582978‑4583124

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 9 ≤` `ATCG/ATCG` `< 23 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`