dear all
     i have ids4 which is subset of id(reads), i want to use ids4
to subset seqs or reads, but still canon work
such is the coding


fastqfile="query.fastq"
library(ShortRead)
reads <- readFastq(fastqfile)
seqs <- sread(reads)
names(seqs)<- id(reads)
ids4<- read.table("ids4.txt")
seqs[ids4]

thx
gaoshan

	[[alternative HTML version deleted]]

