BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	690,106	C→T	9.2%	128 (TAG→TAA)	USA300HOU_RS03325 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	690,106	0	C	T	9.2%	362.2 / 18.7	117	128 (TAG→TAA)	USA300HOU_RS03325	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (54/52); new base T (7/4); total (61/56)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAATATT  >  USA300TCH1516_ALE/690016‑690196

gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                            <  1:17967/91‑1 (MQ=255)

  tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTa                                                                                          >  1:2216619/1‑91 (MQ=255)

  tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTa                                                                                          >  1:3225963/1‑91 (MQ=255)

   ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAACATAACCCCAATAGGGGCTAt                                                                                         <  1:3334394/91‑1 (MQ=255)

   ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                         <  1:1102125/91‑1 (MQ=255)

   ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                         <  1:3438701/91‑1 (MQ=255)

   ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                         <  1:1189600/91‑1 (MQ=255)

   ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAt                                                                                         <  1:2341015/91‑1 (MQ=255)

    tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAtt                                                                                        <  1:3340876/91‑1 (MQ=255)

    tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAtt                                                                                        <  1:3651481/91‑1 (MQ=255)

    tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAtt                                                                                        <  1:2313710/91‑1 (MQ=255)

     cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                                       >  1:837680/1‑91 (MQ=255)

      gTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGCCTATTTa                                                                                      >  1:2085402/1‑91 (MQ=255)

         aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                                   <  1:2098044/91‑1 (MQ=255)

         aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                                   >  1:3496904/1‑91 (MQ=255)

          ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                                  <  1:3972797/91‑1 (MQ=255)

            aaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt                                                                                >  1:679291/1‑91 (MQ=255)

               aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                             <  1:2426814/91‑1 (MQ=255)

               aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                             <  1:2426900/91‑1 (MQ=255)

                  ttATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttt                                                                          <  1:2339649/91‑1 (MQ=255)

                   tATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttttt                                                                         >  1:432966/1‑91 (MQ=255)

                    aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                                        >  1:2634939/1‑91 (MQ=255)

                    aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                                        >  1:2635182/1‑91 (MQ=255)

                    aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                                        >  1:1368435/1‑91 (MQ=255)

                      ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                      >  1:4011673/1‑91 (MQ=255)

                           aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                                 <  1:1463278/91‑1 (MQ=255)

                           aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                                 <  1:560177/91‑1 (MQ=255)

                            aaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGat                                                                <  1:267220/91‑1 (MQ=255)

                            aaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGat                                                                <  1:267071/91‑1 (MQ=255)

                              ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                              >  1:4135983/1‑91 (MQ=255)

                               cTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATAtt                                                             >  1:537640/1‑91 (MQ=255)

                                 gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                           >  1:878465/1‑91 (MQ=255)

                                 gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                           >  1:197550/1‑91 (MQ=255)

                                  cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                          <  1:3768329/91‑1 (MQ=255)

                                    tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                        >  1:362317/1‑91 (MQ=255)

                                    tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                        >  1:362329/1‑91 (MQ=255)

                                    tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                        >  1:363728/1‑91 (MQ=255)

                                      cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                                      >  1:211427/1‑91 (MQ=255)

                                            ttttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTAATTTTTTTATCATGATATTCTACTTCTtcatc                                                >  1:3528615/1‑91 (MQ=255)

                                               tAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATa                                             >  1:2973441/1‑91 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           <  1:3897859/91‑1 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           <  1:3637444/91‑1 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           <  1:74485/91‑1 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           >  1:2769781/1‑91 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           >  1:2770504/1‑91 (MQ=255)

                                                  aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                          <  1:2343374/91‑1 (MQ=255)

                                                    tCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTa                                        >  1:518776/1‑91 (MQ=255)

                                                    tCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTa                                        >  1:519310/1‑91 (MQ=255)

                                                     cATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTaa                                       >  1:1424503/1‑91 (MQ=255)

                                                      aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                      <  1:2407799/91‑1 (MQ=255)

                                                      aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                      <  1:778392/91‑1 (MQ=255)

                                                         tCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTa                                   >  1:1674958/1‑91 (MQ=255)

                                                           aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                                 <  1:3265505/91‑1 (MQ=255)

                                                           aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                                 <  1:946273/91‑1 (MQ=255)

                                                            tCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACt                                >  1:1875606/1‑91 (MQ=255)

                                                            tCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACt                                >  1:1826875/1‑91 (MQ=255)

                                                               cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                             <  1:2944272/91‑1 (MQ=255)

                                                                tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                            <  1:3083575/91‑1 (MQ=255)

                                                                tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                            <  1:1348628/91‑1 (MQ=255)

                                                                 ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                           <  1:238979/91‑1 (MQ=255)

                                                                 ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                           <  1:239098/91‑1 (MQ=255)

                                                                 ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                           <  1:239131/91‑1 (MQ=255)

                                                                 ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                           <  1:1638754/91‑1 (MQ=255)

                                                                  cAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCtt                          <  1:2119100/91‑1 (MQ=255)

                                                                    atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                        >  1:1245626/5‑91 (MQ=37)

                                                                    atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                        >  1:1245771/5‑91 (MQ=37)

                                                                    aTACAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTCTCTTaa                        <  1:3220580/91‑1 (MQ=255)

                                                                     tAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAAc                       <  1:3559417/91‑1 (MQ=255)

                                                                      aaaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACt                      <  1:2595097/91‑1 (MQ=255)

                                                                       aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg                     <  1:2061425/91‑1 (MQ=255)

                                                                        aaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGc                    >  1:2260855/1‑91 (MQ=255)

                                                                         aaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca                   <  1:3501840/91‑1 (MQ=255)

                                                                          aTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTGATCATAACTTAATTTAACTTTTCTTAACTGcac                  >  1:3404598/1‑91 (MQ=255)

                                                                          aTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGcac                  >  1:3404578/1‑91 (MQ=255)

                                                                            agccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt                >  1:933621/3‑91 (MQ=255)

                                                                             gccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt               >  1:1675409/2‑91 (MQ=255)

                                                                             aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt               <  1:1287514/91‑1 (MQ=255)

                                                                             aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt               <  1:1287547/91‑1 (MQ=255)

                                                                             aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt               <  1:1289164/91‑1 (MQ=255)

                                                                             aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt               <  1:2465576/91‑1 (MQ=255)

                                                                              ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa              >  1:1365776/1‑91 (MQ=255)

                                                                              ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa              >  1:3781043/1‑91 (MQ=255)

                                                                              ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa              <  1:2599744/91‑1 (MQ=255)

                                                                              ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa              >  1:1078421/1‑91 (MQ=255)

                                                                               cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa             >  1:3208031/1‑91 (MQ=255)

                                                                               cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa             >  1:3207891/1‑91 (MQ=255)

                                                                                ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa            >  1:2825228/1‑91 (MQ=255)

                                                                                 cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa           >  1:1695619/1‑91 (MQ=255)

                                                                                 cAATAGGGGCTATTTAAATGTTATTTTTTTAGCATGATATTCTAGTTCTTCATTATAGCTTAATTTAACTTTTCTTAACTGCACATTaaaa           <  1:3311904/91‑1 (MQ=255)

                                                                                  tattggggTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaaa          >  1:1607091/5‑91 (MQ=255)

                                                                                    ttggggTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACa        >  1:1429033/3‑91 (MQ=255)

                                                                                     aGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACaa       >  1:3943371/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:1183904/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:2885957/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:41600/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:2560966/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:2560986/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:152230/1‑91 (MQ=255)

                                                                                         gCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAatat   <  1:3256859/91‑1 (MQ=255)

                                                                                          cTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAATAtt  <  1:909771/91‑1 (MQ=255)

GTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAATATT  >  USA300TCH1516_ALE/690016‑690196

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