BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,333,789	T→G	6.6%	Y64D (TAT→GAT)	hflX →	GTPase HflX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,333,789	0	T	G	6.6%	897.4 / 5.5	329	Y64D (TAT→GAT)	hflX	GTPase HflX
Reads supporting (aligned to +/- strand): ref base T (95/202); new base G (27/2); total (123/206)
Fisher's exact test for biased strand distribution p-value = 7.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

TCTTAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAATTATTGAT  >  USA300TCH1516_ALE/1333649‑1333929

  ttAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAatat                                                                                                                                             <  1:1843325/141‑1 (MQ=255)

      tCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATg                                                                                                                                         >  1:2033189/1‑141 (MQ=255)

      tCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATg                                                                                                                                         <  1:1424214/141‑1 (MQ=255)

              cATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAACCCAGAGCTCGTTTAGATCGCAAATATTATCTTggcaaa                                                                                                                                 >  2:3342172/1‑141 (MQ=255)

              cATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaa                                                                                                                                 <  1:2361695/141‑1 (MQ=255)

                tGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaagg                                                                                                                               <  1:107783/141‑1 (MQ=255)

                tGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaagg                                                                                                                               <  1:3143837/141‑1 (MQ=255)

                tGCACAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaagg                                                                                                                               <  2:2570046/141‑1 (MQ=255)

                 gCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTggtaaaggt                                                                                                                              <  1:3498978/141‑1 (MQ=255)

                 gCTCAAGATGAGAAGCAATTTAACTTTGGGTCGACAATGGAGGAATTATCATCTTTCTAAGAGAGGGGCGTACTTGGAGTGTTTGGTCAAATTACTCAAAACAGGGATCGTGTAGAGCGCAAATATTATGTTggtaaaggt                                                                                                                              <  2:1577678/141‑1 (MQ=255)

                      agagGAGAAGCAATTTAATGTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAGCTTTAAGTGTTGGGTCAAATTACTCAAAACAGAGAGCGTGTAGAGCGCAAAGATTATGTTGGTAAAGGTAAAAt                                                                                                                         <  2:520036/141‑1 (MQ=255)

                        atgatAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATCTTGATAAAGCTCAAATTg                                                                                                                       >  2:479653/1‑141 (MQ=255)

                          gatAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAGGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTgaa                                                                                                                     <  2:104907/141‑1 (MQ=255)

                          gatAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTgaa                                                                                                                     <  1:620038/141‑1 (MQ=255)

                          gatAAGCAATTTAAGTTTGAGTATACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTgaa                                                                                                                     <  1:3129806/141‑1 (MQ=255)

                           atAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTgaag                                                                                                                    <  2:35837/141‑1 (MQ=255)

                               gCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAGAGATTACCTTGGCAAAAATAAAATTaaagaact                                                                                                                >  2:1948990/1‑139 (MQ=255)

                               gCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAAt                                                                                                                <  1:2642071/141‑1 (MQ=255)

                                 aatttaatttTGAGTCTAAAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTc                                                                                                              <  2:420255/141‑1 (MQ=255)

                                   tttaatttTGAGTCTACAATGGAGGAATTATCATCTTTATGTGAGACATACCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCaa                                                                                                            <  2:932290/141‑1 (MQ=255)

                                     taatttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGc                                                                                                          <  1:4078806/141‑1 (MQ=255)

                                     taatttTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGc                                                                                                          <  1:2308306/141‑1 (MQ=255)

                                        ttttGAGTATACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAGTTGAAGAAATTCAAGCAtt                                                                                                       <  1:2568986/141‑1 (MQ=255)

                                            gAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGGCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTAtt                                                                                                   <  2:2084487/141‑1 (MQ=255)

                                            gAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCGGAGATTGCGGTAGCATAAGAAAAAATAAAGAAATTCAATAATTTAtt                                                                                                   >  2:2159781/1‑141 (MQ=255)

                                            gAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTAtt                                                                                                   <  1:71496/141‑1 (MQ=255)

                                            gAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTAtt                                                                                                   <  1:2655960/141‑1 (MQ=255)

                                            gAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTAtt                                                                                                   <  1:2769617/141‑1 (MQ=255)

                                              gTCTGCAAGTGAGGTATTATCATCCTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGa                                                                                                 <  2:3081001/141‑1 (MQ=255)

                                                cTACAATGGAGTAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGt                                                                                               <  1:1172979/141‑1 (MQ=255)

                                                cTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGt                                                                                               <  1:2850325/141‑1 (MQ=255)

                                                cTACAATGGAGGAATTATCATCTTTATCAGAGACTAGCCAACTTGAAGTGTTGGGACTAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGt                                                                                               <  2:2910538/141‑1 (MQ=255)

                                                   cAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCa                                                                                            <  1:3776048/141‑1 (MQ=255)

                                                     aTGGAGGAATTATCATCTTTATTACAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCa‑aa                                                                                         <  2:371999/141‑1 (MQ=255)

                                                     aTGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCa‑aa                                                                                         <  1:2124651/141‑1 (MQ=255)

                                                     aTGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCa‑aa                                                                                         <  1:3658536/141‑1 (MQ=255)

                                                      tGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAg                                                                                        <  1:319837/141‑1 (MQ=255)

                                                      tGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAg                                                                                        <  2:4024615/141‑1 (MQ=255)

                                                        gaggaATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGat                                                                                      <  1:4127403/141‑1 (MQ=255)

                                                        gaggaATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGat                                                                                      >  1:216340/1‑141 (MQ=255)

                                                        gaggaATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGat                                                                                      <  1:4022278/141‑1 (MQ=255)

                                                           gaATTAGCATCTTTATCAGAGACTTGCCAACTTGAAGTGATGGGTCAAATTACTCCCCCCTGTCAACGTGGAGGTCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATAtt                                                                                   <  1:1027460/141‑1 (MQ=255)

                                                             aTTATAATTTTTATCAGAGATATGAGAATTTGAAGTGTTGCGTCAAATTACGCAAAACAGAGAACGTGTAGATCGCTAATATTATGTTGGTAAAGGTAAAATTGAAGTAATTCAAGCATTTATTGAGTTCA‑AAGATATTGa                                                                                 <  2:3157263/141‑1 (MQ=255)

                                                               tATCAGCTTAATAAAGGACTTGTCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATg                                                                               <  1:1649922/141‑1 (MQ=255)

                                                               tATCACCCTTATCAGAAACTTGCCAACTTGAAGTGTTGGGTCAAATTAATCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAATGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATg                                                                               <  1:876267/141‑1 (MQ=255)

                                                                atcatcTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGAATATTATCGTAGCAAAAAAAAAATAAAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGt                                                                              >  2:344210/1‑141 (MQ=255)

                                                                 tcatcTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGATCGCCGTAGCATTAAAAAAAAAAAAGAAATTCAAGCATTTATTGAGTTCA‑AAAATATAGATGTa                                                                             >  2:2727576/1‑141 (MQ=255)

                                                                   atcTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGt                                                                           <  2:1679059/141‑1 (MQ=255)

                                                                      tttATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcat                                                                        <  1:4223566/141‑1 (MQ=255)

                                                                       ttATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCGAATATTATGTTGGTAAAGGAAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatc                                                                       >  2:3242212/1‑141 (MQ=255)

                                                                       ttATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatc                                                                       <  2:4136017/141‑1 (MQ=255)

                                                                       ttATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatc                                                                       <  1:4315608/141‑1 (MQ=255)

                                                                       ttATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatc                                                                       <  1:1629875/141‑1 (MQ=255)

                                                                       ttAGCAGAGACTTGCTAACTCGAAGGGGTGGGTCAAATGACTCAAAAATGAGATAGTGTAGATAGCAAAAATTATGTGAGTAAAGAGAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatc                                                                       <  2:1420809/141‑1 (MQ=255)

                                                                        tATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatca                                                                      <  1:3335721/141‑1 (MQ=255)

                                                                        tATCAGAGACGTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGtcatca                                                                      <  1:3834790/141‑1 (MQ=255)

                                                                         aTCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTCCTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATcac                                                                     <  1:3589296/141‑1 (MQ=255)

                                                                          tCAGAGACTTGCCAACTTGAAGTGTTGGGTCAGATTACTCAAAACAGAGATCGTGTAGATCTCGGTGATCGCCGTATCAATAAAAAAATAAAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATcaca                                                                    >  2:3900753/1‑141 (MQ=255)

                                                                          tCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATcaca                                                                    <  1:1926773/141‑1 (MQ=255)

                                                                          tCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATcaca                                                                    <  1:497548/141‑1 (MQ=255)

                                                                            agagACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACaaa                                                                  >  1:2936350/1‑141 (MQ=255)

                                                                            agagACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAAGGTTACCTTGGCATAAGTAAAAATAAAGAAATTCAAGCATTTATTGAGTTCA‑AAAATATTGATGTAGTCACAACaaa                                                                  >  2:1519265/1‑141 (MQ=255)

                                                                            agaaACTTGCCAACTTGAAATTGTGGTTCAAATTACTCAAAAAAGAGATCGTGTAGATCGCAAATATTATGTTTGTAAAGGTAAAATTGAAGAAATTTTAGCAATTATTGAGTTCA‑AAGATATTGATGTAGTCATCACaaa                                                                  <  2:4306450/141‑1 (MQ=255)

                                                                               gaCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTACAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAatga                                                               >  1:300815/1‑141 (MQ=255)

                                                                               gaCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAatga                                                               <  2:2329791/141‑1 (MQ=255)

                                                                               gaCTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAatga                                                               >  1:1519651/1‑141 (MQ=255)

                                                                               gaCTTGCCAACTTGAAGTGATGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAatga                                                               <  1:2462162/141‑1 (MQ=255)

                                                                                   tGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCGGTATCATTAAAAAAAAAAAAGAAATTCAAGCATTTATTGAGTTCA‑AAGAAATTGATCTAGTCATCACAAATGATGaa                                                           >  2:1880058/1‑141 (MQ=255)

                                                                                   tGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCATAGATTATGTTGGCAAAAGAAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGaa                                                           >  2:1135607/1‑141 (MQ=255)

                                                                                    gCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAAt                                                          >  1:2902772/1‑141 (MQ=255)

                                                                                      cAACTAGAAATATTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTa                                                        <  1:987840/141‑1 (MQ=255)

                                                                                       aaCTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTaa                                                       <  1:1131600/141‑1 (MQ=255)

                                                                                         ttttaaGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACg                                                     <  2:4155802/137‑1 (MQ=255)

                                                                                            gAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACt                                                  <  1:1734536/141‑1 (MQ=255)

                                                                                                  ttGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACaa                                            >  1:206497/1‑141 (MQ=255)

                                                                                                   tGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAAt                                           <  1:1269545/141‑1 (MQ=255)

                                                                                                     ggTCAAATTACTCAAAACAGAGATCGTGTAGATCTCGGTGGTCGCCGTAGCATTAAAAAAAATAAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATcc                                         >  2:785131/1‑141 (MQ=255)

                                                                                                     ggTCAAATTACTCAAAACAGAGATCGTGTAGATCGCGGTGATCACCGTAGCAATAAAAAAAATAAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTACCGAATGCACAATcc                                         >  2:2511696/1‑141 (MQ=255)

                                                                                                     ggTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAACTTGCAGAAATTCAAGCATTTATCACCTTGA‑AGGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATcc                                         >  1:4176803/1‑141 (MQ=255)

                                                                                                     ggTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGCTCC‑AAGATATTGCTGTAGTCGTCATAGGTGATGAATTAACGACTGCACAATcc                                         >  1:609834/1‑141 (MQ=255)

                                                                                                       tCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCaa                                       <  1:1987011/141‑1 (MQ=255)

                                                                                                              aCTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTa                                >  1:3850840/1‑141 (MQ=255)

                                                                                                               cgcAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTaa                               <  2:3185048/139‑1 (MQ=255)

                                                                                                               cTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTaa                               <  1:535703/141‑1 (MQ=255)

                                                                                                               cTCAAAACAGAGATCGTGTAAATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTaa                               <  2:3505389/141‑1 (MQ=255)

                                                                                                                  aaaaCAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATg                            <  2:2461601/141‑1 (MQ=255)

                                                                                                                     aaagagATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAg                         <  2:1814728/139‑1 (MQ=255)

                                                                                                                       agagATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCt                       >  1:130919/1‑141 (MQ=255)

                                                                                                                         agaTCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGCGTCCACATGTCCGTGGT‑TGGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCttt                     >  1:3232404/1‑141 (MQ=255)

                                                                                                                         agaTCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCtat                     >  1:2980161/1‑139 (MQ=255)

                                                                                                                                gtAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGCGTTTT‑ATAAGAATGAAGTAGTCATCACAAATTATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTa              <  1:4376266/141‑1 (MQ=255)

                                                                                                                                gtAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTa              <  2:1605632/141‑1 (MQ=255)

                                                                                                                                gtAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTa              <  2:1279616/141‑1 (MQ=255)

                                                                                                                                gtAGATCGCAAATATGATGTTTGTAAAGGTAAGATTGAAGGAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTa              <  2:69720/141‑1 (MQ=255)

                                                                                                                                 gaGATAGCAAATATTGTGTTGGTAAAGGTAAAATTGAAGAAATTAAAGCATGTATTGAGTTCA‑AAGATATTGATGTAGTCATCTCCAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTaa             <  2:1348405/140‑1 (MQ=255)

                                                                                                                                   gATCGTAAATATTATGTTGGTAAATGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTaaaa           <  1:785933/141‑1 (MQ=255)

                                                                                                                                   gATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTaaaa           <  1:1580777/141‑1 (MQ=255)

                                                                                                                                       gCAAATATTATGTTGGTAAAGGTAAAATTAAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAatta       <  1:1374128/141‑1 (MQ=255)

                                                                                                                                           atatTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAATTATTGa   <  1:338991/141‑1 (MQ=255)

                                                                                                                                            tattatGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACACATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAATTATTGAt  >  2:2673953/1‑141 (MQ=255)

TCTTAGTCGGTGTACATGCTCAAGATGATAAGCAATTTAATTTTGAGTCTACAATGGAGGAATTATCATCTTTATCAGAGACTTGCCAACTTGAAGTGTTGGGTCAAATTACTCAAAACAGAGATCGTGTAGATCGCAAATATTATGTTGGTAAAGGTAAAATTGAAGAAATTCAAGCATTTATTGAGTTCA‑AAGATATTGATGTAGTCATCACAAATGATGAATTAACGACTGCACAATCCAAATCACTAAATGAAGCTTTAGGTGTAAAAATTATTGAT  >  USA300TCH1516_ALE/1333649‑1333929

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