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
