BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,021,684	Δ1 bp	coding (542/609 nt)\|coding (1/243 nt)	b4556 → tatD →	hypothetical protein Mg‑dependent DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,021,683	0	C	.	99.2%	590.2 / ‑3.2	130	coding (541/609 nt)	b4556	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base . (68/61); total (68/62)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 137 total aligned reads displayed.

GGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGC  >  NC_000913/4021608‑4021754

ggCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAc                                                                         <  1:203758/76‑1 (MQ=255)

 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcg                                                                        >  1:2322115/1‑75 (MQ=255)

   acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                      >  2:1115011/1‑73 (MQ=255)

   acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                      >  2:1986395/1‑73 (MQ=255)

   acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                      <  1:1467793/76‑4 (MQ=255)

     acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgc                                                                   <  1:4585033/76‑1 (MQ=255)

     acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgc                                                                   <  1:2770837/76‑1 (MQ=255)

      cgcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGAcgcg                                                                  >  2:4141951/1‑76 (MQ=255)

         gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGAc                                                               <  1:1573523/76‑1 (MQ=255)

          aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGCCt                                                              >  2:4372136/1‑76 (MQ=255)

          aaGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACt                                                              >  1:3070096/1‑76 (MQ=255)

          aaGAGAGGCAGGCGTGCGGGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCTGACt                                                              <  2:2865764/76‑1 (MQ=255)

           agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                             >  2:2581991/1‑76 (MQ=255)

           agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTg                                                             <  1:2166417/76‑1 (MQ=255)

            gagaTGCAGGCGTGCGGGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTgg                                                            <  2:958319/76‑1 (MQ=255)

              gaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAg                                                          <  2:2736266/76‑1 (MQ=255)

               aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGc                                                         <  2:3721360/76‑1 (MQ=255)

                tGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCt                                                        >  1:3036524/1‑76 (MQ=255)

                 gCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTg                                                       >  2:3407206/1‑76 (MQ=255)

                  cAGGCGTGCGTGTCGGATGGAATTTATATCGGCATTACCGGTTGTGTTTTCGATGAA‑CGACGCGGACTGGAGCTGc                                                      <  1:2620582/76‑1 (MQ=255)

                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGc                                                      <  1:3747703/76‑1 (MQ=255)

                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGc                                                      <  1:3305661/76‑1 (MQ=255)

                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGc                                                      >  2:492852/1‑76 (MQ=255)

                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCCATGAA‑CGACGCGGACTGGAGCtcc                                                      >  2:757017/1‑74 (MQ=255)

                   aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCg                                                     <  1:2803529/76‑1 (MQ=255)

                   aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCg                                                     >  1:3173258/1‑76 (MQ=255)

                   aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCg                                                     >  1:3854815/1‑76 (MQ=255)

                    ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGTGCTGCgg                                                    >  1:991416/1‑76 (MQ=255)

                    ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCgg                                                    >  2:2250078/1‑76 (MQ=255)

                      cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGa                                                  <  1:921212/76‑1 (MQ=255)

                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCgggca                                                 >  1:2244413/1‑74 (MQ=255)

                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGaa                                                 >  1:2807029/1‑76 (MQ=255)

                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGaa                                                 <  1:2247119/76‑1 (MQ=255)

                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGaa                                                 >  1:3702477/1‑76 (MQ=255)

                        tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAAc                                                >  1:3280238/1‑76 (MQ=255)

                        tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGCAAc                                                >  2:387946/1‑76 (MQ=255)

                         gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACt                                               <  1:3458536/76‑1 (MQ=255)

                          cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtt                                              >  2:4511844/1‑76 (MQ=255)

                            tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACtttt                                            >  2:197724/1‑76 (MQ=255)

                             ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGCA‑CGACGCGGACTGGAGCTGCGGGAACTTTTg                                           >  2:2322601/1‑76 (MQ=255)

                             ggCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACCTTTg                                           <  1:720752/76‑1 (MQ=255)

                              gcgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                          >  1:1419236/1‑76 (MQ=255)

                              gcgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGc                                          >  2:2427814/1‑76 (MQ=255)

                               cgcAGGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGcc                                         <  2:789566/76‑1 (MQ=255)

                                gcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                        >  1:2039957/1‑76 (MQ=255)

                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                      <  2:247433/76‑1 (MQ=255)

                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                      <  1:2099748/76‑1 (MQ=255)

                                   tGGAATTTTTATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTg                                     <  2:1127414/76‑1 (MQ=255)

                                   tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGCCGCGGACTGGAGCTGCGGGAACTTTTGCCGTTg                                     >  2:3942057/1‑76 (MQ=255)

                                   tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTg                                     <  1:4716029/76‑1 (MQ=255)

                                     gAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAt                                   <  1:1482161/76‑1 (MQ=255)

                                      aagttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                  <  2:845511/73‑1 (MQ=255)

                                      aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                  >  2:3239080/1‑76 (MQ=255)

                                      aaTTTATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGAtt                                  >  2:960836/1‑76 (MQ=255)

                                        tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                <  2:2158414/76‑1 (MQ=255)

                                        tttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                >  1:1824628/1‑76 (MQ=255)

                                         ttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                               <  1:772254/76‑1 (MQ=255)

                                         ttATATCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcg                               <  1:1063510/76‑1 (MQ=255)

                                          tataTCGGCATGACCGGCTGGGCTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTGTTGCCGTTGACTTcgg                              >  2:3049479/1‑76 (MQ=255)

                                           atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                             <  2:1741437/76‑1 (MQ=255)

                                           atatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggc                             <  1:846576/76‑1 (MQ=255)

                                            tatCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                            >  1:4308514/1‑76 (MQ=255)

                                             atCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                           >  2:1426222/1‑76 (MQ=255)

                                             atCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                           <  2:1616277/76‑1 (MQ=255)

                                              tCGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                          <  1:4511844/76‑1 (MQ=255)

                                               cGGCATTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaa                         <  1:1050667/76‑1 (MQ=255)

                                                  cATTGCCGGTTGGGATTGCGATGAA‑CGACGCGGACTGGAGCTGCGGTAACTTTTGCCTTTGATTCCGGCGGaaaaa                      <  1:43972/76‑1 (MQ=255)

                                                   aTTACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                     <  1:271080/76‑1 (MQ=255)

                                                    ttACCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                    >  2:1831373/1‑76 (MQ=255)

                                                      aCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                  >  2:1422619/1‑76 (MQ=255)

                                                      aCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                  >  1:560209/1‑76 (MQ=255)

                                                      aCCGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                  >  2:4340095/1‑76 (MQ=255)

                                                       ccGGTTGGGTTTGCGCTGAA‑CGGCGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                 <  1:4438879/76‑1 (MQ=255)

                                                       ccGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                 <  1:1986396/76‑1 (MQ=255)

                                                        cGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGCAAAATTAATg                >  2:3588580/1‑76 (MQ=255)

                                                        cGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                >  1:4172299/1‑76 (MQ=255)

                                                        cGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                >  2:3110375/1‑76 (MQ=255)

                                                        cGGTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTATGCCGTTGATTCCGGCGGAAAAATTACTg                >  1:1112686/1‑76 (MQ=255)

                                                        cGGTTGGGGTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                <  2:3701043/76‑1 (MQ=255)

                                                         ggTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAACATTACTgc               >  2:389370/1‑75 (MQ=255)

                                                         ggTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa               >  1:4146310/1‑76 (MQ=255)

                                                          gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt              <  1:1714372/76‑1 (MQ=255)

                                                          gTTGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt              <  2:4000616/76‑1 (MQ=255)

                                                           ttGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc             <  1:1027926/76‑1 (MQ=255)

                                                            tGGGTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg            >  1:565650/1‑76 (MQ=255)

                                                             gggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGa           >  2:3233435/1‑76 (MQ=255)

                                                              ggTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATAACTGATCGaa          >  2:3503390/1‑76 (MQ=255)

                                                               gTTTGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaaa         <  2:3580502/76‑1 (MQ=255)

                                                                 ttGCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACt       <  1:3464228/76‑1 (MQ=255)

                                                                  tGCGGTGAA‑CGACGCGGACTGGAGCTGCGGGAGCTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTg      <  2:2277136/76‑1 (MQ=255)

                                                                  ggCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTg      <  2:1786180/75‑1 (MQ=255)

                                                                   gCGGTGAA‑CGACGCGGAGTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa     <  2:1101658/76‑1 (MQ=255)

                                                                   gCGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGa     >  2:2482319/1‑76 (MQ=255)

                                                                    cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt    <  2:3896007/76‑1 (MQ=255)

                                                                    cGATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt    >  2:595744/1‑76 (MQ=255)

                                                                     gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg   >  1:3512357/1‑76 (MQ=255)

                                                                     gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg   >  1:731434/1‑76 (MQ=255)

                                                                     gATGAA‑CGACGCGGACTGGAGCTGCGGGAACTCTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg   >  1:1849917/1‑76 (MQ=255)

                                                                     gATGAA‑CGACGCGGACGGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATg   <  2:207091/76‑1 (MQ=255)

                                                                      aTGAA‑CGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc  <  1:2134706/76‑1 (MQ=255)

GGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGC  >  NC_000913/4021608‑4021754

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