BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,774,969	A→G	5.3%	D739D (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,774,969	0	A	G	5.3%	773.8 / 12.8	262	D739D (GAT→GAC)	clfB	Clumping factor B
Reads supporting (aligned to +/- strand): ref base A (115/127); new base G (9/8); total (124/138)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

ATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAG  >  USA300TCH1516_ALE/2774879‑2775059

atcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGaa                                                                                            >  1:4775381/1‑91 (MQ=255)

atcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGaa                                                                                            >  1:4840152/1‑91 (MQ=255)

atcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGaa                                                                                            >  2:3328322/1‑91 (MQ=255)

 tcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAAt                                                                                           >  2:1910719/1‑91 (MQ=255)

 tcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTACCTGAGTCTGACTCGCTGTCTGAGTCTGAAt                                                                                           >  2:171081/1‑91 (MQ=255)

  cGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATc                                                                                          <  2:4924430/91‑1 (MQ=255)

  cGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATc                                                                                          >  1:2680365/1‑91 (MQ=255)

    cTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTCTCTGCATCTGAGTCGCTGTCTGAGTCTGGATCGc                                                                                        >  2:1438119/1‑91 (MQ=255)

    cTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGc                                                                                        <  1:883777/91‑1 (MQ=255)

      aTCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCtg                                                                                      >  2:5260517/1‑90 (MQ=255)

       tCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTat                                                                                     >  1:2544456/1‑91 (MQ=255)

       tCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTat                                                                                     >  1:2587260/1‑91 (MQ=255)

       tCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTat                                                                                     >  1:4957647/1‑91 (MQ=255)

            gTCTGAGCCGCTGTCTGAATCTGAATCGCTGTCTGAGTACGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaa                                                                                >  1:5428560/1‑91 (MQ=255)

             tCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTgtctgaat                                                                               >  2:23748/1‑91 (MQ=255)

             tCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAGTCCGAGTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTgtctgagt                                                                               >  2:5202685/1‑89 (MQ=1)

              cTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatc                                                                              >  1:5321474/1‑91 (MQ=39)

               tGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatct                                                                             <  1:189977/91‑1 (MQ=37)

                  gTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAg                                                                          <  1:1588145/91‑1 (MQ=39)

                  gTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAg                                                                          >  2:4426653/1‑91 (MQ=39)

                   tCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGt                                                                         <  2:2086653/91‑1 (MQ=38)

                    cGCTGTCTGGATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTc                                                                        <  2:5352053/91‑1 (MQ=37)

                    cGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTc                                                                        <  1:3138841/91‑1 (MQ=39)

                    cGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTc                                                                        <  1:4110658/91‑1 (MQ=39)

                      cTGTCTGAATCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGc                                                                      <  2:4201470/91‑1 (MQ=255)

                        gTCTGAATCTGAATCGCTGTCTGAGTCCGAGTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTa                                                                    <  2:119950/91‑1 (MQ=21)

                        gTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTa                                                                    <  1:4056558/91‑1 (MQ=35)

                        gTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTa                                                                    >  1:2620240/1‑91 (MQ=35)

                            gaatctgaatcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctg                                                                <  2:1506570/91‑1 (MQ=33)

                            gaatctgaatcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctg                                                                <  2:157955/91‑1 (MQ=33)

                             aatctgaatcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctga                                                               <  1:4285810/91‑1 (MQ=32)

                              atctgaatcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGtctgag                                                              >  2:4464041/1‑91 (MQ=14)

                                  gactcGCTGTCTGAGTCCGAATCGCTCTCTGAATCTGAGTCGCTGTCTGAGTCTGTATCGCTATCTGAGTCTGAGTCGCTAtctgagtctg                                                          <  2:5509459/88‑1 (MQ=1)

                                    atcGCTGTCTGAGTCCGACTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgag                                                        <  2:4961117/91‑1 (MQ=14)

                                    atcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAAtcgctgtctgagtctgagtcgctgtctgagtctgag                                                        >  2:3182720/1‑91 (MQ=255)

                                    atcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAATCTGAGTCGCTGtctgagtctgag                                                        >  2:2412816/1‑91 (MQ=9)

                                     tcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagt                                                       >  2:183441/1‑91 (MQ=32)

                                      cGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtc                                                      <  1:3965182/91‑1 (MQ=33)

                                      cGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtc                                                      >  1:3187873/1‑91 (MQ=33)

                                         tGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAACCTGAGTCGCTAtctgagtctgagtcgct                                                   >  1:1335161/1‑91 (MQ=18)

                                           tatGAGTCCGAATCGCTATCTGAATATGAGTCGCAGTCTGAGTCTGCATCGCTATACGAATCTGAGTCGCAAtctgagtctgagtcgccgt                                                 >  2:4720090/3‑91 (MQ=11)

                                           tCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgt                                                 >  2:234361/1‑91 (MQ=34)

                                            cTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTCTCTGAATCTGAGTCGCTGtctgagtctgagtcgctgtc                                                >  2:322750/1‑91 (MQ=16)

                                            cTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtc                                                <  1:1239132/91‑1 (MQ=34)

                                            cTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtc                                                <  1:3801410/91‑1 (MQ=34)

                                              gAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtccg                                              >  2:3911488/1‑89 (MQ=14)

                                               aGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGtctgagtctgagtcgctgtctga                                             >  2:3467679/1‑91 (MQ=255)

                                                 tCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTGTCTGAATCTGAGTCGCTAtctgagtctgaatcgctgtctgagt                                           >  2:5000669/1‑91 (MQ=12)

                                                 tCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagt                                           >  1:3494393/1‑91 (MQ=40)

                                                 tCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgcctgagt                                           >  2:2390549/1‑91 (MQ=255)

                                                  ccGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtc                                          <  2:1328747/91‑1 (MQ=40)

                                                   cGACTCGCTATCTGATGGTGAGTCGCTGTCTGTGTCTGAGTCGCTGTCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtct                                         <  2:3067331/91‑1 (MQ=11)

                                                   cGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtct                                         <  1:2990544/91‑1 (MQ=255)

                                                    gAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtctg                                        <  1:1154504/91‑1 (MQ=255)

                                                    gAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtctg                                        >  1:3583572/1‑91 (MQ=255)

                                                    gAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgaatctg                                        >  2:656526/1‑91 (MQ=21)

                                                     aaTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtctga                                       <  1:3931363/91‑1 (MQ=255)

                                                       tCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtctgagt                                     >  1:4960864/1‑91 (MQ=255)

                                                        cGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTAtctgagtctgagtcgctgtctgagtctgagtc                                    >  1:2733883/1‑91 (MQ=255)

                                                             tctgaatctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGACTCTGAGTCGCTATctgagtctgagtcgctgtctgagtctgagtcgctgt                               >  2:5086638/1‑91 (MQ=255)

                                                                gaatctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGagtctgagtcgctgtctgagtctgagtcgctgtctg                            <  1:118797/91‑1 (MQ=255)

                                                                gaatctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGagtctgagtcgctgtctgagtctgagtcgctgtctg                            <  1:2834337/91‑1 (MQ=255)

                                                                gaatctgaGTAGCTATCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGagtctgcgtcgctgtctgactctgagtcgctgtctg                            <  1:5168754/91‑1 (MQ=14)

                                                                 aatctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAgtctgagtcgctgtctgagtctgagtcgctgtctga                           >  2:1034106/1‑91 (MQ=40)

                                                                  gtctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGtctgagtcgctgtctgagtctgagtcgctgtctgag                          <  1:2310150/90‑1 (MQ=39)

                                                                  atctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGtctgagtcgctgtctgagtctgagtcgctgtctgag                          >  2:4095899/1‑91 (MQ=39)

                                                                  atctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGtctgagtcgctgtctgagtctgagtcgctgtctgag                          <  2:3432490/91‑1 (MQ=40)

                                                                   tctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTctgagtcgctgtctgagtctgagtcgctgtctgagt                         >  2:3098467/1‑91 (MQ=39)

                                                                   tctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTctgagtcgctgtctgagtctgagtcgctgtctgagt                         >  2:893127/1‑91 (MQ=40)

                                                                   tctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTctgagtcgctgtctgagtctgagtcgctgtctgagt                         <  1:4556206/91‑1 (MQ=255)

                                                                      gaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGagtcgctgtctgagtctgagtcgctgtctgagtctg                      <  1:3071814/91‑1 (MQ=255)

                                                                      gaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGagtcgctgtctgagtctgagtcgctgtctgagtctg                      <  1:4426721/91‑1 (MQ=255)

                                                                      gaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGagtcgctgtctgagtctgagtcgctgtctgagtctg                      <  1:5546751/91‑1 (MQ=255)

                                                                        gTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGACTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGaa                    >  2:2268375/1‑91 (MQ=255)

                                                                        gTCGCTATCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGaa                    <  1:5644014/91‑1 (MQ=255)

                                                                            cTGTCTGAGTCTGTATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGc                <  2:365970/91‑1 (MQ=255)

                                                                            cTGTCTGAGTCTGAATCGCTATTTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCtgtctgagtctgagtcgctgtctgagtctgagtcgc                >  2:4311285/1‑91 (MQ=255)

                                                                            cTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGc                <  1:4515876/91‑1 (MQ=255)

                                                                            cTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGc                <  1:732113/91‑1 (MQ=255)

                                                                              gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCtg              >  1:244477/1‑90 (MQ=255)

                                                                              gtctgagtctgaATCGCTATCTGAATCTGAATCGCTAACTGAGTCTGAGTCACTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTa              >  2:561268/1‑91 (MQ=255)

                                                                                ctgagtctgagtCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatc            <  2:3253577/91‑1 (MQ=255)

                                                                                ctgagtctgaATCGCTGTCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatc            <  2:4826184/91‑1 (MQ=255)

                                                                                ctgagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatc            <  2:454422/91‑1 (MQ=255)

                                                                                ctgagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatc            >  1:4192466/1‑91 (MQ=255)

                                                                                ctgagtctgaATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAAGCGCtctc            >  2:285321/1‑91 (MQ=11)

                                                                                 ggagtctgaGTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatct           <  2:2307847/90‑1 (MQ=255)

                                                                                  gagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGCCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctg          >  1:2058318/1‑91 (MQ=255)

                                                                                   agtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctga         >  1:2351612/1‑91 (MQ=255)

                                                                                    gtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaa        <  1:669660/91‑1 (MQ=255)

                                                                                     tctgaTTCGCTATCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaat       <  1:3658576/91‑1 (MQ=255)

                                                                                     tctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctggtt       >  2:1655014/1‑88 (MQ=255)

                                                                                      ctgaATCGCTGTCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatc      <  2:785718/91‑1 (MQ=255)

                                                                                      ctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatc      <  1:804104/91‑1 (MQ=255)

                                                                                      ctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatc      >  1:1437152/1‑91 (MQ=255)

                                                                                       tgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTgtctgagtct     >  2:2977242/1‑91 (MQ=255)

                                                                                         aaTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatctga   >  2:2889248/1‑91 (MQ=255)

                                                                                         aaTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaatctga   <  1:549102/91‑1 (MQ=255)

                                                                                          aTCGCTATCTGAATCTGAGTCGCTATCTGAGTTTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTAtctgagtctgag  >  2:1159841/1‑91 (MQ=255)

                                                                                          aTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAg  >  1:1074019/1‑91 (MQ=255)

ATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAG  >  USA300TCH1516_ALE/2774879‑2775059

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: `‑`