BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	348,156	A→C	10.3%	A31A (GCA→GCC)	yezG_5 →	putative antitoxin YezG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	348,156	0	A	C	10.3%	746.1 / 70.9	234	A31A (GCA→GCC)	yezG_5	putative antitoxin YezG
Reads supporting (aligned to +/- strand): ref base A (88/122); new base C (6/18); total (94/140)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 238 total aligned reads displayed.

GACTTTCGAAGAAAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAGTGTGAT  >  USA300TCH1516_ALE/348066‑348245

gACTTTCGAAGAAAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCa                                                                                           <  1:244364/91‑1 (MQ=255)

      cGAAGAAAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCa                                                                                     <  1:1118835/91‑1 (MQ=255)

      cGAAGAAAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCa                                                                                     <  2:2050295/91‑1 (MQ=255)

       gaagaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCCATTGCat                                                                                    <  2:2679006/91‑1 (MQ=255)

       gaagaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCCATTGCat                                                                                    <  2:1477140/91‑1 (MQ=255)

        gagaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCCATTGCata                                                                                   <  2:2283483/90‑1 (MQ=255)

        aagaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCCATTGCata                                                                                   <  2:2930639/91‑1 (MQ=255)

        aagaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCCATTGCata                                                                                   <  1:312964/91‑1 (MQ=255)

         agaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCatat                                                                                  >  1:3665662/1‑91 (MQ=255)

         agaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCatat                                                                                  >  1:3211528/1‑91 (MQ=255)

         agaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCatat                                                                                  <  2:2825594/91‑1 (MQ=255)

         agaaAAATTAAGTGAAATGTAAAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCatat                                                                                  <  1:99027/91‑1 (MQ=255)

          gaaAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCatata                                                                                 >  2:2652779/1‑91 (MQ=255)

               aTTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGAt                                                                            >  2:2395305/1‑91 (MQ=255)

               aTTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGAt                                                                            >  2:1693598/1‑91 (MQ=255)

                ttAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATg                                                                           >  2:1561914/1‑91 (MQ=255)

                ttAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATg                                                                           <  1:790466/91‑1 (MQ=255)

                ttAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATg                                                                           <  2:1272115/91‑1 (MQ=255)

                ttAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATg                                                                           >  1:2659592/1‑91 (MQ=255)

                ttAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATg                                                                           <  2:208926/91‑1 (MQ=255)

                ttAAGTAAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATAAAAGAAAAAAAGAAAAAAGTATATGCAATGGCATATATCGATg                                                                           <  2:51954/91‑1 (MQ=17)

                 tAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATga                                                                          >  1:3543545/1‑91 (MQ=255)

                   aGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATgaag                                                                        <  2:3659270/91‑1 (MQ=255)

                    gTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATgaaga                                                                       <  1:3341530/91‑1 (MQ=255)

                     tGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATgaagaa                                                                      <  1:517057/91‑1 (MQ=255)

                     tGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGCTAGAGTGGGAAAAAGTATATGAAATGGCATATATCGATgaagaa                                                                      >  1:3085398/1‑91 (MQ=255)

                        aaTGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAAtgt                                                                   >  2:3062890/1‑91 (MQ=255)

                         aTGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAAtgtg                                                                  <  2:1796798/91‑1 (MQ=255)

                          tGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTATAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTgg                                                                 <  2:3487890/91‑1 (MQ=255)

                          tGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTgg                                                                 <  1:1561914/91‑1 (MQ=255)

                          tGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTgg                                                                 <  1:606475/91‑1 (MQ=255)

                           gTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGga                                                                <  1:1974626/91‑1 (MQ=255)

                           gTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGga                                                                <  1:1978254/91‑1 (MQ=255)

                            tataGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGgag                                                               >  2:2013067/1‑91 (MQ=255)

                             ataGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGgaga                                                              <  2:1067734/91‑1 (MQ=255)

                              taGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGaa                                                             >  2:756244/1‑91 (MQ=255)

                                gCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAgt                                                           <  1:76715/91‑1 (MQ=255)

                                 cGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAgtg                                                          >  1:2806420/1‑91 (MQ=255)

                                 cGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAgtg                                                          <  1:1720512/91‑1 (MQ=255)

                                  gagaTTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAgtgt                                                         <  2:838614/91‑1 (MQ=255)

                                   agaTTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGtt                                                        <  1:2818454/91‑1 (MQ=255)

                                     aTTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGtttt                                                      >  1:3849891/1‑91 (MQ=255)

                                     aTTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGtttt                                                      <  1:2085941/91‑1 (MQ=255)

                                     aTTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGtttt                                                      >  1:508999/1‑91 (MQ=255)

                                       tGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTAc                                                    <  1:1227278/91‑1 (MQ=255)

                                       tGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTAc                                                    <  1:2721746/91‑1 (MQ=255)

                                        gCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACa                                                   >  1:2036596/1‑91 (MQ=255)

                                        gCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACa                                                   <  1:2698054/91‑1 (MQ=255)

                                         cGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACaa                                                  >  2:2488364/1‑91 (MQ=255)

                                          gAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAAt                                                 <  1:1787192/91‑1 (MQ=255)

                                          gAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAAt                                                 <  1:1478008/91‑1 (MQ=255)

                                              aaGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATtata                                             <  1:400199/91‑1 (MQ=255)

                                              aaGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATtata                                             <  2:3502993/91‑1 (MQ=255)

                                               aGATGATGAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAAGGTGGAGAAGTGTTTTACAATTATAc                                            <  1:167052/91‑1 (MQ=255)

                                                gATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACt                                           <  2:1903726/91‑1 (MQ=255)

                                                  ttattaGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGa                                         >  1:2325400/1‑91 (MQ=255)

                                                   tattaGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGaa                                        <  2:3867653/91‑1 (MQ=255)

                                                   tattaGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGaa                                        <  2:324801/91‑1 (MQ=255)

                                                    attaGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAAc                                       >  2:1413115/1‑91 (MQ=255)

                                                     ttaGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAAcc                                      >  2:767759/1‑91 (MQ=255)

                                                      taGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCa                                     <  1:583723/91‑1 (MQ=255)

                                                         cATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGc                                  >  1:2918933/1‑91 (MQ=255)

                                                         cATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGc                                  >  1:448337/1‑91 (MQ=255)

                                                           tGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAg                                >  1:3534924/1‑91 (MQ=255)

                                                           tGATACCGGGAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAAGGTGGCGAAGTGTTTTACAATTATACTGAACCAAGCAg                                >  2:2336340/1‑91 (MQ=255)

                                                             aTACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGtg                              <  1:3471272/91‑1 (MQ=255)

                                                              tACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGtga                             <  2:3669354/91‑1 (MQ=255)

                                                               aCCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGtgat                            <  1:1283241/91‑1 (MQ=255)

                                                                ccGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGtgatg                           >  2:1065962/1‑91 (MQ=255)

                                                                ccGGTAGAGTGGGAAAAAGTATATGCAATAGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAATCAGtgatg                           >  2:3381170/1‑91 (MQ=255)

                                                                  ggTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGaa                         >  2:3682444/1‑91 (MQ=255)

                                                                   gTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAat                        <  1:1368250/91‑1 (MQ=255)

                                                                    tAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAatt                       <  2:1890234/91‑1 (MQ=255)

                                                                     agagTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAatta                      >  1:3739915/1‑91 (MQ=255)

                                                                     agagTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAatta                      <  1:2925182/91‑1 (MQ=255)

                                                                      gagTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTTCAATTATACTGAACTAAGCAGTGATGAattat                     >  1:3326604/1‑91 (MQ=255)

                                                                       agTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAattatt                    <  1:2906889/91‑1 (MQ=255)

                                                                          gggAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTa                 >  1:3750281/1‑91 (MQ=255)

                                                                          gggAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTa                 >  1:2974049/1‑91 (MQ=255)

                                                                           ggAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGGAGTGTTTTTCAATTATACTGAACCAAGCAGTGATGAATTATTTTAc                <  1:3088831/91‑1 (MQ=255)

                                                                           ggAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACGAAGCCGTGATGAATTATTTTAc                >  2:3127633/1‑91 (MQ=255)

                                                                             aaaaaGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACta              <  1:752978/91‑1 (MQ=255)

                                                                             aaaaaGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACta              <  1:338983/91‑1 (MQ=255)

                                                                               aaaGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACtata            >  2:833509/1‑91 (MQ=255)

                                                                                  gTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGtgatga                          <  2:435446/74‑1 (MQ=255)

                                                                                  gTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACaa         >  2:2918359/1‑91 (MQ=255)

                                                                                  gTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACaa         >  2:1472350/1‑91 (MQ=255)

                                                                                   tataTGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAg        >  1:1157516/1‑91 (MQ=255)

                                                                                   tataTGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAg        <  1:1342672/91‑1 (MQ=255)

                                                                                    atatGCCATTGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgt       >  1:988094/1‑91 (MQ=255)

                                                                                    atatGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgt       >  1:2161744/1‑91 (MQ=255)

                                                                                    atatGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgt       >  2:325568/1‑91 (MQ=255)

                                                                                    atatGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgt       <  2:2388036/91‑1 (MQ=255)

                                                                                    atatGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgt       >  2:3032/1‑91 (MQ=255)

                                                                                     tatGCCATTGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgtg      >  1:2533916/1‑91 (MQ=255)

                                                                                     tatGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAgtg      <  1:483608/91‑1 (MQ=255)

                                                                                       tGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAGtgtg    <  2:1305124/91‑1 (MQ=255)

                                                                                        gCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAGTGTGa   >  2:1050665/1‑91 (MQ=255)

                                                                                        gCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAGTGTGa   <  1:3030532/91‑1 (MQ=255)

                                                                                         cAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAGTGTGAt  <  2:1252574/91‑1 (MQ=255)

GACTTTCGAAGAAAAATTAAGTGAAATGTATAGCGAGATTGCGAATAAGATTATTAGCATGATACCGGTAGAGTGGGAAAAAGTATATGCAATGGCATATATCGATGAAGAATGTGGAGAAGTGTTTTACAATTATACTGAACCAAGCAGTGATGAATTATTTTACTATACAAGTGTGAT  >  USA300TCH1516_ALE/348066‑348245

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