Has anyone assembled 454 data with consed package version 19? I'm having some issues with reading of the .sff files and wondering if anyone has completed an assembly of 454 read data (not using Roche software produced .ace files). I'm using "add454Reads.perl reference.ace sff.fof reference.fa", where the fof specifies the sff files to assembly, but although the script runs, I get an error "doesn't existile /shared/BNFinal/mapping/consed/sff_dir/FPDLD6P02.sff", and the 454 reads are not brought into the assembly; it basically assembles with only the reference sequence. I also tried specifying just the names of reads that I want to bring in to the assembly; these exist in the sff files in the corresponding sff_dir directory to the edit_dir I am running the assembly from. Someone mentioned needing to update phrap, which I will look into, but any other thoughts on this (I'm sure I'm missing something obvious)?
Thanks,
Liz
Thanks,
Liz
Comment