thx i know how to sort reads by ids
rm(list=ls())#
fastqfile="query.fastq"
library(ShortRead)
reads <- readFastq(fastqfile);#
ids<- as.character(id(reads)); #
ids<- sort(ids, method = "sh", index.return=TRUE)
reads<-reads[ids$ix]
[[alternative HTML version deleted]]
_______________________________________________
Bioc-sig-sequencing mailing list
Bioc-sig-sequencing@r-project.org
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
