BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,129,368	A→G	5.7%	intergenic (+201/+37)	hlb_2 → / ← USA300HOU_RS10870	Phospholipase C/putative leukocidin‑like protein 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,129,368	0	A	G	5.7%	406.2 / 8.9	121	intergenic (+201/+37)	hlb_2/USA300HOU_RS10870	Phospholipase C/putative leukocidin‑like protein 1
Reads supporting (aligned to +/- strand): ref base A (54/60); new base G (3/4); total (57/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TCATTAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTT  >  USA300TCH1516_ALE/2129278‑2129458

tCATTAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTa                                                                                            >  1:1206036/1‑91 (MQ=255)

  aTTAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTAcc                                                                                          <  1:1279054/91‑1 (MQ=255)

  aTTAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTAcc                                                                                          <  1:1278741/91‑1 (MQ=255)

  aTTAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAACAAATCTATCTATTCAAATATAAATAAAAGTTAcc                                                                                          <  1:1065551/91‑1 (MQ=255)

   ttAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCCATTCAAATATAAATAAAAGTTACct                                                                                         >  1:1786584/1‑91 (MQ=255)

     aaTACAAGTCTATACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctac                                                                                       <  1:3760613/91‑1 (MQ=255)

      aTACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctact                                                                                      >  1:1095848/1‑91 (MQ=255)

      aTACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctact                                                                                      <  1:3886468/91‑1 (MQ=255)

       tACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctacta                                                                                     <  1:758654/91‑1 (MQ=255)

       tACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACctacta                                                                                     >  1:451132/1‑91 (MQ=255)

        aCAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCtactac                                                                                    <  1:1033574/91‑1 (MQ=255)

          aaGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACAt                                                                                  <  1:1561448/91‑1 (MQ=255)

          aaGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACAt                                                                                  <  1:3408763/91‑1 (MQ=255)

               tttACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTAt                                                                             >  1:126877/1‑91 (MQ=255)

                     ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                       <  1:3550509/91‑1 (MQ=255)

                     ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                       >  1:193177/1‑91 (MQ=255)

                     ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                       >  1:3250919/1‑91 (MQ=255)

                     ctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                       <  1:771787/91‑1 (MQ=255)

                       acacTCATACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                     <  1:2069206/91‑1 (MQ=255)

                       acacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                     <  1:2586562/91‑1 (MQ=255)

                       acacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                     <  1:2425699/91‑1 (MQ=255)

                       acacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAgg                                                                     <  1:600194/91‑1 (MQ=255)

                        cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                    >  1:3505322/1‑91 (MQ=255)

                        cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                    >  1:955271/1‑91 (MQ=255)

                        cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                    >  1:1295139/1‑91 (MQ=255)

                         acTCAGACAAGATTCATACACTGCACGTCATAATAAATCGATCTATTCAAATAGAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                   <  1:3860216/91‑1 (MQ=255)

                         acTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                   <  1:1635355/91‑1 (MQ=255)

                         acTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                   <  1:1483908/91‑1 (MQ=255)

                           tCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAAc                                                                 <  1:1458351/91‑1 (MQ=21)

                           tCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAAc                                                                 <  1:1399893/91‑1 (MQ=255)

                           tCAAACAAGATTCATACACTGCACGGCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAAc                                                                 <  1:3849438/91‑1 (MQ=255)

                           ccAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAAc                                                                 <  1:2167243/90‑1 (MQ=255)

                            cAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACt                                                                >  1:633441/1‑91 (MQ=255)

                               aCAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACtttt                                                             >  1:3549304/1‑91 (MQ=255)

                                 aaGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttat                                                           <  1:289507/91‑1 (MQ=255)

                                    aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTttattta                                                        <  1:1269915/91‑3 (MQ=17)

                                     ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTttatttat                                                       <  1:3883448/91‑4 (MQ=17)

                                     ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                       <  1:2547868/91‑1 (MQ=255)

                                     ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                       <  1:3412227/91‑1 (MQ=255)

                                     ttCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACt                                                       <  1:3355806/91‑1 (MQ=255)

                                      tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                      <  1:944800/91‑1 (MQ=255)

                                      tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                      <  1:1981610/91‑1 (MQ=255)

                                      tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                      <  1:944595/91‑1 (MQ=255)

                                        aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACttct                                                    >  1:2444349/1‑89 (MQ=255)

                                        aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                    >  1:2444166/1‑91 (MQ=255)

                                        aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                    >  1:3842783/1‑91 (MQ=255)

                                        aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                    >  1:3241489/1‑91 (MQ=255)

                                        aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTATATGTAGCAGGCAACTTTTATTACTTAt                                                    >  1:3241251/1‑91 (MQ=255)

                                         tACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAtt                                                   >  1:4054459/1‑91 (MQ=255)

                                                cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                            <  1:1952613/91‑1 (MQ=255)

                                                cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                            <  1:1826547/91‑1 (MQ=255)

                                                cACGTCATAATAAATCTATCTATGCAAATATAAATAAAAGTTACCTACTACATGCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                            <  1:1800818/91‑1 (MQ=255)

                                                  cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                          <  1:2638457/91‑1 (MQ=255)

                                                  cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                          >  1:716604/1‑91 (MQ=255)

                                                  cGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcat                                          >  1:3157006/1‑91 (MQ=255)

                                                      ataataAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatc                                      <  1:3500987/91‑1 (MQ=255)

                                                         ataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatt                                   <  1:1111367/91‑1 (MQ=255)

                                                          taaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatta                                  <  1:1298463/91‑1 (MQ=255)

                                                          taaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatta                                  <  1:3604504/91‑1 (MQ=255)

                                                            aaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAg                                <  1:234332/91‑1 (MQ=255)

                                                            aaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAg                                >  1:719600/1‑91 (MQ=255)

                                                             atctatctATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGt                               >  1:929615/1‑91 (MQ=255)

                                                              tctatctaTTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTa                              <  1:635771/91‑1 (MQ=255)

                                                              tctatctaTTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTa                              >  1:1293146/1‑91 (MQ=255)

                                                                 atctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtt                           <  1:2866113/91‑1 (MQ=255)

                                                                     atTCAAATCTAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTAc                       >  1:637764/1‑91 (MQ=255)

                                                                       tCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAATTCTTTTCATTATCATTAAGTACTTTTACaa                     >  1:3743088/1‑91 (MQ=255)

                                                                        cAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACaaa                    >  1:4105823/1‑91 (MQ=255)

                                                                        cAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACaaa                    <  1:2946657/91‑1 (MQ=255)

                                                                         aaaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAAc                   <  1:676355/91‑1 (MQ=255)

                                                                          aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCGTTATCATTAAGTACTTTTACAAACt                  >  1:3860751/1‑91 (MQ=255)

                                                                          aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACt                  <  1:178745/91‑1 (MQ=255)

                                                                          aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACt                  >  1:789354/1‑91 (MQ=255)

                                                                          aaTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACt                  <  1:178691/91‑1 (MQ=255)

                                                                           atatAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACtt                 >  1:3961623/1‑91 (MQ=255)

                                                                            tataAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTc                <  1:2105378/91‑1 (MQ=255)

                                                                              taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac              >  1:769587/1‑91 (MQ=255)

                                                                              taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac              >  1:770211/1‑91 (MQ=255)

                                                                              gtaataaaAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac              >  1:1777754/3‑91 (MQ=17)

                                                                              gtaataaaAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac              >  1:1751998/3‑91 (MQ=17)

                                                                               aaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcaca             >  1:2913336/1‑91 (MQ=255)

                                                                                aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt            >  1:1283683/1‑91 (MQ=255)

                                                                                aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt            >  1:1320985/1‑91 (MQ=255)

                                                                                aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt            >  1:1287206/1‑91 (MQ=255)

                                                                                aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt            >  1:1283651/1‑91 (MQ=255)

                                                                                 ataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAtt           >  1:1270610/1‑91 (MQ=255)

                                                                                 ataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAtt           <  1:2690648/91‑1 (MQ=255)

                                                                                  taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa          >  1:76410/1‑91 (MQ=255)

                                                                                  taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa          >  1:3341058/1‑91 (MQ=255)

                                                                                  taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa          <  1:3483386/91‑1 (MQ=255)

                                                                                  taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa          >  1:1115612/1‑91 (MQ=255)

                                                                                  taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa          <  1:3800225/91‑1 (MQ=255)

                                                                                  taaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTa          >  1:1329836/1‑91 (MQ=255)

                                                                                    aaaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtg        >  1:3385480/1‑91 (MQ=255)

                                                                                    aaaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtg        >  1:3385513/1‑91 (MQ=255)

                                                                                      aGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACACACTTCACATTAtgtg      <  1:3546681/91‑1 (MQ=255)

                                                                                      aGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgtg      >  1:3602691/1‑91 (MQ=255)

                                                                                       gTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATgtgt     <  1:3645616/91‑1 (MQ=255)

                                                                                       gTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATgtgt     <  1:2536104/91‑1 (MQ=255)

                                                                                         tACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCt   <  1:3129438/91‑1 (MQ=255)

TCATTAATACAAGTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTT  >  USA300TCH1516_ALE/2129278‑2129458

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