SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   solid_denovo_postprocessor.pl removes velvet contig sequence. (http://seqanswers.com/forums/showthread.php?t=3986)

KevinLam 02-09-2010 06:04 PM

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


All times are GMT -8. The time now is 09:37 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.