Unconfigured Ad

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

    #1

    Help with script that seperates reads from fastq file

    hi,
    I have fastq files of reads and i would like to finish this script in order to select the "good" reads from the "bad" reads and after that to save my results in new fastq files(good.fastq bad.fastq) .
    Thanks for the help

    gc()
    rm(list=ls())
    library(ShortRead)
    showMethods(readFastq)
    showMethods(writeFastq)
    library(Biostrings)

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

    #choose case
    num <- 1
    exptPath <- ("Downloads/INDEX1_reads.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)

    # save "good" reads from from fuzzy match in "*.good.fq" file and "bad" reads in *".bad.fq" file


    ?????

Latest Articles

Collapse

  • SEQadmin2
    Beyond CRISPR/Cas9: Understand, Choose, and Use the Right Genome Editing Tool
    by SEQadmin2



    CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).

    Despite this, “CRISPR helped turn genome editing from a specialized technique into
    ...
    07-31-2026, 11:01 AM
  • SEQadmin2
    Proteomic Platforms: How to Choose the Right Analytical Strategy to Improve Detection and Clinical Applications
    by SEQadmin2


    Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.

    The systematic characterization of the human proteome has
    ...
    07-20-2026, 11:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, Yesterday, 10:35 AM
0 responses
9 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-06-2026, 07:41 AM
0 responses
27 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-03-2026, 10:13 AM
0 responses
45 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 07-31-2026, 02:55 AM
0 responses
48 views
0 reactions
Last Post SEQadmin2  
Working...