Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • bio_jit
    Member
    • Jun 2013
    • 13

    #1

    ShortRead streamer

    Hi I am using R in windows 7
    I got a small program to split a large fastq file in 1 million reads/file and name them automatically. Programm running O.k. but I want to change the output folder how I do that?
    my programm

    FastqFile<- "./Sequencing/ERR12730_1.fastq.gz"

    fileBaseName<-sub("fastq$", " ", FastqFile)

    #iterate over Fastqfile

    f<- FastqStreamer(FastqFile, 1000000)
    file_index<-0

    while(length(fq<-yield(f))){
    newName<-paste(fileBaseName, "_", file_index, ".fastq", sep =" ")
    writeFastq( fq, file=newName)
    file_index<- file_index+1
    }
    close(f)

    I want to send the output files in other folder say ./Sequencing/output
    How to do?
    Plz help me out!
  • TiborNagy
    Senior Member
    • Mar 2010
    • 329

    #2
    insert the following code before writeFastq:
    Code:
    newName <- paste('Sequencing/output', newName, sep='/')

    Comment

    • bio_jit
      Member
      • Jun 2013
      • 13

      #3
      not working!

      Hi
      i did as you said but getting this error now

      Error in writeFastq(fq, file = newName) :
      failed to open file '/Sequencing/output./Sequencing/ERR127306_1.fastq.gz_0.fastq'

      any suggestions?

      Comment

      • TiborNagy
        Senior Member
        • Mar 2010
        • 329

        #4
        Sorry, I did not realize the file name contains directory structures.
        Code:
        newName <- paste('Sequencing/output', basename(newName), sep='/')
        The output will be: Sequencing/output/ERR12730_1.fastq.gz_0.fastq
        Last edited by TiborNagy; 03-18-2014, 06:01 AM. Reason: missing characters

        Comment

        Latest Articles

        Collapse

        • SEQadmin2
          16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
          by SEQadmin2





          NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
          10-30-2026, 06:46 AM
        • SEQadmin2
          Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
          by SEQadmin2



          Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

          Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
          ...
          10-14-2026, 09:55 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by SEQadmin2, Today, 07:59 AM
        0 responses
        13 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 11-05-2026, 09:25 AM
        0 responses
        32 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 11-03-2026, 11:54 AM
        0 responses
        45 views
        0 reactions
        Last Post SEQadmin2  
        Started by SEQadmin2, 10-28-2026, 09:57 AM
        0 responses
        170 views
        0 reactions
        Last Post SEQadmin2  
        Working...