I followed the following directions for aligning 454 reads to an assembled rna-seq dataset.
newMapping projdir
addRun [options: -p, -lib, -mcf] [projdir] filedesc
removeRun [projdir] filedesc
setRef [options: -gref, -cref, -random] [projdir] refdesc
runProject [options] [projdir]
I got the following error when I tried to use the -cref option with setRef
my command:
setRef -cref DH1 /fullpathname/grossSilber.fa
error:
Error: File is not a fasta file of sequence: -cref
Is there a particular format for a cref fasta?
It worked when I setRef with no options, but of course runProject won't work with -cdna
Any helpful hints?
newMapping projdir
addRun [options: -p, -lib, -mcf] [projdir] filedesc
removeRun [projdir] filedesc
setRef [options: -gref, -cref, -random] [projdir] refdesc
runProject [options] [projdir]
I got the following error when I tried to use the -cref option with setRef
my command:
setRef -cref DH1 /fullpathname/grossSilber.fa
error:
Error: File is not a fasta file of sequence: -cref
Is there a particular format for a cref fasta?
It worked when I setRef with no options, but of course runProject won't work with -cdna
Any helpful hints?
Comment