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]] _______________________________________________ Bioc-sig-sequencing mailing list Bioc-sig-sequencing@r-project.org https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing