BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	2,455,698	2 bp→AT	coding (1216‑1217/1260 nt)	lyrA →	Lysostaphin resistance protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,455,698	0	C	A	89.9%	391.5 / 37.9	132	coding (1216/1260 nt)	lyrA	Lysostaphin resistance protein A
Reads supporting (aligned to +/- strand): ref base C (0/2); major base A (44/72); minor base . (6/7); total (51/81)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	USA300TCH1516_ALE	2,455,699	0	A	T	97.8%	397.7 / ‑4.1	137	H406L (CAT→CTT)	lyrA	Lysostaphin resistance protein A
Reads supporting (aligned to +/- strand): ref base A (1/2); new base T (50/82); total (51/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

AGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTG  >  USA300TCH1516_ALE/2455559‑2455841

||

agcagAAGCAGTTAAACAACCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGaaa                                                                                                                                        <  2:1092155/148‑1 (MQ=255)

 gcagAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGaaaa                                                                                                                                       >  1:255749/1‑148 (MQ=255)

           ttAAACAAGCACCACAAAGTGACGGTTTGACAAACCAGCCAAATGAAGATGCAATAGAGCAATCCTTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATCATCGAAAAAACTaaagaa                                                                                                                             >  1:337103/1‑148 (MQ=255)

           taaaaCAAGTACCACAAAGGGACGAATTAAAAAAAGATACAAATGAAGATGAAATAGAGGAATCATTAAAAGAACCTGTGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAAAAAAAAAAAAAACTgaagaa                                                                                                                             <  2:742453/146‑1 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACATACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGTA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  2:790247/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCACTCATATAAATAAACTGCGGTGCGCAAATAAGACAGACGTTCATCAGTTGTAAGTTATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  1:1031045/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAAAAATc                                                                                                                        >  1:597604/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaagaaacaac                                                                                                                        >  1:711429/1‑143 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  1:1109203/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  1:1354204/1‑148 (MQ=255)

                  aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCtt                                                                                                                      <  2:1193660/148‑1 (MQ=255)

                  aGCACCACAAAGTGACCATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCtt                                                                                                                      <  2:143317/148‑1 (MQ=255)

                    caccacAAAGTGACGATTTAACAAACGATTCAAATGACGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATCATCGAAAAAACTGAAGAACAATCTTCa                                                                                                                    <  1:916886/148‑1 (MQ=255)

                    caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                                                                                    >  1:1057781/1‑148 (MQ=255)

                     actacaAAGTGACGATTTAACAAGCGATTAAAAAGAAGAAGAATCAGGGAAACCATTAAAATAACCTTCGACTTATAAAGAAGACAGACGTTCAATAGTTGTAATTGATGCAAAAAAA‑TTATCGACAAAACTTAAGAACAATCTTCAg                                                                                                                   <  2:992316/145‑1 (MQ=255)

                     accacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAAGCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                                   <  2:776136/148‑1 (MQ=255)

                      ccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                                  <  2:717151/148‑1 (MQ=255)

                      ccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                                  <  2:488463/148‑1 (MQ=255)

                      ccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                                  <  2:212281/148‑1 (MQ=255)

                       cacaAGGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCatat                                                                                                                 <  2:94089/148‑1 (MQ=255)

                       cacaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                                                                                 <  2:952090/148‑1 (MQ=255)

                         caAAGTGACGATTTAACAAACGATTCAAATCAAGAAGAAATAGGGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaa                                                                                                               <  2:768872/148‑1 (MQ=255)

                           aaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCAGTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                                             <  2:57385/148‑1 (MQ=255)

                           aaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATAATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                                             >  2:882457/1‑148 (MQ=255)

                            aGTGAAGATCTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTCTAGCACAAACTGAAGAACAATCTTCAGATaaaaa                                                                                                            <  2:1362688/148‑1 (MQ=255)

                               gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AATAAATTATCGGAAAAACTAAAGAACAATTCTCAGGTAAGAACaa                                                                                                         >  1:855357/1‑148 (MQ=255)

                               gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCAGCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                                                                                         >  1:737653/1‑148 (MQ=255)

                                aCGATTTAACAAACGATTCAAAGGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaaa                                                                                                        <  2:343432/148‑1 (MQ=255)

                                 cGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                                       >  1:256921/1‑148 (MQ=255)

                                   aTTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataa                                                                                                     <  2:1169366/148‑1 (MQ=255)

                                    tttAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                                    <  2:382573/148‑1 (MQ=255)

                                    tttAACAAACGATTCAAAAGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAACGAAGACAGACGTTCATCATTTGTAATGGATGAAAAAAAAC‑CAACGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                                    <  1:737578/148‑1 (MQ=255)

                                    tttAACAAAAGATTCAAATCAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                                    <  2:1099431/148‑1 (MQ=255)

                                     ttAACCAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                                   <  1:1072474/148‑1 (MQ=255)

                                     ttAACAAAAGAATTAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                                   <  1:735929/148‑1 (MQ=255)

                                             cGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcg                                                                                           >  1:1273676/1‑148 (MQ=255)

                                              gATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgc                                                                                          >  2:1049915/1‑148 (MQ=255)

                                                 tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                                                                       <  2:835086/148‑1 (MQ=255)

                                                      tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGTAACAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  1:271338/148‑1 (MQ=255)

                                                      tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACCTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  1:1190908/148‑1 (MQ=255)

                                                      tGAAGATGAAAGAGAGCAATCATTAACAGAACCTGCAACATATAAAGAAGACAGACGTTTATCAGTTGAAAAAAAAAAA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  2:1398682/148‑1 (MQ=255)

                                                        aaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                                                                                >  2:481991/1‑148 (MQ=255)

                                                         aGATGAAATATAGCAATCATGAAAAGATCCTGCGACTTATAAAGAAGACATACGTTCATCAGTTTCAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAcc                                                                               <  2:60140/148‑1 (MQ=255)

                                                           aTGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTATTACctc                                                                             >  1:680160/1‑148 (MQ=255)

                                                            tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                            <  1:1146270/148‑1 (MQ=255)

                                                             gAAATAGAGTCATCATTAAAAGAACCTGTGACTTAGAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGACAAAACTGAAGAACAATCTTCAGATAAAAACAAATGATAACAAAGCGCTTGCTAGTACCTCtt                                                                           <  2:1083548/148‑1 (MQ=255)

                                                             gAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCtt                                                                           <  1:805428/148‑1 (MQ=255)

                                                             gAAATAGAGCAAACATTAAAAGGACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTTAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCtt                                                                           <  2:912988/148‑1 (MQ=255)

                                                              aaaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                                                                          >  1:868400/1‑148 (MQ=255)

                                                               aagagagCAACCATTAAAAGGACCTGGGACTTATAAAGAAGACGGACGGTGGGGGGTGGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                                                                         <  1:960175/145‑1 (MQ=255)

                                                               aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAACAACAAATAATAACAAAGCGCTTACTAGTACCTCTTaa                                                                         >  1:1110124/1‑148 (MQ=255)

                                                                  agagCACTCATTAAAAGCACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                                      <  1:316006/148‑1 (MQ=255)

                                                                      cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATACAAAAAAAGCGCTTCGTAGTCACTCTTAAAAAgaga                                                                  >  2:348165/1‑146 (MQ=255)

                                                                             tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAg                                                           >  2:1260200/1‑148 (MQ=255)

                                                                             tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGAAAATGCTAg                                                           >  2:1347094/1‑148 (MQ=255)

                                                                               aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCTAGTACCTCTTAAAAAGATGATGCTAGCa                                                         >  2:197109/1‑148 (MQ=255)

                                                                               aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAAATGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATAATGCTAgaa                                                         >  2:660503/1‑146 (MQ=255)

                                                                                 aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAg                                                       >  2:797933/1‑148 (MQ=255)

                                                                                  gAACCCGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgc                                                      <  1:4906/148‑1 (MQ=255)

                                                                                   aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAGTAACAAAGCGATCGATAGTACCTCTTAAAAAGATTATGCTAGCAAgcg                                                     <  1:1037281/148‑1 (MQ=255)

                                                                                    aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                                    <  2:980972/148‑1 (MQ=255)

                                                                                     ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                                   <  2:241143/148‑1 (MQ=255)

                                                                                      cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtt                                                  <  1:512127/148‑1 (MQ=255)

                                                                                      cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtt                                                  <  1:21793/148‑1 (MQ=255)

                                                                                        gCGACTTATAAAGAAGAAAGACGTGCATCAGGTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCGTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtttt                                                <  2:337103/148‑1 (MQ=255)

                                                                                         caaCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  2:1062018/146‑1 (MQ=255)

                                                                                         cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATACTAACAAAGCGCGTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  1:1331429/148‑1 (MQ=255)

                                                                                         cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  1:969764/148‑1 (MQ=255)

                                                                                          gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGACGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                                              <  1:54711/148‑1 (MQ=255)

                                                                                          gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATGATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                                              <  1:1316950/148‑1 (MQ=255)

                                                                                           atttATAACGAAGACAGACGCTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                             <  1:119995/146‑1 (MQ=255)

                                                                                           aCTTAAAAACAAGACAGACGTTCAGCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACATTTAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                             <  1:41646/148‑1 (MQ=255)

                                                                                            cTTATACAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  2:903934/148‑1 (MQ=255)

                                                                                            cTTATAAAGAAGACGGACGTTCAACAGGTGTAATTGATGCA‑AAAAAATGATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAACAATAGCGAGGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  1:46421/148‑1 (MQ=255)

                                                                                            cTTATAAAGAAGACCGACGTTTATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATCACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  1:46801/148‑1 (MQ=255)

                                                                                               ataAAGAAGACAGACCTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                                         <  2:604551/148‑1 (MQ=255)

                                                                                                taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactc                                        >  1:444377/1‑147 (MQ=255)

                                                                                                taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatacta                                        >  2:950846/1‑148 (MQ=255)

                                                                                                   agaagaCAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtagat                                     >  1:97035/1‑148 (MQ=255)

                                                                                                      agacagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtatt                                  >  1:124600/1‑148 (MQ=255)

                                                                                                         cagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAGCTGAAGAACAATATTCAGATAAAAACAAATAATAACCAAGCGCCTGCTAGCACCTCTTAACAAGAAGATCCTACCCAGCGCTTTTCCATCCTACACattttt                               >  2:1264515/1‑148 (MQ=255)

                                                                                                         cagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATattatt                               >  1:509926/1‑148 (MQ=255)

                                                                                                          agacGTTCATCAGTTGTAATTGATGCA‑AAAAAATGATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAttt                              >  2:873506/1‑148 (MQ=255)

                                                                                                           gacTTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGCTAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAtttt                             <  1:949920/148‑1 (MQ=255)

                                                                                                           gacGTTCATCAGTTGGAAAAAAAACA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAtttt                             <  2:1273676/148‑1 (MQ=255)

                                                                                                                  atcaGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaaa                      <  2:256921/148‑1 (MQ=255)

                                                                                                                  atcaGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaaa                      <  2:801241/148‑1 (MQ=255)

                                                                                                                   tcaGTTGTAATTGATGCA‑AAAAAATTATCGAANAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAt                     <  1:159242/148‑1 (MQ=255)

                                                                                                                   aaaGTTGTAATTAATGCAAAAAAAC‑TATCGAAAAAACTAAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATTATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAt                     <  1:999665/146‑1 (MQ=255)

                                                                                                                       ttGTAATTGATG‑AAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTa                 <  1:548658/148‑1 (MQ=255)

                                                                                                                        tGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTaa                >  1:435714/1‑148 (MQ=255)

                                                                                                                        tGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCATCAGATAAAAACAAATAATAACAAAACGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATCCTATATATTATTTATCTAAAAATTTaa                >  1:194381/1‑148 (MQ=255)

                                                                                                                          tAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATAATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTCAAATTGTAcca              >  1:268498/1‑148 (MQ=255)

                                                                                                                            agtGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAcc            <  2:119996/146‑1 (MQ=255)

                                                                                                                                 atgcaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCctct        <  1:1152721/143‑1 (MQ=255)

                                                                                                                                 atgcaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCctct        <  1:453538/143‑1 (MQ=255)

                                                                                                                                      aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGt   <  1:1147573/148‑1 (MQ=255)

                                                                                                                                      aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGt   <  1:1273214/148‑1 (MQ=255)

                                                                                                                                      aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGt   <  1:809975/148‑1 (MQ=255)

                                                                                                                                       aaaaaaGATAGAAAAAACTTAAGCACAATCTTCAGATAAAAACAAATAATACCAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTg  <  2:379618/148‑1 (MQ=255)

||

AGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTG  >  USA300TCH1516_ALE/2455559‑2455841

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