BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	53,905	C→G	8.1%	S281R (AGC→AGG)	USA300HOU_RS00220 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	53,905	0	C	G	8.1%	471.9 / 26.4	137	S281R (AGC→AGG)	USA300HOU_RS00220	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (57/68); new base G (7/4); total (64/73)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 139 total aligned reads displayed.

TGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGTAT  >  USA300TCH1516_ALE/53815‑53995

  gagaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCt                                                                                          >  1:3613840/1‑91 (MQ=255)

   agaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGGCTa                                                                                         <  1:2915187/91‑1 (MQ=255)

   agaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGGCTa                                                                                         <  1:2915940/91‑1 (MQ=255)

   agaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTa                                                                                         >  1:2886556/1‑91 (MQ=255)

   agaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTa                                                                                         >  1:2886858/1‑91 (MQ=255)

   agaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTa                                                                                         >  1:2887140/1‑91 (MQ=255)

     aTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGGCTaaa                                                                                       <  1:599709/91‑1 (MQ=255)

     aTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTaaa                                                                                       <  1:2208616/91‑1 (MQ=255)

     aTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTaaa                                                                                       <  1:2601960/91‑1 (MQ=255)

     aTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTaaa                                                                                       <  1:2620969/91‑1 (MQ=255)

      tttAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAat                                                                                      <  1:3947141/91‑1 (MQ=255)

       ttAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTACatg                                                                                     >  1:3348298/1‑91 (MQ=255)

       ttAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatg                                                                                     >  1:3882624/1‑91 (MQ=255)

       ttAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatg                                                                                     <  1:1084611/91‑1 (MQ=255)

        tAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatga                                                                                    <  1:1203074/91‑1 (MQ=255)

        tAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatga                                                                                    <  1:32015/91‑1 (MQ=255)

         aGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGGCTAAatgat                                                                                   <  1:1686063/91‑1 (MQ=255)

         aGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatgat                                                                                   <  1:2570733/91‑1 (MQ=255)

         aGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatgat                                                                                   <  1:2761037/91‑1 (MQ=255)

           gTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATgg                                                                                 <  1:5977/91‑1 (MQ=255)

           gTAAGGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATgg                                                                                 <  1:3454593/91‑1 (MQ=255)

             aaTGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGtt                                                                               >  1:2750397/1‑91 (MQ=255)

                gAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGa                                                                            >  1:1275156/1‑91 (MQ=255)

                gAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGa                                                                            >  1:1790495/1‑91 (MQ=255)

                 aaGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGAc                                                                           <  1:415103/91‑1 (MQ=255)

                 aaGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGAc                                                                           <  1:415001/91‑1 (MQ=255)

                    aTTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTATGACTAg                                                                        >  1:2069971/1‑91 (MQ=255)

                     tttAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGa                                                                       <  1:439557/91‑1 (MQ=255)

                      ttAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGAc                                                                      <  1:1001625/91‑1 (MQ=255)

                      ttAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGAc                                                                      >  1:1000873/1‑91 (MQ=255)

                      ttAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGAc                                                                      <  1:15098/91‑1 (MQ=255)

                         aaCGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACaaa                                                                   >  1:2325448/1‑91 (MQ=255)

                         aaCGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACaaa                                                                   <  1:382578/91‑1 (MQ=255)

                         aaCGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACaaa                                                                   <  1:381118/91‑1 (MQ=255)

                            gtgtATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAAtt                                                                >  1:1111119/1‑91 (MQ=255)

                            gtgtATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAAtt                                                                >  1:1654986/1‑91 (MQ=255)

                             tgtATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTa                                                               >  1:401714/1‑91 (MQ=255)

                             tgtATAGTTTTTTTGAAAAAGTATTGATAAAATTAACCGATAACATTGAGGCATCATTTAGCCTAGATGATGGTTTGACTAGACACAATTa                                                               >  1:3853791/1‑91 (MQ=255)

                                 taGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCa                                                           <  1:1847251/91‑1 (MQ=255)

                                 taGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCa                                                           >  1:2727033/1‑91 (MQ=255)

                                        tttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAtt                                                    <  1:226547/91‑1 (MQ=255)

                                        tttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAtt                                                    >  1:916097/1‑91 (MQ=255)

                                        tttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAtt                                                    <  1:2818965/91‑1 (MQ=255)

                                        tttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAtt                                                    <  1:1977538/91‑1 (MQ=255)

                                         ttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAttt                                                   <  1:561957/91‑1 (MQ=255)

                                             aaaaGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTaaaa                                               >  1:457381/1‑91 (MQ=255)

                                             aaaaGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTaaaa                                               >  1:3410194/1‑91 (MQ=255)

                                             aaaaGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTaaaa                                               >  1:1690694/1‑91 (MQ=255)

                                               aaGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGt                                             <  1:766196/91‑1 (MQ=255)

                                                aGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTa                                            <  1:1157965/91‑1 (MQ=255)

                                                aGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTa                                            <  1:1157903/91‑1 (MQ=255)

                                                   aTTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCa                                         <  1:653256/91‑1 (MQ=255)

                                                    ttGATAAGATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCaa                                        <  1:2065337/91‑1 (MQ=255)

                                                    ttGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCaa                                        >  1:1024745/1‑91 (MQ=255)

                                                    ttGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCaa                                        >  1:1328087/1‑91 (MQ=255)

                                                     tGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAAt                                       <  1:2618090/91‑1 (MQ=255)

                                                     tGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAAt                                       >  1:1190193/1‑91 (MQ=255)

                                                     tGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAAt                                       >  1:3955810/1‑91 (MQ=255)

                                                      gATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAAtt                                      <  1:3337369/91‑1 (MQ=255)

                                                         aaaaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTcgc                                   <  1:2835344/91‑1 (MQ=255)

                                                         aaaaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTcgc                                   <  1:2522973/91‑1 (MQ=255)

                                                         aaaaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTcgc                                   <  1:3540254/91‑1 (MQ=255)

                                                          aaaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTcgcg                                  <  1:1376434/91‑1 (MQ=255)

                                                           aaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGa                                 <  1:3443991/91‑1 (MQ=255)

                                                           aaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGa                                 <  1:323804/91‑1 (MQ=255)

                                                             ttAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAg                               <  1:3193906/91‑1 (MQ=255)

                                                             ttAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAg                               <  1:3180591/91‑1 (MQ=255)

                                                                aCTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCAc                            >  1:1220799/1‑91 (MQ=255)

                                                                 cTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACt                           >  1:2529465/1‑91 (MQ=255)

                                                                   gATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTgg                         >  1:3770871/1‑91 (MQ=255)

                                                                   gATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTgg                         <  1:2431776/91‑1 (MQ=255)

                                                                   gATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTgg                         >  1:3526525/1‑91 (MQ=255)

                                                                   gATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTgg                         <  1:3909015/91‑1 (MQ=255)

                                                                   gATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTgg                         <  1:2431919/91‑1 (MQ=255)

                                                                    aTAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGt                        >  1:1398025/1‑91 (MQ=255)

                                                                     tAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGtt                       <  1:3586968/91‑1 (MQ=255)

                                                                      aaCATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAACCACTGGTTa                      >  1:3968937/1‑91 (MQ=255)

                                                                       aCATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTaa                     <  1:679073/91‑1 (MQ=255)

                                                                         aTTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATa                   <  1:824428/91‑1 (MQ=255)

                                                                         aTTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATa                   <  1:1507889/91‑1 (MQ=255)

                                                                         aTTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATa                   >  1:1096620/1‑91 (MQ=255)

                                                                         aTTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATa                   >  1:1375146/1‑91 (MQ=255)

                                                                            gAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACAt                >  1:3914391/1‑91 (MQ=255)

                                                                            gAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGAAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACAt                >  1:3296358/1‑91 (MQ=255)

                                                                             aGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACAtt               >  1:244339/1‑91 (MQ=255)

                                                                             aGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACAtt               >  1:3944001/1‑91 (MQ=255)

                                                                             aGGAATCATTTAGCCTAAATGATGGATTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACAtt               <  1:3877765/91‑1 (MQ=255)

                                                                              ggAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTa              >  1:4135676/1‑91 (MQ=255)

                                                                               gAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTaa             <  1:4008817/91‑1 (MQ=255)

                                                                               ccaTCATTTAGGCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTaa             >  1:718681/3‑91 (MQ=255)

                                                                                aaTCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAat            <  1:2699221/91‑1 (MQ=255)

                                                                                   cATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgca         <  1:3470250/91‑1 (MQ=255)

                                                                                   cATTTAGACTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATAGATTAatgca         <  1:2081708/91‑1 (MQ=255)

                                                                                     tttAGGCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgcatg       >  1:3878876/1‑91 (MQ=255)

                                                                                      ttAGGCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgcatgc      >  1:2496849/1‑91 (MQ=255)

                                                                                      ttAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgcatgc      <  1:3657563/91‑1 (MQ=255)

                                                                                      ttAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgcatgc      <  1:2003614/91‑1 (MQ=255)

                                                                                        aGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGt    <  1:872758/91‑1 (MQ=255)

                                                                                         gCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGta   >  1:2838970/1‑91 (MQ=255)

                                                                                          ccTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGtat  >  1:1792987/1‑91 (MQ=255)

TGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGTAT  >  USA300TCH1516_ALE/53815‑53995

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