BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,857,867	C→T	6.4%	L169L (TTG→TTA)	USA300HOU_RS09235 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,857,867	0	C	T	6.4%	637.4 / 15.0	208	L169L (TTG→TTA)	USA300HOU_RS09235	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (92/101); new base T (5/10); total (97/111)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

TTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTT  >  USA300TCH1516_ALE/1857777‑1857957

 tcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCTCTCa                                                                                           >  2:3197597/1‑91 (MQ=255)

 tcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCa                                                                                           <  1:427089/91‑1 (MQ=255)

 tcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCa                                                                                           <  1:4601127/91‑1 (MQ=255)

  cATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCaa                                                                                          >  1:2743393/1‑91 (MQ=255)

  cATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCaa                                                                                          >  1:1809109/1‑91 (MQ=255)

  cATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCaa                                                                                          >  1:2004321/1‑91 (MQ=255)

     ttGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAatt                                                                                       >  2:1935092/1‑91 (MQ=255)

     ttGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAatt                                                                                       >  2:1935728/1‑91 (MQ=255)

      tGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAatta                                                                                      <  1:5108459/91‑1 (MQ=255)

           cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGc                                                                                 >  1:5469454/1‑91 (MQ=255)

           cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGc                                                                                 >  1:766376/1‑91 (MQ=255)

           cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGc                                                                                 >  2:2914693/1‑91 (MQ=255)

            ctcTGCGCTTTTCTTTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCa                                                                                <  2:2182805/89‑1 (MQ=17)

            aGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCa                                                                                >  1:2591724/1‑91 (MQ=255)

             gCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCAc                                                                               >  1:1105345/1‑91 (MQ=255)

              cTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACt                                                                              <  1:4375959/91‑1 (MQ=255)

                gCACTTTTCTCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTg                                                                            <  2:1757223/91‑1 (MQ=255)

                gCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCTTCTTGTGATGTATCACTTAAATTATTTGCACTTg                                                                            <  2:5506738/91‑1 (MQ=255)

                gCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTg                                                                            <  1:5328142/91‑1 (MQ=255)

                gCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTg                                                                            <  2:2285253/91‑1 (MQ=255)

                 cACTTTTTTCTCTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTTAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGc                                                                           <  2:979564/91‑1 (MQ=255)

                   cTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGctt                                                                         >  2:5385889/1‑91 (MQ=255)

                    ttttCTCTGCTTCCGCTTGTGTTTTTTGTTTTGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttc                                                                        <  2:26145/91‑1 (MQ=14)

                    ttttATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCTTCTTGTGATGTATCACTTAAATTATTTGCACTTGcttc                                                                        <  2:2986657/91‑1 (MQ=255)

                    ttttATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttc                                                                        <  1:2821946/91‑1 (MQ=255)

                    ttttATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttc                                                                        >  2:500300/1‑91 (MQ=255)

                        aTCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttct                                                                    >  1:4810137/1‑91 (MQ=255)

                         tCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGcttcttctt                                                                   >  2:1459603/1‑91 (MQ=255)

                         tCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttctt                                                                   >  2:737041/1‑91 (MQ=255)

                         tCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttctt                                                                   >  2:213386/1‑91 (MQ=255)

                         tCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttctt                                                                   >  1:4928538/1‑91 (MQ=255)

                          cTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCttt                                                                  <  2:2070068/91‑1 (MQ=255)

                          cTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCttt                                                                  >  1:4015296/1‑91 (MQ=255)

                            gTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTAGCACTTGCTTCTTCTTTTa                                                                >  1:4336825/1‑91 (MQ=255)

                             tctCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAt                                                               <  2:837583/89‑1 (MQ=38)

                               tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                             <  1:106488/91‑1 (MQ=255)

                               tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                             <  1:732470/91‑1 (MQ=255)

                               tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                             <  2:2681928/91‑1 (MQ=40)

                                ttgcttgctCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgc                                                            <  2:4836072/90‑1 (MQ=255)

                                  gCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                          <  1:4164287/91‑1 (MQ=255)

                                  gCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                          <  1:2196815/91‑1 (MQ=255)

                                  gCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                          <  1:2195649/91‑1 (MQ=40)

                                    ttGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgct                                                        >  1:4958143/1‑91 (MQ=255)

                                      gtgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttg                                                      <  2:4274211/91‑1 (MQ=18)

                                      gtgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttg                                                      <  1:5065981/91‑1 (MQ=40)

                                       tgCTTCGTTCTGAGCTTCTTGTATCTCTTCTGCCTCTTGTGATGTATCACTCACATTTTTTGCATTTGCTTCTTCTTCTGTCGCTGCttga                                                     <  2:2296499/91‑2 (MQ=14)

                                         cttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttgttg                                                   <  1:913764/91‑1 (MQ=255)

                                          ttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgt                                                  <  2:1850496/91‑1 (MQ=40)

                                          ttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTtgttgt                                                  >  2:2793946/1‑91 (MQ=255)

                                           tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTtgtg                                                 >  1:3812409/1‑91 (MQ=255)

                                            cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGc                                                <  2:4498225/91‑1 (MQ=18)

                                            cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGc                                                >  1:1012059/1‑91 (MQ=255)

                                               ttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                         >  1:1728543/1‑79 (MQ=255)

                                               ttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                         <  2:1728543/79‑1 (MQ=40)

                                               ttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGcctt                                             <  2:2004321/91‑4 (MQ=25)

                                                tttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTc                                            <  1:3218360/91‑1 (MQ=255)

                                                    gCTTCTTGAATCTCTTGTGCCTCTTGTTATTTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATg                                        <  2:3077797/91‑1 (MQ=255)

                                                    gCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAatt                                                                                       >  2:1370244/1‑44 (MQ=255)

                                                    gCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATg                                        >  1:3510082/1‑91 (MQ=255)

                                                       tcttGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCa                                     <  2:845768/91‑1 (MQ=18)

                                                       tcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCa                                     >  1:5219149/1‑91 (MQ=255)

                                                        cttGAATCTCTTGTGCCTCTTGTGATTTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCAc                                    <  2:1832688/91‑1 (MQ=37)

                                                        cttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCAc                                    <  1:5564573/91‑1 (MQ=255)

                                                        cttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGCTGTGCTTTCAATGCCAc                                    <  1:2085491/91‑1 (MQ=255)

                                                         ttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACt                                   <  2:586206/91‑1 (MQ=40)

                                                           gAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGc                                 <  1:3044686/91‑1 (MQ=255)

                                                            aaTCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCt                                <  2:164600/91‑1 (MQ=18)

                                                            aaTCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCt                                <  2:3099583/91‑1 (MQ=40)

                                                               ctctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtttt                             >  2:1267000/1‑91 (MQ=38)

                                                                tctTGTTCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtttgg                            <  2:1081487/91‑3 (MQ=255)

                                                                 ctTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGa                           <  2:2022575/91‑1 (MQ=14)

                                                                 ctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgct                                                        >  1:5012944/1‑62 (MQ=255)

                                                                 ctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGa                           <  1:3571360/91‑1 (MQ=40)

                                                                 ctTGTGACTCTTGTCATGTATCACTCAAATTATTTGCAATTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGa                           <  1:4689660/91‑1 (MQ=255)

                                                                  ttGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAt                          >  2:2411658/1‑91 (MQ=255)

                                                                  ttGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAt                          >  2:2891233/1‑91 (MQ=255)

                                                                   tgtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAtt                         >  2:2085346/1‑91 (MQ=255)

                                                                   tgtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAtt                         >  2:3791128/1‑91 (MQ=255)

                                                                   tgtgCCTCTTGTGATATATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAtt                         <  2:87061/91‑1 (MQ=18)

                                                                    gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTc                        <  1:1909580/91‑1 (MQ=255)

                                                                    gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTc                        <  1:2268470/91‑1 (MQ=255)

                                                                    gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTc                        >  2:3440788/1‑91 (MQ=255)

                                                                     tgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttg                                                      >  1:1008586/1‑60 (MQ=255)

                                                                       ccTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCttt                     >  1:65793/1‑91 (MQ=255)

                                                                          ctTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAtt                  <  1:1415820/91‑1 (MQ=255)

                                                                           ttGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTGTAttt                 >  2:3942113/1‑91 (MQ=255)

                                                                              tgATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGAt              <  1:1927269/91‑1 (MQ=255)

                                                                                aTGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGTTTCTTTTTTTGTTAc            >  2:503842/1‑91 (MQ=255)

                                                                                aTGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATAc            >  1:5271460/1‑91 (MQ=255)

                                                                                aTGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATAc            >  1:1388105/1‑91 (MQ=255)

                                                                                aTGTATCACTCAAATAATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATAc            >  1:4321447/1‑91 (MQ=255)

                                                                                 tGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACa           >  1:3411531/1‑91 (MQ=255)

                                                                                  atATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAg          >  2:1773219/2‑91 (MQ=16)

                                                                                  atATCACTTAAATTATTTGCACTTGATTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAg          >  2:5383773/2‑91 (MQ=255)

                                                                                     tCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCca       <  1:2288784/91‑1 (MQ=255)

                                                                                      cACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCcac      <  2:3177448/91‑1 (MQ=255)

                                                                                       aCTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCcaca     <  2:2625585/91‑1 (MQ=255)

                                                                                       aCTCAAATTATTTGCACTTGCTTCTTCTTTTATCCCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCcaca     <  1:3885295/91‑1 (MQ=255)

                                                                                         tCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGTTTCTTTATTTGATACAGCCActct   >  2:4329789/1‑91 (MQ=255)

                                                                                         tCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACt   <  1:4850443/91‑1 (MQ=255)

TTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTT  >  USA300TCH1516_ALE/1857777‑1857957

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