BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	186,459	A→G	8.8%	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	8.8%	641.8 / 41.2	217	Q345R (CAA→CGA)	USA300HOU_RS00855	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (97/101); new base G (16/3); total (113/104)
Fisher's exact test for biased strand distribution p-value = 3.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 224 total aligned reads displayed.

TGCCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACGAGAT  >  USA300TCH1516_ALE/186369‑186549

tGCCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACa                                                                                            >  1:2263129/1‑91 (MQ=255)

tGCCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACa                                                                                            <  2:2693830/91‑1 (MQ=255)

 gCCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACaa                                                                                           >  1:929295/1‑91 (MQ=255)

   cGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATTGTTCAACATGTTGAACGata                                                                                         <  1:2722522/91‑1 (MQ=255)

   cGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATTGTTCAACATGTTGAACGata                                                                                         <  1:2608450/91‑1 (MQ=255)

    gAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAgtat                                                                                        >  1:3129067/1‑91 (MQ=255)

    gAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAatat                                                                                        <  2:728287/91‑1 (MQ=255)

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

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

      gCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatga                                                                                      <  2:334305/91‑1 (MQ=255)

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

        aaGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatgatg                                                                                    >  2:2368378/1‑91 (MQ=255)

        aaGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATatgatg                                                                                    <  1:2079690/91‑1 (MQ=255)

         aGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGc                                                                                   <  1:404985/91‑1 (MQ=255)

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

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

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

             ttATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGt                                                                               >  1:11045/1‑91 (MQ=255)

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

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

                tGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATc                                                                            <  1:2637529/91‑1 (MQ=255)

                 ggCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATcc                                                                           >  2:2265365/1‑91 (MQ=255)

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

                    tATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCaaa                                                                        <  1:832632/91‑1 (MQ=255)

                    tATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCaaa                                                                        <  2:3710149/91‑1 (MQ=255)

                      tCTTGCAGTTAAGAAAGCGTTGGTTCTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAAtt                                                                      <  2:2558087/91‑1 (MQ=255)

                         tGCAGCAAAGAAACCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAAAAATATGATGCGGGTATCCAAATTGAt                                                                   <  1:629184/91‑1 (MQ=255)

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

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

                           cAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATaa                                                                 >  2:2595924/1‑91 (MQ=255)

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

                              cTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATACTgg                                                              >  2:2661019/1‑91 (MQ=255)

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

                                   aaaGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATg                                                         <  2:1811820/91‑1 (MQ=255)

                                   aaaGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATg                                                         >  2:644954/1‑91 (MQ=255)

                                    aaGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGc                                                        <  1:2375275/91‑1 (MQ=255)

                                    aaGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGc                                                        <  1:1530061/91‑1 (MQ=255)

                                     aGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGca                                                       >  2:1092297/1‑91 (MQ=255)

                                     aGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGca                                                       <  1:1041422/91‑1 (MQ=255)

                                       cGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGcaca                                                     <  2:510612/91‑1 (MQ=255)

                                          tGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGc                                                  >  1:2818306/1‑91 (MQ=255)

                                            gTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCa                                                <  1:3197347/91‑1 (MQ=255)

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

                                                tcatcaTGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTta                                            <  1:3779925/91‑1 (MQ=255)

                                                tcatcaTGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTta                                            <  1:858366/91‑1 (MQ=255)

                                                  atcatGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTtata                                          <  1:2510547/91‑1 (MQ=255)

                                                     atGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACt                                       >  1:638909/1‑91 (MQ=255)

                                                     atGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACt                                       >  1:1305644/1‑91 (MQ=255)

                                                      tGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACtt                                      >  1:384387/1‑91 (MQ=255)

                                                      tGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACtt                                      >  1:1362002/1‑91 (MQ=255)

                                                       ggATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtc                                     <  2:3166225/91‑1 (MQ=255)

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

                                                          tAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtcatc                                  >  2:2923081/1‑91 (MQ=255)

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

                                                           aaGCAATCTGATATCGTTCAACATGTTGAACAATATGATGAGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTtcatca                                 >  2:3601147/1‑91 (MQ=255)

                                                            aGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCaa                                <  1:438766/91‑1 (MQ=255)

                                                             gCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAac                               <  1:56578/91‑1 (MQ=255)

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

                                                                aTCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAAGTTCATCAAcact                            <  1:2998600/91‑1 (MQ=255)

                                                                   tGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAAcactcac                         >  2:877346/1‑91 (MQ=255)

                                                                   tGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAAcactcac                         <  1:1399372/91‑1 (MQ=255)

                                                                   tGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAAcactcac                         <  1:1738516/91‑1 (MQ=255)

                                                                   tGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAAcactcac                         <  1:211592/91‑1 (MQ=255)

                                                                    catatTGTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAAcactcact                        >  1:1151865/2‑91 (MQ=255)

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

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

                                                                      tatCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCg                      <  2:2856317/91‑1 (MQ=255)

                                                                        ttgTTCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAg                    >  2:1252930/3‑91 (MQ=255)

                                                                        tCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAg                    >  1:2006728/1‑91 (MQ=255)

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

                                                                         cGTTCAAAATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCCAGt                   <  2:3445606/91‑1 (MQ=255)

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

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

                                                                           ttCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTaa                 >  2:2177234/1‑91 (MQ=255)

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

                                                                            tCAACATTTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                <  1:1936353/91‑1 (MQ=255)

                                                                            tCAACATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAg                >  1:70358/1‑91 (MQ=255)

                                                                             cAACATGTTGAACAATATGATGCGGGTATCCAAATTGATGATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAgg               >  1:2352250/1‑91 (MQ=255)

                                                                              aaCTTGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGa              <  1:3003486/91‑1 (MQ=255)

                                                                              aaCATGTTGAACGATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGa              >  2:1551030/1‑91 (MQ=255)

                                                                              aaCATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGa              <  2:577343/91‑1 (MQ=255)

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

                                                                                cATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAAt            >  2:1454147/1‑91 (MQ=255)

                                                                                cATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAAt            >  1:1804139/1‑91 (MQ=255)

                                                                                 aTGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAAtt           <  1:2368378/91‑1 (MQ=255)

                                                                                  tGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTg          >  1:3651307/1‑91 (MQ=255)

                                                                                  tGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTg          >  1:527735/1‑91 (MQ=255)

                                                                                  tGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTg          >  2:1086138/1‑91 (MQ=255)

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

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

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

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

                                                                                       aaCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACga     >  2:1668509/1‑91 (MQ=255)

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

                                                                                       aaCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACga     >  2:805754/1‑91 (MQ=255)

                                                                                       aaCAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATCTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACga     >  2:3399915/1‑91 (MQ=255)

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

                                                                                         cAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACgaga   >  2:1888200/1‑91 (MQ=255)

                                                                                          aaTATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACGAGAt  >  1:3741997/1‑91 (MQ=255)

                                                                                          aaTATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTACATCAACACTCACTCGAGTAAGGAATTGCACGAGAt  >  2:3836672/1‑91 (MQ=255)

TGCCGAGCAAGAATTATGGCTATCTTGCAGCTAAGAAAGCGTTGGTACTCATCATGGATAAGCAATCTGATATCGTTCAACATGTTGAACAATATGATGCGGGTATCCAAATTGATAATGGCGATGCACATGCCATTTATAACTTCATCAACACTCACTCGAGTAAGGAATTGCACGAGAT  >  USA300TCH1516_ALE/186369‑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: `‑`