BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	341,450	G→C	17.3%	D134H (GAT→CAT)	USA300HOU_RS01525 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	341,450	0	G	C	17.3%	266.3 / 48.4	98	D134H (GAT→CAT)	USA300HOU_RS01525	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (41/40); new base C (9/8); total (50/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGGTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGAA  >  USA300TCH1516_ALE/341360‑341539

gggTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATg                                                                                           >  1:756066/1‑91 (MQ=255)

  gTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGAt                                                                                         >  1:4008896/1‑91 (MQ=255)

   tttCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGAtt                                                                                        <  1:2782456/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                      <  1:210282/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                      <  1:3930199/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                      <  1:1506256/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                      <  1:2618887/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                      <  1:3092495/91‑1 (MQ=255)

      cAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgg                                                                                      <  1:3082584/91‑1 (MQ=255)

       aaaTACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATTgga                                                                                     <  1:438686/91‑2 (MQ=255)

           aCGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACttt                                                                                >  1:3992300/1‑91 (MQ=255)

            cGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATTATATCCTTTATTGATGAGACAAGAAAAGTCCCAATCATT‑GGACtttt                                                                               <  1:4136668/91‑2 (MQ=255)

                ttATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACa                                                                           <  1:2873369/91‑1 (MQ=255)

                ttATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACa                                                                           <  1:2873627/91‑1 (MQ=255)

                  aTCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACaaa                                                                         <  1:1780380/91‑1 (MQ=255)

                   tCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAAt                                                                        >  1:1527259/1‑91 (MQ=255)

                    cAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAAtt                                                                       <  1:3724307/91‑1 (MQ=255)

                    cAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAAtt                                                                       <  1:3020859/91‑1 (MQ=255)

                    cAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAAtt                                                                       <  1:3724426/91‑1 (MQ=255)

                    cAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAAtt                                                                       >  1:3825827/1‑91 (MQ=255)

                      aTAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGa                                                                     >  1:606507/1‑91 (MQ=255)

                        aaaGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACg                                                                   >  1:865039/1‑91 (MQ=255)

                         aaGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAATAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGa                                                                  <  1:3495081/91‑1 (MQ=255)

                         aaGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGa                                                                  <  1:1157208/91‑1 (MQ=255)

                          aGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGAt                                                                 <  1:2752154/91‑1 (MQ=255)

                           gACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATg                                                                <  1:883650/91‑1 (MQ=255)

                           gACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATg                                                                <  1:883608/91‑1 (MQ=255)

                           gACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATg                                                                <  1:665105/91‑1 (MQ=255)

                             cATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGaa                                                              >  1:1162807/1‑91 (MQ=255)

                              aTTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGaaa                                                             >  1:2218923/1‑91 (MQ=255)

                              aTTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGaaa                                                             >  1:2917125/1‑91 (MQ=255)

                              aTTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGaaa                                                             >  1:1590430/1‑91 (MQ=255)

                               ttGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAAc                                                            >  1:1575583/1‑91 (MQ=255)

                                tGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACa                                                           >  1:1539661/1‑91 (MQ=255)

                                 gCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAt                                                          >  1:1488816/1‑91 (MQ=255)

                                  cATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAtt                                                         >  1:1560139/1‑91 (MQ=255)

                                  cATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACAtt                                                         >  1:1560243/1‑91 (MQ=255)

                                   aTGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTg                                                        <  1:2792639/91‑1 (MQ=255)

                                     ggtggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGga                                                      <  1:3956754/91‑1 (MQ=255)

                                       tggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGtaga                                                    >  1:3350788/1‑91 (MQ=255)

                                       tggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGcaga                                                    >  1:3350787/1‑91 (MQ=255)

                                       tggtCGAATGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGGCGATGAAACATTGgaga                                                    <  1:2062760/91‑1 (MQ=255)

                                        ggtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGaa                                                   <  1:3468931/91‑1 (MQ=255)

                                         gtCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAAc                                                  <  1:1529490/91‑1 (MQ=255)

                                          tCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACt                                                 <  1:3659658/91‑1 (MQ=255)

                                           cGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACtt                                                <  1:2854428/91‑1 (MQ=255)

                                            gATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTa                                               >  1:955846/1‑91 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              <  1:1217937/91‑1 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              <  1:3839184/91‑1 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              >  1:3934670/1‑91 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              >  1:3432384/1‑91 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              >  1:3934683/1‑91 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              <  1:396778/91‑1 (MQ=255)

                                             aTTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTaa                                              <  1:3127914/91‑1 (MQ=255)

                                                gTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaata                                           >  1:3404849/1‑91 (MQ=255)

                                                  aCAATCATATCCTTTATTGATGGGACAAGAAAAGT‑CCAGTGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaataat                                         >  1:2754245/1‑91 (MQ=255)

                                                  aCAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAaataat                                         <  1:3969269/91‑1 (MQ=255)

                                                   cAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATg                                        >  1:656368/1‑91 (MQ=255)

                                                      tCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGaaa                                     >  1:3783809/1‑91 (MQ=255)

                                                       cATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAg                                    <  1:3880707/91‑1 (MQ=255)

                                                       cATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAg                                    >  1:3751302/1‑91 (MQ=255)

                                                          atCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTAt                                 >  1:3941392/1‑91 (MQ=255)

                                                          atCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTAt                                 >  1:2012498/1‑91 (MQ=255)

                                                           tCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATc                                <  1:1137567/91‑1 (MQ=255)

                                                             cTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGa                              <  1:221834/91‑1 (MQ=255)

                                                                 aTTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACg                          <  1:1663675/91‑1 (MQ=255)

                                                                   tgatgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGAt                        >  1:118061/1‑91 (MQ=255)

                                                                     atgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTg                      <  1:1358100/91‑1 (MQ=255)

                                                                     atgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTg                      <  1:3163907/91‑1 (MQ=255)

                                                                      tgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTgg                     <  1:3625082/91‑1 (MQ=255)

                                                                      tgaGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTgg                     >  1:1856349/1‑91 (MQ=255)

                                                                       gagaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTggg                    <  1:2706794/91‑1 (MQ=255)

                                                                        agaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGc                   >  1:227523/1‑91 (MQ=255)

                                                                        agaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGc                   <  1:739706/91‑1 (MQ=255)

                                                                        agaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGc                   >  1:4065719/1‑91 (MQ=255)

                                                                         gaCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCg                  >  1:2505787/1‑91 (MQ=255)

                                                                          aCAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCgt                 >  1:3790414/1‑91 (MQ=255)

                                                                             tttaaaGTCCCAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtt              >  1:1291176/4‑91 (MQ=38)

                                                                             aGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtt              >  1:1148761/1‑91 (MQ=255)

                                                                               aaaaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGtttt            <  1:475960/91‑1 (MQ=255)

                                                                                aaaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTa           >  1:455340/1‑91 (MQ=255)

                                                                                 aaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAg          <  1:2393797/91‑1 (MQ=255)

                                                                                 aaGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAg          <  1:1143425/91‑1 (MQ=255)

                                                                                  agtcccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGa         >  1:1910068/5‑91 (MQ=255)

                                                                                  aGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGa         >  1:815670/1‑91 (MQ=255)

                                                                                   gT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGaa        >  1:3094425/1‑91 (MQ=255)

                                                                                   gT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGaa        >  1:2836140/1‑91 (MQ=255)

                                                                                    tcccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGaaa       >  1:4036713/3‑91 (MQ=255)

                                                                                    t‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGaaa       <  1:520165/91‑1 (MQ=255)

                                                                                     cccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAt      >  1:1547098/2‑91 (MQ=255)

                                                                                      ccAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAt      >  1:3507175/1‑91 (MQ=255)

                                                                                      ccAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAt      >  1:889197/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt     >  1:3008931/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt     >  1:594799/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt     >  1:2805764/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt     >  1:1415887/1‑91 (MQ=255)

                                                                                      ccAATCATT‑GGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAAtt     >  1:3220964/1‑91 (MQ=255)

                                                                                          tGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaa  <  1:1706335/91‑1 (MQ=255)

                                                                                          tGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGaa  <  1:1706244/91‑1 (MQ=255)

GGGTTTCAAATACGAATTATCAATAAAGACATTGCATGGTGGTCGATTGTACAATCATATCCTTTATTGATGAGACAAGAAAAGT‑CCAATGATTGGGACTTTAAACAAATTGACGATGAAACATTGGAGAACTTAAATAATGAAAGTATCGATACGATTGGGCGTGTTTTAGAAATTGAA  >  USA300TCH1516_ALE/341360‑341539

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