BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,976,527	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1976527–1977294	1977302	9–776	270 [0]	[0] 273	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1976526	0 (0.000)	238 (0.960)	115/158	0.1	100%	intergenic (‑305/+16)	flhD/insB1	DNA‑binding transcriptional dual regulator with FlhC/IS1 transposase B
*	?	NC_000913	1977303 =	0 (0.000)	238 (0.960)	115/158	0.1	100%	intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

Only 100 of 270 total aligned reads displayed.

CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCG  >  NC_000913/1976434‑1976618

cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGa                                                                                             >  1:3069926‑M1/1‑93 (MQ=255)

  acaCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                           <  1:414353‑M1/93‑4 (MQ=255)

     ttaCACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                        <  1:6064624‑M1/93‑7 (MQ=255)

     cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                        <  1:5737838‑M1/94‑7 (MQ=255)

         tcgCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                    <  1:3495899‑M1/93‑11 (MQ=255)

           tcagatacaACGGAGGGGTGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                  <  1:1453958‑M1/90‑13 (MQ=255)

             gcatACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                <  1:1396510‑M1/93‑15 (MQ=255)

              tatACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                               <  1:5864806‑M1/93‑16 (MQ=255)

              caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                           >  2:6809788‑M1/1‑79 (MQ=255)

              caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                           <  2:587126‑M1/82‑4 (MQ=255)

              caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                           <  2:4746814‑M1/82‑4 (MQ=255)

              aatACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                               <  1:2091006‑M1/93‑16 (MQ=255)

               ttaCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                              <  1:6004554‑M1/93‑17 (MQ=255)

                tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                         >  2:3643459‑M1/1‑77 (MQ=255)

                tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                         >  2:4557754‑M1/1‑77 (MQ=255)

                 gcaACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                            <  1:1504942‑M1/93‑19 (MQ=255)

                    aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                     >  2:2476573‑M1/1‑73 (MQ=255)

                      ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttga                                                                       >  1:6076669‑M1/1‑71 (MQ=255)

                      cgaGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                       <  1:1013505‑M1/93‑24 (MQ=255)

                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                 <  2:3419667‑M1/82‑14 (MQ=255)

                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                 >  2:3865445‑M1/1‑69 (MQ=255)

                         tgggggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                    <  1:5391686‑M1/93‑27 (MQ=255)

                         ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtc                                                                    >  1:5033053‑M1/1‑68 (MQ=255)

                         ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                <  2:1642238‑M1/82‑15 (MQ=255)

                          tggggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                   <  1:3858976‑M1/93‑28 (MQ=255)

                           tgggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                  <  1:250755‑M1/93‑29 (MQ=255)

                            tggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                 <  1:1589472‑M1/93‑30 (MQ=255)

                            tggcTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                 <  1:545751‑M1/93‑30 (MQ=255)

                            gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                             <  2:338075‑M1/82‑18 (MQ=255)

                             ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                            <  2:94808‑M1/82‑19 (MQ=255)

                               cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                          >  2:1227685‑M1/1‑62 (MQ=255)

                                  cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                           <  1:1888372‑M1/94‑36 (MQ=255)

                                   tatTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                          <  1:3541232‑M1/93‑37 (MQ=255)

                                   gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                      <  2:1994572‑M1/82‑25 (MQ=255)

                                    aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcg                                                         >  1:4020872‑M1/1‑57 (MQ=255)

                                     ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctt                                                        >  1:5109635‑M1/1‑56 (MQ=255)

                                     ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                    <  2:6267370‑M1/82‑27 (MQ=255)

                                     ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                    <  2:5376422‑M1/82‑27 (MQ=255)

                                      tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                       <  1:4044644‑M1/94‑40 (MQ=255)

                                      attcAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                       <  1:2333020‑M1/93‑40 (MQ=255)

                                       ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                  <  2:3352472‑M1/82‑29 (MQ=255)

                                         cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                <  2:5580009‑M1/82‑31 (MQ=255)

                                          aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATACATGtataagtcatacttttgttttgggtgtatttcaatctgttaac                                                   >  1:105440‑M1/1‑51 (MQ=255)

                                          aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                               >  2:2117122‑M1/1‑51 (MQ=255)

                                          aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                               >  2:6281784‑M1/1‑51 (MQ=255)

                                           ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                              >  2:3604399‑M1/1‑50 (MQ=255)

                                            taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                             <  2:5297342‑M1/82‑34 (MQ=255)

                                            aaatgCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                 <  1:4272832‑M1/93‑46 (MQ=255)

                                             aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                            <  2:1927866‑M1/82‑35 (MQ=255)

                                               tgagtgaTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                          >  2:7295420‑M1/4‑46 (MQ=255)

                                               tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagta                                              >  1:719989‑M1/1‑46 (MQ=255)

                                               tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                          <  2:1808653‑M1/82‑37 (MQ=255)

                                                 tgtGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                            <  1:5301733‑M1/93‑51 (MQ=255)

                                                  gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttta                                           >  1:6468944‑M1/1‑43 (MQ=255)

                                                  gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                       >  2:6657466‑M1/1‑43 (MQ=255)

                                                    gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                     >  2:3132151‑M1/1‑41 (MQ=255)

                                                      tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                   <  2:2880315‑M1/82‑44 (MQ=255)

                                                      tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                   <  2:3522492‑M1/82‑44 (MQ=255)

                                                        cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctag                                     >  1:2702642‑M1/1‑37 (MQ=255)

                                                          gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgg                                   >  1:245168‑M1/1‑35 (MQ=255)

                                                          gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacga                                   >  1:4241097‑M1/1‑35 (MQ=255)

                                                            tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                             >  2:5131349‑M1/1‑33 (MQ=255)

                                                             cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                <  1:5564579‑M1/94‑63 (MQ=255)

                                                              acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctaga                               >  1:3112766‑M1/1‑31 (MQ=255)

                                                               cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcg                              >  1:905398‑M1/1‑30 (MQ=255)

                                                                acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcag                             >  1:3929889‑M1/1‑29 (MQ=255)

                                                                acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                         >  2:4033963‑M1/1‑29 (MQ=255)

                                                                 caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaat                            >  1:366051‑M1/1‑28 (MQ=255)

                                                                   taacACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                          <  1:6784064‑M1/93‑69 (MQ=255)

                                                                   caacACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctatgctagcaagca                          <  1:6266039‑M1/93‑69 (MQ=255)

                                                                   aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagct                          >  1:3434587‑M1/1‑26 (MQ=255)

                                                                   aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                      <  2:6070448‑M1/82‑57 (MQ=255)

                                                                    aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                         >  1:49481‑M1/1‑25 (MQ=255)

                                                                    aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                     <  2:4829882‑M1/82‑58 (MQ=255)

                                                                      cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaag                       >  1:774253‑M1/1‑23 (MQ=255)

                                                                       cctCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                      <  1:4696487‑M1/93‑73 (MQ=255)

                                                                       acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaag                      >  1:812209‑M1/1‑22 (MQ=255)

                                                                        cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaaa                     >  1:7148622‑M1/1‑21 (MQ=255)

                                                                           gatTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                  <  1:6955772‑M1/93‑77 (MQ=255)

                                                                           aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                              <  2:3990021‑M1/82‑65 (MQ=255)

                                                                              tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                           <  2:2357463‑M1/82‑68 (MQ=255)

                                                                               cctTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag              <  1:448110‑M1/93‑81 (MQ=255)

                                                                               aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacat              >  1:7421923‑M1/1‑14 (MQ=255)

                                                                                 ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggt            >  1:3079921‑M1/1‑12 (MQ=255)

                                                                                  tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                       >  2:382473‑M1/1‑11 (MQ=255)

                                                                                  caacATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa           <  1:671196‑M1/93‑84 (MQ=255)

                                                                                   tacATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat          <  1:6524499‑M1/93‑85 (MQ=255)

                                                                                    aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                     >  2:5933702‑M1/1‑9 (MQ=255)

                                                                                    aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                     >  2:3586446‑M1/1‑9 (MQ=255)

                                                                                    aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                     <  2:6444340‑M1/82‑74 (MQ=255)

                                                                                    aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagctacgctagcaagcaaaaatgaaacaggaatc         >  1:6314590‑M1/1‑9 (MQ=255)

                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat       >  1:1744751‑M1/1‑7 (MQ=255)

                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                   >  2:2912824‑M1/1‑7 (MQ=255)

                                                                                        aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                 <  2:5031936‑M1/82‑78 (MQ=255)

                                                                                         aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga    >  1:3783075‑M1/1‑4 (MQ=255)

                                                                                         aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                >  2:4293376‑M1/1‑4 (MQ=255)

                                                                                         aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                <  2:4200095‑M1/82‑79 (MQ=255)

                                                                                           tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa  <  1:6221152‑M1/94‑93 (MQ=255)

                                                                                            gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg             >  2:2552097‑M1/1‑1 (MQ=255)

                                                                                            gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg             >  2:1352932‑M1/1‑1 (MQ=255)

CGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCG  >  NC_000913/1976434‑1976618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`