BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,329	T→A	8.3%	N133K (AAT→AAA)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,329	0	T	A	8.3%	958.0 / 13.1	322	N133K (AAT→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base T (173/113); new base A (2/34); total (175/147)
Fisher's exact test for biased strand distribution p-value = 1.08e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

CAAATGGCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAA  >  USA300TCH1516_ALE/111189‑111469

cAAATGGCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAAt                                                                                                                                               <  1:588995/141‑1 (MQ=255)

cAAATGGCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAAt                                                                                                                                               >  4:412971/1‑141 (MQ=255)

  aaTGGCAATTCAAAATAAAGGGGAAGCGCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAAGAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑a                                                                                                                                            <  1:27365/141‑1 (MQ=255)

   aTGGCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑aa                                                                                                                                           >  2:271096/1‑141 (MQ=255)

     ggCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑aaaa                                                                                                                                         <  1:1111430/141‑1 (MQ=255)

      gCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑aaaaa                                                                                                                                        >  4:497275/1‑141 (MQ=255)

        aaTTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAGTAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑aaaaaaa                                                                                                                                      >  4:470156/1‑141 (MQ=255)

        aaTTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑aaaaaaa                                                                                                                                      >  1:1080185/1‑141 (MQ=255)

            cAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATGGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAatat                                                                                                                                  >  2:59196/1‑141 (MQ=255)

            cAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAatat                                                                                                                                  >  3:355112/1‑141 (MQ=255)

              aaaTAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATcc                                                                                                                                >  1:642699/1‑141 (MQ=255)

              aaaTAAAGGAGAAGCACAAAAAAAAAAAGAAATACATTTGAAGATAGAAAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTAAAAAAAAAA‑AAAAAAAATATcc                                                                                                                                <  2:772292/141‑1 (MQ=255)

                aTAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAg                                                                                                                              <  4:315911/141‑1 (MQ=255)

                 tAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGt                                                                                                                             >  1:485251/1‑141 (MQ=255)

                 tAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAAAA‑AAAAAAAATATCCAGt                                                                                                                             <  2:39676/141‑1 (MQ=255)

                  aaaGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGtt                                                                                                                            >  2:464177/1‑141 (MQ=255)

                   aaGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTa                                                                                                                           <  1:967212/141‑1 (MQ=255)

                     ggAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTaaa                                                                                                                         <  3:49924/141‑1 (MQ=255)

                      gagaAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTaaaa                                                                                                                        >  3:273612/1‑141 (MQ=255)

                      gagaAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTaaaa                                                                                                                        >  1:677475/1‑141 (MQ=255)

                      gagaAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTaaaa                                                                                                                        <  1:172253/141‑1 (MQ=255)

                      gagaAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGAAAAAAA‑AAAAAAAATATCCAGTTaaaa                                                                                                                        <  2:777338/141‑1 (MQ=255)

                      gagaAGCAATAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTAAAAAAAAAA‑AAAAAAAATATCCAGTTaaaa                                                                                                                        <  2:249344/141‑1 (MQ=255)

                      gagaAGCAAAAAAAAAAAAAAGCATGCTTTTGAAAATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGAAAAAAA‑AAAAAAAATATCCAGTTaaaa                                                                                                                        <  2:700388/141‑1 (MQ=255)

                              aaaaaaaTAAAAGGAATGCTTTTGAAAAAAAATAGAAAAAAAAAAAAAGAAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAAAA‑AAAAAAAATATCCAGTTAAAATGATtaat                                                                                                                <  4:383256/139‑1 (MQ=255)

                               tAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTAAGGATAATC‑AAAAAAGATATCCAGTTAAAATGAGtaata                                                                                                               >  1:941270/1‑141 (MQ=255)

                                 aaaaTAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATtaataat                                                                                                             >  3:90088/1‑141 (MQ=255)

                                  aaaTAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATtaataata                                                                                                            <  1:75725/141‑1 (MQ=255)

                                  aaaTAAAAGGAATGATTTTGAAAATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTTAAAAAAAAA‑AAAAAAAATATCCAGTTAAAATGATtaataata                                                                                                            <  4:172028/141‑1 (MQ=255)

                                       aaaGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTc                                                                                                       >  2:559161/1‑141 (MQ=255)

                                          ggCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCAtt                                                                                                    >  3:404954/1‑141 (MQ=255)

                                            cATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTc                                                                                                  >  1:609102/1‑141 (MQ=255)

                                             aTGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGAATAATAATAAATTCATTTca                                                                                                 >  1:205631/1‑140 (MQ=255)

                                              tccTTTTGAAGATAGATAGAAATACAAGACGTGCAAAAGGCTTCTACTAAACGAATGAAATAAGGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTa                                                                                                >  2:1004819/3‑141 (MQ=255)

                                                cTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACt                                                                                              <  1:266461/141‑1 (MQ=255)

                                                  tttGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGa                                                                                            >  1:1046596/1‑141 (MQ=255)

                                                  tttGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAAAAAAAAA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGa                                                                                            <  2:886230/141‑1 (MQ=255)

                                                    tGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAgg                                                                                          >  3:156000/1‑141 (MQ=255)

                                                    tGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAgg                                                                                          <  2:86659/141‑1 (MQ=255)

                                                     gAAGATAGAAAGAAATACAAGAAGAGCAAAAGGAATTTACTATACTAATGAAAAAAAGACGGAGAAATACGAAGTAGCTAAAAAAAAAA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGa                                                                                         <  2:152494/141‑1 (MQ=255)

                                                            gatagaAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTaaa                                                                                  >  2:1095136/1‑141 (MQ=255)

                                                             atagaAATACAAGAAGTGCAAAAGTATTTTACTATACTAATGAAATAAATACTAAGAAATACTAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaag                                                                                 <  1:410773/141‑1 (MQ=255)

                                                             atagaAATACAAGAAGTGCAAAAGTATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGAATAAAC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaag                                                                                 <  1:502125/141‑1 (MQ=255)

                                                             atagaAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaag                                                                                 <  1:1152654/141‑1 (MQ=255)

                                                             atagaAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaag                                                                                 >  4:276481/1‑141 (MQ=255)

                                                              aagaAATACAAGAAGAGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAAAAAAAAA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaaga                                                                                <  2:808439/140‑1 (MQ=255)

                                                               agaAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaagaa                                                                               >  2:533679/1‑141 (MQ=255)

                                                                gaAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaagaag                                                                              <  1:897549/141‑1 (MQ=255)

                                                                gaAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAaagaag                                                                              <  1:1094709/141‑1 (MQ=255)

                                                                  aaTACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACt                                                                                              >  2:177752/1‑123 (MQ=255)

                                                                     aCAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAAc                                                                         <  1:854273/141‑1 (MQ=255)

                                                                     aCAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAACGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAAc                                                                         <  3:417067/141‑1 (MQ=255)

                                                                      cAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACa                                                                        <  1:568810/141‑1 (MQ=255)

                                                                      cAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATACAGTTAAAATGATTAATAATAAATTCATTTCTAATGAGGAAGTTAAAGAaaaaaaaa                                                                        >  2:10338/1‑139 (MQ=255)

                                                                        agaagTGCAAAAGGATTTTACTATACTAATAAAATAAAGACAGAGAAATACGAAGTAGATCAAGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATa                                                                      <  2:1026335/141‑1 (MQ=255)

                                                                          aagTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATaaa                                                                    >  2:638143/1‑141 (MQ=255)

                                                                          aagTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTAAAGATAAAA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATaaa                                                                    <  4:350245/141‑1 (MQ=255)

                                                                               cAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAg                                                               >  1:1131983/1‑141 (MQ=255)

                                                                                 aaaGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGaa                                                             <  4:201110/141‑1 (MQ=255)

                                                                                 aaaGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGaa                                                             >  1:239035/1‑141 (MQ=255)

                                                                                         ttACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGaaaa                                                     >  2:740837/1‑141 (MQ=255)

                                                                                         ttACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGaaaa                                                     >  1:692982/1‑141 (MQ=255)

                                                                                         ttACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAAAATAAAAAAAGAAATCGaaaa                                                     >  3:458230/1‑141 (MQ=255)

                                                                                          tactatactaATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAAc                                                    >  2:507701/1‑141 (MQ=255)

                                                                                          tacgatactaATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGGTAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCAGTTCTACGGAGTAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGaaaag                                                    >  4:122835/1‑140 (MQ=39)

                                                                                           actatactaATGAAATAAAGACGGATAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACt                                                   <  2:1074478/141‑1 (MQ=255)

                                                                                           actatactaATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACt                                                   >  1:403060/1‑141 (MQ=255)

                                                                                            atatacaaATGAAATAAAGACAGAGAAATACGAAATAGATCAGGATAATA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACtt                                                  <  2:93714/140‑1 (MQ=38)

                                                                                             aatactaATGAAATAAAGACGGAGAAATACGAAGTAGCTCAAAAAAAAA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACttt                                                 <  2:729564/140‑1 (MQ=39)

                                                                                              atactaATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTa                                                <  4:302326/141‑1 (MQ=255)

                                                                                                    aTGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGtttt                                          >  3:215438/1‑141 (MQ=255)

                                                                                                      gAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGtttttt                                        >  1:722445/1‑141 (MQ=255)

                                                                                                      gAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑‑AAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGttttttt                                       >  1:317224/1‑140 (MQ=255)

                                                                                                       aaaTAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTg                                       <  2:79904/141‑1 (MQ=255)

                                                                                                        aaaaaaGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTgc                                      <  3:338351/138‑1 (MQ=255)

                                                                                                        aaTAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTgc                                      >  1:286086/1‑141 (MQ=255)

                                                                                                        aaTAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTgc                                      >  3:302319/1‑141 (MQ=255)

                                                                                                           aaaGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCa                                   >  4:240446/1‑141 (MQ=255)

                                                                                                            aaGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCaa                                  <  2:334501/141‑1 (MQ=255)

                                                                                                             aGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAat                                 >  1:14660/1‑141 (MQ=255)

                                                                                                               aCGGAGAAATACGAAGTAGCTCAGGATAATCAAAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAata                                <  1:73610/141‑1 (MQ=255)

                                                                                                                 ggAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAg                             <  1:1016044/141‑1 (MQ=255)

                                                                                                                   agaAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCa                           >  4:78008/1‑141 (MQ=255)

                                                                                                                         aCGAAGTAGCTCAGGATAATC‑AAAAAAAATATTCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACGTAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAAttttg                     <  3:385544/141‑2 (MQ=255)

                                                                                                                         aCGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTa                     <  1:1004265/141‑1 (MQ=255)

                                                                                                                         aCGAAGTAGCTCAGGATAAAA‑AAAAAAAATATCCAGTTAAAATGATTAAGAAGAAATTCGTTTCTCGTGAGGAAGTTAGGGGGGAAAACATAAAAATAGAAGTCGAAAACTTTAAGTGTTTTGCGCACTAGAGCAATTTTa                     <  4:425817/141‑1 (MQ=255)

                                                                                                                          cGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTTTGCAATATAGCAATTTTaa                    >  2:993979/1‑141 (MQ=255)

                                                                                                                          cGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTaa                    >  2:956170/1‑141 (MQ=255)

                                                                                                                          cGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTaa                    >  2:1148934/1‑141 (MQ=255)

                                                                                                                          cGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTaa                    >  4:332044/1‑141 (MQ=255)

                                                                                                                              gTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGAt                >  1:417115/1‑141 (MQ=255)

                                                                                                                              gTAGCAAAGGATAAAA‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGAt                <  1:429703/141‑1 (MQ=255)

                                                                                                                                aGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGAttt              <  1:582689/141‑1 (MQ=255)

                                                                                                                                 gCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTa             >  1:778354/1‑141 (MQ=255)

                                                                                                                                 gCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTa             >  2:201001/1‑141 (MQ=255)

                                                                                                                                    cAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATg          <  2:1170720/141‑1 (MQ=255)

                                                                                                                                    cAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATg          >  1:382080/1‑141 (MQ=255)

                                                                                                                                     aGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGa         >  2:402660/1‑141 (MQ=255)

                                                                                                                                          aaTC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATtat    <  2:135112/141‑1 (MQ=255)

                                                                                                                                          aaTC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATtat    <  3:118052/141‑1 (MQ=255)

CAAATGGCAATTCAAAATAAAGGAGAAGCTCTAAAAATAAAAGGCATGCTTTTGAAGATAGATAGAAATACAAGAAGTGCAAAAGGATTTTACTATACTAATGAAATAAAGACGGAGAAATACGAAGTAGCTCAGGATAATC‑AAAAAAAATATCCAGTTAAAATGATTAATAATAAATTCATTTCTACTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAA  >  USA300TCH1516_ALE/111189‑111469

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: `‑`