BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	9,716	T→C	23.4%	intergenic (+21/+66)	gyrA → / ← nnrD	DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	9,716	0	T	C	23.4%	330.2 / 110.9	180	intergenic (+21/+66)	gyrA/nnrD	DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand): ref base T (65/73); new base C (19/23); total (84/96)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 187 total aligned reads displayed.

TATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCG  >  USA300TCH1516_ALE/9626‑9806

 aTTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCtata                                                                                             <  2:888624/91‑1 (MQ=255)

 aTTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCtata                                                                                             >  1:1963002/1‑91 (MQ=255)

   tGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATg                                                                                           <  2:954044/91‑1 (MQ=255)

    gAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGt                                                                                          <  1:315108/91‑1 (MQ=255)

     aggTAAGACAAGATGTCATGGAGCGTGTTGAAGAAGAGATACAACAATCATCAGATGAAGAATAATGAAAAATAAGACTCCCCTACAtata                                                                                         <  2:1619734/89‑1 (MQ=37)

     aaGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTa                                                                                         <  2:2170414/91‑1 (MQ=255)

      aGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAg                                                                                        >  1:3073721/1‑91 (MQ=255)

       gTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgg                                                                                       <  1:252902/91‑1 (MQ=38)

        tAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAggg                                                                                      >  1:771646/1‑91 (MQ=255)

        tAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAggg                                                                                      >  1:599559/1‑91 (MQ=255)

          aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                                    <  1:1006918/91‑1 (MQ=255)

          aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGa                                                                                    <  1:134497/91‑1 (MQ=255)

          aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgggaa                                                                                    <  1:1624675/91‑3 (MQ=38)

          aGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATAAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgggaa                                                                                    <  2:1341217/91‑3 (MQ=37)

           gACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgggaag                                                                                   <  1:1208696/91‑4 (MQ=37)

            aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                                  <  1:1873555/91‑1 (MQ=255)

            aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                                  <  1:2102853/91‑1 (MQ=255)

            aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGt                                                                                  <  1:247166/91‑1 (MQ=37)

             cAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTc                                                                                 <  1:1072254/91‑1 (MQ=255)

             cAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTc                                                                                 >  2:2368005/1‑91 (MQ=255)

              aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCt                                                                                <  1:3407026/91‑1 (MQ=255)

              aaGATTTCAGGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAAGAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                                <  2:1905241/91‑1 (MQ=25)

               aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCtt                                                                               >  1:860877/1‑91 (MQ=255)

               aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAAT‑AAAAATAAGACTCCCCTACATATAGGGAAGTCTTa                                                                              <  1:1371060/91‑1 (MQ=25)

                gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTa                                                                              <  1:3476934/91‑1 (MQ=255)

                gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTa                                                                              <  2:2983056/91‑1 (MQ=25)

                gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTa                                                                              <  1:3052801/91‑1 (MQ=25)

                    tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑tttt                                                                         >  1:688297/1‑91 (MQ=255)

                    tCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑tttt                                                                         >  2:3376266/1‑91 (MQ=255)

                     cATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑ttttt                                                                        >  1:2927379/1‑91 (MQ=255)

                        ggATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATg                                                                     <  1:2304744/91‑1 (MQ=255)

                           tCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTa                                                                 <  1:2328598/91‑1 (MQ=255)

                            cGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAg                                                                >  2:3198558/1‑91 (MQ=255)

                               gtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGaaa                                                             >  2:932103/1‑91 (MQ=255)

                               gtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGaaa                                                             >  1:3761316/1‑91 (MQ=255)

                               gtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGaaa                                                             >  1:3496806/1‑91 (MQ=255)

                                 tGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGt                                                           <  2:2033046/91‑1 (MQ=255)

                                 tGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGt                                                           <  1:1146316/91‑1 (MQ=255)

                                 ggaagaaGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGt                                                           <  2:129467/90‑1 (MQ=255)

                                   aagaagATATACATCAATCAGCAGATGAAGAATAATGAAAAATAAGACTTCCCTATATGTAGGGGAGGCTTA‑TTTTTATG‑CTAGAAAGTaa                                                         <  2:3573826/91‑1 (MQ=255)

                                     gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATg                                                       <  1:1716448/91‑1 (MQ=255)

                                      aagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGc                                                      <  2:2013022/91‑1 (MQ=255)

                                       agaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATACGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCt                                                     >  2:410667/1‑91 (MQ=255)

                                        gaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCtt                                                    >  1:3317195/1‑91 (MQ=255)

                                        gaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCtt                                                    <  2:1596880/91‑1 (MQ=255)

                                         atatACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCttt                                                   >  1:2646715/1‑91 (MQ=255)

                                           ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTAc                                                 >  2:1150101/1‑91 (MQ=255)

                                           ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTAc                                                 <  1:2792042/91‑1 (MQ=255)

                                           ataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTAc                                                 >  2:3947671/1‑91 (MQ=255)

                                              caacaaTCATCATATGAAGTATAATAAAAAATAAGACAACCCTATAACTAGGGCAATCTAA‑TTTTTATG‑CTAAAAAATAATGCTCTACtat                                              >  2:2523970/1‑91 (MQ=255)

                                              caacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACtat                                              >  2:2257589/1‑91 (MQ=255)

                                               aacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACtata                                             <  2:1278337/91‑1 (MQ=255)

                                                acaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTatat                                            <  1:1828959/91‑1 (MQ=255)

                                                acaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTatat                                            <  2:2956858/91‑1 (MQ=255)

                                                acaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTatat                                            <  1:1995436/91‑1 (MQ=255)

                                                 caaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATAtt                                           >  1:2014350/1‑91 (MQ=255)

                                                  aaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTc                                          >  1:1722424/1‑91 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATg                                      <  1:2138052/91‑1 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATg                                      <  1:695620/91‑1 (MQ=255)

                                                      atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATg                                      <  1:2215603/91‑1 (MQ=255)

                                                       tcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGATAGTAATGCTTTACTATATTCAATGa                                     >  1:3686255/1‑91 (MQ=255)

                                                        caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGAt                                    >  1:799431/1‑91 (MQ=255)

                                                        caGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGAt                                    >  2:169160/1‑91 (MQ=255)

                                                           aTGAAGAATAATGAAGAATGAGACTTCCCTATGTGTAGGGGAGTCTTA‑TTTGTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAg                                 <  2:1118328/91‑1 (MQ=255)

                                                           aTGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAg                                 >  2:1725377/1‑91 (MQ=255)

                                                            tGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTact                                >  1:1845275/1‑89 (MQ=255)

                                                            tGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGt                                <  2:3840596/91‑1 (MQ=255)

                                                              aagaaTAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTaa                              <  1:2241854/91‑1 (MQ=255)

                                                                 aataatAAAAACTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGa                           <  1:3587163/91‑1 (MQ=255)

                                                                   taataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACt                         >  2:2192943/1‑91 (MQ=255)

                                                                   taataaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACt                         <  1:624699/91‑1 (MQ=255)

                                                                      taaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAAc                      <  1:1757357/91‑1 (MQ=255)

                                                                        aaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACtt                    <  2:3164500/91‑1 (MQ=255)

                                                                         aaaaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTCTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGCTTAGTAATGACTAACttt                   >  2:900351/1‑91 (MQ=255)

                                                                           aaTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATTACTAACtttt                  >  1:1413826/1‑90 (MQ=21)

                                                                           aaTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTc                  >  2:1837227/1‑91 (MQ=21)

                                                                           aaTAAGACTCCCCTACATATAGGGAAGTCTTATTTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTc                  >  1:508142/1‑91 (MQ=21)

                                                                              aaGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              <  2:2637953/91‑1 (MQ=255)

                                                                              aaGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              <  1:318779/91‑1 (MQ=255)

                                                                              aaGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              >  1:493420/1‑91 (MQ=255)

                                                                              aaGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              <  2:2356714/91‑1 (MQ=255)

                                                                              aaGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAt              >  2:2756654/1‑91 (MQ=37)

                                                                              aaGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTAAGACTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTaa               >  2:1253569/1‑91 (MQ=25)

                                                                               aGACTCGCCTACATATAGGGAAGGCTTA‑TTTTTATG‑CTAGTAAGTAATGCTTTACTATATTCTATGATTAGTAATGATTAAGTTTCTAAtt             >  2:267844/1‑91 (MQ=25)

                                                                                gACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTg            >  2:2071565/1‑91 (MQ=255)

                                                                                 actccccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGt           >  1:2876045/5‑91 (MQ=37)

                                                                                 aCTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAACGTAATGCTTTACTATATTCAATGATCAGTAATGACTAACTTTCTAATTGt           >  2:3668800/1‑91 (MQ=255)

                                                                                 aCTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGt           >  2:2407560/1‑91 (MQ=255)

                                                                                 aCTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGt           <  1:3765088/91‑1 (MQ=255)

                                                                                  cTTCCCTATATGTAGGGGCGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATACTCAATGATTAGTAATGACTAACTTTCTAATTGtt          >  2:2122332/1‑91 (MQ=255)

                                                                                  cTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGtt          >  1:423382/1‑91 (MQ=255)

                                                                                   ttCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt         >  2:1520150/1‑91 (MQ=255)

                                                                                   tccccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt         >  1:2839684/3‑91 (MQ=38)

                                                                                   tccccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt         >  2:2394189/3‑91 (MQ=38)

                                                                                    ccccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGCAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTc        >  2:2154192/2‑91 (MQ=38)

                                                                                    ccccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGCAATGCTTTACTATACTCAATGATTAGTAATGACTAACTTTCTAATTGTTTc        >  2:3114012/2‑91 (MQ=37)

                                                                                     cccTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCa       >  2:676808/1‑91 (MQ=255)

                                                                                     cccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTCGTAATGTCTAACTTTCTACTTGTTTCa       >  2:1557054/1‑91 (MQ=37)

                                                                                         atatGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTATACTATATTCAATGATTAGAAATGACTAACTTTCTAATTGTTTCATTGc   <  2:3478531/91‑1 (MQ=255)

                                                                                          tatGTAGGGGAGTCTTA‑TTTTTATG‑CTAGCAAGTAATGCTTTACTATATTCAATGATTCGTAATGACTACCCTTCTAATTGTTTCATTGCg  >  2:567094/1‑91 (MQ=255)

TATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCG  >  USA300TCH1516_ALE/9626‑9806

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