search for: writefasta

It looks like Biostrings function "writeFASTA" overwrites the output file at each run. It seems it does not support the "append" parameter. I have to generate one big file gathering a miRNA identifier and relative sequence followd by a variable number of dara records pertaining such a miRNA target genes transcription. Each...

...write(miRNA.seq,zz, append = TRUE) # geneDesc <- paste (">",gene.id, "|", gene.map[i,"ensembl_transcript_id"], sep="") geneSeq <- gene.seq[i,"3utr"] gene.string <- list(desc=geneDesc, seq=geneSeq) writeFASTA (gene.string, zz) Thank you in advance for your help. Maura tutti i telefonini TIM! [[alternative HTML version deleted]]