BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,313,831	C→A	7.7%	A17E (GCA→GAA)	USA300HOU_RS06490 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,313,831	0	C	A	7.7%	350.1 / 2.2	140	A17E (GCA→GAA)	USA300HOU_RS06490	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (58/64); new base A (2/11); total (61/79)
Fisher's exact test for biased strand distribution p-value = 3.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 145 total aligned reads displayed.

AAAAATAATTCAAAGAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGAAATGGTGTAGAT  >  USA300TCH1516_ALE/1313698‑1313964

aaaaaTAATTCAAAGAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAAATACTCTTTATGAGAGGATGAACAAACCCCACTACCATTCTGTATAGGGGATATCCTACGTAAACTTGCACGACACGc                                                                                                                                       >  2:318353/1‑134 (MQ=255)

aaaaaTAATTCAAAGAATTTAATATTTAATTTAAATAATAGAAAGTATAAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGc                                                                                                                                       <  1:2042068/134‑1 (MQ=255)

aaaaaTAACACCAAAAATATAATATTAATTTAAAATAAAAGAACATATAAAATAAACTAGTATTTATGAGAAGATGAACAAAATATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGc                                                                                                                                       <  1:1792271/134‑1 (MQ=255)

  aaaTAATACAAAAAATTTAATATTTACTTTAAATATTAGAAAATATTAATAAAACTAGAATATAAAAAAGAATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCaa                                                                                                                                     <  1:214898/134‑1 (MQ=255)

    ataatTCAAAAAAATTAATATTTACTTTAAATAATAAAACAAATGAATTAAACTAGTATTTATTAAAGGATAAAAAAAAAATGAAAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAAtt                                                                                                                                   <  1:301484/134‑1 (MQ=255)

      aatTAAAAGAATTTAATATTTATTTTAAATATAAAAACGTATAAATTAAACTGGAATTTATGAGAGGATGAACAAAAAATGAGAAAAATGTTTATAGGGAATATCGAAGTTAAAATTGACAGAGCCGCAAttca                                                                                                                                 >  2:1928814/1‑132 (MQ=255)

       atTCAAAGAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGaa                                                                                                                                >  1:2599265/1‑134 (MQ=255)

         tCAAAGAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAAc                                                                                                                              >  1:2534397/1‑134 (MQ=255)

              gAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAACTGGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACa                                                                                                                         <  2:791489/134‑1 (MQ=255)

              gAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACa                                                                                                                         <  2:723920/134‑1 (MQ=255)

                  ttAATATTTACTTTAAATATTAGAACGTATGAATTAAACTCGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGAGATATCGTAGATAAAATTGGACGAAACGAAATTGAAACGTACATAcc                                                                                                                     >  1:740760/1‑134 (MQ=255)

                    aaTATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTc                                                                                                                   >  2:561598/1‑134 (MQ=255)

                    aaTATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTAGACGAGACGCAATTGAAACGTACATACCTc                                                                                                                   >  1:616381/1‑134 (MQ=255)

                        tttaatttaaATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACt                                                                                                               <  2:1565080/134‑1 (MQ=255)

                          tACTTTAAATATTAGGACGTATGAATTCAACTAGTATTTATGAGAGGCTGAAAAAACAATGAGAATAATATTTCTCGGGGATATCGTAGGTAAAATTGGACGAAATGAACTAGTGACGTACATACCTCAACTGa                                                                                                             <  2:2187410/134‑1 (MQ=255)

                               tAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCaa                                                                                                        <  2:2747413/134‑1 (MQ=255)

                                  atatTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGTATTTAATACGTACATACCTCAACTGAAGCAAAAg                                                                                                     >  2:1286714/1‑134 (MQ=255)

                                     ttAAAACGTCTGAAATAATTTAGTATTAATGAGAAGATGAACAAAACATGAGAAAAATGGTTATAGGGGATATCGTAGGTATAGTTGGACGAGACGCAATTGAAACGTACATACCTCATCTGGATCAAATGttt                                                                                                  <  2:507254/134‑3 (MQ=255)

                                      tAGAATTTATAAATTAAACTGATATTTATGAAAAGAAGAAAAAAAAATGAAAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGtata                                                                                                 <  1:286918/134‑1 (MQ=255)

                                            gTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCaa                                                                                           <  2:2224377/134‑1 (MQ=255)

                                            gTATGAATTAAACTAGTAGTTAGGAGAGGATGAACAGAACATGATAATAATGTTTATAGGGGATATCGTACGTAAAAACGAACTATACGCAATTGAGACGTACATACCTCAACTGAAGCAAAAGTATAAACCaa                                                                                           <  2:2651118/134‑1 (MQ=255)

                                                gAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGt                                                                                       <  2:1991839/134‑1 (MQ=255)

                                                gAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGt                                                                                       >  1:746794/1‑134 (MQ=255)

                                                   ttaaattaaaATTTATAAAGGTATTAACAAAACATAAGAAAAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTAc                                                                                    <  1:1697268/134‑1 (MQ=255)

                                                   ttAAACTAATATAAATGAAAGAATAAACAAAACATAAGAATAATGTTTATAGGGGATATCGTAGGTAAAAATGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTAc                                                                                    <  1:2337684/134‑1 (MQ=255)

                                                    tAAACTAGTATTTATGAGATGATTAACAAAACATAAGAAAAATGATTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACa                                                                                   <  1:837815/134‑1 (MQ=255)

                                                    tAAACTAGGATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATCGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACa                                                                                   <  2:2003925/134‑1 (MQ=255)

                                                     aaaCTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACaa                                                                                  >  1:2817904/1‑134 (MQ=255)

                                                     aaaCTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACaa                                                                                  >  1:251454/1‑134 (MQ=255)

                                                     aaaCTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACaa                                                                                  >  1:1207371/1‑134 (MQ=255)

                                                            tATTTATGAGAGTATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTaaa                                                                           <  2:2465199/134‑1 (MQ=255)

                                                             aTTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAAt                                                                          >  2:2715167/1‑134 (MQ=255)

                                                              tttATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAGTTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATg                                                                         >  2:2649048/1‑134 (MQ=255)

                                                                tATGAGAGTATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAATAGTTACAATTGTAAATGCt                                                                       <  2:1764819/134‑1 (MQ=255)

                                                                tATGAGAGGATGAACAAAACATGAGAATAATGTGTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCt                                                                       <  2:534956/134‑1 (MQ=255)

                                                                 attagagGATAAAAAAAACTTGAGTATATTGATAATAAGGAATATCGTAAGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTg                                                                      <  1:1952763/131‑1 (MQ=255)

                                                                    agagGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGaaa                                                                   <  2:1649507/134‑1 (MQ=255)

                                                                    agagGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGaaa                                                                   <  1:2757746/134‑1 (MQ=255)

                                                                    agagGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGAACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATCCTGcaa                                                                   >  2:2234438/1‑134 (MQ=255)

                                                                       gagtGAAAAAAACATGAGAATAATATATATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATg                                                                <  1:1944400/131‑1 (MQ=255)

                                                                        gATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATgc                                                               <  2:526400/134‑1 (MQ=255)

                                                                         ataaaaaaaaaaTGAAAATAATGTTTATAGGGGATATCGTAGGAAAAATTGGACGAGACGCAATTGGAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGTTGAAAATgca                                                              <  2:2076539/128‑1 (MQ=255)

                                                                         aTGAATAAAATATGAGAATAATGTTTATAGTGGATATATTAGGTAAAATTGGACGAGATGTAATTGAAAAGTACATACCTTAACTGAAGCAAAAGTATAATCAAACAGTTACTATTGTAAATGCTGAAAATgca                                                              <  2:2666393/134‑1 (MQ=255)

                                                                           gAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATgcagc                                                            >  1:2195012/1‑134 (MQ=255)

                                                                           gAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATgcagc                                                            >  2:2295055/1‑134 (MQ=255)

                                                                               aaaaCATCAGAATAATGTTTATAGTGGATATAGTAAGTAAAATTGGACGAGACACAATTGAAACGTACATAACTCAACGGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAAGGCAGCACAt                                                        <  1:204502/134‑1 (MQ=255)

                                                                                taaCATTAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATg                                                       <  2:2589587/133‑1 (MQ=255)

                                                                                aaaaaagagaATAAGATTTATAGGGTATATCGTAGGTAAAATTGGAAGAGACGGGATTGAAACGTACAGACCTCAACTGAAGCAAAAGTATAATCCATCATTTACAATTGTAAATGCTGAAAATGCAGCACATg                                                       <  2:2639283/128‑1 (MQ=255)

                                                                                aaaCATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATg                                                       >  2:1206728/1‑134 (MQ=255)

                                                                                 aaCATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAAAGTACATACATCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATgg                                                      >  1:879093/1‑134 (MQ=255)

                                                                                 aaCATGAAAAAAATGTTAATAGGGGATATCGTAGGTAAAATTAGACGAGACGCAATTGAAACGTACATACATCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATgg                                                      <  1:597559/134‑1 (MQ=255)

                                                                                   cATGAGAAAAATATTTATAGGGGATATCGTAGATAAAATTAAACGAGACGAAATTGAAACGTACATAACTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTa                                                    <  1:2646331/134‑1 (MQ=255)

                                                                                     tGAGAATAATATTTATAGGGGATATCGTAGATAAAATTGGAAAAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTaaa                                                  <  1:2207051/134‑1 (MQ=255)

                                                                                      gagaATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAg                                                 >  1:2389940/1‑134 (MQ=255)

                                                                                       agaATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACAAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAgg                                                >  2:478431/1‑134 (MQ=255)

                                                                                            aatGTTTATAGGGGATATTGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAGAAAAATCAAAAGAATAAATTAAAAATGACGATTGTAAATGCTGAGAATGTTTCACTGGGTAAAGGTTTGa                                           <  2:616381/134‑1 (MQ=255)

                                                                                             atGTTTATAGGGGATATCGTAGGTAAAATTAGACGAGACAAAATTGAAACGTACATACATCAACTGAAGAAAAAGTATAAACCAATAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGAc                                          <  1:2135135/134‑1 (MQ=255)

                                                                                             atGTTTATAGGGGATATAGTAGGTAAAATTGGAAGAGAAGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGAc                                          <  1:2588646/134‑1 (MQ=255)

                                                                                             atGTTTATAGGGGAGATCGTAAGTAAAATTGGACGAAACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGAc                                          <  2:503663/134‑1 (MQ=255)

                                                                                             atGTTTATAGAAAATATAGTAGGTAAAATTGAACTAGACAAAATTGAAAAGTACATATCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGAc                                          <  1:1555085/134‑1 (MQ=255)

                                                                                              tGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACt                                         <  2:1573390/134‑1 (MQ=255)

                                                                                                 ttATAGGTGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAAGCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaa                                      <  2:685364/134‑1 (MQ=255)

                                                                                                    taGGGGATGTCGTAGGTAAAATTGGACGGGATGCGATTGAAACGTACATACCTCAAGTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAGATGCTGAAAGTGCAGCACATGGTAAATGTTTGACTGaaaaa                                   <  2:106256/134‑1 (MQ=255)

                                                                                                    taGGGGATATCGTAGGTAAAATTTGACAAGACGAAATTAAAAAGTACATACCTCAACTGAAGCAAAAGTAAAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaaaaa                                   <  1:1536155/134‑1 (MQ=255)

                                                                                                     aGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAATCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaaaaaa                                  <  2:570225/134‑1 (MQ=255)

                                                                                                     aGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaatata                                  >  1:769852/1‑130 (MQ=255)

                                                                                                     aGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaaaaaa                                  <  2:2101083/134‑1 (MQ=255)

                                                                                                     aGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaaaaaa                                  <  2:548037/134‑1 (MQ=255)

                                                                                                     aGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaaaaaa                                  <  2:1984955/134‑1 (MQ=255)

                                                                                                     aGGGGATATCGTAGGTAAAAATGGACGAGACGAAATTGAAACGTACATACATCAACTAAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACAAGGACCAGGTTCGACTGaacaac                                  >  2:2628529/1‑133 (MQ=255)

                                                                                                     aGGGGATAACGTATGTAAAATTAGACGAAAAGAAAATGAAACGTAAATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGaaaaaa                                  <  1:453066/134‑1 (MQ=255)

                                                                                                         gATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCAAAAGGTAAAGGTTTAACTGAAAAAataga                              >  1:2645292/1‑132 (MQ=255)

                                                                                                           tatCGTAAGTAAAAATGGAAGGGAAGACATTGACAAGAACAAACAAAAACTGATGTAAATCTATATTCTCTCACTTATAATTGTAAATGCTGAAAACGCAGCTATTACAAAAGGATTGACTGAAAAAAtatata                            <  2:2185602/134‑1 (MQ=255)

                                                                                                            atCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGCAAAGGTTTGACTGAAAAAATATATaa                           >  2:1513336/1‑134 (MQ=255)

                                                                                                             tagTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATaaa                          >  2:2824967/3‑134 (MQ=255)

                                                                                                             tCGTATGTAAAATTGGACGAGAAGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAGCAAAACAGTTAAAATTATTGATGGAGAAAGTGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATaaa                          >  2:615002/1‑134 (MQ=255)

                                                                                                             tCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATaaa                          <  2:746794/134‑1 (MQ=255)

                                                                                                              cGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAAc                         <  2:251454/134‑1 (MQ=255)

                                                                                                              cGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAAc                         >  2:375771/1‑134 (MQ=255)

                                                                                                              cGTAGGTAAAATTAGAAGAGACGCAATTGAAACGTACATAAAACAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAAc                         <  1:422817/134‑1 (MQ=255)

                                                                                                               gTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAAACAAAAGTTACAATTGTAAATGCTGAAAATGCCGCACATGTAAAAAGGTTGACTGAAAAAATATATAAACa                        >  1:1478616/1‑134 (MQ=255)

                                                                                                               gTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAAAAAACAGAACCAATTGTAAATGCTGAAAATGCAGCAAATGGTAAAGGTTTGACTGAAAAAATATATAAACa                        >  2:1027741/1‑134 (MQ=255)

                                                                                                                tAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACaa                       >  2:2742209/1‑134 (MQ=255)

                                                                                                                 aGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAACACATGGTAAAGGTTTGACAGAAAAAATATATAAACAAt                      >  1:1095920/1‑134 (MQ=255)

                                                                                                                  ggTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAAtt                     >  1:2838720/1‑134 (MQ=255)

                                                                                                                  ggTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATACAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAAtt                     >  1:1148642/1‑134 (MQ=255)

                                                                                                                   gTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAGAGATTTGACTGAAACAACATATAAACACTTa                    >  2:1674670/1‑134 (MQ=255)

                                                                                                                   gTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTa                    >  1:1417904/1‑134 (MQ=255)

                                                                                                                     aaaaTTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACt                  >  2:137637/1‑134 (MQ=255)

                                                                                                                       aaTTGGAGGAGACGCAATTAGAACATAAATATTTTGAGTGAAGCCAAGATATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTTAAAAAATATATAAACGATTACTaa                <  2:1406884/134‑1 (MQ=255)

                                                                                                                       aaTTGGACGAGATGCAATTGAAACGTACATAACTCAACTTAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGATAAAATATATAAACAATTACTaa                <  2:1362247/134‑1 (MQ=255)

                                                                                                                        aTTGAACGAGAAGAAATTGAAACGTACATAACTAAACAAAAGCAAAAGTATAAACCAACAGTTACAATTGTCAATGCTTAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAg               <  1:1823511/134‑1 (MQ=255)

                                                                                                                         ttGGACGAGACGCAATTGAAACGTACATACATCAACAGAAGCAAAAGTATAAACAAACAGTAACAATTGTAAATGCTGAAAATACAGCACATGGTAAAGAGTTAACTGAAAAAATATATAACAAATTACGAAGa              >  2:2461314/1‑134 (MQ=255)

                                                                                                                         tgagACAAGACGAAACTGAAAATGAAATACTTAAACAGAAGCAAAAGTATAAACCAACAGTTACAATTTTAAATGCTGAAAATTCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGa              <  2:458668/131‑1 (MQ=255)

                                                                                                                           ggACGCGAAGCAATTGATACGTACATACCTCAACTGAAGCAAAAGTATAAACCCACAGTTACAAATGTAAATGCTGAAAATGCAGCACATGGTAAAGGATAGACTGCAAAAAAATATAAACAATCACCAAGaac            >  2:1197688/1‑133 (MQ=255)

                                                                                                                                 gaCGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTCACTGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATCACTAACAAATGctgt      >  2:1805263/1‑134 (MQ=255)

                                                                                                                                 gaCGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGAAATGgtgt      <  2:1907238/134‑1 (MQ=255)

                                                                                                                                  aagCAATAAAAACGTACATACCTCAACTAAAGCAAAAGTAAAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGAAATGGTGTa     <  1:2409191/132‑1 (MQ=255)

                                                                                                                                    gCAATTGAAACGTACATACCTCAACTGAAGAAAAAGTATAAAACAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGAAATGGTGTAGa   <  1:1182623/134‑1 (MQ=255)

                                                                                                                                     cAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGAAATGGTGTAGAt  >  1:1451191/1‑134 (MQ=255)

AAAAATAATTCAAAGAATTTAATATTTACTTTAAATATTAGAACGTATGAATTAAACTAGTATTTATGAGAGGATGAACAAAACATGAGAATAATGTTTATAGGGGATATCGTAGGTAAAATTGGACGAGACGCAATTGAAACGTACATACCTCAACTGAAGCAAAAGTATAAACCAACAGTTACAATTGTAAATGCTGAAAATGCAGCACATGGTAAAGGTTTGACTGAAAAAATATATAAACAATTACTAAGAAATGGTGTAGAT  >  USA300TCH1516_ALE/1313698‑1313964

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