BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,497,638	A→C	F133C (TTC→TGC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,497,638	0	A	C	98.9%	359.8 / ‑4.3	92	F133C (TTC→TGC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base A (0/1); new base C (64/27); total (64/28)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 158 total aligned reads displayed.

CTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGA  >  REL606/3497565‑3497710

cTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcag                                                                        >  1:119532/1‑76 (MQ=255)

 ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagc                                                                       >  2:4351042/1‑76 (MQ=255)

 ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagc                                                                       >  2:3663983/1‑76 (MQ=255)

 ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAAACCTAACGTCGCCAGAGAATCACGGcagc                                                                       >  2:3005793/1‑76 (MQ=255)

  tGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGgcagca                                                                      <  1:1134690/76‑1 (MQ=255)

   gAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaa                                                                     >  2:246010/1‑76 (MQ=255)

   gAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaa                                                                     <  1:3065749/76‑1 (MQ=255)

   gAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaa                                                                     <  1:2721961/76‑1 (MQ=255)

    accATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    <  2:1306408/74‑1 (MQ=255)

    aaCATCCCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    <  1:4881743/76‑1 (MQ=255)

    aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    <  2:823319/76‑1 (MQ=255)

    aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    >  2:383184/1‑76 (MQ=255)

    aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    <  1:4612950/76‑1 (MQ=255)

    aaCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    <  1:4502350/76‑1 (MQ=255)

    aaCATACCGTAGTCATAGCGGATGCCGTAACCCCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaag                                                                    >  2:4376512/1‑76 (MQ=255)

    aaCATACCGTAGTCATAGCGGATGCCGGAACCGCGCCCCCGCAAACCTAACGGCGCCAGAGAATCACGGcaggaaa                                                                    >  1:4947982/1‑75 (MQ=255)

     aCATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGcagcaagc                                                                   <  1:354616/76‑1 (MQ=255)

      cATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgcc                                                                  >  1:3795106/1‑76 (MQ=255)

      cATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgcc                                                                  <  1:1210062/76‑1 (MQ=255)

       aTACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCCAaccg                                                                 >  2:2099108/1‑76 (MQ=255)

       aTACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgccg                                                                 >  2:305034/1‑76 (MQ=255)

        tACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgccgc                                                                >  1:85545/1‑76 (MQ=255)

        tACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgacgc                                                                >  2:2381893/1‑76 (MQ=255)

         aCCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgccgcc                                                               >  2:239856/1‑76 (MQ=255)

          cggTAGTCATAGCGGAGGCGGTACCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCa                                                              <  2:721616/74‑1 (MQ=255)

          ccGTAGTCATAGCGGATGCCGTACCCGCCCCCCGCCAACCCTAACGTCGCCAGAAAATCAAGGCAGCAAGCCGCCa                                                              <  2:4209528/76‑1 (MQ=255)

          ccGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCa                                                              >  2:1346/1‑76 (MQ=255)

          ccGGAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCa                                                              <  1:1475035/76‑1 (MQ=255)

           cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAgccgcccg                                                             >  2:1571326/1‑74 (MQ=255)

           cGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAg                                                             <  1:1187690/76‑1 (MQ=255)

             tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGAc                                                           >  1:601298/1‑76 (MQ=255)

             tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGAc                                                           >  1:3604081/1‑76 (MQ=255)

             tAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGAc                                                           <  1:5045554/76‑1 (MQ=255)

              aGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCCGACg                                                          >  2:3900551/1‑76 (MQ=255)

                 cATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTTACGTCGCCAGAGAATCAAGGCAGGAAGCCGCCCGACGTcc                                                       >  2:4587451/1‑76 (MQ=255)

                  aTAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTccc                                                      >  2:1947239/1‑76 (MQ=255)

                    aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTccccg                                                    >  1:3483559/1‑74 (MQ=255)

                    aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAg                                                    <  1:1338517/76‑1 (MQ=255)

                    aGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAg                                                    >  1:1959779/1‑76 (MQ=255)

                    aGCGGATGCCGAAACCGCGCCCCGGCAACCCTAACGTGGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAg                                                    <  1:1799289/76‑1 (MQ=255)

                     gCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAgg                                                   >  1:4438674/1‑76 (MQ=255)

                      gggATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGc                                                  <  2:1825612/75‑1 (MQ=255)

                        gATGCGGTACCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCa                                                <  2:1826136/76‑1 (MQ=39)

                        gATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCa                                                <  1:3693508/76‑1 (MQ=255)

                        gATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCa                                                <  2:3042993/76‑1 (MQ=255)

                        gATGCCGTAACCGCCCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCa                                                <  1:3479167/76‑1 (MQ=255)

                         aTGCCTTACCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCAc                                               <  1:1690133/76‑1 (MQ=39)

                         aTGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCAc                                               >  1:1048591/1‑76 (MQ=255)

                         aTGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCAc                                               >  1:2404727/1‑76 (MQ=255)

                           gcggTAACCCCTCCCGGGCAACCCTACCGTCCCCAGAGAATAAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCg                                             <  2:1979666/73‑1 (MQ=25)

                            ccGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGt                                            >  2:2309690/1‑76 (MQ=255)

                             cGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGtt                                           >  1:2071455/1‑76 (MQ=255)

                              gTACCCGCGCCCCGCCACCCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                          <  1:2231115/76‑1 (MQ=37)

                              gTAACCGCGCCCCGGGAACCCTAACGTCACCAGAGAAACAAGGCAGCAACACGACAGACGTACCACGCCACTGTTa                                          >  2:4642613/1‑76 (MQ=255)

                              gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                          >  2:2749798/1‑76 (MQ=255)

                              gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                          >  2:2099845/1‑76 (MQ=255)

                              gTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTa                                          <  1:1269765/76‑1 (MQ=255)

                               tccccGCTCCCGGGCACCCCTACCGCCCCAAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCCCCGTTAc                                         <  2:3611578/73‑1 (MQ=16)

                               tAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAc                                         >  2:629925/1‑76 (MQ=255)

                               tAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAc                                         >  1:3192390/1‑76 (MQ=255)

                                aaCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTAcc                                        >  2:2872996/1‑76 (MQ=255)

                                aaCCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTAAcc                                        >  1:1902623/1‑76 (MQ=255)

                                  ccGCGCCCCGGCACCCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGa                                      <  1:4814802/76‑1 (MQ=255)

                                   cgcgCCCCGGCACCCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                     <  1:70685/76‑1 (MQ=255)

                                   cgcgCCCCGGCACCCCGAAGGCCGCCAGAGAACAAGGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAg                                     <  2:2455307/76‑1 (MQ=21)

                                    gcgcCCCGGCACCCCTAACGTCGCCAAAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                    <  2:1223668/76‑1 (MQ=37)

                                    gcgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                    >  2:4234041/1‑76 (MQ=255)

                                     cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGc                                   >  2:4600319/1‑76 (MQ=255)

                                     cgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGc                                   >  2:3387724/1‑76 (MQ=255)

                                      cccccGGCACCCCTAACGTCGCAAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGcc                                  <  2:4838329/75‑1 (MQ=37)

                                       ccccGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGccc                                 >  1:1796244/1‑76 (MQ=255)

                                         ccGGCACCCCTAACGTCGCCAGAAAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCgg                               <  2:4890163/76‑1 (MQ=37)

                                         ccGGCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCgg                               >  2:1416548/1‑76 (MQ=255)

                                           ggCACCCCAAACGTCCCAAGAAAATCAAGGCGCCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGAt                             <  2:4725465/76‑1 (MQ=21)

                                           ggCAACCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGAt                             >  1:2495487/1‑76 (MQ=255)

                                            gCAACCCTAACGCCCCCAGAGAAGCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATc                            <  1:997768/76‑1 (MQ=37)

                                              aaCCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAt                          >  1:4956967/1‑76 (MQ=255)

                                               aCCCTAACGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAtt                         >  1:1498927/1‑76 (MQ=255)

                                                 ccTACCGTCGCCAAAGAACCAGGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCAtttt                       <  1:2236829/76‑1 (MQ=25)

                                                    aaCGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTttctt                    >  1:335435/1‑76 (MQ=255)

                                                    aaCGCCGCCAGAGAACCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTttctt                    <  1:2814384/76‑1 (MQ=37)

                                                      cGTCGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtca                  >  1:2312513/1‑76 (MQ=255)

                                                         cGCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTtcatca               <  1:462009/76‑1 (MQ=255)

                                                          gCCAGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaa              >  1:4122699/1‑76 (MQ=255)

                                                          gCCAGAGAACAAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaa              <  2:44264/76‑1 (MQ=37)

                                                           ccAGAGAACCAGGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaat             <  2:4133901/76‑1 (MQ=38)

                                                           caagaaAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaat             <  2:3530019/74‑1 (MQ=38)

                                                             agagAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           <  2:4586576/76‑1 (MQ=255)

                                                             agagAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           <  2:1157268/76‑1 (MQ=255)

                                                             agagAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           <  1:1662372/76‑1 (MQ=255)

                                                             agaaAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCatcaatca           <  2:3448266/76‑1 (MQ=255)

                                                              gagaATCAGGGCAGAAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAg          <  2:3912327/76‑1 (MQ=38)

                                                              gagaATCAAGCCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAg          <  1:4550374/76‑1 (MQ=255)

                                                                  aTCAAGCCACCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGctct      <  1:3413308/76‑1 (MQ=38)

                                                                   taaaGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt     <  2:4971064/74‑1 (MQ=255)

                                                                   tCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt     >  2:1667488/1‑76 (MQ=255)

                                                                   tCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCtt     >  1:655751/1‑76 (MQ=255)

                                                                     aaGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg   >  2:4860131/1‑76 (MQ=255)

                                                                     aaGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCg   >  2:3778665/1‑76 (MQ=255)

                                                                      aGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCga  >  1:2654508/1‑76 (MQ=255)

CTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCCGGATCATTTTCTTCATCAATCAGCTCTTCGA  >  REL606/3497565‑3497710

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