BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,794,352	A→T	6.1%	Y426N (TAT→AAT)	USA300HOU_RS08960 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,794,352	0	A	T	6.1%	401.4 / 13.3	116	Y426N (TAT→AAT)	USA300HOU_RS08960	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (52/56); new base T (3/4); total (55/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 125 total aligned reads displayed.

AAACGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCG  >  USA300TCH1516_ALE/1794263‑1794442

aaaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATaa                                                                                           <  1:922836/91‑1 (MQ=255)

aaaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATaa                                                                                           >  1:2173696/1‑91 (MQ=255)

aaaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATaa                                                                                           <  1:590287/91‑1 (MQ=255)

 aaCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGCACTCTTTTTTCATGTTCATAAt                                                                                          >  1:3444654/1‑91 (MQ=255)

  aCGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATg                                                                                         <  1:766331/91‑1 (MQ=255)

   cGGCAGTGATTCCACTGTATAAACCACTTGTCGTGGATGGTCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATTATGa                                                                                        <  1:3858679/91‑1 (MQ=255)

   cGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATTATGa                                                                                        <  1:726910/91‑1 (MQ=255)

    ggCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGaa                                                                                       <  1:1672298/91‑1 (MQ=255)

    ggCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGaa                                                                                       <  1:890819/91‑1 (MQ=255)

      cAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACa                                                                                     >  1:2288395/1‑91 (MQ=255)

      cAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACa                                                                                     >  1:2287455/1‑91 (MQ=255)

      cAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACa                                                                                     <  1:3240670/91‑1 (MQ=255)

        aTGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATTATGAACatg                                                                                   <  1:3254761/91‑2 (MQ=255)

         tGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataa                                                                                  <  1:3121781/91‑1 (MQ=255)

         tGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataa                                                                                  <  1:3121066/91‑1 (MQ=255)

          gATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaa                                                                                 <  1:3027346/91‑1 (MQ=255)

          gATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGGTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaa                                                                                 <  1:3285819/91‑1 (MQ=255)

           aTTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCCTAATGAACataaat                                                                                >  1:1302493/1‑91 (MQ=255)

           aTTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACataaat                                                                                >  1:1771031/1‑91 (MQ=255)

              ccACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACAtaaataaa                                                                             <  1:1798946/91‑1 (MQ=255)

              ccACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACAtaaataaa                                                                             <  1:2172837/91‑1 (MQ=255)

               cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            <  1:2934963/91‑1 (MQ=255)

               cACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATaaataaat                                                                            <  1:2935226/91‑1 (MQ=255)

                aCTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATg                                                                           <  1:1966893/91‑1 (MQ=255)

                    tataAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACt                                                                       <  1:1674400/91‑1 (MQ=255)

                      taAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTaa                                                                     >  1:3736091/1‑91 (MQ=255)

                       aaaCCTCATGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAg                                                                    <  1:2074598/91‑1 (MQ=255)

                        aaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGa                                                                   >  1:345184/1‑91 (MQ=255)

                        aaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGa                                                                   >  1:2265243/1‑91 (MQ=255)

                        aaCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGa                                                                   >  1:1439522/1‑91 (MQ=255)

                         aCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGat                                                                  <  1:368924/91‑1 (MQ=255)

                         aCCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGat                                                                  <  1:2010990/91‑1 (MQ=255)

                          ccTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGata                                                                 <  1:3683835/91‑1 (MQ=255)

                            tctTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGatata                                                               <  1:1944561/91‑1 (MQ=255)

                             ctTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGatatat                                                              <  1:1123329/91‑1 (MQ=255)

                             ctTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGatatat                                                              >  1:703932/1‑91 (MQ=255)

                              ttGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATAtctc                                                             >  1:3454389/1‑91 (MQ=255)

                               tGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATca                                                            >  1:3939701/1‑91 (MQ=255)

                                gTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcac                                                           >  1:1579276/1‑91 (MQ=255)

                                 tCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcaca                                                          >  1:2834032/1‑91 (MQ=255)

                                 tCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcaca                                                          >  1:80492/1‑91 (MQ=255)

                                 tCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATcaca                                                          >  1:771652/1‑91 (MQ=255)

                                    tGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCt                                                       >  1:1744468/1‑91 (MQ=255)

                                      gatgatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTaa                                                     <  1:2932512/91‑1 (MQ=255)

                                         gatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAAc                                                  >  1:1415864/1‑91 (MQ=255)

                                         gatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAAc                                                  >  1:1078519/1‑91 (MQ=255)

                                         gatCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAAc                                                  >  1:3412472/1‑91 (MQ=255)

                                           tCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTa                                                <  1:1145703/91‑1 (MQ=255)

                                            cTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAg                                               >  1:2216359/1‑91 (MQ=255)

                                              ttAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTg                                             <  1:2470929/91‑1 (MQ=255)

                                               tAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATGATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGa                                            <  1:2924449/91‑1 (MQ=255)

                                               tAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGa                                            <  1:1440971/91‑1 (MQ=255)

                                               tAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGa                                            >  1:1414680/1‑91 (MQ=255)

                                                aaTTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGAt                                           <  1:1955728/91‑1 (MQ=255)

                                                aaTTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGAt                                           <  1:1669867/91‑1 (MQ=255)

                                                aaTTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGAt                                           <  1:1670014/91‑1 (MQ=255)

                                                 aTTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATg                                          <  1:398871/91‑1 (MQ=255)

                                                     gttgttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCAcc                                      >  1:2367111/1‑91 (MQ=255)

                                                     gttgctTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCAcc                                      >  1:3273259/1‑91 (MQ=255)

                                                      ttgttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCt                                     <  1:1745299/91‑1 (MQ=255)

                                                       tgttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCtt                                    >  1:2425766/1‑91 (MQ=255)

                                                       tgttTTAATAACTGAACTC‑TTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   >  1:1072645/1‑91 (MQ=255)

                                                       tgttTTAATAACTGAACTC‑TTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   >  1:1072757/1‑91 (MQ=255)

                                                        gttTTAATAACTGAACTCTTTTTTCCTGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   <  1:3858492/91‑1 (MQ=255)

                                                        gttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   >  1:1105180/1‑91 (MQ=255)

                                                        gttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   >  1:1105718/1‑91 (MQ=255)

                                                        gttTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCttt                                   >  1:1474037/1‑91 (MQ=255)

                                                          tttAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaa                                 <  1:2758727/91‑1 (MQ=255)

                                                          tttAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaa                                 <  1:1047541/91‑1 (MQ=255)

                                                             aataaCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTtaataa                              >  1:3039754/1‑91 (MQ=255)

                                                                 aCTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAAtt                          <  1:1721621/91‑1 (MQ=255)

                                                                  cTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTa                         >  1:2420864/1‑91 (MQ=255)

                                                                  cTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTa                         >  1:2643186/1‑91 (MQ=255)

                                                                   tGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTaa                        >  1:132105/1‑91 (MQ=255)

                                                                     aaCTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACa                      <  1:361792/91‑1 (MQ=255)

                                                                       ctctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACatt                    <  1:3446004/91‑1 (MQ=255)

                                                                       ctctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACatt                    <  1:3283995/91‑1 (MQ=255)

                                                                       ctctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACatt                    <  1:3168025/91‑1 (MQ=255)

                                                                         ctTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattga                  <  1:2114989/91‑1 (MQ=255)

                                                                          ttttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgat                 >  1:1099117/1‑91 (MQ=255)

                                                                          ttttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgat                 >  1:1508750/1‑91 (MQ=255)

                                                                          ttttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgat                 >  1:1099054/1‑91 (MQ=255)

                                                                            ttttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACattgattg               >  1:3048071/1‑91 (MQ=255)

                                                                             tttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACAttgattga              <  1:1881955/91‑1 (MQ=255)

                                                                             tttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACAttgattga              <  1:1881893/91‑1 (MQ=255)

                                                                              ttCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAc             <  1:2366135/91‑1 (MQ=255)

                                                                                taTGTTCATTATGAACATGAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACtt           >  1:2163856/2‑91 (MQ=255)

                                                                                cATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACtt           >  1:389748/1‑91 (MQ=255)

                                                                                   gTTCATAATGAACATAAATAAATGCACTAAGATATATCACCGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACCtata        >  1:1924906/1‑91 (MQ=255)

                                                                                   gTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtata        >  1:3957148/1‑91 (MQ=255)

                                                                                   gTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTtata        <  1:2888947/91‑1 (MQ=255)

                                                                                     tCATTATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATAtt      >  1:2703097/1‑91 (MQ=255)

                                                                                     tCATAATGATCATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGAATTATAtt      >  1:2079920/1‑91 (MQ=255)

                                                                                     tCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATAtt      <  1:236146/91‑1 (MQ=255)

                                                                                        tAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTcc   <  1:579297/91‑1 (MQ=255)

                                                                                        tAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTcc   <  1:2295955/91‑1 (MQ=255)

                                                                                        tAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTcc   <  1:872261/91‑1 (MQ=255)

                                                                                         aaTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  1:1528705/91‑1 (MQ=255)

                                                                                         aaTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  1:2134070/91‑1 (MQ=255)

                                                                                         aaTGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCg  <  1:2155652/91‑1 (MQ=255)

AAACGGCAATGATTCCACTGTATAAACCTCTTGTCGTGGATGATCTTTAATTTGTTGTTTTAATAACTGAACTCTTTTTTCATGTTCATAATGAACATAAATAAATGCACTAAGATATATCACAGCTAAAACTAGTGATGCACCTTTAATAAAATTAACATTGATTGACTTATATTTCCG  >  USA300TCH1516_ALE/1794263‑1794442

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