Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts

  • cement_head
    replied
    Originally posted by Romualdo View Post
    Well I'm having a hard time with this, Blast2GO-basic remotely blasting just take too long for each sequence, so I got more speed trying it locally with Blast+ Blastx and importing the output xml on Blast2GO for subsequently steps. But Will Nelson is right, is impraticable do this on this computer. Our lab is about to buy a server with 128gb RAM, until then I wanna be more experienced with this, so I made a 100 sequences sample.

    So I got this repeatedly when running the Blastx:

    CFastaReader: Bad gap size at line ***
    CFastaReader: Problem parsing gap mods at line ***

    When "***" are line numbers, this lines matches with sequence id lines, that use this format:

    >?_GroupUn999_2_939_+

    What in this format is generating that error ?
    Make absolutely sure that you buy ECC RAM. Anything less and you will have major problems

    Leave a comment:


  • GenoMax
    replied
    What format are your sequences in?

    That error seems to indicate that there may be a problem with your fasta file. Can you try to replace the "?_" at the beginning of the header? Looks like that may be causing a problem.

    Leave a comment:


  • Romualdo
    replied
    Well I'm having a hard time with this, Blast2GO-basic remotely blasting just take too long for each sequence, so I got more speed trying it locally with Blast+ Blastx and importing the output xml on Blast2GO for subsequently steps. But Will Nelson is right, is impraticable do this on this computer. Our lab is about to buy a server with 128gb RAM, until then I wanna be more experienced with this, so I made a 100 sequences sample.

    So I got this repeatedly when running the Blastx:

    CFastaReader: Bad gap size at line ***
    CFastaReader: Problem parsing gap mods at line ***

    When "***" are line numbers, this lines matches with sequence id lines, that use this format:

    >?_GroupUn999_2_939_+

    What in this format is generating that error ?

    Leave a comment:


  • GenoMax
    replied
    @Will Nelson: curious if you have a roughly equivalent spec computer as the OP. Did you actually time a search?

    Leave a comment:


  • Will Nelson
    replied
    Sorry, my math was all wrong on my last post. Let me try again.

    In reality, it takes at least 5 minutes for blastx to align one transcript to nr. For 16,000 sequences, with 4 threads, that is (16,000x5)/4 = 20,000 minutes, or 13.8 days. Then if you want to get GOs by importing into the blast2GO free version, that takes several more days at least.

    Leave a comment:


  • Will Nelson
    replied
    Blasting against nr is not easy. Even with 4 threads, to blast 16,000 sequences will take around 4,000 minutes, or 66 days. Performing GO assignment is also not easy. Importing 16,000 blastx results into the free version of Blast2GO, and then doing GO assignments, will take many days.

    Leave a comment:


  • cement_head
    replied
    Originally posted by sarvidsson View Post
    No, MPI-BLAST only makes sense on a cluster (if you have access to one). Local multithreading BLAST will be faster than local MPI multiprocessing, and more memory-efficient.

    You might run into trouble with only 8 GB RAM if you want to BLAST the complete nr locally, give it a try but you may get out-of-memory problems. i7-3770 would be 4 cores with hyperthreading, so be prepared for several days of BLASTing...
    Yes, correct - only if you have (access to) a cluster.

    Leave a comment:


  • sarvidsson
    replied
    No, MPI-BLAST only makes sense on a cluster (if you have access to one). Local multithreading BLAST will be faster than local MPI multiprocessing, and more memory-efficient.

    You might run into trouble with only 8 GB RAM if you want to BLAST the complete nr locally, give it a try but you may get out-of-memory problems. i7-3770 would be 4 cores with hyperthreading, so be prepared for several days of BLASTing...

    Leave a comment:


  • Romualdo
    replied
    Originally posted by cement_head View Post
    Do this locally. Download the nr database and use BLASTX against it.

    BUT, it is MUCH faster if you use mpiBLAST on a cluster.

    The command for mpiBLAST (with the correct flags for B2GO) is something like:

    mpiblast -p blastx -d nr -i input.fa -v 20 -b 20 -I T -e 0.001 -m 7 -o output.xml
    mpiBLAST is much faster even in this pc alone ?

    Leave a comment:


  • cement_head
    replied
    Do this locally. Download the nr database and use BLASTX against it.

    BUT, it is MUCH faster if you use mpiBLAST on a cluster.

    The command for mpiBLAST (with the correct flags for B2GO) is something like:

    mpiblast -p blastx -d nr -i input.fa -v 20 -b 20 -I T -e 0.001 -m 7 -o output.xml

    Leave a comment:


  • Romualdo
    started a topic Blast2GO Beginner's Question

    Blast2GO Beginner's Question

    Hi everybody, I'm new in NGS. I used the flow STAR/Cufflinks/Cuffcompare and now i need annotate my transcripts, and I decided to use Blast2GO because it seems more intuitive.
    But not so totaly intuitive for a total beginner like me, and I don't know what fashion of blast I should perform. The basic version of program allow me to use this methods:

    QBlast@NCBI. NCBI o ers a public service that allows searching molecular sequence
    databases with the BLAST algorithm. The main advantages of making use of this service
    are its versatility and that no database maintenance is required. Therefore by selecting
    this option at Blast2GO no additional installations have to be done.

    Remote BLAST. Blast2GO will download the latest BLAST+ executable form NCBI and
    will use it to query NR or other databases remotely.

    Local BLAST against own database. It is possible to use BLAST+ excuteble to query a
    local/own database.

    WWW-BLAST. Alternatively, BLAST can be done locally against a custom database. For
    this, you need to place a copy of your FASTA formatted custom DB plus a WWW-BLAST
    installation on a local BLAST server and indicate Blast2GO their location.

    My fasta have 16450 sequences, and I want to use database NCBI NR Full, I have a i7 3770 8gb ram computer.

    So the question is: with this resources what is the most safe and easy way to Blast ?

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, Today, 12:22 PM
0 responses
9 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-11-2026, 10:35 AM
0 responses
11 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-06-2026, 07:41 AM
0 responses
30 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 08-03-2026, 10:13 AM
0 responses
48 views
0 reactions
Last Post SEQadmin2  
Working...