BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	79,626	A→T	6.9%	intergenic (+235/‑25)	USA300HOU_RS00355 → / → cmoB	hypothetical protein/tRNA U34 carboxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	79,626	0	A	T	6.9%	446.7 / 16.1	130	intergenic (+235/‑25)	USA300HOU_RS00355/cmoB	hypothetical protein/tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand): ref base A (53/68); new base T (7/2); total (60/70)
Fisher's exact test for biased strand distribution p-value = 7.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 134 total aligned reads displayed.

GCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACGAG  >  USA300TCH1516_ALE/79536‑79714

gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                                          >  1:3263004/1‑91 (MQ=255)

gCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGa                                                                                          >  1:3662658/1‑91 (MQ=255)

 cTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGAt                                                                                         <  1:2585559/91‑1 (MQ=255)

     aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                                     >  1:2986162/1‑91 (MQ=255)

     aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                                     <  1:1096923/91‑1 (MQ=255)

     aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                                     >  1:1329708/1‑91 (MQ=255)

     aTTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGa                                                                                     >  1:817842/1‑91 (MQ=255)

       ttttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAAt                                                                                   <  1:2157259/91‑1 (MQ=255)

         ttGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAATCGTAATAATTACGTTTTGAAtag                                                                                 <  1:717849/91‑3 (MQ=37)

           ggTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAg                                                                               >  1:895022/1‑91 (MQ=255)

            gtgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGc                                                                              <  1:929492/91‑1 (MQ=255)

            gtgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGc                                                                              <  1:928613/91‑1 (MQ=255)

             tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                             >  1:1733367/1‑91 (MQ=255)

             tgtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCg                                                                             <  1:1017342/91‑1 (MQ=255)

              gtTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCga                                                                            >  1:343625/1‑91 (MQ=255)

               ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                           >  1:2224117/1‑91 (MQ=255)

               ttcattcaTGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgag                                                                           <  1:3898760/91‑1 (MQ=255)

                tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                                          <  1:1947033/91‑1 (MQ=255)

                tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                                          <  1:1524226/91‑1 (MQ=255)

                tcattcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgagg                                                                          <  1:872869/91‑1 (MQ=255)

                  attcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                                        <  1:2947518/91‑1 (MQ=255)

                  attcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCgaggag                                                                        <  1:2361620/91‑1 (MQ=255)

                   ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                                       >  1:2096469/1‑91 (MQ=255)

                   ttcatGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGgaga                                                                       >  1:2379939/1‑91 (MQ=255)

                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                  >  1:1939048/1‑91 (MQ=255)

                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                  <  1:3945075/91‑1 (MQ=255)

                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                  <  1:2253940/91‑1 (MQ=255)

                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                  >  1:1074656/1‑91 (MQ=255)

                        gACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaa                                                                  >  1:3941401/1‑91 (MQ=255)

                          cTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaat                                                                <  1:2921485/91‑1 (MQ=255)

                           tttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatg                                                               >  1:2306472/1‑91 (MQ=255)

                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                              <  1:2124464/91‑1 (MQ=255)

                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                              <  1:1017839/91‑1 (MQ=255)

                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                              <  1:1483589/91‑1 (MQ=255)

                            ttAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGaatgaaatga                                                              <  1:1483299/91‑1 (MQ=255)

                              aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACt                                                            <  1:382668/91‑1 (MQ=255)

                              aaaaaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAACGAAATGACt                                                            >  1:3720470/1‑91 (MQ=255)

                                 aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                         <  1:1102896/91‑1 (MQ=255)

                                 aaaTCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAg                                                         >  1:3149392/1‑91 (MQ=255)

                                    tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                      <  1:4083511/91‑1 (MQ=255)

                                    tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                      >  1:1346543/1‑91 (MQ=255)

                                    tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                      >  1:4091975/1‑91 (MQ=255)

                                    tCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACt                                                      >  1:505083/1‑91 (MQ=255)

                                       ttATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTAt                                                   <  1:1650748/91‑1 (MQ=255)

                                          ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                                <  1:1290270/91‑1 (MQ=255)

                                          ttGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGAc                                                <  1:3236532/91‑1 (MQ=255)

                                           tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                               <  1:678248/91‑1 (MQ=255)

                                           tGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACg                                               <  1:1134531/91‑1 (MQ=255)

                                             cATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGtt                                             >  1:1062181/1‑91 (MQ=255)

                                              aTAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGttt                                            >  1:2119731/1‑91 (MQ=255)

                                                aaaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtgg                                          <  1:234515/91‑1 (MQ=255)

                                                aaaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtgg                                          >  1:3746356/1‑91 (MQ=255)

                                                 aaTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggt                                         <  1:2692871/91‑1 (MQ=255)

                                                  aTGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTtggtg                                        >  1:3704577/1‑91 (MQ=255)

                                                      aCATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCaa                                    <  1:319387/91‑1 (MQ=255)

                                                       cATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGCTTGGTGGCaaa                                   >  1:3763484/1‑91 (MQ=255)

                                                           gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                               <  1:1742120/91‑1 (MQ=255)

                                                           gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                               <  1:1862287/91‑1 (MQ=255)

                                                           gtgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAgg                               <  1:69524/91‑1 (MQ=255)

                                                            tgtAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGt                              >  1:1720107/1‑91 (MQ=255)

                                                             gtagtaCTATTCATAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                             <  1:2500634/91‑1 (MQ=255)

                                                             gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                             <  1:613265/91‑1 (MQ=255)

                                                             gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                             <  1:2662731/91‑1 (MQ=255)

                                                             gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                             >  1:1440520/1‑91 (MQ=255)

                                                             gtagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTc                             <  1:2500211/91‑1 (MQ=255)

                                                              tagtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCa                            <  1:1261620/91‑1 (MQ=255)

                                                               agtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCaa                           <  1:386136/91‑1 (MQ=255)

                                                                gtaCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAg                          >  1:3568421/1‑91 (MQ=255)

                                                                 taCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGa                         >  1:2978002/1‑91 (MQ=255)

                                                                  aCTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGaa                        <  1:3916320/91‑1 (MQ=255)

                                                                   cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                       >  1:3756864/1‑91 (MQ=255)

                                                                   cTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAAt                       >  1:3757110/1‑91 (MQ=255)

                                                                    tATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATc                      >  1:3265093/1‑91 (MQ=255)

                                                                     aTTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCa                     >  1:728414/1‑91 (MQ=255)

                                                                      ttCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAg                    >  1:1364233/1‑91 (MQ=255)

                                                                      ttCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAg                    >  1:442226/1‑91 (MQ=255)

                                                                       tCAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                   >  1:506817/1‑91 (MQ=38)

                                                                       tCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGCCTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAga                   <  1:3811039/91‑1 (MQ=255)

                                                                        cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat                  >  1:2771319/1‑91 (MQ=39)

                                                                        cAAATCGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat                  >  1:2771320/1‑91 (MQ=39)

                                                                        cAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgat                  <  1:3637718/91‑1 (MQ=255)

                                                                         aaatcGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatg                 >  1:2606121/5‑91 (MQ=39)

                                                                          aatcGTAATAATTACGTTTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga                >  1:1151791/4‑91 (MQ=39)

                                                                          aaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga                >  1:2248792/1‑91 (MQ=255)

                                                                          aaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga                >  1:1880673/1‑91 (MQ=255)

                                                                          aaaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatga                <  1:2939930/91‑1 (MQ=255)

                                                                           aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat               <  1:1268415/91‑1 (MQ=255)

                                                                           aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat               <  1:73481/91‑1 (MQ=255)

                                                                           aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat               <  1:3426163/91‑1 (MQ=255)

                                                                           aaCGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgat               <  1:3033594/91‑1 (MQ=255)

                                                                             cGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatga             <  1:530280/91‑1 (MQ=255)

                                                                              gTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatgat            >  1:558996/1‑91 (MQ=255)

                                                                              gTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAgatgatgat            <  1:1800398/91‑1 (MQ=255)

                                                                               tAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGACAAAAAGGTCAAGAATCAGATGATGata           >  1:2269824/1‑91 (MQ=255)

                                                                                aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat          <  1:2466703/91‑1 (MQ=255)

                                                                                aaTTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGatat          <  1:2191625/91‑1 (MQ=255)

                                                                                 attattACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATg         >  1:3113350/1‑91 (MQ=255)

                                                                                    attaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCa      <  1:1223160/91‑1 (MQ=255)

                                                                                     ttaCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCAc     <  1:3658149/91‑1 (MQ=255)

                                                                                       aCGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACga   >  1:67592/1‑91 (MQ=255)

GCTTAATTTTTGGTGTTCATTCATGACTTTAAAAAATCCTTATTGCATAAATGTACATAGTGTAGTACTATTCAAAACGTAATTATTACGATTTGAATTAAGCGAGGAGAATGAAATGACTAAGACTTATGACGTTTGGTGGCAAAAAGGTCAAGAATCAGATGATGATATGGCACGAG  >  USA300TCH1516_ALE/79536‑79714

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