BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,774,846	A→G	5.4%	S780S (AGT→AGC)	clfB ←	Clumping factor B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,774,846	0	A	G	5.4%	487.9 / 4.7	163	S780S (AGT→AGC)	clfB	Clumping factor B
Reads supporting (aligned to +/- strand): ref base A (73/80); new base G (8/2); total (81/82)
Fisher's exact test for biased strand distribution p-value = 5.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

GCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCG  >  USA300TCH1516_ALE/2774756‑2774936

 cTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCAc                                                                                           <  1:1455653/91‑1 (MQ=255)

  tATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACt                                                                                          <  2:3824648/91‑1 (MQ=255)

   aTCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGCTg                                                                                         >  2:3527913/1‑91 (MQ=255)

   aTCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTg                                                                                         <  2:127108/91‑1 (MQ=255)

   aTCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTg                                                                                         <  2:203994/91‑1 (MQ=255)

         gTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGTATCTGAATCACTGTCTGaa                                                                                   >  2:3126036/1‑91 (MQ=255)

         gTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGaa                                                                                   >  1:283413/1‑91 (MQ=255)

         gTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGaa                                                                                   >  1:1673347/1‑91 (MQ=255)

         gTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGaa                                                                                   >  2:1930378/1‑91 (MQ=255)

          tCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAAt                                                                                  >  1:2555771/1‑91 (MQ=255)

           cTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATc                                                                                 >  2:2885522/1‑91 (MQ=255)

           cTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATc                                                                                 <  1:3580962/91‑1 (MQ=255)

             gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                               <  1:3127739/91‑1 (MQ=255)

             gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                               >  2:2290780/1‑91 (MQ=255)

             gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                               >  2:2290767/1‑91 (MQ=255)

             gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                               <  1:343229/91‑1 (MQ=255)

             gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                               <  1:3453510/91‑1 (MQ=255)

             gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCg                                                                               >  1:631783/1‑91 (MQ=255)

               gTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGag                                                                             >  2:1138296/1‑90 (MQ=255)

                tCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCCGAAt                                                                            >  2:3388048/1‑91 (MQ=255)

                tCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAAt                                                                            >  1:1371535/1‑91 (MQ=255)

                 cGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATc                                                                           >  1:1629552/1‑91 (MQ=255)

                 cGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATc                                                                           >  1:2183808/1‑91 (MQ=255)

                   cTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGc                                                                         <  2:1641729/91‑1 (MQ=255)

                    tATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCt                                                                        <  1:2857739/91‑1 (MQ=255)

                     aTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTa                                                                       <  1:2338305/91‑1 (MQ=255)

                        tGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatct                                                                    <  1:24171/91‑1 (MQ=255)

                        tGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatct                                                                    <  1:1725711/91‑1 (MQ=255)

                        tGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatct                                                                    <  1:1137331/91‑1 (MQ=255)

                         gAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatccg                                                                   >  1:2652227/1‑89 (MQ=255)

                           gTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaa                                                                 <  1:2104651/91‑1 (MQ=255)

                           gTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaa                                                                 >  1:1829526/1‑91 (MQ=255)

                            tCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaat                                                                <  1:2228850/91‑1 (MQ=255)

                             cTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaatc                                                               <  2:1473827/91‑1 (MQ=255)

                                aaTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatctgaatctga                                                            <  2:657710/91‑1 (MQ=255)

                                aaTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaatctga                                                            <  1:1849112/91‑1 (MQ=255)

                                 aTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTAtctgaatctgaa                                                           <  2:2348549/91‑1 (MQ=255)

                                  tCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATctgaatctgaat                                                          >  1:1314639/1‑91 (MQ=255)

                                   cGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAGTCTGAATCACTGTCTGAATCCGAATCGCTGTCtggatctgaatc                                                         >  2:2441366/1‑91 (MQ=255)

                                    gCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCg                                                        >  1:500854/1‑91 (MQ=255)

                                     cTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGc                                                       >  1:1233647/1‑91 (MQ=255)

                                     cTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGc                                                       <  2:1498749/91‑1 (MQ=255)

                                       gtctgagtctgaGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTa                                                     <  1:1025483/91‑1 (MQ=255)

                                       gtctgagtctgaGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTa                                                     >  1:859860/1‑91 (MQ=255)

                                       gtctgagtctgaGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTa                                                     <  1:1428846/91‑1 (MQ=255)

                                       gtctgagtctgaGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTa                                                     <  2:3305005/91‑1 (MQ=255)

                                       gtctgagtctgaGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCAGAATCGCTa                                                     >  1:3482664/1‑91 (MQ=255)

                                       gtctgactctgaGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTa                                                     <  1:1780586/91‑1 (MQ=255)

                                        tctgagtctgagTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTAt                                                    >  1:2548619/1‑91 (MQ=255)

                                         ctgagtctgagtCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATc                                                   <  2:119331/91‑1 (MQ=40)

                                         ctgagtctgagtCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATc                                                   <  2:1107559/91‑1 (MQ=255)

                                          tgagtctgagtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATcc                                                  >  1:3141404/1‑91 (MQ=255)

                                              tctgagtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTATCTGAATCCGGATCGCTATCTGCATCTGGCTCGCTATCCGAGt                                              >  2:904370/1‑91 (MQ=14)

                                               ctgagtcGCTGACCGAGTCTGGGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTc                                             <  2:827660/91‑1 (MQ=255)

                                               ctgagtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAGACACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTc                                             <  2:828020/91‑1 (MQ=255)

                                               ctgagtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTc                                             >  1:2476600/1‑91 (MQ=255)

                                                   gtcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAGTCTGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAg                                         <  2:3055166/91‑1 (MQ=255)

                                                    tcGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCCGAATCGCTATCTGAATCTGAAACGCTATCCGAGTCTGAGt                                        >  2:402929/1‑91 (MQ=255)

                                                     cGCTATCCGAGTCTGCGTCGCTGTCTGAGTCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTc                                       <  2:283413/91‑1 (MQ=14)

                                                     cGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTc                                       >  1:190568/1‑91 (MQ=255)

                                                         aTCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAATCCTAATCGCTATCTGACTCTGCGTCGCTATCTGAGTCTGAGTCGCTg                                   >  2:3475317/1‑91 (MQ=11)

                                                         aTCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTg                                   <  1:2290780/91‑1 (MQ=255)

                                                         aTCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTg                                   <  1:3764169/91‑1 (MQ=255)

                                                         aTCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTg                                   <  1:2290767/91‑1 (MQ=255)

                                                           ccGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtc                                 <  2:204975/91‑1 (MQ=255)

                                                            cGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtct                                >  1:3518369/1‑91 (MQ=255)

                                                             gAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctg                               <  2:3045969/91‑1 (MQ=255)

                                                             gAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctg                               <  2:2143922/91‑1 (MQ=255)

                                                             gAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctg                               >  2:1720529/1‑91 (MQ=255)

                                                               gTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaa                             >  2:2908342/1‑91 (MQ=255)

                                                               gTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaa                             >  1:1393337/1‑91 (MQ=255)

                                                                tCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaat                            <  1:2022328/91‑1 (MQ=255)

                                                                 cTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAGTCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatc                           <  2:2423004/91‑1 (MQ=255)

                                                                 cTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatc                           <  1:2438511/91‑1 (MQ=255)

                                                                  tGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatct                          <  1:3497694/91‑1 (MQ=255)

                                                                   gtgTCGCTGTCCGAATCTGAGTCACTGTCTGAATCCGAATCGCTCTCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatctg                         <  2:244539/89‑1 (MQ=17)

                                                                   gAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatctg                         <  1:2914508/91‑1 (MQ=255)

                                                                   gAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCACTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatctg                         <  2:1950892/91‑1 (MQ=255)

                                                                        gCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCg                    >  1:1034016/1‑91 (MQ=255)

                                                                          tGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCt                  <  1:912386/91‑1 (MQ=255)

                                                                          tGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCt                  >  1:1321763/1‑91 (MQ=255)

                                                                             ctgaatctgaatCACTGTCTGAATCCGAATCGCTTTCTGAATCTGAATCACTTCCTGAGTGCGACGCGCTGGCTGAGTGTGAGTCGCTGtg               >  2:2670349/1‑90 (MQ=2)

                                                                             ctgaatctgaatCACTGTCTGAATCCGAATCGCTAACTGAATCAGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCGGTc               >  2:2993924/1‑91 (MQ=255)

                                                                                aatctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGa            <  1:2080985/91‑1 (MQ=255)

                                                                                 atctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTACCCGAGTCTGAGTCGCCGTCTGAGTCTGAGTCGCTGTCTGaa           >  2:3416135/1‑90 (MQ=14)

                                                                                  tctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGt          <  1:552190/91‑1 (MQ=255)

                                                                                  tctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGt          <  2:609357/91‑1 (MQ=255)

                                                                                   ctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTc         <  1:3165125/91‑1 (MQ=255)

                                                                                     gaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCg       <  1:2908342/91‑1 (MQ=255)

                                                                                     gaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCg       >  1:2222074/1‑91 (MQ=255)

                                                                                       atcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGaa     >  1:3375967/1‑91 (MQ=255)

                                                                                        tcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAAt    >  2:1363419/1‑91 (MQ=255)

                                                                                        tcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAAt    >  1:765357/1‑91 (MQ=255)

                                                                                        tcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAAt    >  2:1318703/1‑91 (MQ=255)

                                                                                        tcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAAt    >  1:2812985/1‑91 (MQ=255)

                                                                                        tcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAAt    >  2:2722093/1‑91 (MQ=255)

                                                                                         cACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGCGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATc   >  1:977335/1‑91 (MQ=255)

                                                                                         cACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATc   >  1:3512202/1‑91 (MQ=255)

                                                                                         cACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATc   <  1:3288943/91‑1 (MQ=255)

                                                                                          aCTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCg  <  2:1315555/91‑1 (MQ=255)

GCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCG  >  USA300TCH1516_ALE/2774756‑2774936

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