BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,309,800	Δ1 bp	15.2%	coding (69/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,800	0	C	.	15.2%	768.0 / 131.7	243	coding (69/1032 nt)	yfiZ_2	putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand): ref base C (105/100); new base . (9/28); total (115/128)
Fisher's exact test for biased strand distribution p-value = 3.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 248 total aligned reads displayed.

GTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTT  >  USA300TCH1516_ALE/2309710‑2309890

gTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATAc                                                                                             >  1:3014506/1‑91 (MQ=255)

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

  gAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGt                                                                                           >  1:2878331/1‑91 (MQ=255)

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

    aTGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGa                                                                                         >  2:1234032/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

          tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATTAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGtgttg                                                                                   <  1:3342151/91‑1 (MQ=255)

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

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

             cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAGAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGc                                                                                <  2:147847/91‑1 (MQ=255)

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

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

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

              tAAAATTAAATTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑g                                                                              <  2:80785/91‑1 (MQ=255)

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

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

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

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

               aaaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACACAC‑TATATACGTGAGTGTTGTGC‑Gt                                                                             >  1:3680517/1‑91 (MQ=255)

                aaaTTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAACACTCACGTATATA‑GTGAGT‑TTGTGcttt                                                                             <  1:570438/91‑4 (MQ=21)

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

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

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

                      attTTCGAAGAACCAATCGCCATATTTAAATTTATACAAATAAAAAGAAAGCACAA‑AATCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttc                                                                      >  1:3510455/1‑91 (MQ=255)

                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttc                                                                      >  1:2636575/1‑91 (MQ=255)

                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttc                                                                      >  1:505542/1‑91 (MQ=255)

                      attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttc                                                                      >  1:547490/1‑91 (MQ=255)

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

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

                           cGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATATGTGAGTGTTGTGC‑GTCTTTTTCTTTTa                                                                 >  2:3078726/1‑91 (MQ=255)

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

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

                               gaaCCAATCGCCATATTTAAATATATACAAAAAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGAt                                                             >  1:3512689/1‑91 (MQ=255)

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

                                   cAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTAt                                                         >  1:3267791/1‑91 (MQ=255)

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

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

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

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

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

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

                                             atTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTgg                                               >  1:2221352/1‑91 (MQ=255)

                                              tttAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGc                                              <  1:158162/91‑1 (MQ=255)

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

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

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

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

                                                 aaaGATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATGCGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                           <  2:1370156/91‑1 (MQ=255)

                                                  aataaataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                          >  1:3123237/1‑91 (MQ=255)

                                                  aaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                          >  1:3601064/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                                          <  2:1126498/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                aaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAGTTCTTTTATAg                            >  1:3584696/1‑91 (MQ=255)

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

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

                                                                  aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                          <  1:1521376/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:3753883/91‑1 (MQ=255)

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

                                                                       agCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTa                     >  1:3853300/1‑91 (MQ=255)

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

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

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

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

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

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

                                                                               aCTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc              >  1:1875147/1‑91 (MQ=255)

                                                                                  cacgtatata‑GTGAGT‑TTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa          >  2:3250788/5‑91 (MQ=255)

                                                                                  cacgtatata‑GTGAGT‑TTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTa          >  1:4036/5‑91 (MQ=255)

                                                                                   acgtatata‑GTGAGT‑TTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTAc         >  1:2563623/4‑91 (MQ=255)

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

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

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

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

                                                                                          taCGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg    <  2:3841111/91‑1 (MQ=255)

                                                                                          taCGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg    >  2:305658/1‑91 (MQ=255)

                                                                                          taCGTGAGTGTTGTGC‑GTCTTTTTCTATTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTg    >  2:2237732/1‑91 (MQ=255)

                                                                                           aCGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGt   <  2:17485/91‑1 (MQ=255)

                                                                                            cGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGtt  >  2:3584155/1‑91 (MQ=255)

GTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTTACTTTTGTT  >  USA300TCH1516_ALE/2309710‑2309890

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