BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB‑5–insA‑5	insB‑5, insA‑5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	117 [0]	[0] 118	insB‑5–insA‑5	insB‑5,insA‑5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	103 (0.990)	69/180	0.3	100%	intergenic (‑305/+16)	flhD/insB‑5	DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB
*	?	NC_000913	1979279 =	0 (0.000)	103 (0.990)	69/180	0.3	100%	intergenic (‑64/‑474)	insA‑5/uspC	IS1 protein InsA/universal stress protein C

Only 100 of 117 total aligned reads displayed.

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTC  >  NC_000913/1978406‑1978599

tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                 >  1:1835115‑M1/1‑97 (MQ=255)

   tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                              >  1:634535‑M1/1‑94 (MQ=255)

    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                             <  2:2654879‑M1/99‑7 (MQ=255)

    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataag                                                                                             <  2:1997243‑M1/99‑7 (MQ=255)

      tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                           >  2:915584‑M1/1‑91 (MQ=255)

         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                        <  1:2952674‑M1/99‑12 (MQ=255)

         cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                        <  1:430191‑M1/99‑12 (MQ=255)

          tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                       <  1:2422092‑M1/99‑13 (MQ=255)

           acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                      <  1:1467911‑M1/99‑14 (MQ=255)

            cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                     <  1:3010332‑M1/99‑15 (MQ=255)

            cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                     <  1:1483960‑M1/99‑15 (MQ=255)

             acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                    <  1:2035087‑M1/99‑16 (MQ=255)

               gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                  <  1:1211365‑M1/99‑18 (MQ=255)

                cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                 <  1:2441341‑M1/99‑19 (MQ=255)

                cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                 <  1:882069‑M1/99‑19 (MQ=255)

                 acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                <  1:2653900‑M1/99‑20 (MQ=255)

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

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

                   aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                              <  2:52741‑M1/99‑22 (MQ=255)

                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                             >  1:821050‑M1/1‑77 (MQ=255)

                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                             >  1:600433‑M1/1‑77 (MQ=255)

                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                             >  1:270422‑M1/1‑77 (MQ=255)

                        aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                         <  1:2288630‑M1/99‑27 (MQ=255)

                        aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                         >  2:764760‑M1/1‑73 (MQ=255)

                          ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                       <  1:1601144‑M1/99‑29 (MQ=255)

                          ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                       >  2:43388‑M1/1‑71 (MQ=255)

                          ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                       <  1:1369596‑M1/99‑29 (MQ=255)

                          ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                       <  1:2963261‑M1/99‑29 (MQ=255)

                           gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                      <  1:385814‑M1/99‑30 (MQ=255)

                           gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                      >  2:663748‑M1/1‑70 (MQ=255)

                           gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                      <  1:1134880‑M1/99‑30 (MQ=255)

                           gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                      >  1:1227143‑M1/1‑70 (MQ=255)

                            aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                     <  2:406886‑M1/99‑31 (MQ=255)

                            aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                     <  2:565509‑M1/99‑31 (MQ=255)

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

                              gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                   <  2:2839769‑M1/99‑33 (MQ=255)

                                gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtttatttcaa                                                                 <  2:623358‑M1/99‑35 (MQ=255)

                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                               <  2:2918956‑M1/99‑37 (MQ=255)

                                  gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                               <  1:1282631‑M1/99‑37 (MQ=255)

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

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

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

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

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

                                           ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                      <  1:2026956‑M1/99‑46 (MQ=255)

                                           ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                      <  2:1257312‑M1/99‑46 (MQ=255)

                                             cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                    <  2:2450184‑M1/99‑48 (MQ=255)

                                              aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                   <  1:2261913‑M1/99‑49 (MQ=255)

                                              aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                   >  1:1541239‑M1/1‑51 (MQ=255)

                                               ataatGCGTGATGCAGATCACACAAAACACTCCATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                  >  2:2270816‑M1/1‑50 (MQ=255)

                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                  <  1:386560‑M1/99‑50 (MQ=255)

                                               ataatGCGTGATGCAAATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                  <  1:1047325‑M1/99‑50 (MQ=255)

                                                taatGCGTGGTGCAGATCACACAAAACCCTCTATTACTTAACATTAATGtatatgtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                 <  2:1156776‑M1/99‑51 (MQ=255)

                                                taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                 >  1:2433289‑M1/1‑49 (MQ=255)

                                                taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                 >  1:1501557‑M1/1‑49 (MQ=255)

                                                 aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                >  1:1441322‑M1/1‑48 (MQ=255)

                                                 aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                <  2:870580‑M1/99‑52 (MQ=255)

                                                   tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcg                                              <  2:2064697‑M1/99‑54 (MQ=255)

                                                    gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                             <  2:1808950‑M1/99‑55 (MQ=255)

                                                      gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                           <  2:383824‑M1/99‑57 (MQ=255)

                                                      gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                           <  2:2787282‑M1/99‑57 (MQ=255)

                                                      gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                           >  2:1883914‑M1/1‑43 (MQ=255)

                                                        gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacg                                         >  2:772598‑M1/1‑41 (MQ=255)

                                                          tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                       >  1:1499879‑M1/1‑39 (MQ=255)

                                                              gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttcgctacgctagcaa                                   <  2:2688411‑M1/99‑65 (MQ=255)

                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                  <  1:539067‑M1/99‑66 (MQ=255)

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

                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                            >  1:3048288‑M1/1‑28 (MQ=255)

                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                            <  2:2519143‑M1/99‑72 (MQ=255)

                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                            >  1:874866‑M1/1‑28 (MQ=255)

                                                                       aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                          >  2:166545‑M1/1‑26 (MQ=255)

                                                                       aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                          >  2:539427‑M1/1‑26 (MQ=255)

                                                                       aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                          >  1:2465803‑M1/1‑26 (MQ=255)

                                                                       aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                          >  1:914504‑M1/1‑26 (MQ=255)

                                                                        aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                         >  2:2343244‑M1/1‑25 (MQ=255)

                                                                          cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                       >  1:2126458‑M1/1‑23 (MQ=255)

                                                                          cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                       >  1:2274623‑M1/1‑23 (MQ=255)

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

                                                                               aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                  >  2:3008762‑M1/1‑18 (MQ=255)

                                                                                aTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                 >  1:3061922‑M1/1‑17 (MQ=255)

                                                                                 ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                >  1:3046594‑M1/1‑16 (MQ=255)

                                                                                  tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat               <  1:2535286‑M1/99‑85 (MQ=255)

                                                                                  tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat               <  2:58712‑M1/99‑85 (MQ=255)

                                                                                   aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata              <  1:2587149‑M1/99‑86 (MQ=255)

                                                                                   aCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata              >  2:183674‑M1/1‑14 (MQ=255)

                                                                                    cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa             >  1:1476900‑M1/1‑13 (MQ=255)

                                                                                    cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa             >  1:1492035‑M1/1‑13 (MQ=255)

                                                                                    cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa             <  1:1264405‑M1/99‑87 (MQ=255)

                                                                                     ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat            >  2:1480590‑M1/1‑12 (MQ=255)

                                                                                     ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat            <  2:596833‑M1/99‑88 (MQ=255)

                                                                                         cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa        >  1:1570733‑M1/1‑8 (MQ=255)

                                                                                         cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa        <  2:1071386‑M1/99‑92 (MQ=255)

                                                                                           tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat      >  1:19086‑M1/1‑6 (MQ=255)

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

                                                                                             aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg    <  2:1489747‑M1/99‑96 (MQ=255)

                                                                                              aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg   >  2:723894‑M1/1‑3 (MQ=255)

                                                                                              aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg   >  2:1442826‑M1/1‑3 (MQ=255)

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

                                                                                               tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga  <  2:2132013‑M1/99‑98 (MQ=255)

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

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTC  >  NC_000913/1978406‑1978599

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