cpockrandt / PhyloCSFpp

PhyloCSF++ computes PhyloCSF tracks for whole-genome multiple sequence alignments, scores single MSA, annotates CDS features in GFF/GTF files with PhyloCSF and confidence scores.
Other
30 stars 4 forks source link

Cannot generate scores for hg38 maf #12

Open asumann opened 1 year ago

asumann commented 1 year ago

Hi,

I cannot run the program with a fasta alignment file, should I only use maf format? Also, if I use publicly available https://hgdownload.soe.ucsc.edu/goldenPath/hg38/multiz100way/maf/chr1_KI270707v1_random.maf.gz as maf and run phylocsf++ build-tracks 58mammals maf it cannot generate scores. It gives long something like this:

OpenBLAS Warning : Detect OpenMP Loop and this application may hang. Please rebuild the library with USE_OPENMP=1 option.
OpenBLAS Warning : Detect OpenMP Loop and this application may hang. Please rebuild the library with USE_OPENMP=1 option.
OpenBLAS Warning : Detect OpenMP Loop and this application may hang. Please rebuild the library with USE_OPENMP=1 option.
Done!
WARNING: gorilla in the model does not occur in alignment file(s). Check --species to select a subset (this affects the power/confidence track).
WARNING: squirrel_monkey in the model does not occur in alignment file(s). Check --species to select a subset (this affects the power/confidence track).
WARNING: bushbaby in the model does not occur in alignment file(s). Check --species to select a subset (this affects the power/confidence track).
WARNING: chinese_tree_shrew in the model does not occur in alignment file(s). Check --species to select a subset (this affects the power/confidence track).

I ignored the first warning and tried to correct the second warning by adding species mapping names, it still gives the same warnings. I wonder if you can replicate the error. Also can you add hg38.100way.maf file into example folder so that I can check file structures. I want to try mammals model with human maf, since running bird example gives no error and I cannot generate scores for human maf file above.

Thanks!