BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	9,710	T→C	19.1%	intergenic (+15/+72)	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,710	0	T	C	19.1%	392.2 / 86.5	184	intergenic (+15/+72)	gyrA/nnrD	DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand): ref base T (68/79); new base C (11/24); total (80/104)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 188 total aligned reads displayed.

TGGACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTC  >  USA300TCH1516_ALE/9620‑9800

tGGACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACtt                                                                                              >  2:1105969/1‑91 (MQ=255)

  gACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTcc                                                                                            >  1:130614/1‑91 (MQ=255)

    ggTATTGAAGTAAGACAAGATTTCATGGAGCGTGTTGAAGAAGATATACAACAATCATTTGATGAAGAATAATAAAAAATAAGACTTCCCt                                                                                          <  2:4805/90‑1 (MQ=255)

    cGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCt                                                                                          >  2:430732/1‑91 (MQ=255)

    cGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCt                                                                                          <  2:2287470/91‑1 (MQ=255)

    cGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCt                                                                                          >  1:2154481/1‑91 (MQ=255)

     gTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCta                                                                                         <  2:747977/91‑1 (MQ=255)

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

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

         tGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCtacata                                                                                     <  1:847926/91‑5 (MQ=255)

          gAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGt                                                                                    >  2:2058002/1‑91 (MQ=255)

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

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

            aGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAg                                                                                  <  1:1985254/91‑1 (MQ=255)

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

             gTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAAAAAAAAAAAATAAGACTCCCCTACATATAgg                                                                                 <  2:2338549/91‑1 (MQ=38)

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

              tAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAggg                                                                                <  1:1739802/91‑1 (MQ=38)

              tAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAggg                                                                                <  1:2220417/91‑1 (MQ=38)

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

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

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

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

                 gACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAgggaag                                                                             <  2:2965929/91‑4 (MQ=37)

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

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

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

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

                    aaGATTTGATGGATCGTGTGGAAGAAGAGATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                          <  2:3586796/91‑1 (MQ=25)

                    aaGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCt                                                                          <  1:2436253/91‑1 (MQ=37)

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

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

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

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

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

                       aTTTCATGGATCGTGTTGCAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCttct                                                                       >  1:1932270/1‑89 (MQ=255)

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

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

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

                                 tCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTa                                                           >  2:1552696/1‑91 (MQ=255)

                                  cGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGCGTCTTA‑TTTTTATG‑CTAg                                                          >  1:1351794/1‑91 (MQ=255)

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

                                     gtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGaaa                                                       <  1:2109316/91‑1 (MQ=255)

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

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

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

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

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

                                           gaagaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAGGTAATg                                                 <  2:217709/91‑1 (MQ=255)

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

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

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

                                             agaTATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCt                                               <  1:623642/91‑1 (MQ=255)

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

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

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

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

                                                tataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTa                                            >  1:1823260/1‑91 (MQ=255)

                                                tataCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTa                                            <  1:436442/91‑1 (MQ=255)

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

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

                                                  taCAACAATCATCAGATGACGAATAATACAAAATAAGACTTCCATATATGTAGGGGAGGCTTA‑CTTATATG‑CTAGAAAGTAATTCGTTACt                                          >  2:2112803/1‑91 (MQ=255)

                                                  taCAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACt                                          <  1:2055838/91‑1 (MQ=255)

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

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

                                                     aacaaTCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACtata                                       >  1:347435/1‑91 (MQ=255)

                                                     aacaaTCATCAGATGAAGAATAATAAAAAATAAGACTGCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACtata                                       <  2:3881612/91‑1 (MQ=255)

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

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

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

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

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

                                                            atcaGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATg                                >  2:2154615/1‑91 (MQ=255)

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

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

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

                                                                  tGAAGAATAATAGAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGt                          <  2:2218375/91‑1 (MQ=255)

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

                                                                       aataatAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGa                     >  1:2325035/1‑91 (MQ=255)

                                                                       aataatAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGa                     <  1:1834260/91‑1 (MQ=255)

                                                                       aataatAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGa                     <  1:598414/91‑1 (MQ=255)

                                                                         taataaAAAATAAGAATTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTACTGCTTTACTATATTCAATGATTAGTAATGACt                   >  1:3373495/1‑91 (MQ=255)

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

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

                                                                             aaaaaaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACt               >  1:502512/1‑91 (MQ=255)

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

                                                                               aaaaTAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACttt             >  2:3671853/1‑91 (MQ=255)

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

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

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

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

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

                                                                                     aGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  1:2279855/1‑91 (MQ=255)

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

                                                                                     aGACTCCCCTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAAtt       >  1:1854221/1‑91 (MQ=37)

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

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

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

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

                                                                                          tCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTc  <  2:2014549/91‑1 (MQ=255)

TGGACGTATTGAAGTAAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTC  >  USA300TCH1516_ALE/9620‑9800

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