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	229 [0]	[0] 225	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)	201 (1.030)	105/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)	201 (1.030)	105/180	0.1	100%	intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

Only 100 of 229 total aligned reads displayed.

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTC  >  NC_000913_3_bme_pgi/1978403‑1978599

 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                   >  1:4273552‑M1/1‑99 (MQ=255)

 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                   >  2:2415451‑M1/1‑99 (MQ=255)

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

     cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                               <  2:2704382‑M1/99‑5 (MQ=255)

        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                           <  1:2041307‑M1/100‑9 (MQ=255)

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

         tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                          <  1:1476376‑M1/100‑10 (MQ=255)

          cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                          >  1:4606990‑M1/1‑90 (MQ=255)

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

              acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                     >  2:285591‑M1/1‑86 (MQ=255)

              acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                     >  2:1551089‑M1/1‑86 (MQ=255)

               cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                    >  2:4364340‑M1/1‑85 (MQ=255)

               cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                    >  1:1325994‑M1/1‑85 (MQ=255)

               cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                      <  2:438682‑M1/98‑14 (MQ=255)

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

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

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

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

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

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

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

                  gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                  <  2:1348693‑M1/99‑18 (MQ=255)

                   cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAGAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                >  2:5125714‑M1/1‑81 (MQ=255)

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

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

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

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAACTGtataagtcatac                                                                                       >  2:1992110‑M1/1‑77 (MQ=255)

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

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                            >  2:3038843‑M1/1‑77 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                            <  2:4153208‑M1/100‑24 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                            >  1:2158238‑M1/1‑77 (MQ=255)

                       tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                            >  1:1790030‑M1/1‑77 (MQ=255)

                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                          >  1:3670199‑M1/1‑75 (MQ=255)

                         cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                           <  2:3263379‑M1/99‑25 (MQ=255)

                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                         >  1:2997665‑M1/1‑74 (MQ=255)

                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                         >  1:1943901‑M1/1‑74 (MQ=255)

                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGAGATGCAGATCACACTAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                          >  1:4087795‑M1/1‑74 (MQ=255)

                            cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                       <  2:1782448‑M1/100‑29 (MQ=255)

                                 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                    >  1:2425830‑M1/1‑67 (MQ=255)

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

                                   gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                <  1:559815‑M1/100‑36 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                               <  2:2928081‑M1/100‑37 (MQ=255)

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                               >  2:5154798‑M1/1‑64 (MQ=255)

                                     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                              <  1:5135414‑M1/100‑38 (MQ=255)

                                     gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                              >  2:1855220‑M1/1‑63 (MQ=255)

                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                             <  1:4554988‑M1/100‑39 (MQ=255)

                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacctttgttttgggtgtatttcaatctgt                                                            >  1:969239‑M1/1‑61 (MQ=255)

                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                            >  1:1646612‑M1/1‑60 (MQ=255)

                                         cGATGTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagttatacttttgttttgggtgtatttcaatctgtt                                                           <  2:4040888‑M1/99‑41 (MQ=255)

                                            ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttctgttttgggtgtatttcaatctgttaaa                                                        >  2:310967‑M1/1‑56 (MQ=255)

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

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

                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                  >  2:1171206‑M1/1‑51 (MQ=255)

                                                  ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACGTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                 >  1:5299644‑M1/1‑50 (MQ=255)

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

                                                     atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                               <  2:3656964‑M1/99‑53 (MQ=255)

                                                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                             <  1:2523916‑M1/100‑55 (MQ=255)

                                                       gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                             >  1:2326356‑M1/1‑45 (MQ=255)

                                                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                           >  1:5055668‑M1/1‑44 (MQ=255)

                                                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                           <  2:3735627‑M1/100‑57 (MQ=255)

                                                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                            <  2:3674304‑M1/99‑56 (MQ=255)

                                                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                             <  1:3347209‑M1/98‑55 (MQ=255)

                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                        >  1:721723‑M1/1‑41 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                      <  1:3138795‑M1/100‑62 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagc                                        >  2:4874432‑M1/1‑39 (MQ=255)

                                                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                    >  1:122817‑M1/1‑36 (MQ=255)

                                                                aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                    >  1:4477913‑M1/1‑36 (MQ=255)

                                                                  aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                 <  1:5240508‑M1/100‑67 (MQ=255)

                                                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                 <  2:3810837‑M1/99‑67 (MQ=255)

                                                                    cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                >  2:4603903‑M1/1‑32 (MQ=255)

                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                               <  1:290742‑M1/99‑69 (MQ=255)

                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                               >  1:3947231‑M1/1‑31 (MQ=255)

                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                          <  2:4507428‑M1/100‑74 (MQ=255)

                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                          >  1:1081284‑M1/1‑27 (MQ=255)

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

                                                                          aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                         <  1:2589627‑M1/100‑75 (MQ=255)

                                                                           aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaatttttttcgcttcgctagcaagcaaaaatga                        <  2:2140277‑M1/100‑76 (MQ=255)

                                                                             cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                      >  1:3175023‑M1/1‑23 (MQ=255)

                                                                             cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                      >  2:1069426‑M1/1‑23 (MQ=255)

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

                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                      >  2:2863251‑M1/1‑22 (MQ=255)

                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                      >  2:2249909‑M1/1‑22 (MQ=255)

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

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

                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgataaaaagaaaaaatgaaaca                    <  2:503920‑M1/100‑80 (MQ=255)

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

                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                  <  2:1701881‑M1/99‑82 (MQ=255)

                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg          <  2:958654‑M1/99‑90 (MQ=255)

                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg          <  2:1934534‑M1/99‑90 (MQ=255)

                                                                                          aaCATAAACGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttcttcgctacgatagcaaaaaaaaatgaaacaggaataatcg          <  2:4208219‑M1/99‑90 (MQ=255)

                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacaataacaaaaaaaaatgaaacaggaataatc           <  2:2564330‑M1/97‑89 (MQ=255)

                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat      >  1:2822083‑M1/1‑7 (MQ=255)

                                                                                             aTAAATGtataagtcatacgtttgttttgggtggatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                     <  2:4305353‑M1/85‑79 (MQ=255)

                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg     >  1:3204222‑M1/1‑6 (MQ=255)

                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagaaagcaaaattgaaa                      >  2:4751203‑M1/1‑6 (MQ=255)

                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg    >  1:3362616‑M1/1‑4 (MQ=255)

                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga  >  2:1617466‑M1/1‑3 (MQ=255)

                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga  >  2:766324‑M1/1‑3 (MQ=255)

                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaaaaaaaaaaaaaaaaatgaaacaggaataatcgaaatggga  <  2:3866244‑M1/100‑98 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga  >  2:913007‑M1/1‑2 (MQ=255)

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTC  >  NC_000913_3_bme_pgi/1978403‑1978599

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