BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,309,793:1	+G	14.7%	coding (76/1032 nt)	yfiZ_2 ←	putative siderophore transport system permease protein YfiZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,309,793	1	.	G	14.7%	845.4 / 158.5	224	coding (76/1032 nt)	yfiZ_2	putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand): ref base . (96/95); new base G (4/29); total (100/124)
Fisher's exact test for biased strand distribution p-value = 3.91e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 254 total aligned reads displayed.

TAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAC  >  USA300TCH1516_ALE/2309703‑2309883

tAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAc                                                                                              <  1:1316361/91‑1 (MQ=255)

tAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAc                                                                                              <  1:2591555/91‑1 (MQ=255)

tAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAc                                                                                              >  1:1604826/1‑91 (MQ=255)

tAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAc                                                                                              >  1:1811840/1‑91 (MQ=255)

tAACAGAGTGAATGATATCGCTGAAATTAATTGTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAc                                                                                              <  2:3196218/91‑1 (MQ=255)

  aCATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑ta                                                                                           <  2:266015/91‑1 (MQ=255)

   cATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCTATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tat                                                                                          >  2:3552542/1‑91 (MQ=255)

   cATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tat                                                                                          >  2:1060194/1‑91 (MQ=255)

   cATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCAGATTTAAATATATACAAATAAAAAGAAAGCACAA‑CCTCAC‑tat                                                                                          >  2:1686061/1‑91 (MQ=255)

     tAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tatat                                                                                        >  1:1702107/1‑91 (MQ=255)

     tAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tatat                                                                                        >  2:2898929/1‑91 (MQ=255)

        tGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCACGtatatag                                                                                      <  1:1491696/91‑2 (MQ=255)

         gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGt                                                                                    <  1:811506/91‑1 (MQ=255)

          aaTGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTg                                                                                   >  2:3221153/1‑91 (MQ=255)

           aTGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGa                                                                                  <  2:500988/91‑1 (MQ=255)

             gATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAg                                                                                 <  2:580314/91‑1 (MQ=38)

              atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAgt                                                                                <  1:3417125/91‑1 (MQ=38)

              atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAgt                                                                                <  1:385080/91‑1 (MQ=38)

               tatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAgtt                                                                               <  2:3786821/91‑2 (MQ=37)

               tatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAgtt                                                                               <  1:2818114/91‑2 (MQ=37)

                atCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAgttt                                                                              <  2:428964/91‑3 (MQ=37)

                 tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAgtttg                                                                             <  1:1359591/91‑4 (MQ=37)

                 tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGtgttg                                                                            >  2:3656763/1‑91 (MQ=255)

                   gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑Ttgtg                                                                          <  1:520259/91‑1 (MQ=25)

                    cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGc                                                                         <  1:625506/91‑1 (MQ=25)

                    cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGct                                                                        <  2:1520606/91‑2 (MQ=17)

                     tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGct                                                                        <  2:754630/91‑2 (MQ=25)

                     tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctt                                                                       <  2:2117023/91‑3 (MQ=17)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctt                                                                       <  1:2333842/91‑3 (MQ=21)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                      <  1:3736597/91‑1 (MQ=255)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                      <  2:2783803/91‑1 (MQ=255)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                      <  2:2139140/91‑1 (MQ=255)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                      <  1:2524164/91‑1 (MQ=255)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGcttt                                                                      <  1:3636025/91‑4 (MQ=12)

                       aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtc                                                                     <  1:275243/91‑1 (MQ=255)

                       aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctttc                                                                     <  2:557892/91‑5 (MQ=9)

                       aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctttc                                                                     <  1:3487516/91‑5 (MQ=9)

                        aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctttc                                                                     <  2:2423895/91‑5 (MQ=21)

                        aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtct                                                                    <  1:3720102/91‑1 (MQ=255)

                         attaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctt                                                                   <  1:353000/91‑1 (MQ=255)

                            aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtcttttt                                                                <  2:375904/91‑1 (MQ=255)

                                 tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttctttt                                                           <  1:874771/91‑1 (MQ=255)

                                   gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTAt                                                         <  2:2831337/91‑1 (MQ=255)

                                    aagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATg                                                        <  1:56633/91‑1 (MQ=255)

                                     agaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGa                                                       <  2:2759803/91‑1 (MQ=255)

                                      gaaTCAATCGCCATATTTAAATATATACAAAAAAAAAGAAAGCACAA‑ACTCAC‑TATAGACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGAt                                                      >  1:3512607/1‑91 (MQ=255)

                                       aaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATc                                                     >  1:3890997/1‑91 (MQ=255)

                                       aaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATc                                                     <  1:3793708/91‑1 (MQ=255)

                                         ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTa                                                   >  1:2639021/1‑91 (MQ=255)

                                         ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTa                                                   >  1:308611/1‑91 (MQ=255)

                                            aTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATgg                                                >  2:1585982/1‑91 (MQ=255)

                                             tCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGc                                               <  1:433185/91‑1 (MQ=255)

                                             tCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGc                                               >  1:358001/1‑91 (MQ=255)

                                                 cATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAt                                           <  2:3826783/91‑1 (MQ=255)

                                                 cATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAt                                           <  1:2229855/91‑1 (MQ=255)

                                                  atatTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAtt                                          >  1:197711/1‑91 (MQ=255)

                                                  atatTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAtt                                          <  1:688681/91‑1 (MQ=255)

                                                  atatTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGCGTTGTGC‑GTCTTTTTCTATTATGATCTATGGCAGAtt                                          >  1:1043593/1‑91 (MQ=255)

                                                   tatTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTg                                         <  2:2719851/91‑1 (MQ=255)

                                                      ttAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCt                                      <  2:73462/91‑1 (MQ=255)

                                                      ttAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCt                                      >  2:1486773/1‑91 (MQ=255)

                                                       tAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTa                                     <  2:1301158/91‑1 (MQ=255)

                                                       tAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTa                                     <  1:1209811/91‑1 (MQ=255)

                                                        aaaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                    >  2:470797/1‑91 (MQ=255)

                                                         aaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                   >  1:205486/1‑91 (MQ=255)

                                                         aaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                   <  2:3866893/91‑1 (MQ=255)

                                                         aaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                   >  1:3600779/1‑91 (MQ=255)

                                                         aaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                   <  2:1126498/91‑1 (MQ=255)

                                                          atatatACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACtt                                  <  1:2381061/91‑1 (MQ=255)

                                                          atatatACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACtt                                  <  1:431294/91‑1 (MQ=255)

                                                           tatataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                 >  1:2092059/1‑91 (MQ=255)

                                                           tatataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                 >  1:2939283/1‑91 (MQ=255)

                                                            atataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAt                                <  1:277234/91‑1 (MQ=255)

                                                              ataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                              <  2:426272/91‑1 (MQ=255)

                                                                 cAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtt                           >  1:2666543/1‑91 (MQ=255)

                                                                  aaaTAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCttt                          >  2:3141143/1‑91 (MQ=255)

                                                                  aaaTAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCttt                          >  1:1132027/1‑91 (MQ=255)

                                                                       aaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAg                     >  1:2427783/1‑91 (MQ=255)

                                                                        aaagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCCGATTGGCTACTTATTTCTTTTATAGc                    >  1:3314589/1‑91 (MQ=255)

                                                                        aaagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGc                    <  1:2210745/91‑1 (MQ=255)

                                                                        aaagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATCTCTTTCATAGc                    >  2:1832044/1‑91 (MQ=255)

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                   <  1:3335573/91‑1 (MQ=255)

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGca                   <  2:1401614/91‑2 (MQ=255)

                                                                           gaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAt                 >  2:2969251/1‑91 (MQ=255)

                                                                            aaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                <  1:1064512/91‑1 (MQ=255)

                                                                            aaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                <  1:3228388/91‑1 (MQ=255)

                                                                             acgCACAACACTCACGTATATA‑GTGAGT‑TTGTGCTTTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCact                >  2:1982210/3‑89 (MQ=9)

                                                                             aagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt               >  2:253248/1‑91 (MQ=255)

                                                                               gCACAACACTCACGTATATA‑GTGAGT‑TTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc             >  1:2539699/1‑91 (MQ=25)

                                                                               gCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc             >  2:1286443/1‑91 (MQ=255)

                                                                               gCACAA‑ACTCAC‑TATATACGTGAGCGTTGTGC‑GTCTTTTTCTTTTATGATATATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc             >  1:245271/1‑91 (MQ=255)

                                                                                 acaA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt           >  2:3067490/1‑91 (MQ=255)

                                                                                  caA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc          <  2:2060397/91‑1 (MQ=255)

                                                                                      aCTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc       <  1:1527325/91‑1 (MQ=255)

                                                                                          aC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGTAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa   <  2:3492072/91‑1 (MQ=255)

                                                                                          aC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa   >  2:1370446/1‑91 (MQ=255)

                                                                                           c‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATTGCAGATTGTCTACTTATTTCTTTTATAGCCATTTACTTTCACCtttc  >  2:820275/1‑89 (MQ=255)

                                                                                           c‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc  >  2:2601669/1‑91 (MQ=255)

                                                                                           c‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc  >  2:3866851/1‑91 (MQ=255)

                                                                                           c‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc  >  1:1604736/1‑91 (MQ=255)

TAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAC  >  USA300TCH1516_ALE/2309703‑2309883

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