BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,375,937	G→A	5.5%	Y29Y (TAC→TAT)	topB ←	DNA topoisomerase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,375,937	0	G	A	5.5%	784.1 / 20.2	234	Y29Y (TAC→TAT)	topB	DNA topoisomerase 3
Reads supporting (aligned to +/- strand): ref base G (123/98); new base A (2/11); total (125/109)
Fisher's exact test for biased strand distribution p-value = 7.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 240 total aligned reads displayed.

TTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCATAAAT  >  USA300TCH1516_ALE/2375847‑2376027

ttCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAg                                                                                            >  2:3344855/1‑91 (MQ=255)

  ccTTTAAGTTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                          <  1:1936560/91‑1 (MQ=255)

  ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCCCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                          >  1:3364107/1‑91 (MQ=255)

  ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                          >  1:3034842/1‑91 (MQ=255)

  ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                          >  2:1937242/1‑91 (MQ=255)

  ccTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTa                                                                                          >  2:3547424/1‑91 (MQ=255)

   cTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTaa                                                                                         <  2:2142631/91‑1 (MQ=255)

   cTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATaa                                                                                         <  1:3171673/91‑1 (MQ=255)

    tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAAc                                                                                        <  1:1702428/91‑1 (MQ=255)

    tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAAc                                                                                        >  2:3742610/1‑91 (MQ=255)

    tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAAc                                                                                        <  2:3437080/91‑1 (MQ=255)

    tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAAc                                                                                        <  1:1742584/91‑1 (MQ=255)

    tttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAAc                                                                                        <  2:1847877/91‑1 (MQ=255)

     ttAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAAcc                                                                                       <  2:2814731/91‑1 (MQ=255)

      tAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAACCa                                                                                      <  2:1030466/91‑1 (MQ=255)

      tAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTACTTTGAAAATAACCa                                                                                      <  1:1260377/91‑1 (MQ=255)

        aaTTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCAtt                                                                                    >  1:680332/1‑91 (MQ=255)

          tttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTAc                                                                                  <  2:2099635/91‑1 (MQ=255)

          tttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTAc                                                                                  <  1:112541/91‑1 (MQ=255)

           ttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACg                                                                                 >  2:1635351/1‑91 (MQ=255)

           ttttATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACg                                                                                 >  2:2912579/1‑91 (MQ=255)

             ttATCGTATTGTTCAGGTGTCACATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCt                                                                               <  2:2777448/91‑1 (MQ=255)

              tATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCtt                                                                              >  2:2485818/1‑91 (MQ=255)

               aTCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTc                                                                             <  1:84183/91‑1 (MQ=255)

                tCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCt                                                                            <  1:2567679/91‑1 (MQ=255)

                tCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCt                                                                            <  2:3746524/91‑1 (MQ=255)

                 cGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTg                                                                           >  1:3653295/1‑91 (MQ=255)

                  gTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGa                                                                          >  1:425452/1‑91 (MQ=255)

                  gTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGa                                                                          >  2:1562375/1‑91 (MQ=255)

                    aTTGTTCTGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGatt                                                                        <  1:1942664/91‑1 (MQ=255)

                     ttGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattt                                                                       <  2:2267295/91‑1 (MQ=255)

                     ttGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattt                                                                       <  1:3788629/91‑1 (MQ=255)

                      tGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGattta                                                                      <  2:1195562/91‑1 (MQ=255)

                           aGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGt                                                                 >  1:2681263/1‑91 (MQ=255)

                           aGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGt                                                                 >  2:2699832/1‑91 (MQ=255)

                           aGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTGATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGt                                                                 <  1:3269309/91‑1 (MQ=255)

                            ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                                >  1:2402919/1‑91 (MQ=255)

                            ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                                <  2:1700276/91‑1 (MQ=255)

                            ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                                >  1:2076223/1‑91 (MQ=255)

                            ggTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTa                                                                <  1:2122919/91‑1 (MQ=255)

                                   cATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATc                                                         >  2:3106704/1‑91 (MQ=255)

                                    aTTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCa                                                        >  2:1526335/1‑91 (MQ=255)

                                    aTTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCa                                                        >  1:904451/1‑91 (MQ=255)

                                      ttGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGc                                                      >  2:959213/1‑91 (MQ=255)

                                        gTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCaa                                                    <  1:3250882/91‑1 (MQ=255)

                                         tCACTAGATGACCTAACGCCCACGTGACAGTGTATTGGTTATTTTCAAAGTAAACATTACGCTTCTGATTTATTTGTAAAGCATCAGCAat                                                   <  1:3103969/91‑1 (MQ=255)

                                         tCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAat                                                   >  1:1714192/1‑91 (MQ=255)

                                         tCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAat                                                   <  1:1552067/91‑1 (MQ=255)

                                           aCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAatat                                                 >  2:2719195/1‑91 (MQ=255)

                                           aCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAatat                                                 >  2:2412158/1‑91 (MQ=255)

                                           aCTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAatat                                                 >  2:2087208/1‑91 (MQ=255)

                                            cTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATAtc                                                <  1:1841570/91‑1 (MQ=255)

                                            cTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATAtc                                                >  1:296675/1‑91 (MQ=255)

                                               gATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATctct                                             <  1:2341733/91‑1 (MQ=255)

                                                aTGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCtt                                            <  1:971949/91‑1 (MQ=255)

                                                aTGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCtt                                            <  2:3671903/91‑1 (MQ=255)

                                                 tGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTg                                           >  2:1931982/1‑91 (MQ=255)

                                                     cTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGAc                                       >  2:2041453/1‑91 (MQ=255)

                                                       aaCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTg                                     <  2:715755/91‑1 (MQ=255)

                                                        aCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGa                                    >  2:1835989/1‑91 (MQ=255)

                                                        aCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGa                                    >  2:3682195/1‑91 (MQ=255)

                                                        aCGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGa                                    >  2:1359868/1‑91 (MQ=255)

                                                         cGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGAt                                   >  2:3524545/1‑91 (MQ=255)

                                                          gCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATg                                  >  2:3239336/1‑91 (MQ=255)

                                                          gCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATg                                  <  1:2764028/91‑1 (MQ=255)

                                                          gCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATg                                  >  1:2506761/1‑91 (MQ=255)

                                                           cccACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATgg                                 <  2:1823143/91‑1 (MQ=255)

                                                           cccACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATgg                                 <  1:1526335/91‑1 (MQ=255)

                                                             cACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtt                               >  2:2302362/1‑91 (MQ=255)

                                                             cACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtt                               <  1:3452672/91‑1 (MQ=255)

                                                             cACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtt                               >  1:3213809/1‑91 (MQ=255)

                                                             cACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtt                               <  1:3201388/91‑1 (MQ=255)

                                                               cGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtttt                             >  2:2627894/1‑91 (MQ=255)

                                                               cGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGtttt                             <  1:3630783/91‑1 (MQ=255)

                                                                 tGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTc                           >  2:472510/1‑91 (MQ=255)

                                                                 tGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTc                           >  1:1001967/1‑91 (MQ=255)

                                                                  gACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCa                          >  1:197993/1‑91 (MQ=255)

                                                                   aCAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAg                         >  2:229115/1‑91 (MQ=255)

                                                                    cAATATAATGGTTATTGTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGc                        <  2:3705411/91‑1 (MQ=255)

                                                                     aaTATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCt                       <  2:3112609/91‑1 (MQ=255)

                                                                     aaTATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCt                       <  1:1988492/91‑1 (MQ=255)

                                                                       tatatTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTaa                     >  2:2932728/1‑91 (MQ=255)

                                                                       tatatTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTaa                     >  2:2931880/1‑91 (MQ=255)

                                                                         tatTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAata                   >  2:2146173/1‑91 (MQ=255)

                                                                          atTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAatat                  <  1:1628708/91‑1 (MQ=255)

                                                                          atTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAatat                  >  2:301060/1‑91 (MQ=255)

                                                                          atTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAatac                  >  2:304192/1‑90 (MQ=255)

                                                                           ttGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTCATAtt                 >  2:3144213/1‑91 (MQ=255)

                                                                            tGGTTACTTTGAAAATAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTa                >  1:2285982/1‑91 (MQ=255)

                                                                               ttATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAg             >  1:1990049/1‑91 (MQ=255)

                                                                                  ttttCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGAtt          >  2:3422549/1‑91 (MQ=255)

                                                                                    ttCAAAGTAACCATTACGCTTCTGTTTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTc        <  2:3616199/91‑1 (MQ=255)

                                                                                    ttCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTc        <  1:3433863/91‑1 (MQ=255)

                                                                                    ttCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTc        <  1:458608/91‑1 (MQ=255)

                                                                                      cAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCAt      >  1:2086076/1‑91 (MQ=255)

                                                                                      cAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCAt      >  1:3512551/1‑91 (MQ=255)

                                                                                       aaaGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCATa     <  1:57296/91‑1 (MQ=255)

                                                                                       aaaGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACCGATGGTTTTTCAGCTAATATTAAAGATTTCATa     >  1:3174953/1‑91 (MQ=255)

                                                                                          gTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCATAAAt  <  2:533863/91‑1 (MQ=255)

                                                                                          gTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCATAAAt  <  2:586392/91‑1 (MQ=255)

TTCCTTTAAATTTTTATCGTATTGTTCAGGTGTCGCATTTGTCACTAGATGACCTAACGCCCACGTGACAATATATTGGTTATTTTCAAAGTAACCATTACGCTTCTGATTTATTTGTAAAGCATCAGCAATATCTCTTGCGACTGATGGTTTTTCAGCTAATATTAAAGATTTCATAAAT  >  USA300TCH1516_ALE/2375847‑2376027

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