Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • felvis56
    replied
    Thank you for the fast reply.

    Apologies if I am being basic. If there is a single base with Q<25 but the following bases are ok will the read be cut at that point or is there a set number of bases needed to be below Q<20 resulting in the read being cut?

    When I look at my fastq files on FastQC the error bars do dip below 20 but I was thinking this was due to a small number of bases over multiple reads.

    Thanks

    Leave a comment:


  • fkrueger
    replied
    Hi Fiona,

    quality trimming removes the portion of the read where the qualities become minimal, but does not remove then entire read (pair) completely. This is taken from Cutadapt --help:

    Code:
     --quality-base=QUALITY_BASE
                            Assume that quality values are encoded as
                            ascii(quality + QUALITY_BASE). The default (33) is
                            usually correct, except for reads produced by some
                            versions of the Illumina pipeline, where this should
                            be set to 64. (Default: 33)

    Leave a comment:


  • felvis56
    replied
    I wanted to ask a question about the quality trimming of sequences. I have Illumina reads and use Trim Galore to remove adapters and primers with good success.

    Is the quality trimming based on the average Phred score of the read or if I use a cutoff of 25 and have any/some bases below 25, will the whole read be removed?

    Thanks,

    Fiona

    Leave a comment:


  • fkrueger
    replied
    Great, enjoy!

    Leave a comment:


  • padmoo
    replied
    It is still running!!!!! What a stupid mistake! Thanks for finding it and shortening my command!

    Leave a comment:


  • fkrueger
    replied
    no that shouldn't matter... Could it be that your queuing system does not take you to the same directory (option like -cwd)? What happens if you just run the command without using the queue?

    By the way you can shorten the command like this:

    Code:
    /gpfs/scratch/cbh12wsu/trim/trim_galore --paired --path_to_cutadapt /gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt [COLOR="Red"]-o adaptrim.fastq[/COLOR] -length 50 1_L.fastq 2_R.fastq
    By the way -o is meant to be an output directory, not a file name. This might be the reason why it can't create the output files in the directory called adaptrim.fastq

    Leave a comment:


  • padmoo
    replied
    these are the folder permissions: drwxr-xr-x
    I also just did a fastqc run and that adds the files to the folder, so it should be fine.

    could it be a problem that the to write file is called .fastqxx.txt?

    Leave a comment:


  • fkrueger
    replied
    Ok, Cutadapt seems to be working now. Do you have write permission in the folder where you keep your FastQ files?

    Leave a comment:


  • padmoo
    replied
    I've installed it and when I run felix's command I get:
    cutadapt version 1.8.1
    ...

    But I still get this error:
    Path to Cutadapt set as: '/gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt' (user defined)
    1.8.1
    Cutadapt seems to be working fine (tested command '/gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt --version')
    Failed to write to file '1_L.fastq_trimming_report.txt': No such file or directory

    Used this command:
    bsub -N -q short -o cut3.txt "/gpfs/scratch/cbh12wsu/trim/trim_galore -phred33 --illumina --paired -phred33 --path_to_cutadapt /gpfs/scratch/cbh12wsu/trim/cutadapt/cutadapt/scripts/cutadapt -o adaptrim.fastq -length 50 1_L.fastq 2_R.fastq"

    More ideas?

    Leave a comment:


  • dpryan
    replied
    Just install it as a user (pip install --user cutadapt) or use virtualenv.

    Leave a comment:


  • fkrueger
    replied
    Hmm, maybe you need to contact your sys-admins then, you obviously need to get Cutadapt to be working before you can start using it. You should even be able to install it in your home directory and then point there?

    Leave a comment:


  • padmoo
    replied
    I put it back where it belongs and tried again. Same error as before.

    With normal installation you mean the pip command? I tried that before but I'm not allowed to do that on the cluster. I can load a cutadapt module but I was hoping making it an executable would solve the problem.

    Leave a comment:


  • fkrueger
    replied
    Cutadapt probably needs its internal structure to be intact, can't you just do a normal installation and then point to the cutadapt executable (without renaming it to .py?)

    Leave a comment:


  • padmoo
    replied
    Hi Felix,

    thanks for replying so quick!

    I get the following error when I try your command:
    File "/gpfs/scratch/cbh12wsu/trim/cutadapt.py", line 99
    print(rest, match.read.name, file=self.file)
    ^
    SyntaxError: invalid syntax

    I think I should mention that I made cutadapt.py an executable. And copied it from it's original folder into the trim folder.

    Leave a comment:


  • fkrueger
    replied
    Hmm, when you just type "/gpfs/scratch/cbh12wsu/trim/cutadapt.py" on the command line do you see something like:

    Code:
    cutadapt version 1.8.1
    Copyright (C) 2010-2015 Marcel Martin <[email protected]>
    
    cutadapt removes adapter sequences from high-throughput sequencing reads.
    
    Usage:
        cutadapt -a ADAPTER [options] [-o output.fastq] input.fastq
    ?

    Leave a comment:

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 02:06 PM
0 responses
7 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
27 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
47 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
59 views
0 likes
Last Post seqadmin  
Working...
X