Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • salamay
    Member
    • May 2014
    • 20

    #1

    GNU Parallel Losing Records

    I am using GNU parallel to speed up BLAST and I am using this command:

    cat test.fna | parallel --gnu --block 100k --recstart '>' --pipe blastn -db database -query - > test.fna.blastn

    I am noticing that errors arise such as:

    "Error: (1431.1) FASTA-Reader: Warning: FASTA-Reader: No residues given
    BLAST engine error: Warning: Sequence contains no data"

    As well, if I investigate the BLAST output, the first query name is not accurate (it is just a small string from the whole query name), and that some reads are lost in the process (e.g. 4979 reads in BLAST output, but 5000 reads in the original reads file).

    I have checked to make sure there isn't another '>' in the header names so that recstart doesn't split incorrectly. Any other ideas?

    Thanks!

    edit: the pattern I am finding is that for each block that parallel makes, it loses the first ~7-8 reads in that block, and for the first read that actually gets put into the block, the query name is actually read in as the sequence and so it gives an error. Not sure if that helps in solving this issue or not.
    Last edited by salamay; 12-09-2014, 12:10 PM.

Latest Articles

Collapse

  • kinpict.com
    Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by kinpict.com
    The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
    Today, 06:37 PM
  • 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

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 11-12-2026, 07:59 AM
0 responses
27 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-05-2026, 09:25 AM
0 responses
35 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-03-2026, 11:54 AM
2 responses
56 views
0 reactions
Last Post kinpict.com  
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
175 views
0 reactions
Last Post SEQadmin2  
Working...