BRESEQ :: Evidence

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

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

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

Only 100 of 169 total aligned reads displayed.

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601

cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGTTGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     >  1:3517142‑M1/1‑100 (MQ=255)

cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  2:262987‑M1/100‑1 (MQ=255)

cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                     <  1:3401699‑M1/100‑1 (MQ=255)

 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                    >  1:3401920‑M1/1‑99 (MQ=255)

  gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                   <  2:1741893‑M1/100‑3 (MQ=255)

  gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                   >  1:1026864‑M1/1‑98 (MQ=255)

     cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                >  2:3769652‑M1/1‑95 (MQ=255)

     cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                >  2:4149197‑M1/1‑95 (MQ=255)

       cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              <  2:2550516‑M1/100‑8 (MQ=255)

         tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                            <  2:3685563‑M1/100‑10 (MQ=255)

          cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                           >  1:357845‑M1/1‑90 (MQ=255)

          cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGagatcggaag                                                                                           >  2:3044880‑M1/1‑90 (MQ=255)

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

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

            cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcacac                                                                                         >  1:1488271‑M1/1‑88 (MQ=255)

             tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                        >  2:3234597‑M1/1‑87 (MQ=255)

             tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                        <  2:1483072‑M1/100‑14 (MQ=255)

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

               cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                      <  1:478314‑M1/100‑16 (MQ=255)

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

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

                acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                     >  2:3385151‑M1/1‑84 (MQ=255)

                  gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                   <  2:67880‑M1/100‑19 (MQ=255)

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

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

                        aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:137368‑M1/100‑25 (MQ=255)

                        aCAACGGAGGGGGGCTGCGATTTTAAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  1:450416‑M1/100‑25 (MQ=255)

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

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

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

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

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

                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                           <  2:4025387‑M1/100‑27 (MQ=255)

                             ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCGGATCACTCAAAACACTCAATTACTTAACATAAATGtataagtaatacttcagtttgggatgtat                                                                        >  1:3269346‑M1/1‑71 (MQ=255)

                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  2:2374556‑M1/100‑33 (MQ=255)

                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  1:1781622‑M1/100‑33 (MQ=255)

                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  2:3341247‑M1/100‑33 (MQ=255)

                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                     <  2:284241‑M1/100‑33 (MQ=255)

                                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                      >  1:3544547‑M1/1‑68 (MQ=255)

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

                                    ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACCTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 >  1:2422604‑M1/1‑64 (MQ=255)

                                     gCTGCGATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTTCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 <  1:2100829‑M1/99‑37 (MQ=255)

                                     gCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAAGGtataagacatacttttgttttgggtgtatttcaat                                                                  <  1:1623596‑M1/98‑36 (MQ=255)

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

                                      cTGCGATTTTCAATAATGCGTGATGAAGATCACACAAATCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                   <  1:3094644‑M1/96‑35 (MQ=255)

                                       tGCGATTTTCATTAATGCGTGATGCAGATCACACAAAACGCTCTATGACTTAACATAAATGgataagtcatacttttgttttgggtgtatttcaatc                                                                 <  1:1515319‑M1/97‑37 (MQ=255)

                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              <  2:3933023‑M1/100‑40 (MQ=255)

                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                              >  2:2660981‑M1/1‑61 (MQ=255)

                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            <  1:743646‑M1/100‑42 (MQ=255)

                                           aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                          >  2:1764030‑M1/1‑57 (MQ=255)

                                             tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                        <  2:144451‑M1/100‑46 (MQ=255)

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

                                               tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                      >  2:2036757‑M1/1‑53 (MQ=255)

                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttttttttttttttatttcaatctgtgaaaaagttt                                                    >  1:1019479‑M1/1‑51 (MQ=255)

                                                   taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                  <  2:4195505‑M1/100‑52 (MQ=255)

                                                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                 <  2:1018812‑M1/100‑53 (MQ=255)

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

                                                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                               >  1:2544418‑M1/1‑46 (MQ=255)

                                                      tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                               <  2:2851535‑M1/100‑55 (MQ=255)

                                                         gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                            <  2:1620529‑M1/100‑58 (MQ=255)

                                                            aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         <  2:1344373‑M1/100‑61 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        >  1:4275631‑M1/1‑39 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        >  1:2488339‑M1/1‑39 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                        >  2:4230223‑M1/1‑39 (MQ=255)

                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         <  2:959292‑M1/99‑61 (MQ=255)

                                                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                       <  1:2998799‑M1/100‑63 (MQ=255)

                                                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                       <  1:795021‑M1/100‑63 (MQ=255)

                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcga                                    >  2:2279130‑M1/1‑35 (MQ=255)

                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                    <  2:2869309‑M1/100‑66 (MQ=255)

                                                                  ntcACACAAAACACTCAATTACTTAACATAAATGtataactcatacttttgtattgggtgtatatcaatctgttaaaaagtttttcgctacgctagcaac                                   >  1:3043925‑M1/2‑34 (MQ=255)

                                                                  aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaan                                   <  2:4335950‑M1/100‑67 (MQ=255)

                                                                      cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttgttttttttcaatctgttaaaaagtttttcgctacactagcaatca                                 >  1:2384468‑M1/1‑30 (MQ=255)

                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctatagcaaagtg                              >  1:2610218‑M1/1‑29 (MQ=255)

                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatggtaaaaagttttttgctacgctagcaagcaaa                               >  1:3100153‑M1/1‑29 (MQ=255)

                                                                       acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgataaaaagtttttcgctacgctagtaagcaaa                               >  1:3879230‑M1/1‑29 (MQ=255)

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

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

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

                                                                           aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                          >  2:2074637‑M1/1‑25 (MQ=255)

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

                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                      >  2:2389822‑M1/1‑21 (MQ=255)

                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  1:931703‑M1/99‑79 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    <  2:2687903‑M1/100‑82 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    >  2:4056701‑M1/1‑19 (MQ=255)

                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    >  1:3097384‑M1/1‑19 (MQ=255)

                                                                                   aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  1:3332587‑M1/95‑79 (MQ=255)

                                                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa               <  2:1463227‑M1/100‑87 (MQ=255)

                                                                                      aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa               >  2:806512‑M1/1‑14 (MQ=255)

                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat              <  2:4243056‑M1/100‑88 (MQ=255)

                                                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg            <  2:3779354‑M1/100‑90 (MQ=255)

                                                                                         tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg            <  2:3459377‑M1/100‑90 (MQ=255)

                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat        <  2:2617968‑M1/100‑94 (MQ=255)

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

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

                                                                                              tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       <  2:3028919‑M1/100‑95 (MQ=255)

                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaataagttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     <  1:1374575‑M1/100‑97 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgcaatgggat   <  2:1919595‑M1/100‑99 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat   <  1:3231406‑M1/100‑99 (MQ=255)

                                                                                                  tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat   >  1:943163‑M1/1‑2 (MQ=255)

                                                                                                   gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttttgctacgctagcaatcaaacatgaaaaaggattaatcgaaatgggatg  >  1:466580‑M1/1‑1 (MQ=255)

CTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978403‑1978601

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