BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,309,793:1	+G	9.1%	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	9.1%	877.6 / 82.5	208	coding (76/1032 nt)	yfiZ_2	putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand): ref base . (82/107); new base G (9/10); total (91/117)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 228 total aligned reads displayed.

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

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

  aCATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑ta                                                                                           >  2:2403041/1‑91 (MQ=255)

  aCATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑ta                                                                                           >  1:95159/1‑91 (MQ=255)

   cATAGTGAATGATATCGCTAAAATTAATTTTCGCAGATCCAATCGCGATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tat                                                                                          >  1:3091692/1‑91 (MQ=255)

    aTAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tata                                                                                         >  2:3781922/1‑91 (MQ=255)

    aTAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tata                                                                                         >  1:548736/1‑91 (MQ=255)

      aGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tatata                                                                                       >  1:1115658/1‑91 (MQ=255)

      aGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tatata                                                                                       >  1:1116100/1‑91 (MQ=255)

      aGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tatata                                                                                       >  2:870311/1‑91 (MQ=255)

        tGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGtatata                                                                                       >  2:2685960/1‑91 (MQ=255)

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

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

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

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

             gATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAgt                                                                                <  1:3564470/91‑1 (MQ=255)

               tatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAgtgt                                                                              <  1:101563/91‑1 (MQ=255)

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

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

                   gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTtgtg                                                                          >  2:2735840/1‑91 (MQ=255)

                    cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGc                                                                         >  1:2802798/1‑91 (MQ=255)

                      aaacTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                      <  1:3007614/91‑1 (MQ=255)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctt                                                                       <  2:2677729/91‑3 (MQ=21)

                      aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                      >  2:2043154/1‑91 (MQ=255)

                       aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGcttt                                                                      <  2:3481820/91‑4 (MQ=21)

                       aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtc                                                                     >  2:3636335/1‑91 (MQ=255)

                       aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtc                                                                     >  2:65289/1‑91 (MQ=255)

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

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

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

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

                         attaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCGC‑TATATACGTGAGTGTTGTGC‑Gtctt                                                                   <  2:1917247/91‑1 (MQ=255)

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

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

                         attaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctt                                                                   >  1:3421418/1‑91 (MQ=255)

                           taattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttt                                                                 >  1:1543407/1‑91 (MQ=255)

                            aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtcttttt                                                                >  2:636361/1‑91 (MQ=255)

                             attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttc                                                               <  2:1821138/91‑1 (MQ=255)

                              ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtatttttct                                                              >  2:3449437/1‑91 (MQ=255)

                                ttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttcttt                                                            <  2:1364575/91‑1 (MQ=255)

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

                                 tCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttctttt                                                           >  1:463576/1‑91 (MQ=255)

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

                                  cGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTa                                                          <  1:2323962/91‑1 (MQ=255)

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

                                   gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTAt                                                         >  2:129717/1‑91 (MQ=255)

                                      gaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGAt                                                      >  1:3445510/1‑91 (MQ=255)

                                        aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCt                                                    <  1:1290585/91‑1 (MQ=255)

                                        aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCt                                                    >  1:3158777/1‑91 (MQ=255)

                                               gCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAg                                             <  1:2391662/91‑1 (MQ=255)

                                                ccATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGa                                            >  1:2352187/1‑91 (MQ=255)

                                                ccATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGa                                            <  2:60796/91‑1 (MQ=255)

                                                 cATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAt                                           >  2:2587981/1‑91 (MQ=255)

                                                     tttAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGc                                       >  2:2854692/1‑91 (MQ=255)

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

                                                        aaaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                    <  2:2796514/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                           tatataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                                 <  1:2590653/91‑1 (MQ=255)

                                                             tataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAtt                               <  2:2985974/91‑1 (MQ=255)

                                                               taCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTc                             <  1:392165/91‑1 (MQ=255)

                                                                aCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCt                            <  2:286350/91‑1 (MQ=255)

                                                                  aaaTAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCttt                          <  2:3556939/91‑1 (MQ=255)

                                                                     tAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtat                       >  1:1721638/1‑91 (MQ=255)

                                                                     tAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtat                       >  1:1722154/1‑91 (MQ=255)

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

                                                                        aaagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGc                    >  2:2060508/1‑91 (MQ=255)

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

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                   >  1:2554350/1‑91 (MQ=255)

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                   >  2:2124849/1‑91 (MQ=255)

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                   >  2:1176642/1‑91 (MQ=255)

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

                                                                           agacgCACAACACTCACGTATATA‑GTGAGT‑TTGTGCTTTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                  >  2:3820226/5‑91 (MQ=9)

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

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

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

                                                                            aaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt                >  2:734826/1‑91 (MQ=255)

                                                                             aagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt               <  1:1100816/91‑1 (MQ=255)

                                                                              agCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTa              <  1:2834920/91‑1 (MQ=255)

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

                                                                                cacaA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACt            >  1:2252478/1‑91 (MQ=255)

                                                                                cacaA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACt            >  2:482770/1‑91 (MQ=255)

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

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

                                                                                  caA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc          >  1:1897715/1‑91 (MQ=255)

                                                                                  caA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc          >  1:2944218/1‑91 (MQ=255)

                                                                                   aa‑aCTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTcc         <  2:927100/91‑1 (MQ=255)

                                                                                   aa‑aCTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTcc         <  2:1710964/91‑1 (MQ=255)

                                                                                      aCTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCat       >  2:3850508/1‑90 (MQ=255)

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

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

                                                                                       cTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc      <  1:482770/91‑1 (MQ=255)

                                                                                        tCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt     <  2:128463/91‑1 (MQ=255)

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

                                                                                           c‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc  >  1:180944/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` `< 32 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`