BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	267,908	T→C	6.9%	Y537Y (TAT→TAC)	USA300HOU_RS01175 →	Staphylocoagulase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	267,908	0	T	C	6.9%	376.4 / 4.6	143	Y537Y (TAT→TAC)	USA300HOU_RS01175	Staphylocoagulase
Reads supporting (aligned to +/- strand): ref base T (62/66); new base C (9/2); total (73/70)
Fisher's exact test for biased strand distribution p-value = 5.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

TTTGGATATGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACAAGAA  >  USA300TCH1516_ALE/267774‑268048

tttGGATATGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATCTAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTa                                                                                                                                        >  1:1061860/1‑141 (MQ=255)

  tGGATATGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTaac                                                                                                                                      <  2:878253/141‑1 (MQ=255)

    gATATGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGCAAAATGGTCAAGTAACATACGGAGCTCGCCCAACACACAACAAGCCAAGC‑AAAACAAACGCATACAACGTaaaa                                                                                                                                     >  4:239065/1‑139 (MQ=21)

     atatGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTaacaac                                                                                                                                   >  3:126291/1‑141 (MQ=255)

     atatGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTaacaac                                                                                                                                   <  1:409405/141‑1 (MQ=255)

        tGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAAcaca                                                                                                                                >  1:646175/1‑141 (MQ=255)

              gagaCCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGaaa                                                                                                                          <  1:712968/141‑1 (MQ=255)

              gagaCCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGaaa                                                                                                                          >  3:5277/1‑141 (MQ=255)

               agaCCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAAc                                                                                                                         >  3:144486/1‑141 (MQ=255)

                aaCCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACg                                                                                                                        <  2:735063/140‑1 (MQ=255)

                 aCCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACgg                                                                                                                       <  1:777950/141‑1 (MQ=255)

                   cAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGcc                                                                                                                     >  2:49430/1‑141 (MQ=255)

                   cAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGcc                                                                                                                     >  1:170827/1‑141 (MQ=255)

                     aGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCaa                                                                                                                   >  1:184060/1‑141 (MQ=255)

                      gATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAg                                                                                                                  <  3:329342/141‑1 (MQ=255)

                       aTTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGt                                                                                                                 >  1:359352/1‑141 (MQ=255)

                        ttCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACATACAAGAAGCCAAGC‑GAAACGAATGCATACAATGTAACAACACATGCAAACGGCCAAGTa                                                                                                                >  1:622156/1‑141 (MQ=14)

                        ttCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTa                                                                                                                <  1:865302/141‑1 (MQ=255)

                          cAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATACAACGTAACAACACATGGAAACGGCCAAGAATc                                                                                                              >  2:422631/1‑141 (MQ=255)

                            aTAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCat                                                                                                            <  1:46653/141‑1 (MQ=255)

                            aTAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCat                                                                                                            >  2:547052/1‑141 (MQ=255)

                            aTAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATACAACGTAACAACACATGCAAACGGCCAAGTATCaa                                                                                                            >  4:84400/1‑140 (MQ=255)

                            aTAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTAAAGTATCATACGGAGCACGTCCGACACAAAACAAGCCAAGC‑AAAACAAACACATAAAACGAAACAACACATGGAAAAGGAAAAGTATCaa                                                                                                            >  2:438658/1‑140 (MQ=21)

                             tAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGTTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCata                                                                                                           <  4:274857/141‑1 (MQ=255)

                                 ccATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGCAAACGGCCAAGTATCATATggg                                                                                                       >  4:307318/1‑140 (MQ=255)

                                  cATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGgcg                                                                                                      <  1:630945/141‑1 (MQ=255)

                                   aTCAGAAACAAATGCATATAACGTAACTACTTATTC‑ATGTTGCGAAGAACCATGCGGTCCTCGCCCGACGCAAAACAAGCCCAGT‑ACACGAAACTCATACAACGTAACTACCCCCAAAGACGGCCAATTATCATATGgctc                                                                                                     <  4:342298/141‑3 (MQ=255)

                                     cAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGCCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATACAACGTAACAACACATGCAAACGGCCAAGTATCATATGGCCCTc                                                                                                   >  2:530897/1‑141 (MQ=255)

                                      aGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATACAACGTAACAACACATGGAAACGGCCAAGTAACATATAGCGCACg                                                                                                  >  4:14700/1‑141 (MQ=255)

                                      aGAAACAAATGCATAGAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCg                                                                                                  <  1:557424/141‑1 (MQ=255)

                                            aaaTGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAac                                                                                            <  3:202379/141‑1 (MQ=255)

                                             aaggCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAaca                                                                                           <  2:258093/138‑1 (MQ=255)

                                                gCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACaa                                                                                        >  3:97208/1‑141 (MQ=255)

                                                 cATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTAACATATGGCGCTCGCCCAACACaaa                                                                                       >  2:427951/1‑141 (MQ=255)

                                                 cATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATACAACGTAACAACACATGGAAACGGCCAAGTAACATATGGCGCTCGCCCAACACaaa                                                                                       >  2:967128/1‑141 (MQ=255)

                                                     taACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAAcaa                                                                                   >  1:410174/1‑141 (MQ=255)

                                                       aCGTAACAACACATGC‑AAATGGTCAAGTATCATACGGTGCTCGCCCAACACAAAAAAAGCCAAGC‑AAAACAAATGCATACAACGTAACAACACATGCAAATGGTCAAGTATCATATGGCGCTCGCCCGACACAAAAaaagc                                                                                 <  3:277876/141‑1 (MQ=11)

                                                       aCGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAAcaagc                                                                                 >  1:467551/1‑141 (MQ=255)

                                                       aCGTAACAACACATGC‑AAATGATTAAGTATCATACGGTGCTAGCATAAAAAAAAAAAAGTAAAGC‑GAAATAAATGAATAAAACGTAAAAATACATGAAAATGGTAAGGTGTCATATGGCGATAGAAAAACAAAAAAaaagc                                                                                 >  4:106448/1‑141 (MQ=255)

                                                           aacaacACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAAcaagccaag                                                                             <  3:40253/141‑1 (MQ=255)

                                                           aacaacACAAGC‑AAATGGTCAAGTATAATACGAAGCTCGACCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGACCAATTATCATATGTCACCCGCCCAACAAAAAAccagtcaaa                                                                             >  2:634518/1‑140 (MQ=21)

                                                            acgacacATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGACCAAGTATCATATGGCGCTCGCCCAACACAAAAcaagccaagc                                                                            <  4:177808/138‑1 (MQ=255)

                                                              aacaCATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCaa                                                                          <  3:76772/141‑1 (MQ=255)

                                                              aacaCATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCaa                                                                          <  3:117106/141‑1 (MQ=255)

                                                                cacaTGC‑AAATGGTCAAGTATCATACGGAGCTCGCCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAATGTAACAACACATGGAAACGGCCAAGTGACATATGGCGCTCGACCAACACAAAACAAGCCAAACtaaa                                                                        >  1:959204/1‑141 (MQ=21)

                                                                  caTGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACa                                                                      >  1:533247/1‑141 (MQ=255)

                                                                    tGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACaaa                                                                    <  1:804236/141‑1 (MQ=255)

                                                                     gC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAAt                                                                   <  2:825451/141‑1 (MQ=255)

                                                                        aaaTGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGc                                                                 >  4:274627/1‑141 (MQ=255)

                                                                         aaTGGTCAAGTATCATACGGAGCCCGACCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAAAGGAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAATCAAAACAAATGCa                                                                <  2:112166/141‑1 (MQ=25)

                                                                           tGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGGTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATa                                                              >  3:386996/1‑141 (MQ=255)

                                                                           tGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATa                                                              <  3:127792/141‑1 (MQ=255)

                                                                           tGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATa                                                              >  2:787266/1‑141 (MQ=255)

                                                                           tGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGACGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATa                                                              >  4:292618/1‑141 (MQ=255)

                                                                             gTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACaaa                                                                                       >  2:366538/1‑114 (MQ=255)

                                                                                  gTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTa                                                       <  2:1003896/141‑1 (MQ=255)

                                                                                       aTACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTaacaac                                                  >  1:206667/1‑141 (MQ=255)

                                                                                         aCGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAacac                                                <  1:1039778/141‑1 (MQ=255)

                                                                                            gAGCCCTTCGGACACAAAACAAGCCAAAC‑AAAACAAACGCATAAAGCGTAACAACACAAGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATg                                             <  2:999769/141‑1 (MQ=21)

                                                                                              gCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCa                                           >  2:234416/1‑141 (MQ=255)

                                                                                              gCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCa                                           >  1:449896/1‑141 (MQ=255)

                                                                                                tCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCaaa                                         >  2:358754/1‑141 (MQ=255)

                                                                                                 cGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAAc                                        <  1:287918/141‑1 (MQ=255)

                                                                                                   tacGACACAAAACAAGCCAAGA‑AAAACAAACGAATATAACGTAACAACACATGAAAACGGCCAAGTATCATATGGAGCGCGACCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACgg                                      <  2:541342/139‑1 (MQ=18)

                                                                                                   tCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACgg                                      >  1:634708/1‑141 (MQ=255)

                                                                                                    ccGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGt                                     <  1:372838/141‑1 (MQ=255)

                                                                                                    ccGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGt                                     <  1:458234/141‑1 (MQ=255)

                                                                                                    ccGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGt                                     <  1:148977/141‑1 (MQ=255)

                                                                                                        cacaAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAg                                 <  2:1085168/141‑1 (MQ=255)

                                                                                                         acaAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAgt                                >  2:88582/1‑141 (MQ=255)

                                                                                                          caAAACAAGCCAAGC‑AAAACAAACGCACATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCACGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAgtg                               <  4:143357/141‑1 (MQ=255)

                                                                                                            aaaCAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTc                             <  1:64650/141‑1 (MQ=255)

                                                                                                            aaaCAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTc                             >  1:983784/1‑141 (MQ=255)

                                                                                                                aagccaagc‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGACATAc                         >  2:151464/1‑141 (MQ=255)

                                                                                                                aagccaagc‑AAAACAAACGAATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATAc                         <  4:19192/141‑1 (MQ=255)

                                                                                                                  gccaagc‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATACGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTACAGCGACATACgg                       >  2:225602/1‑141 (MQ=18)

                                                                                                                    caagcaAAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCGCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGa                      <  2:661816/141‑1 (MQ=255)

                                                                                                                    caagc‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAg                     >  3:267446/1‑141 (MQ=255)

                                                                                                                    caagc‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAg                     >  2:462357/1‑141 (MQ=255)

                                                                                                                     aagc‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGc                    >  1:883220/1‑141 (MQ=255)

                                                                                                                       gc‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTc                  >  1:949800/1‑141 (MQ=255)

                                                                                                                        c‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCg                 <  2:853242/141‑1 (MQ=255)

                                                                                                                          aaaaCAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGc                <  1:105798/141‑1 (MQ=255)

                                                                                                                           aaacaaacGCATATAACGTAAAAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGcc               <  4:260864/141‑1 (MQ=255)

                                                                                                                             aaaaacGGATAGGACGTGACGACACATGGAGCCGGCCCAGTAGGATAGGGGGCGGGCCCAACACAAAAGAAGCCACGCAAAACAAAGGCATACAACGGAACAACACATGCAAACGGTCAAGTGCCATACGGGGCGCGCCCg             <  4:55206/139‑1 (MQ=255)

                                                                                                                               aaacGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACATCACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGAc           <  2:935905/141‑1 (MQ=255)

                                                                                                                               aaacGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGAc           >  3:413424/1‑141 (MQ=255)

                                                                                                                               aaacGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGCCATACGGAGCTCGCCCGAc           >  4:368437/1‑141 (MQ=255)

                                                                                                                                aacGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACa          >  1:269816/1‑141 (MQ=255)

                                                                                                                                 acGCATATAACGTAACAACACATGGAAACGGCCACGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACAt         <  1:1004639/141‑1 (MQ=255)

                                                                                                                                 acGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACAt         >  1:517516/1‑141 (MQ=255)

                                                                                                                                 acGCATATAACGTAACAACACATGGAAACGGCCAACTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACAt         <  1:593512/141‑1 (MQ=255)

                                                                                                                                  cGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGACAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATa        <  2:429358/141‑1 (MQ=255)

                                                                                                                                   gCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATAc       >  1:231597/1‑141 (MQ=255)

                                                                                                                                   gCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATAc       <  1:365583/141‑1 (MQ=255)

                                                                                                                                    tatatAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACa      >  1:1031399/2‑141 (MQ=255)

                                                                                                                                    cATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACa      <  1:27739/141‑1 (MQ=255)

                                                                                                                                    cATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACa      <  2:900947/141‑1 (MQ=255)

                                                                                                                                     atatAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACaa     >  1:60962/1‑141 (MQ=255)

                                                                                                                                        taACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACaagaa  <  3:261700/141‑1 (MQ=255)

TTTGGATATGAAGCGAGACCAAGATTCAATAAGCCATCAGAAACAAATGCATATAACGTAACAACACATGC‑AAATGGTCAAGTATCATACGGAGCTCGTCCGACACAAAACAAGCCAAGC‑AAAACAAACGCATATAACGTAACAACACATGGAAACGGCCAAGTATCATATGGCGCTCGCCCAACACAAAACAAGCCAAGCAAAACAAATGCATACAACGTAACAACACATGCAAACGGTCAAGTGTCATACGGAGCTCGCCCGACATACAAGAA  >  USA300TCH1516_ALE/267774‑268048

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