Hi Hrishi, load the probe package for your array type, which contains the sequence of the 25mers. http://www.bioconductor.org/data/probes/Packages/hgu95av2probe_1.0.zip or http://www.bioconductor.org/data/probes/Packages/hgu95av2probe_1.0.tar .gz library(hgu95av2probe) data(hgu95av2probe) summary(hgu95av2probe) my.out.seq <- hgu95av2probe$sequence[out.of.bounds] my.in.seq <- hgu95av2probe$sequence[in.bounds] Justin ----- Justin Borevitz Ecology and Evolution (CLSC 915E) University of Chicago 1101 E. 57th St. Chicago, IL 60637 USA http://natural-systems.org ___________________________________ Message: 15 Date: Tue, 28 Dec 2004 12:41:22 -0800 (PST) From: hrishikesh deshmukh <d_hrishikesh@yahoo.com> Subject: [BioC] Sequence data and filtering To: bioconductor@stat.math.ethz.ch Message-ID: <20041228204122.92550.qmail@web50301.mail.yahoo.com> Content-Type: text/plain; charset=us-ascii Dear Bioconductorians, My sincere apologies if this email is duplicated.I have searched the mailing list and could not find a way to solve my problem: Find out the probe intensity (PM adn MM) for HG_U95Av2 which fall in (200 to 20000 units) range and extract their respective PM sequences. I have been able to filter the probes but how do i extract the sequence data?! I am including the small script which i have written for filtering, now how do i extract sequence information. Thanks in advance. Hrishi