Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • zisis86
    Member
    • Mar 2013
    • 10

    Usage of writeFastq function from ShortRead

    Hi,
    I am totaly new and i have a question about the syntax of writeFastq function.
    I have this script that search for primer sequences in a file and saves good reads from fuzzy in .good.fq file and bad reads in .bad.fq file> i want to use writeFastq functions from ShortRead package to save "good" reads from fuzzy match in "*.good.fq" file and "bad" reads in *".bad.fq" file.Can you help me with the code that i have to fill in in order to do this .
    here is my code until now.

    gc()
    rm(list=ls())
    library(ShortRead)
    #library(Biostrings)

    #read bait sequences
    baitseq <- readLines("/home/users/bait_seq.txt")
    baitseq

    #choose case
    num <- 1
    exptPath <- "/home/users/INDEX1.fastq"
    exptPath

    #read data
    fq <- readFastq(dirPath=exptPath)
    sread(fq)

    #count ocurrences of bait
    pos <- grepl(baitseq[num],sread(fq))
    pos[1:25]
    # number of reads
    length(pos)
    # number of hits
    sum(pos)

    # fuzzy match - 2 substitutions
    fpos <- agrep(baitseq[num],sread(fq), max = list(sub = 2, ins=0, del=0))
    fpos[1:25]
    # number of fuzzy hits
    length(fpos)

    # writeFastq functions in ShortRead package to save "good" reads from fuzzy match in "*.good.fq" file and "bad" reads in *".bad.fq" file # ???????
    ............


    #tables
    #table <- tables(sread(fq), n=5)
    #table$top
    Last edited by zisis86; 10-20-2013, 11:32 PM.

Latest Articles

Collapse

  • SEQadmin2
    Nine Things a Sample Prep Scientist Thinks About Before Sequencing
    by SEQadmin2


    I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.


    Here are nine questions we think about, in roughly the order they matter, before...
    06-18-2026, 07:11 AM
  • SEQadmin2
    From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data
    by SEQadmin2


    Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.


    The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
    ...
    06-02-2026, 10:05 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 06-17-2026, 06:09 AM
0 responses
24 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-09-2026, 11:58 AM
0 responses
42 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-05-2026, 10:09 AM
0 responses
48 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-04-2026, 08:59 AM
0 responses
49 views
0 reactions
Last Post SEQadmin2  
Working...