BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	676,374	Δ11,477 bp	IS5‑mediated	[ybeQ]–gltI	[ybeQ], ybeR, djlB, ybeT, ybeU, djlC, hscC, rihA, gltL, gltK, gltJ, sroC, gltI

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	676374	689014–687846	11473–12641	1082 [0]	[638] 652	[ybeQ]–insH3	[ybeQ],ybeR,djlB,ybeT,ybeU,djlC,hscC,rihA,gltL,gltK,gltJ,sroC,gltI,insH3

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	273955 =	NA (NA)	1071 (1.210)	195/200	0.0	100%	noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 676373	0 (0.000)	1071 (1.210)	195/200	0.0	100%	coding (175/978 nt)	ybeQ	Sel1 repeat‑containing protein YbeQ

Only 100 of 652 total aligned reads displayed.

TCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCC  >  NC_000913/688952‑689115

tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTa                                                                 <  2:18799673/101‑1 (MQ=11)

 cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTaa                                                                >  1:17240669/1‑101 (MQ=31)

 cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTaa                                                                <  2:8816264/101‑1 (MQ=31)

 cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTaa                                                                <  1:15147450/101‑1 (MQ=31)

  ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAAc                                                               >  1:10834789/1‑101 (MQ=32)

   gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACa                                                              >  1:2855177/1‑101 (MQ=32)

    cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAt                                                             <  1:4657662/101‑1 (MQ=33)

    cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAt                                                             <  2:3601845/101‑1 (MQ=33)

     gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                            >  2:15997240/1‑101 (MQ=35)

     gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                            >  2:5303638/1‑101 (MQ=35)

     gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAtt                                                            <  1:8830852/101‑1 (MQ=35)

      aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAttt                                                           <  2:17020579/101‑1 (MQ=38)

      aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAttt                                                           <  1:18339027/101‑1 (MQ=38)

      aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAttt                                                           >  2:7222238/1‑101 (MQ=35)

      aaGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACAttt                                                           >  1:16157904/1‑101 (MQ=35)

       aGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTa                                                          <  2:14216661/101‑1 (MQ=38)

       aGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTa                                                          >  1:8532662/1‑101 (MQ=38)

         gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGt                                                        <  2:14677627/101‑1 (MQ=40)

         gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGt                                                        <  2:18430590/101‑1 (MQ=40)

          tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGtt                                                       <  1:381095/101‑1 (MQ=40)

            aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTg                                                     <  1:7651981/101‑1 (MQ=255)

            aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTg                                                     >  2:941045/1‑101 (MQ=255)

            aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTg                                                     <  2:6751011/101‑1 (MQ=255)

             gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTgg                                                    <  1:17729701/101‑1 (MQ=255)

             gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTgg                                                    <  2:12885603/101‑1 (MQ=255)

              ttGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGc                                                   >  2:17943888/1‑101 (MQ=255)

               tgatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCt                                                  >  1:6581313/1‑101 (MQ=255)

                gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCTCACCTTCCTTAGGTAACATTTAGTTTGGCTa                                                 <  1:10006868/101‑1 (MQ=255)

                gatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTa                                                 >  1:4310838/1‑101 (MQ=255)

                   gaCTCATGTTGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAAt                                              >  2:5647243/1‑101 (MQ=255)

                   gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAAt                                              >  2:16800271/1‑101 (MQ=255)

                   gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAAt                                              >  1:6965237/1‑101 (MQ=255)

                    aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATg                                             >  2:6023335/1‑101 (MQ=255)

                    aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATg                                             >  2:1010136/1‑101 (MQ=255)

                    aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATg                                             >  2:9014883/1‑101 (MQ=255)

                    aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAAGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATg                                             <  1:3583373/101‑1 (MQ=255)

                     cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGt                                            >  1:5225184/1‑101 (MQ=255)

                     cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGt                                            <  2:8513827/101‑1 (MQ=255)

                        atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTaaa                                         <  1:2162683/101‑1 (MQ=255)

                         tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAg                                        <  1:19967436/101‑1 (MQ=255)

                         tgatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAg                                        >  2:12149124/1‑101 (MQ=255)

                          gatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGa                                       >  1:7367133/1‑101 (MQ=255)

                          gatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGa                                       >  1:3135379/1‑101 (MQ=255)

                            tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGata                                     <  1:16814217/101‑1 (MQ=255)

                            tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGata                                     <  1:20158897/101‑1 (MQ=255)

                             gaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGatat                                    <  1:12861349/101‑1 (MQ=255)

                              aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATAtt                                   <  2:19619804/101‑1 (MQ=255)

                              aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATAtt                                   >  1:13666980/1‑101 (MQ=255)

                              aaCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCATTAGGTAACATTTAGTTTGGCTAAATGTAAAGATAtt                                   >  2:20282734/1‑101 (MQ=255)

                               aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTg                                  >  1:1013578/1‑101 (MQ=255)

                                cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGc                                 <  1:18820596/101‑1 (MQ=255)

                                cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGc                                 <  1:19058732/101‑1 (MQ=255)

                                cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGc                                 <  1:5879601/101‑1 (MQ=255)

                                cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGc                                 >  1:19747057/1‑101 (MQ=255)

                                 ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCt                                >  1:11030100/1‑101 (MQ=255)

                                  cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                               <  2:3047921/101‑1 (MQ=255)

                                  cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                               <  1:12632922/101‑1 (MQ=255)

                                  cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                               <  1:6131135/101‑1 (MQ=255)

                                  cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTg                               >  2:4195403/1‑101 (MQ=255)

                                   tGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGt                              <  2:11722014/101‑1 (MQ=255)

                                    gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtt                             <  2:17805416/101‑1 (MQ=255)

                                     ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGttt                            <  1:15430005/101‑1 (MQ=255)

                                      tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGtttt                           >  1:11673635/1‑101 (MQ=255)

                                       cTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTa                          >  1:13057276/1‑101 (MQ=255)

                                        tATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAt                         <  2:17433971/101‑1 (MQ=255)

                                         aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                        <  2:20640780/101‑1 (MQ=255)

                                         aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                        >  1:5525241/1‑101 (MQ=255)

                                         aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                        >  1:13984926/1‑101 (MQ=255)

                                         aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAtt                        <  1:11210729/101‑1 (MQ=255)

                                            gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtt                     <  2:7971924/101‑1 (MQ=255)

                                            gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtt                     <  2:3176765/101‑1 (MQ=255)

                                            gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtt                     <  1:18135319/101‑1 (MQ=255)

                                            gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtt                     <  1:16435707/101‑1 (MQ=255)

                                              tCCAGATGACAAACATGATCTCATATCAGGGCCTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtttg                   <  2:19938854/101‑1 (MQ=255)

                                              tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtttg                   <  2:4521556/101‑1 (MQ=255)

                                               ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTAttgtttgt                  <  1:12194708/101‑1 (MQ=255)

                                                cAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATtgtttgtt                 >  1:14659854/1‑101 (MQ=255)

                                                 aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTgtttgttt                <  2:13671791/101‑1 (MQ=255)

                                                  gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGtttt               <  2:9715814/101‑1 (MQ=255)

                                                  gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGtttt               >  1:2355733/1‑101 (MQ=255)

                                                   aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGttttt              <  1:15714354/101‑1 (MQ=255)

                                                   aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGttttt              <  2:18207478/101‑1 (MQ=255)

                                                    tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTATGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGttttgc             <  2:778751/101‑3 (MQ=255)

                                                    tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTg             >  2:11343895/1‑101 (MQ=255)

                                                    tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTg             <  1:9159383/101‑1 (MQ=255)

                                                    tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTg             >  2:6975129/1‑101 (MQ=255)

                                                      aCAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCg           <  2:2669395/101‑1 (MQ=255)

                                                       cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCga          <  1:13500611/101‑1 (MQ=255)

                                                        aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCgag         <  2:12395159/101‑1 (MQ=255)

                                                         aaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCgaga        <  2:13321896/101‑1 (MQ=255)

                                                          aCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGGTATTGCTGTTTTATTGTTTGTTTTTGCGAGAt       <  2:16385669/101‑1 (MQ=255)

                                                           cATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATg      <  2:15152713/101‑1 (MQ=255)

                                                           cATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATg      <  1:20434561/101‑1 (MQ=255)

                                                           cATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATg      <  2:11059074/101‑1 (MQ=255)

                                                            aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATgc     >  2:10372063/1‑101 (MQ=255)

                                                            aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATgc     <  2:5012647/101‑1 (MQ=255)

                                                             tGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATgcg    >  2:6544273/1‑101 (MQ=255)

                                                             tGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATgcg    <  2:10700882/101‑1 (MQ=255)

                                                              gATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATgcgc   >  2:14519979/1‑101 (MQ=255)

                                                               aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGcgcc  >  1:9513763/1‑101 (MQ=255)

TCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTAACATTTAGTTTGGCTAAATGTAAAGATATTGCTGTTTTATTGTTTGTTTTTGCGAGATGCGCC  >  NC_000913/688952‑689115

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