Pour utiliser tous les modules de metawrap, il faut toutes les banques et indexes qu'ils indiquent dans les instructions (c'est-à-dire CheckM, Kraken1, Kraken2, NCBI_nt, NCBI_tax, Indexed hg38)
Nous utilisons principalement CheckM et Indexed hg38 s'il n'est pas possible de tout installer
J'ai rencontré un problème avec checkm, est-ce que checkm data setRoot a été exécuté ? Malheureusement je ne semble pas avoir la permission de le faire moi-même
Selon le tutoriel :
Now you need to tell CheckM where to find this data before running anything:
checkm data setRoot # CheckM will prompt to to chose your storage location
On newer versions of CheckM, you would run:
checkm data setRoot /path/to/your/dir/MY_CHECKM_FOLDER
J'ai malheureusement toujours le même problème (checkm ne semble pas pouvoir accéder aux profils hmm, qui sont cependant installés sous /shared/bank/checkm/2015-01-16)
Quelque chose ne semble pas aller avec checkm lui-même, j'obtiens le message d'erreur suivant:
Il y a un problème avec la base de données Kraken:
message d'erreur:
The folder /shared/home/emendes/KRAKEN_DB doesnt exist. Please consult the metaWRAP *****
***** database guite to download and build the KRAKEN database
Il n'y a pas encore de base de données installées pour Kraken.
Il convient donc pour l'instant de télécharger ou construire la base qui vous convient dans votre home (~/KRAKEN_DB).
il semble qu'il y ait aussi un problème avec la banque hg38
/shared/home/rpopall/BMTAGGER_DB/hg38.bitmask file doesnt exist. Please configure your *****bmtagger genome index *****
Selon le tutoriel :
Making host genome index for bmtagger
If you want to remove human (see end of instrutions for non-human hosts) reads from tour sequencing in the READ_QC module, you will need to dowlnoad and index the human genome. See the official bmtagger manual for detailed instructions: https://www.hmpdacc.org/hmp/doc/HumanSequenceRemoval_SOP.pdf
First, lets download and merge the human genome hg38:
Now lets index the human genome. Note that the file names of the indeces must be exactly as specified for metaWRAP to recognize them! Also note that indexing will take considerable memory and time (here I pass 100GB of RAM as a -M parameter).
Note: metaWRAP looks for files hg38.bitmask and hg38.srprism - make sure they are names exactly like this. Done! Now dont forget to specify the BMTAGGER_DB variable in the config-metawrap file! Run which config-metawrap to find it.
Idem, il n'y a pas de base pour bmtagger.
Il convient donc pour l'instant de la créer (ou simplement indexer une base existante type /shared/bank/homo_sapiens/hg38 avec bmtool) dans votre home (~/BMTAGGER_DB).
Merci ! J'ai crée la base de données dans ~/BMTAGGER_INDEX, mais je n'arrive pas à accéder au fichier config-metawrap pour spécifier BMTAGGER_DB=/path/to/your/index/BMTAGGER_INDEX