BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,313,113	A→T	9.6%	intergenic (+17/+280)	rny_1 → / ← USA300HOU_RS06485	Ribonuclease Y/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,313,113	0	A	T	9.6%	414.6 / 26.9	126	intergenic (+17/+280)	rny_1/USA300HOU_RS06485	Ribonuclease Y/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (54/59); new base T (3/9); total (58/68)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 131 total aligned reads displayed.

ATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTG  >  USA300TCH1516_ALE/1313023‑1313202

aTGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTcaca                                                                                           <  1:3807616/91‑1 (MQ=255)

  gAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACaaa                                                                                         >  1:798691/1‑91 (MQ=255)

   aaTTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAt                                                                                        <  1:3875917/91‑1 (MQ=255)

   aaTTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAt                                                                                        <  1:2695101/91‑1 (MQ=255)

     ttACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTa                                                                                      <  1:3192168/91‑1 (MQ=255)

     ttACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTa                                                                                      <  1:212539/91‑1 (MQ=255)

     ttACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTa                                                                                      >  1:3826029/1‑91 (MQ=255)

      tACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtttg                                                                                     <  1:2421318/91‑3 (MQ=255)

      tACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAg                                                                                     <  1:185188/91‑1 (MQ=255)

      tACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAg                                                                                     <  1:182824/91‑1 (MQ=255)

          atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                 <  1:2048087/91‑1 (MQ=255)

          atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                 >  1:1776098/1‑91 (MQ=255)

          atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                 >  1:3841279/1‑91 (MQ=255)

          atatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAg                                                                                 >  1:3839831/1‑91 (MQ=255)

           tatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAgg                                                                                <  1:1548120/91‑1 (MQ=255)

           tatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAgg                                                                                >  1:1444372/1‑91 (MQ=255)

             tCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGa                                                                              >  1:3882617/1‑91 (MQ=255)

              ccTGGTCATATCAAGGTCACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAg                                                                             <  1:3551049/91‑1 (MQ=255)

               cTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGc                                                                            <  1:1621930/91‑1 (MQ=255)

               cTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGc                                                                            <  1:179279/91‑1 (MQ=255)

                tGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCt                                                                           >  1:2910140/1‑91 (MQ=255)

                 ggTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAATTTGTGAGGGAgaca                                                                          <  1:3951028/91‑4 (MQ=39)

                  gTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCttt                                                                         >  1:2017538/1‑91 (MQ=255)

                    cATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCttttt                                                                       <  1:3777419/91‑1 (MQ=255)

                     atatCAAGGTGACAGTTGTCCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtttttt                                                                      >  1:519586/1‑91 (MQ=255)

                      tatgaaGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                     <  1:3734007/87‑1 (MQ=255)

                      tatCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                     <  1:266885/91‑1 (MQ=255)

                      tatCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                     >  1:2055401/1‑91 (MQ=255)

                      tatCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTa                                                                     <  1:1446104/91‑1 (MQ=255)

                        tCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAg                                                                   >  1:3772672/1‑91 (MQ=255)

                        tCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAg                                                                   <  1:3870980/91‑1 (MQ=255)

                         cAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAAAAATTTTTGTCTCCCTCACCAATTAGTGAGGGAGCGTTTTTAAGt                                                                  >  1:894728/1‑91 (MQ=255)

                           aGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTg                                                                <  1:1899412/91‑1 (MQ=255)

                           aGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTg                                                                <  1:4077957/91‑1 (MQ=255)

                            ggTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGt                                                               >  1:1810944/1‑91 (MQ=255)

                            ggTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGt                                                               >  1:3893292/1‑91 (MQ=255)

                             gTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTa                                                              <  1:3998692/91‑1 (MQ=255)

                              tGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAg                                                             >  1:2589468/1‑91 (MQ=255)

                              tGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAg                                                             >  1:2589789/1‑91 (MQ=255)

                               gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                            <  1:1158724/91‑1 (MQ=255)

                               gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                            <  1:1159103/91‑1 (MQ=255)

                               gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                            <  1:1159440/91‑1 (MQ=255)

                               gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                            <  1:2748008/91‑1 (MQ=255)

                                aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                           <  1:1901417/91‑1 (MQ=255)

                                aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                           >  1:1697130/1‑91 (MQ=255)

                                 cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                          >  1:2694803/1‑91 (MQ=255)

                                 cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                          >  1:4132924/1‑91 (MQ=255)

                                 cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                          <  1:1494068/91‑1 (MQ=255)

                                   gttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTa                                                        >  1:2262429/1‑91 (MQ=255)

                                    ttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTaa                                                       >  1:1700107/1‑91 (MQ=255)

                                    ttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTaa                                                       <  1:3349987/91‑1 (MQ=255)

                                     tgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAAt                                                      <  1:1704807/91‑1 (MQ=255)

                                      gttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATc                                                     >  1:2849231/1‑91 (MQ=255)

                                          gagagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGt                                                 <  1:3815003/91‑1 (MQ=255)

                                           agaggCAAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGtt                                                <  1:1118116/91‑1 (MQ=255)

                                           agagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGtt                                                >  1:1092491/1‑91 (MQ=255)

                                            gagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTa                                               <  1:1904725/91‑1 (MQ=255)

                                            gagaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTa                                               >  1:1205890/1‑91 (MQ=255)

                                               aCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGAc                                            >  1:2961545/1‑91 (MQ=255)

                                               aCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGAc                                            <  1:2279412/91‑1 (MQ=255)

                                                cTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACa                                           >  1:1029218/1‑91 (MQ=255)

                                                 tAGAGGAGTAGAATATGCGAAATAACTTTTGTCTCCCTCGCAAATTAGTGAGGCAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAg                                          <  1:3861137/91‑1 (MQ=255)

                                                 tAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAg                                          <  1:2667775/91‑1 (MQ=255)

                                                  agagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGc                                         <  1:1611415/91‑1 (MQ=255)

                                                  agagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGc                                         >  1:3683727/1‑91 (MQ=255)

                                                   gagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCa                                        >  1:3020246/1‑91 (MQ=255)

                                                    agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                       >  1:2778924/1‑91 (MQ=255)

                                                    agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                       <  1:30869/91‑1 (MQ=255)

                                                        gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTa                                   <  1:1037916/91‑1 (MQ=255)

                                                        gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTa                                   <  1:1041038/91‑1 (MQ=255)

                                                         tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                  <  1:812345/91‑1 (MQ=255)

                                                         tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                  <  1:1773338/91‑1 (MQ=255)

                                                         tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                  <  1:2831519/91‑1 (MQ=255)

                                                         tAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTAt                                  <  1:812199/91‑1 (MQ=255)

                                                               atGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaa                            >  1:4047042/1‑91 (MQ=255)

                                                               atGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAAGGAGCTTTTTTAAGTTGTAGTCTTAATTTAGTTAGACAGCACTTTATAGGtaa                            >  1:3865950/1‑91 (MQ=255)

                                                                tGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaat                           <  1:3472032/91‑1 (MQ=255)

                                                                 gCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaata                          <  1:1757036/91‑1 (MQ=255)

                                                                  cGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGtaataa                         >  1:760198/1‑91 (MQ=255)

                                                                   gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTATTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                        <  1:3497040/91‑1 (MQ=255)

                                                                   gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                        <  1:1604337/91‑1 (MQ=255)

                                                                   gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                        <  1:1317382/91‑1 (MQ=255)

                                                                        aatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTatat                   <  1:3380785/91‑1 (MQ=255)

                                                                        aatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTatat                   <  1:3379494/91‑1 (MQ=255)

                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  1:460041/1‑91 (MQ=255)

                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  1:24584/1‑91 (MQ=255)

                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  1:24361/1‑91 (MQ=255)

                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATAGTaa                >  1:3216588/1‑91 (MQ=255)

                                                                            tttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaaa               >  1:280270/1‑91 (MQ=255)

                                                                                gctcCCTCACTAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATGTTAAACagt           >  1:2641694/2‑91 (MQ=39)

                                                                                 ctcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagta          <  1:3814026/91‑1 (MQ=255)

                                                                                  tcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtag         >  1:4018261/1‑91 (MQ=255)

                                                                                  tcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtag         >  1:109314/1‑91 (MQ=255)

                                                                                   cccTCACAAATTAGTGGGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagt        >  1:1209444/1‑91 (MQ=255)

                                                                                   cccTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagt        >  1:3239860/1‑91 (MQ=255)

                                                                                    ccTCACTAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGtt       >  1:1001339/1‑91 (MQ=255)

                                                                                    ccTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGtt       <  1:3169748/91‑1 (MQ=255)

                                                                                     cTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTCAATCTAATTAGACAGCACTTTATCGGTAATAACTATAGTAAACAGTGGTTa      >  1:3203182/1‑91 (MQ=255)

                                                                                       cactAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTAtt    >  1:760715/1‑91 (MQ=255)

                                                                                         caAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTg  >  1:2889343/1‑91 (MQ=255)

ATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTG  >  USA300TCH1516_ALE/1313023‑1313202

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