BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,512	A→G	32.1%	S194S (TCA→TCG)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,512	0	A	G	32.1%	711.1 / inf	466	S194S (TCA→TCG)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand): ref base A (144/172); new base G (52/98); total (196/270)
Fisher's exact test for biased strand distribution p-value = 2.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 475 total aligned reads displayed.

TAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTG  >  USA300TCH1516_ALE/111422‑111602

tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCa                                                                                            <  2:4862096/91‑1 (MQ=255)

 aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAg                                                                                           >  1:2435180/1‑91 (MQ=255)

      ttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCtaaa                                                                                      <  1:5582611/91‑4 (MQ=255)

       tttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAat                                                                                     >  2:3896610/1‑91 (MQ=255)

        ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                    <  1:2841709/91‑1 (MQ=255)

        ttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAata                                                                                    <  1:344304/91‑1 (MQ=255)

         tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAatat                                                                                   >  1:519792/1‑91 (MQ=255)

         tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAatat                                                                                   >  2:2443886/1‑91 (MQ=255)

          gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATc                                                                                  >  2:3805008/1‑91 (MQ=255)

          gcgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATc                                                                                  >  2:3499354/1‑91 (MQ=255)

           cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                 <  2:4519357/91‑1 (MQ=255)

           cgcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                 >  1:1359526/1‑91 (MQ=255)

           cgcAATATAGCAATTTTAAAGATTTAATGAATTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                 <  2:1964693/91‑1 (MQ=255)

            gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                <  2:2741902/91‑1 (MQ=255)

            gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                <  1:2883531/91‑1 (MQ=255)

            gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCaa                                                                                <  1:1524775/91‑1 (MQ=255)

            gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCaa                                                                                >  1:175893/1‑91 (MQ=255)

             ttatatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                               >  2:4026390/3‑91 (MQ=255)

             aaaTATAGCAATTTGAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                               <  2:709727/90‑1 (MQ=255)

              aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAAtt                                                                              <  1:3231739/91‑1 (MQ=255)

              aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAtt                                                                              <  2:1110244/91‑1 (MQ=255)

              aaTATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAAAAAAGAGGTGACGAGTTATTCGGCTAAATATCAAtt                                                                              <  2:5448271/91‑1 (MQ=17)

               atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                             <  2:3711363/91‑1 (MQ=255)

               atatAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                             >  2:3326417/1‑91 (MQ=255)

                ttttgCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGTTCCCGAGTTATTCATCTCAATATCAATTaa                                                                            >  2:4391308/5‑91 (MQ=21)

                    gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                        >  1:3117819/1‑91 (MQ=255)

                    gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTaa                                                                        >  2:4750492/1‑91 (MQ=255)

                    gCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCa                                                                                 >  1:3512932/1‑82 (MQ=255)

                     cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                       <  1:4472702/91‑1 (MQ=255)

                     cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                       <  1:2495663/91‑1 (MQ=255)

                     cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                       <  2:1654509/91‑1 (MQ=255)

                     cAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                       >  1:4026734/1‑91 (MQ=255)

                     cAATTTTAAAGAATTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                       <  2:5720077/91‑1 (MQ=255)

                      aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                                                      >  1:3801146/1‑91 (MQ=255)

                      aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                                                      <  2:5375192/91‑1 (MQ=255)

                      aaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatg                                                                      <  1:1640561/91‑1 (MQ=255)

                       aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                     <  2:5634356/91‑1 (MQ=255)

                       aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                     >  2:2927119/1‑91 (MQ=255)

                       aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatga                                                                     >  2:5154481/1‑91 (MQ=255)

                       aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatga                                                                     >  2:1392594/1‑91 (MQ=255)

                        ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                    <  2:827261/91‑1 (MQ=255)

                        ttttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                    <  1:5317804/91‑1 (MQ=255)

                           tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATgatgat                                                                 >  2:1483282/1‑91 (MQ=255)

                            aaaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGAtt                                                                <  1:4980914/91‑1 (MQ=255)

                             aaGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATta                                                               <  2:2283423/91‑1 (MQ=255)

                              aGATTTAATGCATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                              <  2:4663597/91‑1 (MQ=14)

                              aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                              <  1:3145264/91‑1 (MQ=255)

                              aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                              >  1:2977778/1‑91 (MQ=255)

                              aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                              <  2:4045164/91‑1 (MQ=255)

                              aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATtat                                                              >  2:5240738/1‑91 (MQ=255)

                              aGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                              <  1:2668606/91‑1 (MQ=255)

                                aTTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATaa                                                            <  2:5459526/91‑1 (MQ=255)

                                 tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                           <  1:5858140/91‑1 (MQ=255)

                                 tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAAt                                                           >  1:1469239/1‑91 (MQ=255)

                                 tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                           <  2:5312359/91‑1 (MQ=255)

                                 tttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                           >  1:5475538/1‑91 (MQ=255)

                                  ttAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATg                                                          >  1:522834/1‑91 (MQ=255)

                                   tAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGt                                                         >  1:3928820/1‑91 (MQ=255)

                                    aaTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTa                                                        >  2:3638937/1‑91 (MQ=255)

                                     aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                       >  1:2143599/1‑91 (MQ=255)

                                     aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTaa                                                       <  1:3028756/91‑1 (MQ=255)

                                      tGAATTATAAAGATTGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCACAATATCAATTAACTAATGATGATTATAATGTaaa                                                      <  2:1034465/91‑1 (MQ=255)

                                      tGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTaaa                                                      <  1:707607/91‑1 (MQ=255)

                                         aTTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                   >  1:2074598/1‑91 (MQ=255)

                                              aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                              <  2:175893/91‑1 (MQ=255)

                                              aaaGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAc                                              <  1:1870907/91‑1 (MQ=255)

                                                aGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGt                                            >  1:4693042/1‑91 (MQ=255)

                                                 gATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAGGATGATTATAATGTAAAACAATTACGTa                                           <  2:5782196/91‑1 (MQ=255)

                                                  aTGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTaa                                          <  2:2593295/91‑1 (MQ=255)

                                                       gaTGTATCATATAATCCAGAGGTGCCGAGTTCTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                     <  1:5437858/91‑1 (MQ=255)

                                                       gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                     <  1:1264521/91‑1 (MQ=255)

                                                       gaTATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGat                                     >  1:407926/1‑91 (MQ=255)

                                                             tCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACa                               >  1:894049/1‑91 (MQ=255)

                                                              cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACAt                              <  1:2787991/91‑1 (MQ=255)

                                                              cATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAAGGATGATTATAATGTAAAACAATTACGTAAAAGATATGACAt                              <  1:5798387/91‑1 (MQ=255)

                                                                   aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaa                         <  2:4318011/91‑1 (MQ=255)

                                                                   aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaa                         >  2:4137267/1‑91 (MQ=255)

                                                                   aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaa                         <  1:5664840/91‑1 (MQ=255)

                                                                   aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaa                         <  1:3369182/91‑1 (MQ=255)

                                                                   aaTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaa                         >  2:2069815/1‑91 (MQ=255)

                                                                    aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                        <  2:3104204/91‑1 (MQ=255)

                                                                    aTCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaac                        <  2:18548/91‑1 (MQ=255)

                                                                      ccAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                      <  2:2182418/91‑1 (MQ=255)

                                                                        agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAAt                    >  1:4194483/1‑91 (MQ=255)

                                                                          agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaa                  >  2:505770/1‑91 (MQ=255)

                                                                           ggTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTACAACAATTACGTAAAAGATATGACGTACCAACAAATaaa                 >  1:2239663/1‑91 (MQ=255)

                                                                            gTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAg                >  1:5111452/1‑91 (MQ=255)

                                                                              gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgcg              <  2:1295682/91‑1 (MQ=255)

                                                                                cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc            >  1:461964/1‑91 (MQ=255)

                                                                                cGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGcc            >  2:2798101/1‑91 (MQ=255)

                                                                                 gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTGTAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg           >  2:3437142/1‑91 (MQ=255)

                                                                                 gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg           <  1:4427390/91‑1 (MQ=255)

                                                                                 gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCg           <  1:1152164/91‑1 (MQ=255)

                                                                                   gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa         <  1:2580718/91‑1 (MQ=255)

                                                                                   gTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGaa         <  1:1027133/91‑1 (MQ=255)

                                                                                      aTTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCt      >  2:1826289/1‑91 (MQ=255)

                                                                                        tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt    <  1:4771769/91‑1 (MQ=255)

                                                                                        tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt    <  2:1217931/91‑1 (MQ=255)

                                                                                        tCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgt    >  1:3676381/1‑91 (MQ=255)

                                                                                          aGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg  <  1:1932763/91‑1 (MQ=255)

TAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTG  >  USA300TCH1516_ALE/111422‑111602

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