BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,309,787:1	+C	12.7%	coding (82/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,787	1	.	C	12.7%	983.4 / 152.8	251	coding (82/1032 nt)	yfiZ_2	putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand): ref base . (113/106); new base C (4/28); total (117/134)
Fisher's exact test for biased strand distribution p-value = 1.94e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 254 total aligned reads displayed.

GACCAGTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCA  >  USA300TCH1516_ALE/2309697‑2309877

  ccAGTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑Ac                                                                                            <  2:3471439/91‑1 (MQ=255)

  ccAGTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑Ac                                                                                            <  1:2956616/91‑1 (MQ=255)

    aGTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTc                                                                                          <  1:3305453/91‑1 (MQ=255)

    aGTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTc                                                                                          <  2:2061398/91‑1 (MQ=255)

     gTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCa                                                                                         >  1:982014/1‑91 (MQ=255)

     gTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCa                                                                                         >  2:119036/1‑91 (MQ=255)

     gTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCa                                                                                         >  2:2322340/1‑91 (MQ=255)

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

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

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

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

        aCATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑ta                                                                                     <  2:25525/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:1062258/1‑91 (MQ=255)

          aTAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tata                                                                                   <  1:1721995/91‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                          cTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGc                                                                   <  1:27761/91‑1 (MQ=255)

                           tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctt                                                                 <  1:2468624/91‑3 (MQ=17)

                           tAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAG‑ACGCACAACACTCACGTATATA‑GTGAGT‑TTGTGctt                                                                 <  1:252818/91‑3 (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                                                                 <  2:1603282/91‑3 (MQ=21)

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

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

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

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

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

                                ttaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtcttt                                                            <  2:404474/91‑1 (MQ=255)

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

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

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

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

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

                                    ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑Gtctttttct                                                        >  1:34867/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‑TATATACGTGAGTGTTGTGC‑Gtctttttcttaaa                                                    <  1:2410189/91‑4 (MQ=255)

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

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

                                           agaaCCAATCGCCATATTTAAATATATACGAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGa                                                 <  2:511500/91‑1 (MQ=255)

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

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

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

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

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

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

                                                    cGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCa                                        <  2:1802792/91‑1 (MQ=255)

                                                     gCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGAGTTGTGC‑GTCTATTTCTTTTATGATCTATCGCAg                                       >  2:1854686/1‑91 (MQ=255)

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

                                                       cATATTTAAATATATACAAATAAAAAGAAAGCACAA‑AGTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAt                                     >  1:3326580/1‑91 (MQ=255)

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

                                                       cATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGAt                                     <  2:3826783/91‑1 (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)

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

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

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

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

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

                                                               aaTATATACAAATAAGTAGAAAGCACAA‑AGTCAC‑TATATTCGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                             <  2:3866868/91‑1 (MQ=255)

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

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

                                                                    ataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                        >  2:2780124/1‑91 (MQ=255)

                                                                        aaaTAAAAAGATAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCttt                    <  1:3585941/91‑1 (MQ=255)

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

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

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

                                                                            aaaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtata                <  1:2672384/91‑1 (MQ=255)

                                                                             aaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAg               >  1:2427784/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:3228388/91‑1 (MQ=255)

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

                                                                                   acgCACAACACTCACGTATATA‑GTGAGT‑TTGTGCTTTCTTTTTCTTTTATGATCTATGGCAGATTCGCTACTTATTTCTTTTATAGCCAtt          >  1:3696805/3‑91 (MQ=9)

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

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

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

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

                                                                                       acaA‑ACTCAC‑TATATACGCGAGTGTAGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTCACtt     >  2:1449397/1‑91 (MQ=255)

                                                                                          a‑aCTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCa  >  1:2167452/1‑91 (MQ=255)

GACCAGTAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGC‑GTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCA  >  USA300TCH1516_ALE/2309697‑2309877

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