BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_bme_pgi	1,978,503	Δ776 bp		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	1978503–1979270	1979278	9–776	178 [0]	[0] 179	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1978502	0 (0.000)	167 (0.950)	97/180	0.1	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
*	?	NC_000913_3_bme_pgi	1979279 =	0 (0.000)	167 (0.950)	97/180	0.1	100%	intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

Only 100 of 178 total aligned reads displayed.

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCC  >  NC_000913_3_bme_pgi/1978403‑1978600

    aCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                 >  2:3136738‑M1/1‑96 (MQ=255)

     cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCATATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                <  2:4556181‑M1/99‑5 (MQ=255)

     cTCGTCACTACACGCACATACAACGGAGAAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                <  2:715005‑M1/99‑5 (MQ=255)

       cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                             >  1:339616‑M1/1‑93 (MQ=255)

        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                             <  2:4459248‑M1/99‑8 (MQ=255)

        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                             >  1:2833057‑M1/1‑92 (MQ=255)

         tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                           >  1:3633007‑M1/1‑91 (MQ=255)

            cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                        >  2:750165‑M1/1‑88 (MQ=255)

            cTACACGCACATACAACGGAGGGGGACTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                        >  1:1114126‑M1/1‑88 (MQ=255)

             tACACGCACATACAACGGAGGCGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                       >  2:4738510‑M1/1‑87 (MQ=255)

                acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                    <  2:956474‑M1/100‑17 (MQ=255)

                 cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                   >  1:2381621‑M1/1‑83 (MQ=255)

                 cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatac                                                                                        >  2:4054602‑M1/1‑83 (MQ=255)

                  gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                  <  1:1200981‑M1/100‑19 (MQ=255)

                   cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                 <  2:1750705‑M1/100‑20 (MQ=255)

                    acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                <  1:537737‑M1/100‑21 (MQ=255)

                     caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                >  2:3692760‑M1/1‑79 (MQ=255)

                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                              >  1:1213232‑M1/1‑78 (MQ=255)

                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                               >  1:3318528‑M1/1‑78 (MQ=255)

                      aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttctt                                                                              >  2:852013‑M1/1‑78 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                             <  1:3953467‑M1/100‑24 (MQ=255)

                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                           <  1:4132093‑M1/100‑26 (MQ=255)

                            cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                        <  1:65238‑M1/100‑29 (MQ=255)

                              gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                      <  1:3707248‑M1/100‑31 (MQ=255)

                                 gggggCTGCGATTTTCAATATTTCGTGATGCAGATCACACAAAACACTCATTTACTTAACATAAATGtataagtcatactttttttttgggtgtatttc                                                                    >  2:3529171‑M1/1‑67 (MQ=255)

                                 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                    <  1:1372268‑M1/99‑33 (MQ=255)

                                 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                    >  1:2154835‑M1/1‑67 (MQ=255)

                                 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACAC‑AAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                  >  1:1420414‑M1/1‑66 (MQ=255)

                                  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATATATGtataagtcatacttttgttttgggtgtatttcaa                                                                  <  2:730953‑M1/100‑35 (MQ=255)

                                  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                  <  1:1440704‑M1/100‑35 (MQ=255)

                                  ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                  <  2:4281782‑M1/100‑35 (MQ=255)

                                   gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                 >  2:1996120‑M1/1‑65 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                >  1:2570859‑M1/1‑64 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                <  1:2302514‑M1/100‑37 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgcatttcaatc                                                                >  2:414374‑M1/1‑64 (MQ=255)

                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                              >  2:2527287‑M1/1‑62 (MQ=255)

                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                              <  2:116030‑M1/100‑39 (MQ=255)

                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                >  2:229061‑M1/1‑62 (MQ=255)

                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                             >  2:1531266‑M1/1‑61 (MQ=255)

                                           aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                          >  1:180246‑M1/1‑57 (MQ=255)

                                           aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaa                                                          <  2:2088815‑M1/99‑43 (MQ=255)

                                             tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                         >  2:1068015‑M1/1‑55 (MQ=255)

                                              ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                      >  2:1697970‑M1/1‑54 (MQ=255)

                                              ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                      >  2:3194391‑M1/1‑54 (MQ=255)

                                              ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                      >  2:2544808‑M1/1‑54 (MQ=255)

                                              ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                      >  1:2862643‑M1/1‑54 (MQ=255)

                                              gtcAGTAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtattactcatacttttgttttgggtgtatttcaatctgttaaaaa                                                       <  1:4015308‑M1/98‑46 (MQ=255)

                                                cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                    >  1:4732665‑M1/1‑52 (MQ=255)

                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                   <  2:4244478‑M1/100‑50 (MQ=255)

                                                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                >  1:1804951‑M1/1‑48 (MQ=255)

                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                               <  2:3782697‑M1/100‑54 (MQ=255)

                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                               <  2:2688240‑M1/100‑54 (MQ=255)

                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatccgttaaaaagtttttcg                                               <  2:3123330‑M1/100‑54 (MQ=255)

                                                        cGTGATGAAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                             <  2:4525698‑M1/99‑56 (MQ=255)

                                                          tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                          <  1:3546427‑M1/100‑59 (MQ=255)

                                                          tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcacacttttgttttgggtgtatttcagt                                                                 >  2:4412929‑M1/1‑42 (MQ=255)

                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                <  1:4431405‑M1/77‑37 (MQ=255)

                                                            aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                         >  1:124831‑M1/1‑40 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                         >  1:1567175‑M1/1‑39 (MQ=255)

                                                               cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                      >  2:2778237‑M1/1‑37 (MQ=255)

                                                               cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgccag                                      >  1:1998602‑M1/1‑37 (MQ=255)

                                                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatctcaatctgttaaaaagtttttcgctacgctagc                                     >  1:3323169‑M1/1‑36 (MQ=255)

                                                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                  >  1:2628318‑M1/1‑33 (MQ=255)

                                                                    cacacaAAACACTCAATTACTTAACATAAATGtatcagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                 >  1:4510674‑M1/1‑32 (MQ=255)

                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                               <  1:3466652‑M1/100‑70 (MQ=255)

                                                                      cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagccagcaa                               >  1:3808531‑M1/1‑30 (MQ=255)

                                                                      cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                               >  1:1021168‑M1/1‑30 (MQ=255)

                                                                      cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                               >  2:1941372‑M1/1‑30 (MQ=255)

                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                              >  1:3323548‑M1/1‑29 (MQ=255)

                                                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                            >  2:1830224‑M1/1‑28 (MQ=255)

                                                                        caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                            <  1:3230201‑M1/100‑73 (MQ=255)

                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                           <  1:971289‑M1/100‑74 (MQ=255)

                                                                          aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                          <  2:2923202‑M1/100‑75 (MQ=255)

                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagtaaaaatgaaac                      <  1:739179‑M1/100‑79 (MQ=255)

                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                      <  1:3013543‑M1/100‑79 (MQ=255)

                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                      <  2:3436465‑M1/100‑79 (MQ=255)

                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                       >  2:1911554‑M1/1‑21 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgccagcaagcaaaaatgaaacag                    >  2:4690172‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                    >  1:1002539‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                    >  1:165667‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                    >  1:1632074‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                    >  1:1491058‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                    >  2:3089020‑M1/1‑19 (MQ=255)

                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                   <  1:4237896‑M1/99‑82 (MQ=255)

                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                   <  1:3776773‑M1/99‑82 (MQ=255)

                                                                                    ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                <  1:3836478‑M1/100‑85 (MQ=255)

                                                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtattacaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                  >  2:2849179‑M1/1‑14 (MQ=255)

                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat             >  1:3512681‑M1/1‑13 (MQ=255)

                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                 <  1:2859079‑M1/96‑84 (MQ=255)

                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaaaaaaaaaaaaaaaaatgaaacaggaat                <  2:1858287‑M1/97‑85 (MQ=255)

                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc            <  2:1239261‑M1/98‑89 (MQ=255)

                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc            <  1:1549373‑M1/98‑89 (MQ=255)

                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga          <  2:190795‑M1/99‑91 (MQ=255)

                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgataaaaaaaaaaaatgaaacaggaataatcg           <  2:2373926‑M1/98‑90 (MQ=255)

                                                                                            cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa        >  1:3005190‑M1/1‑8 (MQ=255)

                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat       <  1:1071146‑M1/100‑94 (MQ=255)

                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg      <  1:4089541‑M1/100‑95 (MQ=255)

                                                                                                attgtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg    <  1:1963683‑M1/98‑97 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat  <  2:1265446‑M1/100‑99 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat  <  2:680711‑M1/100‑99 (MQ=255)

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCC  >  NC_000913_3_bme_pgi/1978403‑1978600

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