BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,219,574	T→A	7.5%	Q299L (CAG→CTG)	USA300HOU_RS06070 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,219,574	0	T	A	7.5%	655.9 / 28.8	197	Q299L (CAG→CTG)	USA300HOU_RS06070	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (94/88); new base A (5/10); total (99/98)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

TCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTA  >  USA300TCH1516_ALE/1219484‑1219664

tCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAActgct                                                                                            >  2:15935/1‑91 (MQ=255)

  aCCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGt                                                                                          <  2:3351483/91‑1 (MQ=255)

   ccATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGtt                                                                                         >  1:814309/1‑91 (MQ=255)

   ccATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGtt                                                                                         <  1:2368801/91‑1 (MQ=255)

    cATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTg                                                                                        >  1:910492/1‑91 (MQ=255)

     atatAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGa                                                                                       <  1:40811/91‑1 (MQ=255)

      tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGaa                                                                                      <  2:1247087/91‑1 (MQ=255)

       ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                     <  2:3799783/91‑1 (MQ=255)

       ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                     <  2:1848312/91‑1 (MQ=255)

        taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACa                                                                                    <  1:3484712/91‑1 (MQ=255)

        taACTGTTCAGTATCTTTA‑TTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaa                                                                                   <  2:3014097/91‑1 (MQ=37)

          aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaaa                                                                                  <  1:306551/91‑1 (MQ=38)

          aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaaa                                                                                  >  1:2716200/1‑91 (MQ=255)

           cTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACAAAt                                                                                 <  1:773377/91‑1 (MQ=38)

           cTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAAt                                                                                 >  2:2157075/1‑91 (MQ=255)

            tGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATc                                                                                <  1:1134049/91‑1 (MQ=255)

              ttCAGGATCTGTATTTTGAGACTGCTCATATTCTTCAAGCAACTTCGCTAATTTATTTTGATATGTGTACAACGGCTGTTGAACAAATCGa                                                                              <  1:3718486/91‑1 (MQ=255)

               tCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGaa                                                                             <  2:1579223/91‑1 (MQ=255)

                cAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAAc                                                                            >  1:659242/1‑91 (MQ=255)

                  gTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTa                                                                          <  2:2417007/91‑1 (MQ=255)

                  gTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTa                                                                          >  1:2313676/1‑91 (MQ=255)

                   tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                         >  2:2596546/1‑91 (MQ=255)

                   tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                         <  1:2497743/91‑1 (MQ=255)

                        ttATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCAt                                                                    <  1:2577310/91‑1 (MQ=255)

                           tttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCg                                                                 <  2:842500/91‑1 (MQ=255)

                            ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                                <  2:3464622/91‑1 (MQ=255)

                            ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                                >  2:1453356/1‑91 (MQ=255)

                            ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                                >  1:1390064/1‑91 (MQ=255)

                            ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                                >  2:2764319/1‑91 (MQ=255)

                             tttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCa                                                               >  1:3797933/1‑91 (MQ=255)

                             tttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCa                                                               <  1:3466227/91‑1 (MQ=255)

                               tAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACg                                                             <  2:1661031/91‑1 (MQ=255)

                               tAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACg                                                             <  1:3213334/91‑1 (MQ=255)

                               tAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACg                                                             <  2:2313676/91‑1 (MQ=255)

                                aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                            >  2:1039290/1‑91 (MQ=255)

                                aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                            >  1:522693/1‑91 (MQ=255)

                                 gACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgtt                                                           <  1:1159562/91‑1 (MQ=255)

                                    tGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgttgtt                                                        <  1:3085844/91‑1 (MQ=255)

                                     gCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTTAACAAATCGAACTAAATCATTCGCACGTTGttt                                                       >  2:1338072/1‑91 (MQ=255)

                                     gCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGttt                                                       <  1:2882449/91‑1 (MQ=255)

                                      cTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTa                                                      >  2:3277879/1‑91 (MQ=255)

                                       tCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTaa                                                     <  2:575944/91‑1 (MQ=255)

                                          tatTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgc                                                  >  1:1877268/1‑91 (MQ=255)

                                           atTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgcg                                                 >  1:22391/1‑91 (MQ=255)

                                           atTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgcg                                                 >  1:3858709/1‑91 (MQ=255)

                                             tcttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGtt                                               >  1:447437/1‑91 (MQ=255)

                                                  aatcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGt                                          >  2:137715/1‑91 (MQ=255)

                                                  aatcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTCTTTAACGCGTTCACGt                                          <  1:3570325/91‑1 (MQ=255)

                                                   atcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGtt                                         <  2:1009140/91‑1 (MQ=255)

                                                    tcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACCAAATCATTCGCACGTTGTTTAACGCGTTCACGTTc                                        >  2:213513/1‑91 (MQ=255)

                                                       aCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGcc                                     <  1:1330553/91‑1 (MQ=255)

                                                        cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGccc                                    >  1:875797/1‑90 (MQ=255)

                                                        cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                    <  2:897581/91‑1 (MQ=255)

                                                        cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                    <  2:369259/91‑1 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCCcgc                                 >  1:1870785/1‑91 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                 >  2:1098357/1‑91 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                 >  1:1524520/1‑91 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                 <  1:2139045/91‑1 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                 >  1:2595480/1‑91 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                 <  1:2570214/91‑1 (MQ=255)

                                                            gCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgcg                                >  2:448965/1‑91 (MQ=255)

                                                             cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                               >  1:1037236/1‑91 (MQ=255)

                                                             cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                               >  2:3558455/1‑91 (MQ=255)

                                                              tAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGca                              >  2:3622855/1‑91 (MQ=255)

                                                               aaTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcat                             >  2:2051834/1‑91 (MQ=255)

                                                               aaTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcat                             <  1:2090486/91‑1 (MQ=255)

                                                                atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                            >  2:1430066/1‑91 (MQ=255)

                                                                atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                            >  1:3145926/1‑91 (MQ=255)

                                                                 tttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatca                           <  2:2509596/91‑1 (MQ=255)

                                                                  ttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatcat                          >  2:150657/1‑91 (MQ=255)

                                                                   tatttTTATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATa                         <  1:3501417/91‑1 (MQ=255)

                                                                   tatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATa                         <  2:67768/91‑1 (MQ=255)

                                                                     ttttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaa                       <  1:1614319/91‑1 (MQ=255)

                                                                      tttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaa                      >  2:402569/1‑91 (MQ=255)

                                                                      tttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaa                      >  2:500635/1‑91 (MQ=255)

                                                                      tttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaa                      <  2:2966721/91‑1 (MQ=255)

                                                                       ttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaaa                     >  2:3746888/1‑91 (MQ=255)

                                                                       ttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaaa                     <  2:2808435/91‑1 (MQ=255)

                                                                         gATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCCTCATAAAAACg                   >  2:172881/1‑91 (MQ=255)

                                                                         gATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACg                   <  1:138019/91‑1 (MQ=255)

                                                                             tttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCa               <  2:3170058/91‑1 (MQ=255)

                                                                              ttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCaa              >  2:304365/1‑91 (MQ=255)

                                                                              ttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCaa              >  1:3336368/1‑91 (MQ=255)

                                                                               tgttcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAg             >  2:1597116/5‑91 (MQ=255)

                                                                               tGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGAATCATAAAAACGATCAAg             <  2:2939998/91‑1 (MQ=255)

                                                                                gTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGc            >  1:2809030/1‑91 (MQ=255)

                                                                                 tACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCa           <  1:1281929/91‑1 (MQ=255)

                                                                                  aCAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaa          >  1:1250374/1‑91 (MQ=255)

                                                                                  aCAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaa          <  1:2074435/91‑1 (MQ=255)

                                                                                   cAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaaa         <  2:3590057/91‑1 (MQ=255)

                                                                                   cAACTGCTGTTAAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaaa         <  2:3536587/91‑1 (MQ=255)

                                                                                     aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc       >  2:1743893/1‑91 (MQ=255)

                                                                                     aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc       >  1:2642532/1‑91 (MQ=255)

                                                                                      ctgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCa      >  1:3380084/1‑91 (MQ=255)

                                                                                      ctgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCa      >  1:2405420/1‑91 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  >  1:2265156/1‑91 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  <  2:2075995/91‑1 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  <  2:747453/91‑1 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  <  2:268827/91‑1 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  <  1:1502007/91‑1 (MQ=255)

TCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTA  >  USA300TCH1516_ALE/1219484‑1219664

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