BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,538,743	A→T	6.0%	intergenic (‑257/+70)	gpmA_2 ← / ← fieF	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase/Ferrous‑iron efflux pump FieF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,538,743	0	A	T	6.0%	805.4 / 27.3	238	intergenic (‑257/+70)	gpmA_2/fieF	2,3‑bisphosphoglycerate‑dependent phosphoglycerate mutase/Ferrous‑iron efflux pump FieF
Reads supporting (aligned to +/- strand): ref base A (89/134); new base T (4/11); total (93/145)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

TCATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCT  >  USA300TCH1516_ALE/2538654‑2538831

tCATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTaa                                                                                         <  2:2253373/91‑1 (MQ=255)

tCATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTaa                                                                                         >  1:2094860/1‑91 (MQ=255)

 cATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTaaa                                                                                        <  1:1104873/91‑1 (MQ=255)

 cATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTaaa                                                                                        <  1:868046/91‑1 (MQ=255)

  aTAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAAc                                                                                       <  2:3682504/91‑1 (MQ=255)

    aaaTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCt                                                                                     <  2:3210220/91‑1 (MQ=255)

     aaTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTa                                                                                    <  2:1514765/91‑1 (MQ=255)

      aTGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaa                                                                                   >  2:3844693/1‑91 (MQ=255)

       tGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaa                                                                                  <  2:1245405/91‑1 (MQ=255)

        taaTCAATAGCCTCCTAATAACGATTCATGTGATTCTCCTTTGTGCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTaaaa                                                                                 <  2:1282495/90‑1 (MQ=255)

        gAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaaaa                                                                                 >  1:1479061/1‑91 (MQ=255)

         aaTCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTaaaaa                                                                                >  1:2113141/1‑91 (MQ=255)

           tCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTAAAAATa                                                                              <  1:2258863/91‑1 (MQ=255)

           tCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATa                                                                              >  2:547834/1‑91 (MQ=255)

           tCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATa                                                                              >  1:314472/1‑91 (MQ=255)

            cAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTAAAAATat                                                                             <  1:1682631/91‑2 (MQ=255)

             aaTAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTATCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGa                                                                            <  2:829154/91‑1 (MQ=255)

               tAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTTTACCTAAAAATatttt                                                                          <  1:1593107/91‑5 (MQ=38)

               tAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATg                                                                          <  1:1798219/91‑1 (MQ=255)

                aaCCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGt                                                                         <  1:1350596/91‑1 (MQ=255)

                   caccTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGtaat                                                                      <  1:451111/91‑1 (MQ=255)

                    accTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGtaata                                                                     >  2:3037279/1‑91 (MQ=255)

                    accTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGtaata                                                                     <  1:839172/91‑1 (MQ=255)

                      cTAATAACAAATCATATTATACACGTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATaaa                                                                   <  2:3618529/91‑1 (MQ=255)

                       taataaCAAATCATATTATACACCTTTGTCCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATATAATaaaa                                                                  >  2:3380589/1‑91 (MQ=255)

                        aataaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATaaaaa                                                                 <  1:3085102/91‑1 (MQ=255)

                         ataaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAAc                                                                <  2:963407/91‑1 (MQ=255)

                         ataaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAAc                                                                <  1:2924306/91‑1 (MQ=255)

                          taaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACg                                                               <  1:940899/91‑1 (MQ=255)

                           aaCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGc                                                              <  2:1583715/91‑1 (MQ=255)

                            aCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcc                                                             <  2:361051/91‑1 (MQ=255)

                            aCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcc                                                             >  1:1628738/1‑91 (MQ=255)

                            aCAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGAAGTAATAAAAACGcc                                                             <  2:691087/91‑1 (MQ=255)

                             cAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcct                                                            >  2:484550/1‑91 (MQ=255)

                              aaaTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctc                                                           >  2:3865583/1‑91 (MQ=255)

                                aTCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGcctcct                                                         <  1:323029/91‑1 (MQ=255)

                                 tCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCctcctc                                                        >  2:1908447/1‑91 (MQ=255)

                                 tCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCctcctc                                                        >  1:1708914/1‑91 (MQ=255)

                                 tCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCctcctc                                                        <  1:2084976/91‑1 (MQ=255)

                                    tattatACATCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGa                                                     <  1:1254570/91‑1 (MQ=255)

                                    tattatACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGa                                                     >  2:1000357/1‑91 (MQ=255)

                                     attatACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGat                                                    >  1:595357/1‑91 (MQ=255)

                                     attatACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGat                                                    >  1:522389/1‑91 (MQ=255)

                                      ttatACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGata                                                   <  1:1004686/91‑1 (MQ=255)

                                       tataCACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGatat                                                  <  2:1300279/91‑1 (MQ=255)

                                       tataCACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGatat                                                  >  2:2480376/1‑91 (MQ=255)

                                         taCACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATAttt                                                >  2:2654423/1‑91 (MQ=255)

                                           cacCTTTGTTCGCTATTTTACTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGAtatataa                                              <  2:2603159/91‑2 (MQ=255)

                                            acCTTTGTTCGCTATTTTTATAAAGATTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTtata                                             <  2:2245466/91‑1 (MQ=255)

                                                ttGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATc                                         <  2:1478028/91‑1 (MQ=255)

                                                ttGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATc                                         <  2:920150/91‑1 (MQ=255)

                                                 tGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCt                                        <  1:3543415/91‑1 (MQ=255)

                                                   ttCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTAt                                      <  2:3713303/91‑1 (MQ=255)

                                                    tCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATg                                     <  1:2760502/91‑1 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAg                                  >  2:2900262/1‑91 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAg                                  >  1:2632937/1‑91 (MQ=255)

                                                       cTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAg                                  >  2:1184951/1‑91 (MQ=255)

                                                        tATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGa                                 >  1:2690182/1‑91 (MQ=255)

                                                          tttttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACg                               >  2:1017765/1‑91 (MQ=255)

                                                          tttttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACg                               <  1:2782030/91‑1 (MQ=255)

                                                          tttttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACg                               >  2:1668141/1‑91 (MQ=255)

                                                           ttttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGt                              <  1:2735345/91‑1 (MQ=255)

                                                            tttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGtt                             >  1:3546546/1‑91 (MQ=255)

                                                            tttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGtt                             <  1:154299/91‑1 (MQ=255)

                                                            tttCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGtt                             >  2:2924456/1‑91 (MQ=255)

                                                              tCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTa                           >  1:2865224/1‑91 (MQ=255)

                                                               cTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTaa                          <  2:2426183/91‑1 (MQ=255)

                                                               cTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTaa                          <  1:3840674/91‑1 (MQ=255)

                                                                tAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTaaa                         >  1:649884/1‑91 (MQ=255)

                                                                 aaGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAAt                        >  2:1901037/1‑91 (MQ=255)

                                                                  aGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATa                       >  2:1647197/1‑91 (MQ=255)

                                                                   ggTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTta                                               <  2:307289/66‑1 (MQ=255)

                                                                   ggTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTta                                               >  1:307289/1‑66 (MQ=255)

                                                                    ttttAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACa                     >  2:111331/2‑91 (MQ=255)

                                                                     tttAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACAt                    <  1:218363/91‑1 (MQ=255)

                                                                      ttAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACAtt                   <  2:2967865/91‑1 (MQ=255)

                                                                       tAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTaaa                         <  2:1267068/84‑1 (MQ=255)

                                                                        aaaaaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATtat                 <  2:2172522/91‑1 (MQ=255)

                                                                          aaaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAg               >  2:3572060/1‑91 (MQ=255)

                                                                          aaaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAg               <  1:1507123/91‑1 (MQ=255)

                                                                           aaaTATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGa              <  1:523936/91‑1 (MQ=255)

                                                                             atatTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATg            >  2:2734875/1‑91 (MQ=255)

                                                                              tatTTTTCGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATgg           <  1:1945216/91‑1 (MQ=255)

                                                                               atTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGt          >  1:2946190/1‑91 (MQ=255)

                                                                               atTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGt          <  1:41343/91‑1 (MQ=255)

                                                                                tttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTc         <  1:484550/91‑1 (MQ=255)

                                                                                tttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTc         <  1:2019865/91‑1 (MQ=255)

                                                                                 ttttAGGTTTACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCt        >  1:2603177/1‑91 (MQ=255)

                                                                                 ttttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCt        >  1:777084/1‑91 (MQ=255)

                                                                                  tttAGGTTTACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTg       >  2:3461886/1‑91 (MQ=255)

                                                                                  tttAGGTTTACCTAAAAATAGATGTAATAAAAACGCCACCTCAGATATTTATATATCTATGAAGATGTTTAAATACATTATAGATCGTCTg       >  1:3616079/1‑91 (MQ=255)

                                                                                   ttAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTgg      <  2:2135370/91‑1 (MQ=255)

                                                                                     aGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGtt    >  1:1200872/1‑91 (MQ=255)

                                                                                     aGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGtt    <  2:277361/91‑1 (MQ=255)

                                                                                     aGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGtt    >  2:25484/1‑91 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc   >  1:708935/1‑91 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc   <  2:104302/91‑1 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc   <  1:3004843/91‑1 (MQ=255)

                                                                                      ggTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTc   <  1:2480376/91‑1 (MQ=255)

                                                                                       gTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCt  <  1:3459933/91‑1 (MQ=255)

TCATAAATGAATCAATAACCACCTAATAACAAATCATATTATACACCTTTGTTCGCTATTTTTCTAAGGTTTAAAAAATATTTTTAGGTAAACCTAAAAATAGATGTAATAAAAACGCCTCCTCAGATATTTATATATCTATGAAGACGTTTAAATACATTATAGATGGTCTGGTTCT  >  USA300TCH1516_ALE/2538654‑2538831

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