BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,530,791	C→T	5.9%	intergenic (‑42/+231)	yefM ← / ← bdbD	Antitoxin YefM/Disulfide bond formation protein D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,530,791	0	C	T	5.9%	761.1 / 21.8	237	intergenic (‑42/+231)	yefM/bdbD	Antitoxin YefM/Disulfide bond formation protein D
Reads supporting (aligned to +/- strand): ref base C (112/110); new base T (4/10); total (116/121)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 238 total aligned reads displayed.

TAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCCTA  >  USA300TCH1516_ALE/2530701‑2530880

 aaGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCt                                                                                          <  1:2750682/91‑1 (MQ=255)

  aGTGCTTTTAAATTCTGTCGCTCGTTTGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTCCAAAGATCTg                                                                                         >  2:3643415/1‑91 (MQ=255)

  aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTg                                                                                         >  2:2412768/1‑91 (MQ=255)

  aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTg                                                                                         <  1:729795/91‑1 (MQ=255)

  aGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAACATCTg                                                                                         >  2:3644510/1‑91 (MQ=255)

  aGTACCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAGTATCTg                                                                                         <  1:3070425/91‑1 (MQ=255)

   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                        <  2:892968/91‑1 (MQ=255)

   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                        <  2:1870552/91‑1 (MQ=255)

   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                        <  1:3877546/91‑1 (MQ=255)

   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGt                                                                                        <  2:1087995/91‑1 (MQ=255)

   gTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGt                                                                                        <  1:1395799/91‑1 (MQ=255)

    tGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTa                                                                                       <  1:3293825/91‑1 (MQ=255)

     gCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGTAc                                                                                      <  2:2364920/91‑1 (MQ=255)

      ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                     >  1:3768824/1‑91 (MQ=255)

      ccTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACa                                                                                     >  2:1943050/1‑91 (MQ=255)

        tttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAAt                                                                                   >  1:544941/1‑91 (MQ=255)

         ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAGATATTTGTACAAtt                                                                                  <  1:464952/91‑2 (MQ=255)

         ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATa                                                                                  >  1:1450556/1‑91 (MQ=255)

         ttAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATa                                                                                  <  1:2355682/91‑1 (MQ=255)

          tAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATaa                                                                                 <  2:2583349/91‑1 (MQ=255)

             aTTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTCATCGCCTCCCTTTAACACACCTCAATTGTACACATATCTGTACAATAAAAc                                                                              >  1:257298/1‑91 (MQ=255)

                cTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAg                                                                           >  1:2314998/1‑91 (MQ=255)

                  gTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCg                                                                         >  1:2870972/1‑91 (MQ=255)

                   tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                        >  2:2459885/1‑91 (MQ=255)

                   tCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGa                                                                        >  2:1782227/1‑91 (MQ=255)

                    cGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGac                                                                       >  1:408363/1‑91 (MQ=255)

                      ggCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGacac                                                                     <  2:1156292/90‑1 (MQ=255)

                       gCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGAcaca                                                                    <  2:1808692/91‑1 (MQ=255)

                         gTATGAATAATTTTCAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAAAATCTGCACAATAAAACAAGTGACACAGa                                                                  >  2:3355653/1‑91 (MQ=255)

                          tATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaa                                                                 <  2:2664314/91‑1 (MQ=255)

                          tATGAATAATTTTTAATAATCATATTTATCGCCCCCTATTAACACACCTCAATTGTACAAATAGCTGTACAATAAAACAAGCGACACAGaa                                                                 <  1:124504/91‑1 (MQ=255)

                            tGAATAAATTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGgaaa                                                               <  2:3501465/91‑1 (MQ=255)

                             gAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaa                                                              <  1:1017877/91‑1 (MQ=255)

                             gAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaa                                                              <  1:2070746/91‑1 (MQ=255)

                               ataatTTTTAATAATCATATTTATCGCCTCCTTTTAACACAGCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaaaa                                                            >  1:3430132/1‑91 (MQ=255)

                               ataatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGaaaaaaa                                                            >  1:225334/1‑91 (MQ=255)

                                taatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAac                                                           >  1:899745/1‑91 (MQ=255)

                                 aatTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAact                                                          <  2:632540/91‑1 (MQ=255)

                                   tttatAATAAACATATTTATCGCGTCCTTTAATCACAACTCAAGTGTACAGATATCTGGACAATAAAACAAGCGACACAGAAAAAAactac                                                        <  2:2239666/91‑1 (MQ=255)

                                    ttttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAactact                                                       >  1:3173353/1‑91 (MQ=255)

                                      ttAGTAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTaa                                                     <  1:3066567/91‑1 (MQ=255)

                                      ttAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTaa                                                     <  1:2052121/91‑1 (MQ=255)

                                       taataaTCATATTTATCGCCGCCTTTTAACACACCTCAATTGTACAAAGATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAAt                                                    <  2:1565714/91‑1 (MQ=255)

                                        aataatCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATc                                                   <  1:3543244/91‑1 (MQ=255)

                                            atCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGa                                               >  1:2393167/1‑91 (MQ=255)

                                            atCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGa                                               >  2:407818/1‑91 (MQ=255)

                                             tCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGAt                                              >  1:3717026/1‑91 (MQ=255)

                                                 atTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCg                                          <  1:994290/91‑1 (MQ=255)

                                                     aTCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACg                                      <  2:1417082/91‑1 (MQ=255)

                                                     aTCGCCTCCTTTTAACACACCGCAATTGTACAAATATCTGTACAAGAAAACAAGCGACACAGAAAAAAACTGCTAATCAAGATTTCGCACg                                      <  2:3460444/91‑1 (MQ=255)

                                                      tCGCCTCCTTTTAACACTCCTCACTTGTACAAATATCTGTACCATAAAACACGCCACACAGAAAAAAACTACTAATCAAGATTTCacactg                                     >  2:229444/1‑89 (MQ=255)

                                                      tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAATATTTCGCACgg                                     >  1:3587729/1‑91 (MQ=255)

                                                      tCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgg                                     <  2:3311073/91‑1 (MQ=255)

                                                       cGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACggg                                    <  1:2966705/91‑1 (MQ=255)

                                                        gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAATTACTAATCAAGATTTCGCACgggg                                   >  2:3688706/1‑91 (MQ=255)

                                                        gCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACgggg                                   <  1:2459210/91‑1 (MQ=255)

                                                         cctcctTTTAACACACCTCAATTGTACAAATATGTGTACAATAAAACAAGCGACACAGATAAAAACTACTAATCAAGATTTCGCACGGGga                                  <  1:3511798/91‑1 (MQ=255)

                                                         cctcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGga                                  <  2:3332461/91‑1 (MQ=255)

                                                          ctcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaa                                 >  1:1285676/1‑91 (MQ=255)

                                                           tcctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaag                                <  2:3065731/91‑1 (MQ=255)

                                                            cctTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGgaaga                               <  1:2696022/91‑1 (MQ=255)

                                                                 tAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCaa                          >  1:2987352/1‑91 (MQ=255)

                                                                   acacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATa                        <  2:753133/91‑1 (MQ=255)

                                                                    cacacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAc                       <  2:2118291/91‑1 (MQ=255)

                                                                     acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                      <  1:933351/91‑1 (MQ=255)

                                                                     acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                      <  1:2435016/91‑1 (MQ=255)

                                                                     acacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATAcc                      <  1:304023/91‑1 (MQ=255)

                                                                      cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                     >  1:1074021/1‑91 (MQ=255)

                                                                      cacCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCt                     <  1:1679602/91‑1 (MQ=255)

                                                                       acCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTa                    <  2:3842878/91‑1 (MQ=255)

                                                                        ccTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTAt                   >  1:1805618/1‑91 (MQ=255)

                                                                        ccTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTAt                   >  2:2942133/1‑91 (MQ=255)

                                                                         cTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATc                  >  1:2853541/1‑91 (MQ=255)

                                                                            aaTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTg               >  2:271497/1‑91 (MQ=255)

                                                                            aaTTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTg               <  1:2937238/91‑1 (MQ=255)

                                                                             aTTGTACAAATCTCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATCTCGCACGGGGAAGAATCAATACCTATCCTGt              <  2:3024394/91‑1 (MQ=255)

                                                                              ttGTACAGATATTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTa             >  1:388447/1‑91 (MQ=255)

                                                                               tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc            >  2:2903447/1‑91 (MQ=255)

                                                                               tGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc            <  1:3579266/91‑1 (MQ=255)

                                                                               tGTACAAATAGCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTAc            <  2:2443298/91‑1 (MQ=255)

                                                                                gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt           <  1:3141609/91‑1 (MQ=255)

                                                                                gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt           <  2:1050812/91‑1 (MQ=255)

                                                                                gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt           <  2:2073449/91‑1 (MQ=255)

                                                                                gTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACt           >  2:3827345/1‑91 (MQ=255)

                                                                                  acagatatTTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa         >  2:3403627/5‑91 (MQ=255)

                                                                                  aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa         >  2:3235651/1‑91 (MQ=255)

                                                                                  aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa         >  1:2623417/1‑91 (MQ=255)

                                                                                  aCAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGa         >  1:3625915/1‑91 (MQ=255)

                                                                                   cAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGAt        >  2:3647465/1‑91 (MQ=255)

                                                                                    aaaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGAtt       >  1:2599025/1‑91 (MQ=255)

                                                                                     aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc      <  1:2114140/91‑1 (MQ=255)

                                                                                     aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc      >  2:3899704/1‑91 (MQ=255)

                                                                                     aaTATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc      <  2:461052/91‑1 (MQ=255)

                                                                                     aaTAGCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTc      <  2:415839/91‑1 (MQ=255)

                                                                                      atatCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTcc     >  1:1686159/1‑91 (MQ=255)

                                                                                        atCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCct   <  2:2560041/91‑1 (MQ=255)

                                                                                        atCTGTACAATAAAACAAACGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCct   <  2:1021925/91‑1 (MQ=255)

                                                                                         tCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAATATTTCGCACGGGGAAGACTCAATACCGATCCTGTACTGGTTCCctc  >  2:1764248/1‑90 (MQ=255)

                                                                                         tCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCcta  >  1:1159655/1‑91 (MQ=255)

                                                                                         tCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATGTCGCACGGGGACGAATCAATACCTATCCTGTACTGATTCCcta  >  2:536742/1‑91 (MQ=255)

TAAGTGCCTTTAAATTCTGTCGAGCGTATGAATAATTTTTAATAATCATATTTATCGCCTCCTTTTAACACACCTCAATTGTACAAATATCTGTACAATAAAACAAGCGACACAGAAAAAAACTACTAATCAAGATTTCGCACGGGGAAGAATCAATACCTATCCTGTACTGATTCCCTA  >  USA300TCH1516_ALE/2530701‑2530880

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