BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,309,793:1	+G	8.6%	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	8.6%	485.2 / 39.6	116	coding (76/1032 nt)	yfiZ_2	putative siderophore transport system permease protein YfiZ
Reads supporting (aligned to +/- strand): ref base . (53/53); new base G (2/8); total (55/61)
Fisher's exact test for biased strand distribution p-value = 9.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 130 total aligned reads displayed.

TAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTT  >  USA300TCH1516_ALE/2309703‑2309881

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

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

 aaCATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑t                                                                                         >  1:1472717/1‑91 (MQ=255)

 aaCATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑t                                                                                         >  1:429740/1‑91 (MQ=255)

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

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

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

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

     tAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATCTTTAAATATATACAAATAAAAAGAAAGCAGAA‑ACTCAC‑tatat                                                                                     >  1:4144994/1‑91 (MQ=255)

     tAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑tatat                                                                                     <  1:2142914/91‑1 (MQ=255)

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

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

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

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

        tGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACg                                                                                  >  1:3602010/1‑91 (MQ=255)

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

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

           aTGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCGC‑TATATACGTGa                                                                               <  1:3198828/91‑1 (MQ=255)

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

            tGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAg                                                                              >  1:2152469/1‑91 (MQ=255)

            tGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAg                                                                              >  1:2149014/1‑91 (MQ=255)

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

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

              atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAgtg                                                                            <  1:190869/91‑1 (MQ=255)

              atatCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAgtg                                                                            <  1:1296775/91‑1 (MQ=255)

                atCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGtgtt                                                                          >  1:1634348/1‑91 (MQ=255)

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

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

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

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

                 tCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGtgttg                                                                         <  1:1631083/91‑1 (MQ=255)

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

                  cGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGtgttgt                                                                        <  1:549229/91‑1 (MQ=255)

                   gCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTtgtg                                                                       <  1:1676460/91‑1 (MQ=255)

                   gCTAAAATTAATCTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTtgtg                                                                       >  1:3199174/1‑91 (MQ=255)

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

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

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

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

                        aattaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtct                                                                  <  1:1963773/91‑1 (MQ=255)

                         attaattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctt                                                                 <  1:3360773/91‑1 (MQ=255)

                            aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtcttttt                                                              >  1:2859769/1‑91 (MQ=255)

                            aattTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtcttttt                                                              >  1:1333460/1‑91 (MQ=255)

                             attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctttttc                                                             <  1:661767/91‑1 (MQ=255)

                             attTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctttttc                                                             >  1:1437042/1‑91 (MQ=255)

                              ttttCGAAGAACCAATCGCCATATTTAAATATATATAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctttttct                                                            >  1:2652277/1‑91 (MQ=255)

                              ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctttttct                                                            >  1:2653359/1‑91 (MQ=255)

                              ttttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctttttct                                                            >  1:1767664/1‑91 (MQ=255)

                                ttCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGtctttttcttt                                                          >  1:2171759/1‑91 (MQ=255)

                                   gaagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTCAGTGTTGTGCGTCTTTTTCTTTTAt                                                       >  1:3870321/1‑91 (MQ=255)

                                    aagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATg                                                      <  1:3810432/91‑1 (MQ=255)

                                    aagaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATg                                                      <  1:3104959/91‑1 (MQ=255)

                                      gaaCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGAt                                                    <  1:1840373/91‑1 (MQ=255)

                                        aCCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCt                                                  >  1:491508/1‑91 (MQ=255)

                                         ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                 >  1:1346758/1‑91 (MQ=255)

                                         ccAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTa                                                 >  1:1347201/1‑91 (MQ=255)

                                          cAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTAt                                                >  1:3947165/1‑91 (MQ=255)

                                            aTCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATgg                                              >  1:3470902/1‑91 (MQ=255)

                                              cGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCa                                            >  1:2158551/1‑91 (MQ=255)

                                                ccATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGa                                          <  1:1792116/91‑1 (MQ=255)

                                                 cATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGAt                                         <  1:1889061/91‑1 (MQ=255)

                                                   tatTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGCGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTg                                       >  1:253254/1‑91 (MQ=255)

                                                     tttAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGc                                     <  1:1392145/91‑1 (MQ=255)

                                                      ttAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCt                                    <  1:3442777/91‑1 (MQ=255)

                                                        aaaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                  >  1:1110139/1‑91 (MQ=255)

                                                        aaaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                  >  1:4004141/1‑91 (MQ=255)

                                                        aaaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGCGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTAc                                  <  1:1445249/91‑1 (MQ=255)

                                                         aaTATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACt                                 >  1:3348928/1‑91 (MQ=255)

                                                           tatataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTa                               >  1:784605/1‑91 (MQ=255)

                                                             tataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAtt                             <  1:1952725/91‑1 (MQ=255)

                                                              ataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                            <  1:2768699/91‑1 (MQ=255)

                                                              ataCAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTAttt                            <  1:3741249/91‑1 (MQ=255)

                                                                 cAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCtt                         >  1:24540/1‑91 (MQ=255)

                                                                     tAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtat                     <  1:519779/91‑1 (MQ=255)

                                                                      aaaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtata                    <  1:1025995/91‑1 (MQ=255)

                                                                      aaaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTtata                    <  1:3170257/91‑1 (MQ=255)

                                                                       aaaaGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAg                   >  1:1261987/1‑91 (MQ=255)

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                 >  1:819221/1‑91 (MQ=255)

                                                                         aagaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGcc                 >  1:3692775/1‑91 (MQ=255)

                                                                          agaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                <  1:745256/91‑1 (MQ=255)

                                                                          agaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                <  1:745243/91‑1 (MQ=255)

                                                                          agaaagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCa                <  1:2753651/91‑1 (MQ=255)

                                                                            gacgCACAACACTCACGTATATA‑GTGAGT‑TTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAtt              >  1:3892070/4‑91 (MQ=21)

                                                                            aaatCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGACAtt              <  1:3272931/91‑1 (MQ=255)

                                                                             aagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCGAttt             >  1:3456101/1‑91 (MQ=255)

                                                                             aagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt             >  1:3642670/1‑91 (MQ=255)

                                                                             aagCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCAttt             >  1:2921087/1‑91 (MQ=255)

                                                                               gCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc           >  1:3259370/1‑91 (MQ=255)

                                                                               gCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTAc           >  1:4050544/1‑91 (MQ=255)

                                                                                 acaA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACtt         <  1:2152100/91‑1 (MQ=255)

                                                                                  caA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc        <  1:2132608/91‑1 (MQ=255)

                                                                                  caA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc        <  1:2132592/91‑1 (MQ=255)

                                                                                  caA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTc        >  1:1948769/1‑91 (MQ=255)

                                                                                    a‑aCTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCa      >  1:3982294/1‑91 (MQ=255)

                                                                                      aCTCACGTATATA‑GTGAGT‑TTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc    >  1:2709528/1‑91 (MQ=38)

                                                                                      aCTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAc     >  1:2322500/1‑91 (MQ=255)

                                                                                       cTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCAcc    >  1:2709898/1‑91 (MQ=255)

                                                                                        tCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt   <  1:3455929/91‑1 (MQ=255)

                                                                                        tCAC‑GATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCt   <  1:3457231/91‑1 (MQ=255)

                                                                                         cAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCtt  <  1:1405369/91‑1 (MQ=255)

TAACATAGTGAATGATATCGCTAAAATTAATTTTCGAAGAACCAATCGCCATATTTAAATATATACAAATAAAAAGAAAGCACAA‑ACTCAC‑TATATACGTGAGTGTTGTGCGTCTTTTTCTTTTATGATCTATGGCAGATTGGCTACTTATTTCTTTTATAGCCATTTACTTCCACCTT  >  USA300TCH1516_ALE/2309703‑2309881

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