BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,141,234	T→A	6.6%	S285T (TCA→ACA)	USA300HOU_RS05655 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,141,234	0	T	A	6.6%	781.3 / 10.2	317	S285T (TCA→ACA)	USA300HOU_RS05655	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (165/119); new base A (2/25); total (168/149)
Fisher's exact test for biased strand distribution p-value = 1.80e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 337 total aligned reads displayed.

AAGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCAATATATG  >  USA300TCH1516_ALE/1141094‑1141374

aaGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAGATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAAt                                                                                                                                              >  1:2370210/1‑141 (MQ=255)

aaGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAAt                                                                                                                                              >  1:1178667/1‑141 (MQ=255)

aaGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAAt                                                                                                                                              >  1:3636598/1‑141 (MQ=255)

 aGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACATTAAAAATATAAGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAAta                                                                                                                                             >  2:2654092/1‑140 (MQ=255)

  gTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCGTCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCa                                                                                                                                            >  2:2265786/1‑141 (MQ=255)

  gTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAACTCa                                                                                                                                            >  1:576073/1‑141 (MQ=255)

  gTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCa                                                                                                                                            <  1:4184182/141‑1 (MQ=255)

  gTGCAAAAAAACAATAATATAAATATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCa                                                                                                                                            <  1:2752843/141‑1 (MQ=255)

      aagaaGAAAGGAAGAAAAAGATAGAAAATGAGCCGATCTTAA‑CAGGAGGTTCATGAATAGGTGAAAATAGCCCTCCAACAAGAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCaaaaa                                                                                                                                        <  2:1786492/138‑1 (MQ=255)

         aaGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAgg                                                                                                                                     >  1:3479175/1‑141 (MQ=255)

         aaGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATAACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAgg                                                                                                                                     >  1:3289934/1‑141 (MQ=255)

           gACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaa                                                                                                                                   >  2:1627084/1‑141 (MQ=255)

           gACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaa                                                                                                                                   <  1:3767405/141‑1 (MQ=255)

            aCAGTAATAAAAAGAAAGAAAATGAG‑CGACCTAAAACAACAAGTACATTAAAAATTGAAAATAGAACTCAAACAATAAAAAAGGAAGAAAATAATCAAAAAGAAAAACCAAAACATAAAGATGAAAAATCAAAAAAGgaaa                                                                                                                                  <  2:2497650/141‑1 (MQ=255)

              aataaaaaAAAGAATGAAAATGAA‑AAACCTAAAGAATCAGATACATTAAAAGTAGAAAAAAGACCACCAACAATAAAAAAGGATAAAAATAAAAAAAAAGAGAAAACGAAACATAAAGATGAAAAATAAAAAAAGgaaaag                                                                                                                                <  2:2820318/139‑1 (MQ=255)

               gTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaaga                                                                                                                               >  1:4231623/1‑141 (MQ=255)

               gTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaaga                                                                                                                               <  1:2629028/141‑1 (MQ=255)

                 aataaAAAGATTGAAAATGAG‑CGACTTAAAGCATCAGGTAAATTAAAAGATGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAGAGCGCAACCGAAACAAAAAGATGAAAAATCAAAAAAGgaaaagaaa                                                                                                                             >  2:929267/1‑141 (MQ=255)

                 aataaAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGgaaaagaaa                                                                                                                             <  1:879272/141‑1 (MQ=255)

                         gATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAagtagt                                                                                                                     >  2:2081574/1‑141 (MQ=255)

                             gAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAACGGCACAGAAAAAAGTAGTTgaa                                                                                                                 >  1:323110/1‑141 (MQ=255)

                              aaaaTGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACAAAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAATTgaaa                                                                                                                >  2:185759/1‑141 (MQ=255)

                              aaaaAGAA‑CGACCCACAGCATCAGGTCCATTAAAAGTTCAAAATAGCCCTCCAAAAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAAATGAAAAATCAAAAAATGAAAAGAAAAAAGTAGTTgaaa                                                                                                                <  1:302960/141‑1 (MQ=255)

                               aaaTGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaa                                                                                                               >  2:2369203/1‑141 (MQ=255)

                                 aTGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaag                                                                                                             <  1:4287942/141‑1 (MQ=255)

                                    aG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaagaaa                                                                                                          <  1:1613824/141‑1 (MQ=255)

                                       cGACCTAAAGCATCAGGTACATTAAAAGTTGACAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaagaaaaa                                                                                                        >  2:3921890/1‑141 (MQ=255)

                                       cGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaagaaaaa                                                                                                        >  2:2934045/1‑141 (MQ=255)

                                       cGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTgaaaaagaaaaa                                                                                                        >  1:2017835/1‑141 (MQ=255)

                                            caaaGCATCACGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGtt                                                                                                   <  1:1889895/140‑1 (MQ=255)

                                             aaaGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAAAAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTa                                                                                                  <  1:1285496/141‑1 (MQ=255)

                                                gCATAAGGTAAATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCa                                                                                               <  1:2421929/141‑1 (MQ=255)

                                                acATCAGGTAAATAAAAAGATGAAAATAGCCATACAAAAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATAAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCa                                                                                               <  2:1717581/140‑1 (MQ=255)

                                                  aTCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGc                                                                                             >  1:3703512/1‑141 (MQ=255)

                                                      ggTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCtttt                                                                                         >  1:4087336/1‑141 (MQ=255)

                                                      ggTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCtttt                                                                                         <  1:133199/141‑1 (MQ=255)

                                                       gTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTa                                                                                        >  1:4118771/1‑141 (MQ=255)

                                                        tACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGAAGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTaa                                                                                       <  2:131452/141‑1 (MQ=255)

                                                         aCATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAAt                                                                                      >  1:2246364/1‑141 (MQ=255)

                                                           aTTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGCAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGAAGTTGAAAAAGAAAAAGCGTAACCAGCTTTAAATaa                                                                                    >  2:126836/1‑140 (MQ=255)

                                                            ttAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATaga                                                                                   >  1:791816/1‑141 (MQ=255)

                                                            ttAAAAGGTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAAAATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGAAGTTGAAAAAGAAAAAGCATTACCACAATTTAATaga                                                                                   >  2:2777243/1‑141 (MQ=255)

                                                             aaaaaATAGAAAATAACCATACAACAAAAAAAAAGGTTAAAAATAATCACAAAAAGCAACAAAAAAATAAAGAAAAAAAAACAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATagag                                                                                  <  2:4356225/140‑1 (MQ=255)

                                                              aaaaGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAgaga                                                                                 >  1:1973856/1‑141 (MQ=255)

                                                                 aGTTGAAAACAGCCCTCCAACAATAAAAAAGGTTGAAAATAAGCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAGAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatga                                                                              <  1:82975/141‑1 (MQ=255)

                                                                     gAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAgatgatgat                                                                          >  1:1132349/1‑141 (MQ=255)

                                                                       aaaTAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAg                                                                        <  1:681812/141‑1 (MQ=255)

                                                                          tAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCaa                                                                     <  1:2858721/141‑1 (MQ=255)

                                                                           aGCCCTCCAACAATAAAAAAGTTTGATAATAATCACAAAGAGCAACCGAAACATAAAGCTGCAAAATCAGAAATGGAAAAGAAAAAAGTAGTTTAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCaaa                                                                    >  2:3585088/1‑140 (MQ=255)

                                                                           aGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAg                                                                    >  1:2114889/1‑141 (MQ=255)

                                                                            aacccacaacaaTAAAAAAGGTTGAAAATAATCAAAAAGAGCAAACGAAACATAAAGAAGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGa                                                                   <  2:1178667/135‑1 (MQ=255)

                                                                               cTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATGGCAAGAGta                                                                >  2:1531911/1‑141 (MQ=255)

                                                                                  aaacaaAAAAAAAGATTAAAAATAATGACAAAGAACAAACTAAACATAAAGATGAAAAAGCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagta                                                             <  2:3116286/140‑1 (MQ=255)

                                                                                   aacaaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGGAAGAAtagtag                                                            >  1:1759434/1‑141 (MQ=255)

                                                                                   aacaaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagtag                                                            >  2:3356617/1‑141 (MQ=255)

                                                                                   aacaaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagtag                                                            <  1:2188615/141‑1 (MQ=255)

                                                                                   aaaaagaaaaaaGGTTAAAAATAATCAAAAAGAGCAACAGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAAtagtag                                                            <  2:196043/135‑1 (MQ=255)

                                                                                    acaaTAAAAAAGGTTGCAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATagtagt                                                           >  1:92089/1‑141 (MQ=255)

                                                                                      aaTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGAGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCa                                                         >  2:616043/1‑141 (MQ=255)

                                                                                       aaaaaaaaAGTAAAAAAAAAACAAAAAGAGAAAAAAAAACAGAAAGAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCaa                                                        <  2:1587425/139‑1 (MQ=255)

                                                                                       aTAAAAAATGTAAAAAATAATAACAAAAAGCAACCGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCaa                                                        <  2:3479175/141‑1 (MQ=255)

                                                                                       aTAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCaa                                                        >  1:509505/1‑141 (MQ=255)

                                                                                        aaaaaaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAATAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAAt                                                       <  1:1464621/140‑1 (MQ=255)

                                                                                        aaaaaaaGGTAGAAAAAAATCAAAAAGAGCAACAGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAAt                                                       <  2:2185231/140‑1 (MQ=255)

                                                                                         aaaaaaGGTTGAAAATAATCACAAATAACAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAAtt                                                      >  1:1065403/1‑141 (MQ=255)

                                                                                         aaaaaaAGTAAAAAAAAAAAAAAAAGAGAAACAGAAAAATAAAAAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAAtt                                                      <  2:2736706/141‑1 (MQ=255)

                                                                                           aaaaGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTAt                                                    >  1:172014/1‑141 (MQ=255)

                                                                                             aaGGTTGAAATTAATCACCAAGAGCAAACGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCt                                                  <  1:288922/141‑1 (MQ=255)

                                                                                                  tGAAAATAAACAAAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGa                                             <  1:3738076/141‑1 (MQ=255)

                                                                                                    aaaaTAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGata                                           >  1:3970320/1‑141 (MQ=255)

                                                                                                         aaaaacaAAGAGAAAACGAAAAATAAAGAAGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTaaa                                      <  2:3639739/137‑1 (MQ=255)

                                                                                                          atCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAg                                     >  2:554704/1‑141 (MQ=255)

                                                                                                          atCACAAAGAGCAACCGAAACATAAAGATGAAAAAACAAAGAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTGTAATAGAGATGATGATAGCAAGAATAGTAGTAAATTATCTAGTGTTATTaacg                                     >  2:3696267/1‑139 (MQ=255)

                                                                                                            cacaAAGAGCAACCGAAACATAAAGATGAAAAATCAAACAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGaa                                   >  1:1369985/1‑141 (MQ=255)

                                                                                                            cacaAAGAGCAACCGAAACATAAAAATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGaa                                   <  1:264436/141‑1 (MQ=255)

                                                                                                            aacaAAGAGAAACCGAAACATAAAAAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGaa                                   <  2:1876254/140‑1 (MQ=255)

                                                                                                            aaaaaaGAGAAAAAGAAACATAAAGATGAAAAATAAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGaa                                   <  2:3738148/138‑1 (MQ=255)

                                                                                                                 agagCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtga                              <  1:1822000/141‑1 (MQ=255)

                                                                                                                  gagCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATCGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgat                             >  1:3448721/1‑141 (MQ=255)

                                                                                                                   agCAAACGAAAAATAAAGAAGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatg                            <  2:4260356/141‑1 (MQ=255)

                                                                                                                    gCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatga                           <  1:539707/141‑1 (MQ=255)

                                                                                                                    gCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTtgatga                           <  2:4231623/141‑1 (MQ=255)

                                                                                                                      aaCCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAAc                         >  1:1567492/1‑141 (MQ=255)

                                                                                                                          aaaaCAAAAAGATGAAAAAAAAAAAAAGGAAAAGAAAAAAGTAAAAAAAAAAAAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCaaa                     <  2:2592920/140‑1 (MQ=255)

                                                                                                                           aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAAt                    >  2:646323/1‑141 (MQ=255)

                                                                                                                           aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAAt                    >  1:2954932/1‑141 (MQ=255)

                                                                                                                           aaaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAAt                    <  1:2513037/141‑1 (MQ=255)

                                                                                                                            aaCATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATc                   <  2:3970320/141‑1 (MQ=255)

                                                                                                                              cATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCAt                 <  1:3563044/141‑1 (MQ=255)

                                                                                                                                tAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTATTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGTTATTAAAGAACTTGATGAACCAAATCATaa               >  2:825913/1‑141 (MQ=255)

                                                                                                                                aaaaGATAAAAAAAAAAAAAAGAAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaa               <  2:2388066/140‑1 (MQ=255)

                                                                                                                                aaaaGAAAAAAAAAAAAAAAAGGAAAAGAAAAAAGTAGTTGATAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATaa               <  2:4187224/140‑1 (MQ=255)

                                                                                                                                 aaaGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACAAGATATTAATAGAGATGATGATAGAACGCACAGCTATCAAGCATATAGTGAGACAGAAGAACTTGATGAACCAAATCATaaa              >  1:2062988/1‑141 (MQ=255)

                                                                                                                                  aaGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTCTAATAGGTTTGGTGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACAAAATCATaaaa             >  1:4390551/1‑141 (MQ=255)

                                                                                                                                    aaagAAAAAACAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAg           <  1:337367/138‑1 (MQ=255)

                                                                                                                                       gAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCaa        >  1:2995447/1‑141 (MQ=255)

                                                                                                                                       gAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGACGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCaa        >  1:3158764/1‑141 (MQ=255)

                                                                                                                                       aaaaaaaCAAAAAAAGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCaa        <  2:7397/140‑1 (MQ=255)

                                                                                                                                       aaaaaaTCAAAAAATGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCTGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCaa        <  2:2899847/140‑1 (MQ=255)

                                                                                                                                           aaTCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAACAAGCAatata    >  2:107908/1‑141 (MQ=255)

AAGTGCAACAAGACAGTAATAAAAAGATTGAAAATGAG‑CGACCTAAAGCATCAGGTACATTAAAAGTTGAAAATAGCCCTCCAACAATAAAAAAGGTTGAAAATAATCACAAAGAGCAACCGAAACATAAAGATGAAAAATCAAAAAAGGAAAAGAAAAAAGTAGTTGAAAAAGAAAAAGCGTTACCAGCTTTTAATAGAGATGATGATAGCAAGAATAGTAGTCAATTATCTAGTGATATTAAAGAACTTGATGAACCAAATCATAAAAAGCAATATATG  >  USA300TCH1516_ALE/1141094‑1141374

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