SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics


Reply
 
Thread Tools
Old 02-09-2010, 05:04 PM   #1
KevinLam
Senior Member
 
Join Date: Nov 2009
Location: SEA
Posts: 114
Default solid_denovo_postprocessor.pl removes velvet contig sequence.

This question might be better directed to SOLID developers
But I have a problem where my final basespace output I have just a fasta header and no sequence. Did the solid scripts remove it because of low quality?

These are the last few steps.

velvet_0.7.55/velvetg_de ./S4output_directory/ -read_trkg yes -amos_file yes

solid_denovo_postprocessor.pl --afgfile ./S4output_directory/velvet_asm.afg --csfasta ./S4-output/colorspace_input.csfasta --run fragment --output ./sample4_input

echo 'Running denovoadp'
time denovoadp ./sample4_input 30000 > ./sample4_output


so after step 1
I get this in contigs.fa
>NODE_780_length_22_cov_13.590909
CCCCCCCCCCCCCCCCCACATGGCCGCGCCCCCCCCCCCCCC

from step 2 onwards the contig has empty sequence i.e.
>NODE_780-0
>NODE_1882-0
KevinLam is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:55 PM.


Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2010, Jelsoft Enterprises Ltd.