I have used metavelvet before and it produced output from velvetg velveth and metavelvetg, however with the script below it only produced the meta-velvetg.contigs.fa file
That was the first issue and secondly I used metavelvet to assemble part of a metagenomic dataset which had been aligned to a group of genomes. The reads that aligned were about 500Mb, however the assembly is only 35Mb. This seems far too small or is that acceptable?
Sorry if these are silly questions.
Thanks,
Tom
Code:
velveth out-dir 51 -fasta NZ_bowtie_aligned.done /cm/shared/apps/velvet/1.2.10/velvetg out-dir -exp_cov auto /ibers/ernie/home/thh32/bin/meta-velvetg out-dir -amos_file yes | tee logfile
Sorry if these are silly questions.
Thanks,
Tom
Comment