BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,775,053	A→G	5.5%	D711D (GAT→GAC)	clfB ←	Clumping factor B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,775,053	0	A	G	5.5%	474.5 / 2.4	169	D711D (GAT→GAC)	clfB	Clumping factor B
Reads supporting (aligned to +/- strand): ref base A (87/70); new base G (6/6); total (93/76)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

GTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGA  >  USA300TCH1516_ALE/2774963‑2775142

gtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGCGTCTGACTCGCTgtctgaa                                                                                           >  2:2745282/1‑91 (MQ=255)

gtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaa                                                                                           >  1:3024900/1‑91 (MQ=255)

 tctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaat                                                                                          <  1:393373/91‑1 (MQ=255)

    gactCGCTATCTGCATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatctg                                                                                       <  2:3518369/88‑1 (MQ=255)

    gaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatctg                                                                                       <  2:3518332/91‑1 (MQ=255)

      aTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAg                                                                                     <  1:1332017/91‑1 (MQ=255)

       tCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGt                                                                                    >  1:2392377/1‑91 (MQ=255)

        cGCTATCTGAACCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGGGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTc                                                                                   <  1:2884830/91‑1 (MQ=255)

         gCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAGTCg                                                                                  >  2:2394997/1‑91 (MQ=255)

           tATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCt                                                                                <  2:1810380/91‑1 (MQ=255)

            gtctgagtctgagTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTg                                                                               <  2:260923/90‑1 (MQ=255)

            atctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTg                                                                               <  1:1364816/91‑1 (MQ=255)

            atctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTg                                                                               <  1:3405397/91‑1 (MQ=255)

             tctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAATCTGAGTCGCTgt                                                                              >  2:151503/1‑91 (MQ=255)

             tctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgt                                                                              >  2:2856204/1‑91 (MQ=255)

             tctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgt                                                                              >  2:3435433/1‑91 (MQ=255)

              ctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtc                                                                             >  2:3039493/1‑91 (MQ=255)

                gaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctg                                                                           >  1:2704451/1‑91 (MQ=255)

                  atctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgag                                                                         <  1:1921838/91‑1 (MQ=255)

                  atctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgag                                                                         <  1:2483427/91‑1 (MQ=255)

                  atctgaGTCGCTATCTGAGTCTGAGACGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgag                                                                         <  1:3601987/91‑1 (MQ=255)

                   tctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgagt                                                                        <  2:2431820/91‑1 (MQ=255)

                   tctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgagt                                                                        >  1:653713/1‑91 (MQ=255)

                    ctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgagtc                                                                       <  1:2395314/91‑1 (MQ=255)

                     tgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgagtct                                                                      >  1:2100993/1‑91 (MQ=255)

                      gagtcgctatctgagtctgagtcgctatctgagtctGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgagtctg                                                                     <  2:3831556/91‑1 (MQ=255)

                       cctCGCTATCTGAGTCCGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTATCTGAATCTGAGTCGCTgtctgagtctga                                                                    <  2:1359308/89‑1 (MQ=11)

                        gTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGtctgagtctgag                                                                   >  1:592472/1‑91 (MQ=255)

                         tCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTctgagtctgagt                                                                  >  1:3533641/1‑91 (MQ=255)

                         tCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTctgagtctgagt                                                                  >  2:74513/1‑91 (MQ=255)

                         tCGCTATCTGAGACTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAGTCGCTGCgtgggtctgcgt                                                                  >  1:1426615/1‑91 (MQ=14)

                          cGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCtgagtctgagtc                                                                 <  1:2535381/91‑1 (MQ=255)

                           gCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcg                                                                <  1:904370/91‑1 (MQ=255)

                            cTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgc                                                               <  1:350229/91‑1 (MQ=255)

                            cTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgc                                                               >  2:567332/1‑91 (MQ=255)

                            cTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgc                                                               >  1:3685296/1‑91 (MQ=255)

                              ttctgagtctgagtcgctgtctgagtctgagtcgctgTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctg                                                             <  2:269948/90‑1 (MQ=255)

                              aTCTGAGTCTGAGTAGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTAtctgagtctgagtcgctgtctgagtctgagtcgctg                                                             >  2:3297476/1‑91 (MQ=255)

                               tctgagtctgagtcgctgtctgagtctgagtcgctgTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgt                                                            >  2:2013450/1‑91 (MQ=255)

                                ctgagtctgagtcgctgtctgagtctgagtcgctgtCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtc                                                           <  2:1094770/91‑1 (MQ=255)

                                  gagtctgagtcgctgtctgagtctgagtcgctgtctGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctg                                                         >  1:2012220/1‑91 (MQ=255)

                                   agtctgagtcgctatctgaatctgagtcgctgtctgAATCTGAATCGCTATCTGAATCTGAATCGCTGTCTGAAtctgagtcgctgtctga                                                        <  1:2830436/91‑1 (MQ=11)

                                    gtctgagtcgctgtctgagtctgagtcgctgtctgaGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaa                                                       >  2:2935115/1‑91 (MQ=255)

                                    atctgggtcgctgtctgagtctgagtcgctgtctgagTCTGAGTCGCTATCTGAgtctgagtcgctgtctgagtctgagtcgctgtctgag                                                       <  2:1925356/90‑2 (MQ=11)

                                      ctgagtcgctgtctgagtctgagtcgctgtctgagtCTGAATCGTTATCTGAAGCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatc                                                     <  1:3100153/91‑1 (MQ=255)

                                      ctgagtcgctgtctgagtctgagtcgctgtctgagtCTGAATCGCTATCTGAATCtgagtcgctgtctgagtctgagtcgctgtctgagtc                                                     >  2:850143/1‑91 (MQ=255)

                                        gagtcgctgtctgagtctgagtcgctgtctgagtctGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctg                                                   <  2:1393337/91‑1 (MQ=255)

                                         agtcgctgtctgagtctgagtcgctgtctgagtctgAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctga                                                  >  2:2386431/1‑91 (MQ=255)

                                           tcgctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGTATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagt                                                <  1:402929/91‑1 (MQ=255)

                                           tcgctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagt                                                >  1:2858258/1‑91 (MQ=255)

                                           tcgctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagt                                                >  1:1962602/1‑91 (MQ=255)

                                             gctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcg                                              <  1:1318703/91‑1 (MQ=255)

                                             gctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcg                                              >  1:2955977/1‑91 (MQ=255)

                                              ctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgc                                             <  1:728981/91‑1 (MQ=255)

                                              ctgtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgc                                             >  1:2872807/1‑91 (MQ=255)

                                                gtctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctg                                           >  1:3762641/1‑91 (MQ=255)

                                                 tctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctgt                                          <  2:1049197/91‑1 (MQ=255)

                                                  ctgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctatctgagtctgagtcgctgtc                                         >  2:1467887/1‑91 (MQ=255)

                                                   tgagtctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGtctgagtcgctgtctgaatctgaatcactgtct                                        <  1:1148280/91‑1 (MQ=11)

                                                        ctgagtcgctgtctgagtctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCtgagtcgctgtctgaatctgagtcgctgtctgaatc                                   <  2:3472736/91‑1 (MQ=255)

                                                        ctgagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCtgagtcgctgtctgaatctgagtcgctgtctgaatc                                   <  1:2944659/91‑1 (MQ=255)

                                                          gagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGagtcgctgtctgaatctgagtcgctgtctgaatctg                                 >  1:1017317/1‑91 (MQ=255)

                                                          gagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGagtcgctgtctgaatctgagtcgctgtctgaatctg                                 <  2:3018279/91‑1 (MQ=255)

                                                          gagtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGagtcgctgtctgaatctgagtcgctgtctgaatctg                                 <  1:3679249/91‑1 (MQ=255)

                                                            gtcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGtctgaatctgaa                               >  2:1593271/1‑91 (MQ=255)

                                                             tcgctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTctgaatctgaat                              <  1:1065597/91‑1 (MQ=255)

                                                               gctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCg                            <  1:3449208/91‑1 (MQ=255)

                                                               gctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCACTGTCTGATTCTGAATCg                            >  2:646161/1‑91 (MQ=255)

                                                                ctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGc                           <  1:1124289/91‑1 (MQ=255)

                                                                ctgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGc                           <  1:2433607/91‑1 (MQ=255)

                                                                 tgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctg                                                   <  2:3332833/66‑1 (MQ=255)

                                                                 tgtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctg                                                   <  2:3054067/66‑1 (MQ=255)

                                                                  gtctgagtctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTg                         <  2:2519979/91‑1 (MQ=255)

                                                                  gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTg                         <  1:152705/91‑1 (MQ=255)

                                                                   tctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgt                        >  2:2315013/1‑91 (MQ=255)

                                                                   tctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCCGAATCTGAATCGCTgt                        >  2:3709371/1‑91 (MQ=255)

                                                                    ctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtc                       <  2:1655637/91‑1 (MQ=255)

                                                                     tgagtcagaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtct                      >  1:3658040/1‑91 (MQ=255)

                                                                      gagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctg                     >  1:2926591/1‑91 (MQ=255)

                                                                        gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgag                   >  2:3673877/1‑91 (MQ=255)

                                                                        gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgag                   >  1:551428/1‑91 (MQ=255)

                                                                        gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgag                   >  2:1221800/1‑91 (MQ=255)

                                                                        gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgag                   <  1:791436/91‑1 (MQ=255)

                                                                         tctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagt                  >  1:34958/1‑91 (MQ=255)

                                                                         tctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagt                  >  1:2611853/1‑91 (MQ=255)

                                                                         tctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagt                  >  1:34557/1‑91 (MQ=255)

                                                                            gaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagtctg               <  1:1526024/91‑1 (MQ=255)

                                                                             aaTCGCTCTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagtctga              <  2:3077430/91‑1 (MQ=255)

                                                                             aaTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGtctgaatctgaa                               >  1:1563384/1‑74 (MQ=255)

                                                                             aaTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGtctgaatctgaa                               <  2:1563384/74‑1 (MQ=255)

                                                                             aaTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagtctgg              >  2:2106367/1‑90 (MQ=255)

                                                                              gtCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaa             <  2:967191/90‑1 (MQ=255)

                                                                                cGCTATCTGAGTCTGGGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatc           <  2:1034016/91‑1 (MQ=255)

                                                                                 gCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgtatcg          >  1:1550468/1‑91 (MQ=255)

                                                                                  cTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgc         >  1:1718055/1‑91 (MQ=255)

                                                                                    atctgaatctgaGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgcta       >  1:1880727/1‑91 (MQ=255)

                                                                                    atctgaatctgaGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgcta       >  2:133711/1‑91 (MQ=255)

                                                                                     tctgagtctgagtcgctgtctgagtctgagtcgctgTCTGAATCTGAATCGCTGTCTGAATCTGAATCGCTGtctgattctgatccgctga      >  2:1718058/1‑89 (MQ=1)

                                                                                          atctgaGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctgtc                                         <  2:3207646/51‑1 (MQ=38)

                                                                                          atctgaGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctgtc                                         >  1:3207646/1‑51 (MQ=38)

GTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGA  >  USA300TCH1516_ALE/2774963‑2775142

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