BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	186,459	A→G	9.0%	Q345R (CAA→CGA)	USA300HOU_RS00855 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	186,459	0	A	G	9.0%	337.6 / 19.4	111	Q345R (CAA→CGA)	USA300HOU_RS00855	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (38/63); new base G (8/2); total (46/65)
Fisher's exact test for biased strand distribution p-value = 1.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 112 total aligned reads displayed.

CCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACGAGAT  >  USA300TCH1516_ALE/186371‑186549

ccGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAat                                                                                          <  1:203314/91‑1 (MQ=255)

ccGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAat                                                                                          <  1:3391956/91‑1 (MQ=255)

ccGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAat                                                                                          <  1:520909/91‑1 (MQ=255)

ccGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAat                                                                                          <  1:3667855/91‑1 (MQ=255)

   aGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatg                                                                                       <  1:1464302/91‑1 (MQ=255)

    gCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatga                                                                                      >  1:1778359/1‑91 (MQ=255)

    gCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatga                                                                                      >  1:4055088/1‑91 (MQ=255)

     cAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatgat                                                                                     <  1:2448685/91‑1 (MQ=255)

        gAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCg                                                                                  >  1:3112517/1‑91 (MQ=255)

         aaTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAATAATATGATGCgg                                                                                 <  1:3200101/91‑1 (MQ=255)

         aaTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCgg                                                                                 <  1:3832986/91‑1 (MQ=255)

         aaTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCgg                                                                                 >  1:3409518/1‑91 (MQ=255)

          aTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCggg                                                                                <  1:2414885/91‑1 (MQ=255)

          aTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCggg                                                                                <  1:2315365/91‑1 (MQ=255)

          aTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCggg                                                                                >  1:3957594/1‑91 (MQ=255)

          aTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCggg                                                                                >  1:2668795/1‑91 (MQ=255)

          aTTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCggg                                                                                <  1:3171104/91‑1 (MQ=255)

             aTGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTAt                                                                             >  1:3836855/1‑91 (MQ=255)

             aTGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTAt                                                                             <  1:753641/91‑1 (MQ=255)

              tGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATc                                                                            >  1:177190/1‑91 (MQ=255)

               ggCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATcc                                                                           <  1:1098433/91‑1 (MQ=255)

                gCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCa                                                                          <  1:2482578/91‑1 (MQ=255)

                gCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCa                                                                          <  1:2482085/91‑1 (MQ=255)

                gCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCa                                                                          <  1:2482015/91‑1 (MQ=255)

                gCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCa                                                                          <  1:906592/91‑1 (MQ=255)

                 cTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCaa                                                                         >  1:2554067/1‑91 (MQ=255)

                   aTCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAAt                                                                       <  1:508073/91‑1 (MQ=255)

                    tCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAAtt                                                                      <  1:1868014/91‑1 (MQ=255)

                      ttGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGa                                                                    <  1:549842/91‑1 (MQ=255)

                        gCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATa                                                                  <  1:1603205/91‑1 (MQ=255)

                         cAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATaa                                                                 <  1:4015979/91‑1 (MQ=255)

                         cAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATaa                                                                 <  1:2279778/91‑1 (MQ=255)

                           gCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATg                                                               >  1:807414/1‑91 (MQ=255)

                           gCTAAGAAAGCGATGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATg                                                               >  1:3178180/1‑91 (MQ=255)

                             tAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGc                                                             <  1:3529936/91‑1 (MQ=255)

                              aaGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCg                                                            <  1:1525988/91‑1 (MQ=255)

                              aaGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCg                                                            <  1:3441039/91‑1 (MQ=255)

                               aGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGa                                                           >  1:1904795/1‑91 (MQ=255)

                                  aaGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGc                                                        >  1:1225439/1‑91 (MQ=255)

                                    gCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATGATGGCGATGcac                                                      <  1:3184907/91‑1 (MQ=255)

                                    gCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGcac                                                      <  1:1265071/91‑1 (MQ=255)

                                      gTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATTTGATGCGGGTATCCAAATTGATAATGGCGATGCACAt                                                    <  1:3511111/91‑1 (MQ=255)

                                      gTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACAt                                                    <  1:636316/91‑1 (MQ=255)

                                        tGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGc                                                  <  1:185570/91‑1 (MQ=255)

                                         ggTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGcc                                                 <  1:368625/91‑1 (MQ=255)

                                         ggTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGcc                                                 >  1:1124758/1‑91 (MQ=255)

                                         ggTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCAATGCACATGcc                                                 <  1:1191762/91‑1 (MQ=255)

                                          gTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCa                                                >  1:718940/1‑91 (MQ=255)

                                            aCTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCAtt                                              <  1:245306/91‑1 (MQ=255)

                                            aCTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCAtt                                              <  1:2513826/91‑1 (MQ=255)

                                            aCTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCAtt                                              <  1:2513747/91‑1 (MQ=255)

                                                  catGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAAc                                        <  1:4038958/91‑1 (MQ=255)

                                                  catGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAAc                                        <  1:404163/91‑1 (MQ=255)

                                                       aTAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtcat                                   >  1:2898605/1‑91 (MQ=255)

                                                       aTAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtcat                                   >  1:1168122/1‑91 (MQ=255)

                                                        tAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtcatc                                  <  1:3198855/91‑1 (MQ=255)

                                                         aaGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtcatca                                 >  1:1683475/1‑91 (MQ=255)

                                                          aGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCaa                                >  1:2117796/1‑91 (MQ=255)

                                                          aGCAATCTGATATCGCTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCaa                                >  1:4100660/1‑91 (MQ=255)

                                                            cAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAaca                              <  1:1369597/91‑1 (MQ=255)

                                                             aaTCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAacac                             <  1:2872228/91‑1 (MQ=255)

                                                             aaTCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAacac                             >  1:760569/1‑91 (MQ=255)

                                                                 tGATATCGTTCAACATGTTGAACAATATGATGCTGGTATCCAAATTGATAATGGCGATGGACATGCCATTTATAACTTCATCAAcactcac                         <  1:3686230/91‑1 (MQ=255)

                                                                    tatCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCg                      >  1:709652/1‑91 (MQ=255)

                                                                     atCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCgc                     >  1:1622531/1‑90 (MQ=255)

                                                                     atCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGa                     >  1:2675832/1‑91 (MQ=255)

                                                                      tCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAg                    <  1:808631/91‑1 (MQ=255)

                                                                       tgTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGt                   >  1:186160/2‑91 (MQ=255)

                                                                        gTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTTACTCGAGTa                  >  1:3980758/1‑91 (MQ=255)

                                                                        gTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTa                  >  1:3449508/1‑91 (MQ=255)

                                                                        gTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTa                  >  1:4006552/1‑91 (MQ=255)

                                                                        gTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTa                  >  1:3764401/1‑91 (MQ=255)

                                                                        gTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTa                  <  1:3110043/91‑1 (MQ=255)

                                                                         ttCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTaa                 <  1:1334701/91‑1 (MQ=255)

                                                                         ttCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTaa                 <  1:1434297/91‑1 (MQ=255)

                                                                         ttCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTaa                 <  1:3110642/91‑1 (MQ=255)

                                                                          tCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                >  1:957180/1‑91 (MQ=255)

                                                                          tCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                <  1:686728/91‑1 (MQ=255)

                                                                          tCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                <  1:4130218/91‑1 (MQ=255)

                                                                          tCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                >  1:2130855/1‑91 (MQ=255)

                                                                          tCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                >  1:2131341/1‑91 (MQ=255)

                                                                           cAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAgg               >  1:1637928/1‑91 (MQ=255)

                                                                           cAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAgg               >  1:1939089/1‑91 (MQ=255)

                                                                           cAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAgg               <  1:2267862/91‑1 (MQ=255)

                                                                           cAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAgg               >  1:3720513/1‑91 (MQ=255)

                                                                            aaCATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGa              >  1:1272755/1‑91 (MQ=255)

                                                                             aCATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGaa             >  1:1529821/1‑91 (MQ=255)

                                                                              cATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAAt            <  1:470258/91‑1 (MQ=255)

                                                                              cATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAAt            <  1:470173/91‑1 (MQ=255)

                                                                                 gTTTAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGc         <  1:3216569/91‑1 (MQ=255)

                                                                                 gTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGc         >  1:1894716/1‑91 (MQ=255)

                                                                                  ttGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCa        >  1:2619971/1‑91 (MQ=255)

                                                                                  ttGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCa        >  1:2829413/1‑91 (MQ=255)

                                                                                   tGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCAc       <  1:2762173/91‑1 (MQ=255)

                                                                                     aaCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACga     <  1:3635410/91‑1 (MQ=255)

                                                                                      aCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACgag    >  1:1995762/1‑91 (MQ=255)

                                                                                      aCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACgag    >  1:1993464/1‑91 (MQ=255)

                                                                                      aCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACgag    <  1:3381795/91‑1 (MQ=255)

                                                                                       cAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACgaga   <  1:3338500/91‑1 (MQ=255)

                                                                                        aaTATGATGCGGGTATCCAAATTGATAAGGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACGAGAt  <  1:3926257/91‑1 (MQ=255)

CCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACGAGAT  >  USA300TCH1516_ALE/186371‑186549

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