Hi,
I'm using the following function to extract a subselection of reads based on
their ID from a FastQ file:
extract.reads.fq<-function(dir,fq.file,id.list)
{
rfq <- readFastq(dir, pattern=fq.file)
tmp=extract.id(rfq)
sv.reads=rfq[which(id.list%in%tmp)]
writeFastq(sv.reads, paste("sv_",fq.file,sep=""))
}

The function works well but for one file i got the following error :
Error in .local(dirPath, pattern, ...) :
  ShortRead internal: too many 'snap' entries
Tried to google it but nothing is coming out, anyone encoutred this error ?

Thanks in advance for your help.




----------------------------------------------------------------
Kind regards / Met vriendelijke groet,
Mohamed-Ramzi Temanni,
----------------------------------------------------------------

	[[alternative HTML version deleted]]

