[BioC] Sequence output formatting

Tobias Straub tstraub at med.uni-muenchen.de
Mon Jul 9 17:37:14 CEST 2007


sorry... forgot something:

I'd recommend using sth like

con <- file("fastaseq.txt", open="wt")
for (i in 1:length(acc))
{
	writeLines(paste(">",acc[i],sep=""), con, sep="\n")
	writeLines(getSEQ(acc[i]), con, sep="\n")
}
close(con)

On Jul 9, 2007, at 5:26 PM, Tobias Straub wrote:

> I'd recommend using sth like
>
> con <- file("fastaseq.txt", open="wt")
> for (i in 1:length(acc))
> {
> 	writeLines(paste(">",acc[i],sep=""))
> 	writeLines(getSEQ(acc[i]))
> }
> close(con)

======================================================================
Dr. Tobias Straub         Adolf-Butenandt-Institute, Molecular Biology
tel: +49-89-2180 75 439         Schillerstr. 44, 80336 Munich, Germany



More information about the Bioconductor mailing list