BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,857,867	C→T	5.6%	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	5.6%	786.6 / 22.0	236	L169L (TTG→TTA)	USA300HOU_RS09235	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (111/111); new base T (8/6); total (119/117)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

TTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACT  >  USA300TCH1516_ALE/1857777‑1857956

ttcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTc                                                                                           <  1:2899643/91‑1 (MQ=255)

ttcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGCATCACTc                                                                                           <  2:510144/91‑1 (MQ=255)

      tGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAatta                                                                                     >  1:1054586/1‑91 (MQ=255)

       tttACTGCTGCACTTTTATCCGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCTTCTTGTGATGTATCACTTAAattat                                                                                    <  2:71809/89‑1 (MQ=1)

        ttACTGCTGCACTTTTATCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAattatt                                                                                   <  2:208406/90‑1 (MQ=255)

        aTACAGCTGCACTTTTATCTGTATCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAattatt                                                                                   >  2:4452353/1‑91 (MQ=255)

        aTACAGCTGAACTTTTCTCTGTTTCTGCTCGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAattatt                                                                                   >  1:5376410/1‑91 (MQ=255)

         tACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTAttt                                                                                  >  2:3573440/1‑91 (MQ=255)

         tACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTAttt                                                                                  >  1:2496723/1‑91 (MQ=255)

         tACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTAttt                                                                                  >  1:2610351/1‑91 (MQ=255)

         tACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTATGCCTCTTGTGATGTATCACTCAAATTAttt                                                                                  >  1:5432065/1‑91 (MQ=255)

           cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGc                                                                                <  1:1374655/91‑1 (MQ=255)

           cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGc                                                                                <  1:3648148/91‑1 (MQ=255)

            aGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGCATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCa                                                                               >  1:5700043/1‑91 (MQ=255)

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

             gttGCAGTTTTATCTGTTACTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCAc                                                                              <  1:4913797/89‑1 (MQ=255)

             gCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCAc                                                                              <  1:3782817/91‑1 (MQ=255)

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

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

               tGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACtt                                                                            <  1:5084148/91‑1 (MQ=255)

                gCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTg                                                                           >  1:5238319/1‑91 (MQ=255)

                 cACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGc                                                                          >  1:420550/1‑91 (MQ=255)

                   cTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGctt                                                                        >  2:3202079/1‑91 (MQ=255)

                     tttATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttct                                                                      <  1:4417647/91‑1 (MQ=255)

                      ttATCTGTTTCTGCTCGTGCTTCCTTTTTTGCTTCTTAAATCTCTTGTGCGTTCTGTGATGTATCACTCAAATTATTTGCACTTGcttctt                                                                     <  2:3729501/91‑1 (MQ=14)

                       tATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttc                                                                    >  2:5349707/1‑91 (MQ=255)

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

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

                             tttCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAt                                                              >  2:115659/1‑91 (MQ=255)

                               tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                            >  2:4145754/1‑91 (MQ=18)

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

                               tCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                            >  1:1129723/1‑91 (MQ=40)

                                 tGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgct                                                          <  1:4452353/91‑1 (MQ=255)

                                  gCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                         >  2:3161895/1‑91 (MQ=18)

                                  gCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                         >  2:1669977/1‑91 (MQ=255)

                                  gCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                         >  1:4967328/1‑91 (MQ=255)

                                   cTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCCCTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                        >  2:2679288/1‑91 (MQ=255)

                                   cTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                        <  1:952553/91‑1 (MQ=255)

                                   cTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                        >  2:3968113/1‑91 (MQ=255)

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

                                      gtgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttg                                                     <  2:3239330/91‑1 (MQ=255)

                                       tgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttgt                                                    <  1:4653296/91‑1 (MQ=255)

                                          ttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGc                                                                          >  2:3494975/1‑66 (MQ=255)

                                          ttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGc                                                                          >  2:1396772/1‑66 (MQ=255)

                                           tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGc                                                                          <  1:5168184/65‑1 (MQ=255)

                                           tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGc                                                                          >  2:5168184/1‑65 (MQ=255)

                                           tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                        >  1:1830057/1‑83 (MQ=255)

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

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

                                            cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                            <  1:1561805/78‑1 (MQ=255)

                                              tttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCtt                                             >  2:1972181/1‑91 (MQ=40)

                                               ttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCttt                                            >  1:3486593/1‑91 (MQ=255)

                                                 ttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCa                                          <  1:57781/91‑1 (MQ=255)

                                                  tAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCaa                                         <  2:5238319/91‑1 (MQ=255)

                                                   aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAAt                                        >  1:2456877/1‑91 (MQ=255)

                                                   aGCTTCTTGAATCTCTTGTGCCTCTTGTGAGGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgct                                                          <  2:2809526/73‑1 (MQ=255)

                                                    gCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGcttcttct                                                                   <  1:2814577/63‑1 (MQ=255)

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

                                                     cttcttGATTTTCTTGTGCCTCTTGTGTTGTATCACTCGACTTTTTTGCGCTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGc                                      <  2:3056878/91‑1 (MQ=14)

                                                       tcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGc                                               <  1:479275/80‑1 (MQ=40)

                                                       tcttGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCa                                    <  2:3690194/91‑1 (MQ=255)

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

                                                        cttGAATCTCTTGTGTCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCAc                                   <  1:4143882/91‑1 (MQ=255)

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

                                                           gAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGc                                <  2:558855/91‑1 (MQ=40)

                                                             aTCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtt                              >  1:1086131/1‑91 (MQ=255)

                                                              tctcTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                        <  1:778566/64‑1 (MQ=255)

                                                              tctcTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCttt                             >  2:146182/1‑91 (MQ=40)

                                                              tccctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGTACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCTCTGCttt                             >  1:2018309/4‑91 (MQ=255)

                                                               ctctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTTCtttt                            >  2:1602173/1‑91 (MQ=255)

                                                               ctctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtttt                            >  1:2703965/1‑91 (MQ=39)

                                                                 ctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGa                          >  1:671187/1‑91 (MQ=40)

                                                                 ctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGa                          >  2:554878/1‑91 (MQ=40)

                                                                  ttGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAt                                                              <  2:5719155/54‑1 (MQ=255)

                                                                  ttGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCGCTGCTTTTGAt                         >  2:584462/1‑91 (MQ=255)

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

                                                                     tgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCt                      >  1:5841186/1‑91 (MQ=255)

                                                                      gCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCtt                     <  1:285695/91‑1 (MQ=255)

                                                                        ctctTGTGATATATCACTCAAATTATTTGCACTTGc                                                                          >  1:4003168/1‑36 (MQ=255)

                                                                         tctTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAt                  >  2:2913701/1‑91 (MQ=255)

                                                                         tctTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAt                  <  2:5376410/91‑1 (MQ=255)

                                                                         tctTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAt                  <  1:984581/91‑1 (MQ=255)

                                                                          ctTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAAt                                        <  1:1041977/68‑1 (MQ=255)

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

                                                                           ttGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATg                                       <  1:1178977/68‑1 (MQ=255)

                                                                           ttGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAttt                >  2:783763/1‑91 (MQ=255)

                                                                           ttGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAttt                >  2:5271858/1‑91 (MQ=255)

                                                                           ttGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAttt                >  1:4058840/1‑91 (MQ=255)

                                                                              tgATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCg                                                            <  1:1233422/44‑1 (MQ=40)

                                                                              tgATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGAt             >  2:4051270/1‑91 (MQ=255)

                                                                                  gTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCt                      <  2:390144/78‑1 (MQ=255)

                                                                                  atATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAg         >  1:3269569/2‑91 (MQ=9)

                                                                                   tATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGc        >  1:5627946/1‑91 (MQ=255)

                                                                                    aTCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGcc       <  1:1909158/91‑1 (MQ=255)

                                                                                    aTCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTAAATGCCACTGCTTTTGATTCTTTTTTTGATACAGcc       >  2:1749299/1‑91 (MQ=255)

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

                                                                                      cACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTTCTTTTGTTTCTTTATTTGTTACAGCcac     >  2:1055243/1‑91 (MQ=255)

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

                                                                                         tCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACt  >  1:2545907/1‑91 (MQ=255)

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

TTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACT  >  USA300TCH1516_ALE/1857777‑1857956

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