BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	164,586	T→A	14.4%	Q30L (CAA→CTA)	phnD ←	Phosphate‑import protein PhnD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	164,586	0	T	A	14.4%	341.7 / 44.7	117	Q30L (CAA→CTA)	phnD	Phosphate‑import protein PhnD
Reads supporting (aligned to +/- strand): ref base T (50/50); new base A (2/15); total (52/65)
Fisher's exact test for biased strand distribution p-value = 3.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

CCAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCAT  >  USA300TCH1516_ALE/164496‑164674

ccAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCtttt                                                                                          >  1:913308/1‑91 (MQ=255)

ccAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCtttt                                                                                          <  1:1663544/91‑1 (MQ=255)

 cAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTg                                                                                         <  1:347554/91‑1 (MQ=255)

 cAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTg                                                                                         <  1:312609/91‑1 (MQ=255)

 cAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTGTTg                                                                                         >  1:3155419/1‑91 (MQ=255)

  aGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGa                                                                                        <  1:2047410/91‑1 (MQ=255)

   gCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAt                                                                                       >  1:245512/1‑91 (MQ=255)

    cATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAtt                                                                                      >  1:813364/1‑91 (MQ=255)

    cATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAtt                                                                                      >  1:2026255/1‑91 (MQ=255)

    cATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGAtt                                                                                      >  1:1018112/1‑91 (MQ=255)

     aTTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTa                                                                                     <  1:3013009/91‑1 (MQ=255)

      ttttGCGAAGGTTCAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTAt                                                                                    >  1:3175332/1‑91 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATc                                                                                   >  1:728370/1‑91 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATc                                                                                   >  1:1321779/1‑91 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                   <  1:2657031/91‑1 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                   <  1:2657096/91‑1 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                   <  1:1836110/91‑1 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                   <  1:3013890/91‑1 (MQ=255)

       tttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATc                                                                                   <  1:3013912/91‑1 (MQ=255)

        ttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCt                                                                                  <  1:3508515/91‑1 (MQ=255)

        ttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATca                                                                                  <  1:1103303/91‑2 (MQ=255)

        ttGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATca                                                                                  <  1:3930495/91‑2 (MQ=255)

         tGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATcaa                                                                                 <  1:3798191/91‑3 (MQ=255)

         tGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATcaa                                                                                 <  1:1107046/91‑3 (MQ=255)

          gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                >  1:1420922/1‑91 (MQ=255)

          gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                >  1:2945371/1‑91 (MQ=255)

          gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                <  1:487499/91‑1 (MQ=255)

          gCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaa                                                                                <  1:2806186/91‑1 (MQ=255)

           cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                               >  1:599396/1‑91 (MQ=255)

           cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                               >  1:739621/1‑91 (MQ=255)

           cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                               >  1:739819/1‑91 (MQ=255)

           cGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTaaa                                                                               >  1:607218/1‑91 (MQ=255)

            gAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAg                                                                              >  1:3221357/1‑91 (MQ=255)

            gAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAg                                                                              <  1:2339115/91‑1 (MQ=255)

            gAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAg                                                                              <  1:2339091/91‑1 (MQ=255)

             aaGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGa                                                                             >  1:1846697/1‑91 (MQ=255)

              aGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTAGATAATCAAAAGaa                                                                            <  1:3364400/91‑1 (MQ=255)

               ggTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAAc                                                                           <  1:1450706/91‑1 (MQ=255)

               ggTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAAc                                                                           >  1:1986358/1‑91 (MQ=255)

                gTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACt                                                                          <  1:1456064/91‑1 (MQ=255)

                 tACAAATTGGACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                         <  1:2875607/91‑1 (MQ=255)

                 tACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                         <  1:3712655/91‑1 (MQ=255)

                 tACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACtt                                                                         <  1:1664647/91‑1 (MQ=255)

                  aCAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTg                                                                        <  1:1510858/91‑1 (MQ=255)

                   cAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGa                                                                       <  1:3236129/91‑1 (MQ=255)

                     aaTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGt                                                                     >  1:259500/1‑91 (MQ=255)

                     aaTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGt                                                                     >  1:3054614/1‑91 (MQ=255)

                     aaTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGt                                                                     <  1:746945/91‑1 (MQ=255)

                      aTTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGtt                                                                    <  1:283273/91‑1 (MQ=255)

                        tGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTc                                                                  <  1:3323396/91‑1 (MQ=255)

                        tGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTc                                                                  <  1:214317/91‑1 (MQ=255)

                          aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                <  1:1861410/91‑1 (MQ=255)

                          aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                >  1:1429593/1‑91 (MQ=255)

                          aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                >  1:1429232/1‑91 (MQ=255)

                          aaCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                <  1:510807/91‑1 (MQ=255)

                          aaCGGTTAAGTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCca                                                                <  1:3965829/91‑1 (MQ=255)

                           aCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcac                                                               <  1:2787220/91‑1 (MQ=255)

                           aCGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcac                                                               <  1:1246320/91‑1 (MQ=255)

                            cGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcaca                                                              >  1:1075159/1‑91 (MQ=255)

                            cGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCcaca                                                              >  1:1075507/1‑91 (MQ=255)

                               ttAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgc                                                           >  1:3009779/1‑91 (MQ=255)

                                 aaTTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgctg                                                         <  1:1606907/91‑1 (MQ=255)

                                    tCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgctgctg                                                      <  1:2700596/91‑1 (MQ=255)

                                    tCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACAtgctgctg                                                      <  1:2700490/91‑1 (MQ=255)

                                       ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                   >  1:4013811/1‑91 (MQ=255)

                                       ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                   >  1:3979227/1‑91 (MQ=255)

                                       ttAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaa                                                   >  1:3000367/1‑91 (MQ=255)

                                        tAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaaa                                                  <  1:3893392/91‑1 (MQ=255)

                                        tAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCaaa                                                  >  1:3191511/1‑91 (MQ=255)

                                           gTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAAc                                               <  1:449233/91‑1 (MQ=255)

                                           gTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAAc                                               >  1:3580704/1‑91 (MQ=255)

                                            tttGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACa                                              <  1:1304222/91‑1 (MQ=255)

                                                     tcctGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCt                                     <  1:3673715/91‑1 (MQ=255)

                                                         gATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACa                                 >  1:544701/1‑91 (MQ=255)

                                                           tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                               <  1:1368626/91‑1 (MQ=255)

                                                           tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                               <  1:412072/91‑1 (MQ=255)

                                                           tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTCCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                               <  1:411750/91‑1 (MQ=255)

                                                           tttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTCCCACATGCTGCTGCAAAAACAATGACTGCTAACAtt                               <  1:412426/91‑1 (MQ=255)

                                                            ttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTa                              >  1:3040954/1‑91 (MQ=255)

                                                            ttAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTa                              >  1:3041214/1‑91 (MQ=255)

                                                              aGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAAt                            >  1:1903004/1‑91 (MQ=255)

                                                              aGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAAt                            >  1:3467227/1‑91 (MQ=255)

                                                                aaTCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATAc                          <  1:963613/91‑1 (MQ=255)

                                                                aaTCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATAc                          >  1:4743/1‑91 (MQ=255)

                                                                   ccGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACaaa                       >  1:619403/1‑91 (MQ=255)

                                                                     gAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATa                     >  1:2779163/1‑91 (MQ=255)

                                                                          aTTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAacac                <  1:875653/91‑1 (MQ=255)

                                                                           ttACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACt               >  1:2135500/1‑91 (MQ=255)

                                                                            tACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACtt              <  1:3665124/91‑1 (MQ=255)

                                                                            tACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACtt              <  1:1623026/91‑1 (MQ=255)

                                                                              cTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaa            >  1:3022351/1‑91 (MQ=255)

                                                                              cTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaa            >  1:3022330/1‑91 (MQ=255)

                                                                                aGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaaaa          <  1:60709/91‑1 (MQ=255)

                                                                                aGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTaaaa          <  1:3098957/91‑1 (MQ=255)

                                                                                   gTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAAttt       >  1:2635703/1‑91 (MQ=255)

                                                                                    ttCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAAtttt      >  1:72515/1‑91 (MQ=255)

                                                                                     tCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAAttttt     <  1:348425/91‑1 (MQ=255)

                                                                                       ttttGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCa   >  1:2907115/1‑91 (MQ=255)

                                                                                       ttttGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCa   >  1:2511069/1‑91 (MQ=255)

                                                                                        tttGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCAt  >  1:3358869/1‑91 (MQ=255)

CCAGCATTTTGCGAAGGTACAAATTGAACGGTTAATTCTTTAGGTTTGTATCCTCCTGATTTAGAATCCGAATCATTACTAGCGTTCTTTTGATTATCTAAAGAACTTGAGTTTCCACATGCTGCTGCAAAAACAATGACTGCTAACATTAATACAAATAAACACTTAAAATTTTTCAT  >  USA300TCH1516_ALE/164496‑164674

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