harvest icon indicating copy to clipboard operation
harvest copied to clipboard

parsnp failed to generate phylogeny

Open haol04 opened this issue 7 years ago • 2 comments

I am trying to run parsnp for 330 bacteria genome (genome size ~2mb). The alignment seems to work properly. but it failed in the last step, the error message is pasted below. Your help would be greatly appreciated! Li

-bash-4.1$ ./parsnp -g ./ref/2603V-R_NC004116.gb -d GBS_WGS_2017/ -c -C -p 32

<<Parsnp started>>

-->Reading Genome (asm, fasta) files from GBS_WGS_2017/.. |->[OK] -->Reading Genbank file(s) for reference (.gbk) ./ref/2603V-R_NC004116.gb.. |->[OK] -->Running Parsnp multi-MUM search and libMUSCLE aligner.. |->[OK] -->Running PhiPack on LCBs to detect recombination.. |->[SKIP] -->Reconstructing core genome phylogeny.. |->[OK] -->Creating Gingr input file.. ERROR The following command failed:

/tmp/_MEI2r348L/harvest --midpoint-reroot -u -q -i /afs/grid.pfizer.com/vaccine/ngs/Tools/Parsnp-Linux64-v1.2/P_2017_03_16_161557720869/parsnp.ggr -o /afs/grid.pfizer.com/vaccine/ngs/Tools/Parsnp-Linux64-v1.2/P_2017_03_16_161557720869/parsnp.ggr -n /afs/grid.pfizer.com/vaccine/ngs/Tools/Parsnp-Linux64-v1.2/P_2017_03_16_161557720869/parsnp.tree Please veryify input data and restart Parsnp. If the problem persists please contact the Parsnp development team.

haol04 avatar Mar 16 '17 20:03 haol04