BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,774,846	A→G	13.0%	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	13.0%	316.4 / 9.7	172	S780S (AGT→AGC)	clfB	Clumping factor B
Reads supporting (aligned to +/- strand): ref base A (70/70); new base G (23/3); total (94/78)
Fisher's exact test for biased strand distribution p-value = 1.98e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

TTGAATCGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAAT  >  USA300TCH1516_ALE/2774707‑2774982

ttGAATCGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTAGCTATCCGAGTCTGAGTCGCTGTCT‑GGATCTGAATCAc                                                                                                                                         >  1:159081/1‑141 (MQ=255)

 tGAATCGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAGTCTGAATCACt                                                                                                                                        >  2:450507/1‑141 (MQ=255)

  gAATCGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTg                                                                                                                                       >  2:212634/1‑141 (MQ=255)

   aagcGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGt                                                                                                                                      <  1:351893/138‑1 (MQ=255)

     tCGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTCTCCAAGCCCGAGTGGCTGTCT‑GAATCTGAATCAcgtgcc                                                                                                                                    >  2:464789/1‑136 (MQ=18)

      gggAGGCACTATCGGGGTCTGACTCGCTATCTGAATCCGAGTCGCTATCCGGGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTg                                                                                                                                   <  2:560116/140‑1 (MQ=17)

      cGGAGTCACTATCGGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATTCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTg                                                                                                                                   <  2:1016898/141‑1 (MQ=255)

         aGTCACTGTCTGAGTGCGAGTCGCTAGCTGCAGCCGAGTCGCTATGCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAAt                                                                                                                                <  2:756656/141‑1 (MQ=18)

         aGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCC‑GAATCTGAATCACTGTCTGCAt                                                                                                                                >  1:175848/1‑141 (MQ=255)

          gTCATTTTCTGAGTCGGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATc                                                                                                                               <  2:84194/141‑1 (MQ=255)

            caatATCTGGGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGGCTGAGTCGCTATCTGAGGCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATcct                                                                                                                             <  2:207415/138‑2 (MQ=255)

             aCTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGAATCCGa                                                                                                                            >  1:480561/1‑141 (MQ=255)

             aCTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGa                                                                                                                            >  1:397488/1‑141 (MQ=255)

             aCTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGa                                                                                                                            >  1:253270/1‑141 (MQ=255)

              cTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGGATCCGaa                                                                                                                           >  1:743326/1‑141 (MQ=255)

                        ctgggtcTCTATTTGAGTCCGGGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTatctg                                                                                                                 <  2:187666/141‑1 (MQ=18)

                          gagtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaa                                                                                                               >  1:217195/1‑141 (MQ=255)

                           agtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAGTCTGAATCGCTGTCTGAATCCGTATCGCTgtctgaat                                                                                                              >  2:48287/1‑141 (MQ=255)

                            gtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAGTCGCTatctgaaac                                                                                                             >  2:414069/1‑139 (MQ=255)

                            gtcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAAGCGCTATCTGAGTCTGAGTCGCTGTCT‑GGATCTGAATCACAGTCAGGATCTGAGTGGCAatctggatg                                                                                                             >  2:939069/1‑140 (MQ=14)

                             tcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GACTCTGAGTCGCTGTCTGAGTCGGAATCGCTatctggacct                                                                                                            >  2:614497/1‑136 (MQ=14)

                             tcGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GACTCTGAATCACTGTCTGAATCCGAATCGCTatctgcatct                                                                                                            >  1:288077/1‑141 (MQ=255)

                              cGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTatctgaatctg                                                                                                           <  2:678615/141‑1 (MQ=255)

                                cTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTGTGAATCGCTGTCTGTGTCTGAGTCGCTATCTGAGTCTGAGTCGCTTTCT‑GAGTCTGATTCGCTATCTGAATCCGAATCCCTAtctgaatcggaa                                                                                                         >  2:764035/1‑141 (MQ=14)

                                cTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAGTCGGAATCGCGatctgagtcagca                                                                                                         >  2:1197313/1‑137 (MQ=17)

                                cTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTAtctgaatctgaa                                                                                                         >  1:841874/1‑141 (MQ=255)

                                cTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAATCTGAGTCGTTGTCT‑GAATCTGAATCGCTGTCGGAATCAGAATCGCTGtcagaatctaaa                                                                                                         >  2:817252/1‑141 (MQ=17)

                                 tATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATctgaatctgaat                                                                                                        <  1:949960/141‑1 (MQ=255)

                                  aTCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCtgaatctgaatc                                                                                                       <  1:490465/141‑1 (MQ=255)

                                   tCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCT‑GAGTCTGAGTGGCTGTCTGAATCCGAATCGCGATACGGATCTGAATCg                                                                                                      >  2:405633/1‑141 (MQ=14)

                                    cTGATTCCGAGTCGCTTTCCGGGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGc                                                                                                     <  2:71664/141‑1 (MQ=255)

                                    cTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTTAGACTTTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGc                                                                                                     <  1:1178012/141‑1 (MQ=255)

                                    cTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGATTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAAGCTGAATCACTGTCTTAATCCGAATCCCTCTTTGAGTCTGAATTGc                                                                                                     >  1:229574/1‑141 (MQ=17)

                                     tGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCt                                                                                                    >  1:972236/1‑141 (MQ=255)

                                     tGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCt                                                                                                    >  1:1053018/1‑141 (MQ=255)

                                      gAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGCATCGCTATCTGAATCTGAATCGCtg                                                                                                   >  1:154488/1‑140 (MQ=255)

                                        aTCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATc                                                                                                 >  1:424696/1‑141 (MQ=255)

                                         tCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGAATCTGAATCGCTGTCTGAATCTGAATCGCTGTcc                                                                                                >  2:335053/1‑141 (MQ=255)

                                          ccGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCg                                                                                               <  2:1050292/141‑1 (MQ=255)

                                            gAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAg                                                                                             >  1:26552/1‑141 (MQ=255)

                                               tCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCt                                                                                          >  2:396415/1‑141 (MQ=255)

                                                     tCCGAGTCTGGGTCGCTATCTGAGGCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCg                                                                                    <  2:907994/141‑1 (MQ=255)

                                                        gAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAGTCGCTGTCTGAATCCGAATCGCGATCTGAATCTGAATCGGTGTCCGAGTCTGAGCAGCTg                                                                                 >  2:90304/1‑141 (MQ=14)

                                                        gAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTg                                                                                 <  2:1119293/141‑1 (MQ=255)

                                                          gTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATGCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtc                                                                               >  1:792975/1‑141 (MQ=255)

                                                          gTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAATCTGAGTCGCTGTCT‑GAATCTGCATCGCTGGATGAGTCGGCATCGCTGGCGGCATCTGAATCGCTCGCCGAGCCAGAGTCGCTGtc                                                                               >  2:63801/1‑141 (MQ=14)

                                                          gTCTGAGTCGCTATCTGAGTCTGAATCGCTGGCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGCATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtc                                                                               <  1:232347/141‑1 (MQ=255)

                                                            gtGATTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctg                                                                             <  2:766118/140‑1 (MQ=255)

                                                            cTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGGATCAGAATCGCGATCTGGATCTGAGACGCTATCAGAAGCTGAGACGCGGtctg                                                                             >  2:971676/1‑141 (MQ=14)

                                                            cTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctg                                                                             <  1:208619/141‑1 (MQ=255)

                                                              gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGAGTCGCTGtctgaa                                                                           >  1:767953/1‑141 (MQ=255)

                                                              gAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGGGTCGCTGtctgat                                                                           >  1:338903/1‑140 (MQ=255)

                                                               aGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGAATCCGAATCGCTGTCTGAATCTGAATCGCTATCTGAGTCTCCGTCGCTGtctgtgt                                                                          >  1:210901/1‑138 (MQ=17)

                                                                 tCGCTTTCTGGGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAGTCACTGTCTGAGTCCGAGTCGCTATCTGAATCTGATTCGCTATCCGAGTCTGAGTCGCTGtctgaatct                                                                        <  2:1092974/141‑1 (MQ=14)

                                                                  cGCTATCTGAGTCTGAATCGCTGTCTGAGGCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCGGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGtctgaatctg                                                                       <  1:711762/141‑1 (MQ=255)

                                                                      aTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAGTCTGAATCACTGTCTGAATCCGAATCGCTGTCTGAGTCTGAGTCGCTATCGGAGTCTGAGGCGCTgtctgagtctgactg                                                                   >  2:610134/1‑138 (MQ=14)

                                                                      aTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGACTCACTGTCTGAATCCGAATCGCTGTCTGAGTCGGAGTCGCTGTCTGAGTCCGCGGCACTGTCtgagtctgaatc                                                                   >  2:279876/1‑141 (MQ=14)

                                                                      aTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCtgaatctgaatc                                                                   >  2:1193141/1‑141 (MQ=255)

                                                                      aTCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCtgaatctgaatc                                                                   <  1:805778/141‑1 (MQ=255)

                                                                          gAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTg                                                               >  2:526703/1‑141 (MQ=255)

                                                                              cgtattCGCTGTCTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTGTCT‑GAATCTGAGTCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTg                                                           <  2:633966/136‑1 (MQ=11)

                                                                               tGAATCCCTGTCTGAGTCGGAGTCGCTATCCGAGTCTGAGTCGCCGTATAGGAT‑GCAGCCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGa                                                          <  1:929039/141‑1 (MQ=17)

                                                                                 aaTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGt                                                        <  1:87260/141‑1 (MQ=255)

                                                                                  aTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCGCTGTCTGAGTCCGAATCGCTGTCTGAATCTGTATCGCTTTGTGAGGCTGTGTGGCTGGCTGAATCTGAGTCGCTTGCTGATTc                                                       >  2:286085/1‑141 (MQ=11)

                                                                                  aTCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTc                                                       <  1:196269/141‑1 (MQ=255)

                                                                                          ctgagtctgagtCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGGGTCGCTGTCTGAACCTGAATCGCTGTCTGAGGCCGGATCGc                                               >  2:629649/1‑141 (MQ=255)

                                                                                          ctgagtctgagtCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTGTGAGTCCGAATCGc                                               >  2:12192/1‑141 (MQ=255)

                                                                                          ctgagtcggcggTGCTATCCGAGTTTGAGTCTCTTTCT‑GAATCTGAATCTCTCTCTGACTCCGAATCGTTATCTGAATCTGAATCGCTATCCGAGTCGGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGc                                               <  2:717750/141‑1 (MQ=17)

                                                                                               tcggagtcGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTTTCTGAATCCGAATCGCTATCTGAATCTGACTCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatct                                          <  2:387459/138‑1 (MQ=255)

                                                                                                     tcGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCGGAGTCTGAGTCGCTATCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTttctgaatct                                    >  1:278326/1‑141 (MQ=255)

                                                                                                     tcGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatctgagtct                                    >  1:68754/1‑141 (MQ=255)

                                                                                                       gCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatctgaatctga                                  <  1:560056/141‑1 (MQ=255)

                                                                                                       gCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATAACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCTGAGTCTGATTCGCTGTCTGATTCTGAATCGCTGTCTGAGTCTGAATCGCTatctgaatctga                                  >  2:573441/1‑141 (MQ=14)

                                                                                                       ccTATCCGTGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTatctgaatctga                                  <  1:135494/140‑1 (MQ=255)

                                                                                                        cTATCCGAGTCTGAGTCGCTGTCT‑GAATGGGAGTCAGTGTCTGAGGCAGAATGGCTATCTGATGCTGAATAGCTAGGTGGGTGTGGGACGATTTCTGAAACAGAATCGATGGCTGAGACCGAATAGCTATCTGAATCTGAg                                 >  2:780203/1‑141 (MQ=12)

                                                                                                        cTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCt                                                                                                                                    >  2:804945/1‑42 (MQ=40)

                                                                                                        cTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAg                                 >  1:1094855/1‑141 (MQ=255)

                                                                                                          aTCCGAGTCTGAGTCGCAGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTCTCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGATTCGCTGTCTGAGTCCGAATGGCTATCTGGATCTGAGtt                               >  1:110323/1‑140 (MQ=17)

                                                                                                            ccGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGc                             <  2:1061573/141‑1 (MQ=255)

                                                                                                              gAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTg                           <  1:450507/141‑1 (MQ=255)

                                                                                                                 tCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtct                        <  1:1007254/141‑1 (MQ=255)

                                                                                                                 tCTGAGTCGCTGTCC‑GAATCTGACTCCCTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAACCGCTATCTGAATCTGAGTCGCTgtct                        <  2:743326/141‑1 (MQ=255)

                                                                                                                  cTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAAGCTGATTCGCTgactg                       >  1:1860/1‑141 (MQ=255)

                                                                                                                   tGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCGGAGTCGCTgtctga                      >  2:157563/1‑141 (MQ=255)

                                                                                                                   tGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAAAAGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctga                      <  1:630688/141‑1 (MQ=255)

                                                                                                                   tGAGTCGCTGTCT‑GAATCTGAATCACTATCTGAATCCGAATCGCTATCTGAATCTGAGTCGCTGTCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGGATCGCGAGCTGAAGCTGAGGCGCGttctgg                      >  2:50246/1‑140 (MQ=14)

                                                                                                                     ggTCGCTGTCT‑GAATCTGAATCACTCTCTGAATCCGAATCGCTATCTGAATCTGACTCGCTATCCGAGTCTGAGTCGCTTTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagt                    <  2:876849/140‑1 (MQ=255)

                                                                                                                      gTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAGTCTGAATCGCTATCTGAGTCTGAGTCGCTATCTGAATCTGAATCGCTAACTGAGTCTGAGTCGCTATCTGAGTCTGAGTCGCTagctgagtc                   >  2:529817/1‑141 (MQ=11)

                                                                                                                      gTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgggtc                   >  1:756506/1‑141 (MQ=255)

                                                                                                                      gTCGCTGTCG‑GCATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagtc                   <  1:1034827/141‑1 (MQ=255)

                                                                                                                       tCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAGTCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCTGAATCGCTGTCTGTATGTGAGTCGCTgtctgcgtcg                  >  2:104454/1‑140 (MQ=17)

                                                                                                                        cGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagtctg                 <  2:22306/141‑1 (MQ=255)

                                                                                                                         gCTTTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTgtctgagtctga                <  1:170681/141‑1 (MQ=255)

                                                                                                                              ct‑gaatctgaatCACTGTCTGAATCCGAATCGCTATCTGAAACTGAACTGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGc           <  1:620820/141‑1 (MQ=255)

                                                                                                                               t‑gaatctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGCATCGCTGTCTGAGTCCGAATCGCCATCTGCATCTGGGTCGCTGTCTGAGTCTGGACCGCt          >  1:911538/1‑141 (MQ=17)

                                                                                                                               ccgaatctgactcACTGTCTGAACCCGAATCGCTATCTGATTCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGc           <  2:1054520/139‑1 (MQ=255)

                                                                                                                                 gaatctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTa         <  1:396415/141‑1 (MQ=255)

                                                                                                                                   atctgaatcACTGTCGGAATCCGAATCGCTATCTGAACCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatc       <  1:414069/141‑1 (MQ=255)

                                                                                                                                     ctgaatcACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctg     >  1:202800/1‑141 (MQ=255)

                                                                                                                                        agtcACTGTCTGAGGCCGAGTCGCTATCTGAATCTGAATCGCTCTCTGAGTAAAAATCGCTCTCTGAATCTGCGTCGCTTTCTGAGTCCGAATCGCTAGCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTatctgaat  <  2:1149866/139‑1 (MQ=11)

TTGAATCGGAGTCACTATCTGAGTCTGAGTCGCTATCTGAATCCGAGTCGCTATCCGAGTCTGAGTCGCTATCTGAGTCTGAATCGCTGTCTGAGTCTGAGTCGCTATCCGAGTCTGAGTCGCTGTCT‑GAATCTGAATCACTGTCTGAATCCGAATCGCTATCTGAATCTGAATCGCTATCCGAGTCTGAGTCGCTGTCTGAATCTGAATCGCTGTCTGAGTCCGAATCGCTATCTGAATCTGAGTCGCTGTCTGAGTCTGAATCGCTATCTGAAT  >  USA300TCH1516_ALE/2774707‑2774982

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