BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	53,905	C→G	5.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	5.1%	859.9 / 21.0	246	S281R (AGC→AGG)	USA300HOU_RS00220	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (110/123); new base G (2/11); total (112/134)
Fisher's exact test for biased strand distribution p-value = 4.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

TGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGT  >  USA300TCH1516_ALE/53815‑53993

tGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGTTTTTGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGc                                                                                          >  2:3394729/1‑91 (MQ=255)

tGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGc                                                                                          >  2:2082824/1‑91 (MQ=255)

tGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGc                                                                                          <  1:2214864/91‑1 (MQ=255)

tGGAGATTTAGGTAATGAAGATTTAAACGTGTACAGTCTTTTTGAACAAGTATTGATAAGATTAACTGATAACATTGAGGAATCATTTAGc                                                                                          >  1:3092603/1‑91 (MQ=255)

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

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

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

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

    gaTTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGATAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTaa                                                                                      >  2:3553408/1‑91 (MQ=255)

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

      tttAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAat                                                                                    >  2:2444066/1‑91 (MQ=255)

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

      tttAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAat                                                                                    >  1:2486541/1‑91 (MQ=255)

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

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

       ttAGGTAATGAAGATTTAAACGTGTATAAAAAAAATAAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatg                                                                                   <  2:3862346/91‑1 (MQ=255)

        tAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatga                                                                                  >  1:2923610/1‑91 (MQ=255)

        tAGGTAATGAAGATCTAAACGTGTTTAGTTTTTTTGAAAAAGTCTTGATAGAATTAACTGATAACATTGAGGAATCATTTAGGTTAAatga                                                                                  <  2:3610902/91‑1 (MQ=255)

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

          ggTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatgatg                                                                                >  2:176098/1‑91 (MQ=255)

          ggTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatgatg                                                                                <  2:2383963/91‑1 (MQ=255)

          ggTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAatgatg                                                                                <  2:1741546/91‑1 (MQ=255)

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

           gTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATgg                                                                               >  2:636801/1‑91 (MQ=255)

           gTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATgg                                                                               >  1:1595176/1‑91 (MQ=255)

            tAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGt                                                                              >  1:3134048/1‑91 (MQ=255)

             aaTGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGtt                                                                             <  1:2539063/91‑1 (MQ=255)

             aaTGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGtt                                                                             <  1:2662244/91‑1 (MQ=255)

             aaTGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGtt                                                                             <  1:1834577/91‑1 (MQ=255)

              aTGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGttt                                                                            <  2:2134675/91‑1 (MQ=255)

              aTGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGttt                                                                            >  1:106633/1‑91 (MQ=255)

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

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

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

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

                  aGATTTAAACGTGTATAGTTTTTTTTTTTTTTTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACt                                                                        >  2:3760888/1‑91 (MQ=255)

                  aGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATGATTTAGCCTAAATGATGGTTTGACt                                                                        >  1:3602899/1‑91 (MQ=255)

                   gATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTa                                                                       <  1:1771968/91‑1 (MQ=255)

                    aTTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAg                                                                      >  2:3672758/1‑91 (MQ=255)

                    aTTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAg                                                                      >  1:1486154/1‑91 (MQ=255)

                       tAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACa                                                                   <  2:3583865/91‑1 (MQ=255)

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

                          aCGTGTATGGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACaaaa                                                                >  1:3377765/1‑91 (MQ=255)

                           cGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAAt                                                               <  2:2822770/91‑1 (MQ=255)

                           cGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAAt                                                               <  1:1381780/91‑1 (MQ=255)

                                ataGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGc                                                          <  1:3006053/91‑1 (MQ=255)

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

                                  aGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCaa                                                        <  2:2742587/91‑1 (MQ=255)

                                  aGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCaa                                                        <  1:2332734/91‑1 (MQ=255)

                                   gTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAag                                                       >  1:622015/1‑91 (MQ=255)

                                   gTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAag                                                       <  2:2892090/91‑1 (MQ=255)

                                    tttttttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAaga                                                      >  1:2194150/1‑91 (MQ=255)

                                    tttttttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAaga                                                      <  2:1533219/91‑1 (MQ=255)

                                      tttttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAgaga                                                    <  2:1950951/91‑1 (MQ=255)

                                       ttttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAt                                                   <  2:3277093/91‑1 (MQ=255)

                                       ttttGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCAAAATGATGGTTTGACTAGACAAAATTATGCAAGAGAt                                                   <  1:3278486/91‑1 (MQ=255)

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

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

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

                                           gAAAAAGTATTGATAAAATTAACTGAGAACATTGAGGAATCATTTAGCCTGAATGATGGTTTGACTGGACAAAATTATGCAAGAGATTTaa                                               <  2:250635/91‑1 (MQ=255)

                                            aaaaaGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTaaa                                              <  2:2306377/91‑1 (MQ=255)

                                             aaaaGTATTGATAAAATTAACTGATAACATTGAGGAATCATATAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTaaaa                                             >  1:3366644/1‑91 (MQ=255)

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

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

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

                                                 gTATTGATAAAATTAACTGATAACATTGCGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAACGTAg                                         >  1:468846/1‑91 (MQ=255)

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

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

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

                                                       aTAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTc                                   <  1:3351345/91‑1 (MQ=255)

                                                       aTAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTc                                   >  2:709419/1‑91 (MQ=255)

                                                         aaaaTTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTcgc                                 >  2:3034434/1‑91 (MQ=255)

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

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

                                                              tAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGc                            <  1:553051/91‑1 (MQ=255)

                                                               aaCTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCa                           <  2:1812682/91‑1 (MQ=255)

                                                                aCTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTcg                                  <  2:2538971/83‑1 (MQ=255)

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

                                                                aCTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCAc                          <  2:1128454/91‑1 (MQ=255)

                                                                aCTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAAGTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCAc                          <  1:3625845/91‑1 (MQ=255)

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

                                                                  tGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTg                        >  1:3230636/1‑91 (MQ=255)

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

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

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

                                                                       ataTTGATGAATCATGTAGCCGAAATGAGGGTTTGACGAGACAAAATTATGCAAGATATTTAAAAGTAGCAATTCGCGAAGCACTGGTTaa                   <  1:2523896/89‑1 (MQ=255)

                                                                           tGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACa               <  2:741966/91‑1 (MQ=255)

                                                                           tGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACa               <  1:383753/91‑1 (MQ=255)

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

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

                                                                                caTCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAat          >  1:1336092/2‑91 (MQ=255)

                                                                                 aTCATTTAGGCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatg         >  1:2923966/1‑91 (MQ=255)

                                                                                 aTCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatg         <  2:3842056/91‑1 (MQ=255)

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

                                                                                      ttAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAgtgcatgc    >  2:3435735/1‑91 (MQ=255)

                                                                                      ttAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgcatgt    <  2:3134048/91‑2 (MQ=255)

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

                                                                                      ttAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAatgcatgc    >  2:3208717/1‑91 (MQ=255)

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

                                                                                        aGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGt  >  1:492321/1‑91 (MQ=255)

TGGAGATTTAGGTAATGAAGATTTAAACGTGTATAGTTTTTTTGAAAAAGTATTGATAAAATTAACTGATAACATTGAGGAATCATTTAGCCTAAATGATGGTTTGACTAGACAAAATTATGCAAGAGATTTAAAAGTAGCAATTCGCGAAGCACTGGTTAATACATTAATGCATGCGT  >  USA300TCH1516_ALE/53815‑53993

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