BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	9,725	G→A	17.1%	intergenic (+30/+57)	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,725	0	G	A	17.1%	428.5 / 59.9	175	intergenic (+30/+57)	gyrA/nnrD	DNA gyrase subunit A/ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase
Reads supporting (aligned to +/- strand): ref base G (66/77); new base A (19/11); total (87/88)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

AAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGTATT  >  USA300TCH1516_ALE/9635‑9815

aaGAAAAGATTGCAAGGATCGTGTTGAAGAAGATATACAACAATCAAAAGATGAAAAAAAAAAAAAAATAAGACTTCCCTATATGTAgggg                                                                                              <  2:2671438/91‑1 (MQ=255)

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

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

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

  gACAAGAGTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCATATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAg                                                                                            <  2:1489241/91‑1 (MQ=37)

   aCAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGt                                                                                           <  1:2771188/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)

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

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

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

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

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

      aGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCtt                                                                                        <  2:445495/91‑1 (MQ=37)

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

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

       gATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTCCCCTACATATAGGGAAGTCTTa                                                                                       <  1:3052676/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)

              tGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATCTGTAGGGGAGTCTTA‑TTTTTAt                                                                               >  2:2603028/1‑91 (MQ=255)

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

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

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

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

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

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

                      gtTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGaaa                                                                      >  1:3761326/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:129524/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                  aTGAAGAATAATGAAGAATGAGACTTCCCTATGTGTAGGGGAGTCTTA‑TTTGTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAg                                          <  2:1118328/91‑1 (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)

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

                                                             taaAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAAc                               <  1:2811162/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‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTc                           >  1:508142/1‑91 (MQ=21)

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

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

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

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

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

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

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

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

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

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

                                                                          ttCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGttt                  <  1:3820514/91‑1 (MQ=255)

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

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

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

                                                                            cccTACATATAGGGAAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCa                >  1:2191038/1‑91 (MQ=38)

                                                                             ccTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCAt               <  2:3493138/91‑1 (MQ=255)

                                                                                atatGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGc            >  1:3107716/1‑91 (MQ=255)

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

                                                                                     tAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAg       <  2:1889923/91‑1 (MQ=255)

                                                                                      aGGGGAGTCTTA‑TTTTTATG‑CTAGAGAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAgg      <  2:2133549/91‑1 (MQ=255)

                                                                                      aGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATTCTTTACTATATTCAATGATTAGTAATGACTAACGTTCTAATTGTTTCATTGCGTAAgg      <  2:2642720/91‑1 (MQ=255)

                                                                                      aGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAgg      >  2:3145204/1‑91 (MQ=255)

                                                                                      aGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAgg      >  2:1490549/1‑91 (MQ=255)

                                                                                      aGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAgg      >  2:2474360/1‑91 (MQ=255)

                                                                                      aGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATGAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAgg      <  2:377767/91‑1 (MQ=255)

                                                                                       ggggAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGt     >  2:951936/1‑91 (MQ=255)

                                                                                        gggAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGTa    <  1:615250/91‑1 (MQ=255)

                                                                                         ggAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGTAt   >  2:2867473/1‑91 (MQ=255)

                                                                                         ggAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGTAt   >  2:1813198/1‑91 (MQ=255)

                                                                                         ggAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGTAt   >  2:1812502/1‑91 (MQ=255)

AAGACAAGATTTCATGGATCGTGTTGAAGAAGATATACAACAATCATCAGATGAAGAATAATAAAAAATAAGACTTCCCTATATGTAGGGGAGTCTTA‑TTTTTATG‑CTAGAAAGTAATGCTTTACTATATTCAATGATTAGTAATGACTAACTTTCTAATTGTTTCATTGCGTAAGGTATT  >  USA300TCH1516_ALE/9635‑9815

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