BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,668,042	A→G	5.3%	Y50Y (TAT→TAC)	USA300HOU_RS13735 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,668,042	0	A	G	5.3%	508.4 / 11.6	165	Y50Y (TAT→TAC)	USA300HOU_RS13735	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (63/92); new base G (6/3); total (70/95)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 165 total aligned reads displayed.

TAAACCAACAATGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTACACCTAACA  >  USA300TCH1516_ALE/2667904‑2668181

tAAACCAACAATGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATa                                                                                                                                           >  2:1046899/1‑141 (MQ=255)

 aaaCCAACAATGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATaa                                                                                                                                          <  1:528838/141‑1 (MQ=255)

   aCCAACAATGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACt                                                                                                                                        <  2:46342/141‑1 (MQ=255)

    ccAACAATGACCGCAACTACTGATGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACtt                                                                                                                                       <  4:317205/141‑1 (MQ=255)

       acaaTGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAg                                                                                                                                    <  2:198433/141‑1 (MQ=255)

        caaTGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGa                                                                                                                                   <  1:799142/141‑1 (MQ=255)

         aaTGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGaa                                                                                                                                  >  3:275290/1‑141 (MQ=255)

            gACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGc                                                                                                                               >  1:16854/1‑141 (MQ=255)

            gACCGCAAATACTGGTGCAAAGTCTAATTTTACGTACGTTAGAAGACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATACCTTGTAGATCTTTAATATCAATGACGCCACCAATAGCGTAACTTAATAATGc                                                                                                                               <  4:200908/141‑1 (MQ=255)

            aaCCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCGTAACTTAAGAATGc                                                                                                                               >  4:368322/2‑141 (MQ=255)

            aaCCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCGTAACTTAAGAATGc                                                                                                                               >  2:29784/2‑141 (MQ=255)

             aCCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGCAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCt                                                                                                                              <  2:326410/141‑1 (MQ=255)

             aCCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCt                                                                                                                              <  1:135987/141‑1 (MQ=255)

                   aaTACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGt                                                                                                                        >  3:233706/1‑141 (MQ=255)

                    aTACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTa                                                                                                                       >  2:995597/1‑141 (MQ=255)

                    aTACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTa                                                                                                                       <  2:581693/141‑1 (MQ=255)

                     tACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAg                                                                                                                      <  1:25478/141‑1 (MQ=255)

                      tcTGGTGCAACTCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAATCTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGa                                                                                                                     <  2:786767/140‑1 (MQ=255)

                        tGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGaaa                                                                                                                   >  2:632359/1‑141 (MQ=255)

                          gTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATg                                                                                                                 >  2:959511/1‑141 (MQ=255)

                                cGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGcac                                                                                                           <  2:829540/141‑1 (MQ=255)

                                cGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGcac                                                                                                           <  3:389022/141‑1 (MQ=255)

                                 gCCTAAGTTTAAGCACGTTAGAATACCTGCTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGcaca                                                                                                          <  3:56996/141‑1 (MQ=255)

                                 gCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGcaca                                                                                                          <  2:977681/141‑1 (MQ=255)

                                   cTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAg                                                                                                        >  4:171014/1‑141 (MQ=255)

                                   cTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAg                                                                                                        >  2:12833/1‑141 (MQ=255)

                                     aaGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCt                                                                                                      >  2:372330/1‑141 (MQ=255)

                                        tttAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTAGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCa                                                                                                   >  2:1041282/1‑141 (MQ=255)

                                          tAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCAcc                                                                                                 <  1:321389/141‑1 (MQ=255)

                                          tAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCAcc                                                                                                 <  1:831732/141‑1 (MQ=255)

                                            aGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCGTAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCActgg                                                                                               <  2:703807/141‑4 (MQ=255)

                                            aGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCgg                                                                                               >  1:142888/1‑141 (MQ=255)

                                            aGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCgg                                                                                               <  1:275305/141‑1 (MQ=255)

                                             gCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGa                                                                                              <  2:195955/141‑1 (MQ=255)

                                             gCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTAGTAGAAATGTCGCACAAGCTGCACCGGa                                                                                              >  2:787211/1‑141 (MQ=255)

                                              cACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGAt                                                                                             <  1:54240/141‑1 (MQ=255)

                                                cGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATaa                                                                                           >  1:726347/1‑141 (MQ=255)

                                                cGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATaa                                                                                           <  2:327690/141‑1 (MQ=255)

                                                    aGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCt                                                                                       <  1:628300/141‑1 (MQ=255)

                                                     gAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCtt                                                                                      <  3:339277/141‑1 (MQ=255)

                                                      aaTACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTg                                                                                     <  3:389136/141‑1 (MQ=255)

                                                       aTACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGt                                                                                    <  3:15265/141‑1 (MQ=255)

                                                       aTACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGt                                                                                    <  1:929587/141‑1 (MQ=255)

                                                       aTACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGt                                                                                    <  1:802154/141‑1 (MQ=255)

                                                        tACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTa                                                                                   <  1:935641/141‑1 (MQ=255)

                                                        tACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTa                                                                                   >  4:35903/1‑141 (MQ=255)

                                                          ccTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGt                                                                                 <  1:593439/141‑1 (MQ=255)

                                                          ccTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGt                                                                                 >  3:286071/1‑141 (MQ=255)

                                                             aCTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCt                                                                              >  2:741553/1‑141 (MQ=255)

                                                              cTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCtt                                                                             >  3:114493/1‑141 (MQ=255)

                                                              cTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCtt                                                                             >  2:331828/1‑141 (MQ=255)

                                                              cTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCtt                                                                             >  4:173803/1‑141 (MQ=255)

                                                                   aTACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTa                                                                        >  3:347140/1‑141 (MQ=255)

                                                                   aTACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTa                                                                        <  3:427847/141‑1 (MQ=255)

                                                                   aTACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTa                                                                        <  4:286071/141‑1 (MQ=255)

                                                                   aTACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTa                                                                        <  2:831486/141‑1 (MQ=255)

                                                                      ccGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCg                                                                     >  1:64458/1‑141 (MQ=255)

                                                                      ccGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCg                                                                     <  1:307077/141‑1 (MQ=255)

                                                                            aaaCTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGt                                                               <  1:720669/141‑1 (MQ=255)

                                                                             aaCTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTg                                                              >  1:268161/1‑141 (MQ=255)

                                                                              aCTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGc                                                             <  3:226358/141‑1 (MQ=255)

                                                                               cTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCACATACGCCCCCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCa                                                            >  4:173701/1‑141 (MQ=255)

                                                                               cTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCa                                                            <  2:415846/141‑1 (MQ=255)

                                                                                 tCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATa                                                          >  4:166515/1‑141 (MQ=255)

                                                                                  ccTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATAc                                                         <  1:778990/141‑1 (MQ=255)

                                                                                   cTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCTACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACt                                                        <  4:42979/141‑1 (MQ=255)

                                                                                   cTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACt                                                        <  3:332292/141‑1 (MQ=255)

                                                                                   cTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACt                                                        <  1:292764/141‑1 (MQ=255)

                                                                                    tGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACtt                                                       <  1:753199/141‑1 (MQ=255)

                                                                                     gCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCGTAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACttt                                                      >  1:516051/1‑141 (MQ=255)

                                                                                     gCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCGTAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACttt                                                      >  1:828558/1‑141 (MQ=255)

                                                                                      cTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTa                                                     >  3:13550/1‑141 (MQ=255)

                                                                                             cccAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAg                                              >  2:288770/1‑141 (MQ=255)

                                                                                              ccAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGc                                             <  4:267400/141‑1 (MQ=255)

                                                                                              ccAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGc                                             <  2:637049/141‑1 (MQ=255)

                                                                                               cAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCg                                            <  1:86051/141‑1 (MQ=255)

                                                                                                aaTATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGt                                           <  2:904089/141‑1 (MQ=255)

                                                                                                aaTATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGt                                           <  1:995259/141‑1 (MQ=255)

                                                                                                 atatCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTc                                          >  4:49248/1‑141 (MQ=255)

                                                                                                 atatCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTc                                          >  2:240478/1‑141 (MQ=255)

                                                                                                      ttGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGAt                                     <  1:597858/141‑1 (MQ=255)

                                                                                                      ttGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGAt                                     >  3:155644/1‑141 (MQ=255)

                                                                                                         tAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCAc                                  <  1:495706/141‑1 (MQ=255)

                                                                                                          aGAAATTTAATATCAAATCCGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACt                                 <  4:142841/141‑1 (MQ=255)

                                                                                                                  aaTATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCaa                         >  2:268773/1‑141 (MQ=255)

                                                                                                                  aaTATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCaa                         >  2:524213/1‑141 (MQ=255)

                                                                                                                  aaTATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCaa                         <  1:694379/141‑1 (MQ=255)

                                                                                                                   ataaCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAg                        <  4:162621/141‑1 (MQ=255)

                                                                                                                     atCAAATACGCCACCAATAGCGTAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGAc                      >  1:443054/1‑141 (MQ=255)

                                                                                                                     atCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGAc                      >  1:658127/1‑141 (MQ=255)

                                                                                                                     atCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGAc                      >  1:685142/1‑141 (MQ=255)

                                                                                                                     atCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGAc                      >  2:645590/1‑141 (MQ=255)

                                                                                                                       cAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTa                    <  1:574484/141‑1 (MQ=255)

                                                                                                                             cGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGc              >  1:1058230/1‑141 (MQ=255)

                                                                                                                                caccaATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTa           >  2:256286/1‑141 (MQ=255)

                                                                                                                                 accaATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATTCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTac          >  4:255558/1‑141 (MQ=255)

                                                                                                                                  ccaATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCGCAAGCCGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTaca         <  1:288416/141‑1 (MQ=255)

                                                                                                                                    aaTAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTACAcc       <  1:378462/141‑1 (MQ=255)

                                                                                                                                       aGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTACACCTaa    >  4:339225/1‑141 (MQ=255)

                                                                                                                                         cATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTACACCTAACa  <  2:949757/141‑1 (MQ=255)

TAAACCAACAATGACCGCAAATACTGGTGCAACGCCTAAGTTTAAGCACGTTAGAATACCTACTGCGATACCGCTTAAACTTCCTGCTAAATCCCCAATATCTTGTAGAAATTTAATATCAAATACGCCACCAATAGCATAACTTAAGAATGCTTGTGGTAGAAATGTCGCACAAGCTGCACCGGATAATGCTTGTAGTCCTTGTTTACCGTACGGTGCATACTTTAAAAATAGCGTCATGATCACTAAAACCAAGACTAATGTGCCTACACCTAACA  >  USA300TCH1516_ALE/2667904‑2668181

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