BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,775,053	A→G	6.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	6.5%	469.5 / 4.4	183	D711D (GAT→GAC)	clfB	Clumping factor B
Reads supporting (aligned to +/- strand): ref base A (82/83); new base G (9/5); total (93/90)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

GAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGAGTCACTGT  >  USA300TCH1516_ALE/2774913‑2775192

gaatcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTatctgaa                                                                                                                                             >  1:1105958/1‑141 (MQ=255)

  atcGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTatctgaatc                                                                                                                                           <  1:1505620/141‑1 (MQ=255)

    cGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTatctgaatctg                                                                                                                                         >  3:659103/1‑141 (MQ=255)

    cGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCACTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCACTatcttaatctg                                                                                                                                         >  4:664868/1‑141 (MQ=255)

     gCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAATCTGAATCGCTgtctgagtctga                                                                                                                                        >  1:601987/1‑141 (MQ=255)

      cTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAg                                                                                                                                       <  3:232371/141‑1 (MQ=255)

      cTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAg                                                                                                                                       <  1:20854/141‑1 (MQ=255)

      cTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATAGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAg                                                                                                                                       >  1:1620473/1‑141 (MQ=255)

        gTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTc                                                                                                                                     >  3:605851/1‑141 (MQ=255)

         tCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCg                                                                                                                                    >  2:466303/1‑141 (MQ=255)

         tCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCg                                                                                                                                    >  1:15219/1‑141 (MQ=255)

          cTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGc                                                                                                                                   >  1:1474313/1‑141 (MQ=255)

            gAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGACTCTGAACCGCTATCTGATTCTGCGTCGCTg                                                                                                                                 >  1:769642/1‑141 (MQ=255)

              gTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtc                                                                                                                               >  3:504106/1‑141 (MQ=255)

               tCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtct                                                                                                                              <  3:190990/141‑1 (MQ=255)

                    aTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTatctgaatctg                                                                                                                                         >  3:210684/1‑125 (MQ=255)

                    aTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTatctgaatctg                                                                                                                                         <  4:210684/125‑1 (MQ=255)

                         tATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTctgagtctgagt                                                                                                                    <  1:1058145/141‑1 (MQ=255)

                           tctgaatctgaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcg                                                                                                                  >  4:490823/1‑141 (MQ=255)

                                 tctgaGTGGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtct                                                                                                            >  1:1312274/1‑141 (MQ=255)

                                    gaGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaa                                                                                                         >  1:799309/1‑141 (MQ=255)

                                      gTCGCTCTCTGAGTCTGACTCGCTATCTGAATCGGAGTCGCTATCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatc                                                                                                       <  2:1225388/141‑1 (MQ=14)

                                      atCGCTGTCCGAATCTGAGTCGCTATCTGAATCTGAGTCACTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatc                                                                                                       >  1:333459/2‑141 (MQ=11)

                                       tCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGACTCGCTGTCTGAGTCTGAATCGCTGTCTGAGtctgagtcgctgtctgaaact                                                                                                      >  2:574423/1‑141 (MQ=11)

                                        cGCTGTCTGAGTCTGAATCGCTGTCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctg                                                                                                     <  2:1151827/141‑1 (MQ=255)

                                        cGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctg                                                                                                     >  1:1391297/1‑141 (MQ=255)

                                         gCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctga                                                                                                    >  1:889620/1‑141 (MQ=255)

                                          cTGTCTGAGTCTGAATCGCTGTCTGAATCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgag                                                                                                   >  2:654806/1‑141 (MQ=255)

                                            ttctgagtctgagTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtc                                                                                                 <  2:1457486/140‑1 (MQ=255)

                                             tctgagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTGTCTGAATCTGAGTCGctgtctgagtctgagtcgctgtctgagtctgagtcg                                                                                                >  2:745646/1‑141 (MQ=255)

                                             tctgagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAGTCTGAGTCGctgtctgagtctgagtcgctgtctgagtctgagtcg                                                                                                >  4:463147/1‑141 (MQ=255)

                                                gagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctg                                                                                             <  1:1167372/141‑1 (MQ=255)

                                                gagtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctg                                                                                             <  1:97161/141‑1 (MQ=255)

                                                    ctgaGTCGCTATCTGAATCTGAGTCACTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctgtctg                                                                                         <  2:333459/141‑1 (MQ=255)

                                                    ctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctgtctg                                                                                         <  2:1039931/141‑1 (MQ=255)

                                                      gaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGtctgagtcgctgtctgaatctgagtcgctgtctgaa                                                                                       >  4:234950/1‑141 (MQ=255)

                                                        aTCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCtgagtcgctgtctgaatctgagtcgctgtctgaatc                                                                                     <  2:987637/141‑1 (MQ=255)

                                                           gCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAgtcgctgtctgaatctgagtcgctgtctgaatctga                                                                                  <  3:7353/141‑1 (MQ=255)

                                                            cTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGtctgaatctgaa                                                                                 <  2:799309/141‑1 (MQ=255)

                                                              gtctgagtctgagTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAGTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCtgaatctgaatc                                                                               <  2:1615834/140‑1 (MQ=255)

                                                              atctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCtgaatctgaatc                                                                               <  2:154475/141‑1 (MQ=255)

                                                               tctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCg                                                                              >  4:191384/1‑141 (MQ=255)

                                                                ctgaatctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCTCTGTCTGAATCTGTATCGCTGTCTGAGTCTGCATCGc                                                                             >  4:576023/1‑141 (MQ=14)

                                                                    atctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtc                                                                         <  2:42305/141‑1 (MQ=255)

                                                                      ctgaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTgtctgagtc                                                                                                                         >  3:417129/1‑91 (MQ=255)

                                                                        gaGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTAAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgag                                                                     >  2:898379/1‑141 (MQ=255)

                                                                            cGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTgtctgagtctg                                                                 >  3:522994/1‑141 (MQ=255)

                                                                                aTCGGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatc                                                             <  2:63813/141‑1 (MQ=255)

                                                                                   tgagtctgagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtcTGAGTCTGAATCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAATCGTTGTCTGCATCTGAATCACTATCTGCATCTGAATCGCTCtctgcctctgtatcgct                                                          >  2:1575421/1‑141 (MQ=14)

                                                                                   tgagtctgagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtcTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTATCTGAATCTGAATCGCTGtctgagtctgaatcact                                                          >  2:1256099/1‑141 (MQ=255)

                                                                                    gagtctgagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctGAGTCTGAATCGCTATCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGATTCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgattctgaatcgcta                                                         >  2:919468/1‑141 (MQ=14)

                                                                                      gtctgagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgaGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatc                                                       >  1:377791/1‑141 (MQ=255)

                                                                                        ctgagtccctgtctgagtctgattcgct‑‑‑‑‑‑gtctgagtCTGAGTCTCTATCTTCATCTGAGTCGCTGTCTGAGTCTGAGTCGCTATCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctg                                                     <  4:35059/141‑1 (MQ=14)

                                                                                         tgagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtcGGAGTCGCTATCTGAATCTGAGACTCTGTCTGAGTCTGAGTCGCTGTCTGCATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctgtctga                                                    >  1:1448838/1‑141 (MQ=14)

                                                                                          gagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctGAATCGCTGTCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgag                                                   <  3:333313/141‑1 (MQ=255)

                                                                                          gagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgag                                                   <  1:176323/141‑1 (MQ=255)

                                                                                          gagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagcctGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgag                                                   >  1:1331811/1‑141 (MQ=255)

                                                                                          gagtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagcctGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgag                                                   <  2:1331811/141‑1 (MQ=255)

                                                                                           actctctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagt                                                  <  2:140888/136‑1 (MQ=255)

                                                                                            gtcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATGTGAATCGCTTgctgagcctgaatctctagctgattg                                                 >  4:147141/1‑138 (MQ=17)

                                                                                             tcgctgtctgagtctgagtcgctgtctgagtctgagTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatct                                                      >  2:962085/1‑141 (MQ=17)

                                                                                             tcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGTGTCTGAGTCGCTGTCTGCATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtct                                                >  4:535253/1‑141 (MQ=255)

                                                                                             tcgctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtc                                                                                                             >  1:345489/1‑80 (MQ=255)

                                                                                             tcgctgtcggattctgattcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtct                                                <  2:751054/141‑1 (MQ=255)

                                                                                               gctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctga                                              <  2:91260/141‑1 (MQ=255)

                                                                                                ctgtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctgaa                                             >  4:438985/1‑141 (MQ=255)

                                                                                                  gtctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctgaatc                                           <  1:701856/141‑1 (MQ=255)

                                                                                                   tctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctgaatcg                                          <  2:740892/141‑1 (MQ=255)

                                                                                                   tctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctgaatcg                                          >  2:192323/1‑141 (MQ=255)

                                                                                                    ctgagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaCTCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctgaatcgc                                         <  1:474350/141‑1 (MQ=255)

                                                                                                      gagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTGTCTGAATCTGAGTCTCTGTCTGCGTCTTAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGCGTCGCTGtctgagtctgggtcactgtctgggtccgagtcgaga                                       >  4:112002/1‑138 (MQ=11)

                                                                                                      gagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgagtcgctgtc                                                                                           <  2:693257/89‑1 (MQ=255)

                                                                                                      gagtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatcgctatctgagtctgaatcgcta                                       >  1:1639540/1‑141 (MQ=255)

                                                                                                       agtctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTctgagtctgaatcgctatctgagtctgaatcgctat                                      >  3:114832/1‑141 (MQ=255)

                                                                                                         tctgagtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTgagtctgaatcgctatctgagtctgaatcgctatct                                    <  3:362784/141‑1 (MQ=255)

                                                                                                          ctgactcgct‑‑‑‑‑‑gtctgattctgaATCGCTATCTGATTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTATCTGAATCTGAATCGCTATCTGagtctgaatcgctatctgagtctgaatcgctatctg                                   <  4:369894/141‑1 (MQ=14)

                                                                                                              gtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTGTCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCtgaatcgctatctgagtctgaatcgctatctgagtc                               <  1:55474/141‑1 (MQ=255)

                                                                                                              gtcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCtgaatcgctatctgagtctgaatcgctatctgagtc                               >  4:274804/1‑141 (MQ=255)

                                                                                                              gtcgct‑‑‑‑‑‑gtctgagtccgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTCAGTCGCTGACTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCtgaatcgctatctgagtctgaatcgctatctgagtc                               <  3:577926/141‑1 (MQ=255)

                                                                                                               tcgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTATGAATCTGCATCGCTGTCTGAGTCTgcatagctatccgagaatgcatcgatatatgagact                              >  1:908695/1‑141 (MQ=17)

                                                                                                                cgct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGaatcgctatctgagtctgaatcgctatctgagtctg                             >  2:128327/1‑141 (MQ=255)

                                                                                                                 gct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgag                                                                                                   >  2:618270/1‑70 (MQ=255)

                                                                                                                 gct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGtctgagtcgctgtctgaatctgag                                                                                                   <  1:618270/70‑1 (MQ=255)

                                                                                                                  ct‑‑‑‑‑‑gtctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAAtcgctatctgagtctgaatcgctatctgagtctgaa                           <  1:1016060/141‑1 (MQ=255)

                                                                                                                           tctgagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAATCGCTATCTGACTCTGAATCa                        >  2:664652/1‑141 (MQ=255)

                                                                                                                           tctgagtctgaATCGCTATCTGAATCTGAGTCGCTGCCTGAGTCTGAGTCGCAGTATGCATGGGAGTGGCTGTCTTACACCGAATCGCCGCCTGAAGGAGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCa                        >  1:745135/1‑141 (MQ=17)

                                                                                                                              gagtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTg                     <  1:48692/141‑1 (MQ=255)

                                                                                                                                gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTc                   >  3:446044/1‑141 (MQ=255)

                                                                                                                                gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTc                   <  3:148517/141‑1 (MQ=255)

                                                                                                                                gtctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTc                   <  2:544336/141‑1 (MQ=255)

                                                                                                                                gtctgaATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAATCGCTGTCCGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTc                   <  2:519885/141‑1 (MQ=255)

                                                                                                                                  ctgaATTGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTTGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTg                 <  2:1583949/141‑1 (MQ=255)

                                                                                                                                  ctgaATCGCTATCTGAATCTGAGTCGCTGTCTGATTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTg                 <  1:1307193/141‑1 (MQ=255)

                                                                                                                                  ctgaATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCCCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTg                 <  1:1303050/141‑1 (MQ=255)

                                                                                                                                         gCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGa          <  1:677066/141‑1 (MQ=255)

                                                                                                                                             tctgaatctgaGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGAGTCa      >  4:409282/1‑141 (MQ=255)

                                                                                                                                                gaatctgaGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTGTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGAGTCACTg   >  2:626015/1‑141 (MQ=255)

                                                                                                                                                gaatctgaGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGAGTCACTg   >  2:406161/1‑141 (MQ=255)

                                                                                                                                                 aatctgaGTCGCTGTCNGAGTCTGAGTCGNTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGAGTCCCTGt  >  3:458128/1‑141 (MQ=255)

                                                                                                                                                  atctgaGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGtctgagtctgaatc                                                             >  2:1500849/1‑81 (MQ=255)

GAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCTGAGTCTGAGTCGCT‑‑‑‑‑‑GTCTGAGTCTGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAGTCGCTGTCTGAATCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAATCACTGTCTGAGTCCGAGTCACTGT  >  USA300TCH1516_ALE/2774913‑2775192

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