Bonjour,
Je travaille sur le séquençage d'individus du moustique Anopheles gambiae. Mon jeu de données prolonge une précédente étude du consortium 1000 génomes moustiques (AG1000G) dont le SNP calling a été réalisé à l'aide de GATK 3 entre 2014 et 2017.
De manière à me servir du jeu de données du 1000 génomes comme référence, je souhaiterais pouvoir utiliser GATK 3.8.
Malheureusement, ce dernier n'est pas disponible en module, ni en environnement conda partagé.
J'ai tenté de l'installer à partir d'une installation de miniconda3 que j'ai réalisée dans mon home/ (conda install -c bioconda gatk) mais ça ne fonctionne pas.
Serait-il possible d'installer GATK 3.8 en module ? Si nécessaire, j'ai à ma disposition un jar (package-archive-gatk-GenomeAnalysisTK-3.8-0-ge9d806836.tar.bz2) dans mon répertoire aedes_amplicon.
En vous remerciant par avance pour l'aide précieuse apportée,
Cordialement,
Bonjour,
Demande d'installation en cours:
En revanche, même si j'entends la problématique de prolonger l'étude précédente, je ne peux que vous inviter à utiliser une version plus récente de GATK. De nombreuses améliorations et corrections sont apportés.
gatk3 n'est d'ailleurs plus maintenu et l'équipe recommande d’utiliser gatk4.
Bonjour,
Nous utilisons d'habitude HaplotypeCaller et sommes bien conscients des améliorations apportées par celui-ci, mais cette étude nécessite de découvrir les variants selon la même méthode que celle utilisée en 2014 de manière à ne pas introduire de biais. Il s'agit d'un important ajout de populations de moustiques au projet AG1000G, sur des régions qui n'avaient pas encore été étudiées et qui sont extrêmement intéressantes.
En tous cas, merci pour votre réponse.
Bonjour,
Le module est maintenant disponible. Comme d'habitude:
module load gatk/3.8
Bonne journée
Bonjour.
J'ai aussi un besoin momentanné de faire tourner certains outils de gatk3. Dans mon cas, après 'module load gatk/3.8', 'srun gatk DepthOfCoverage' me donne la sortie suivante : 'slurmstepd: error: execve(): gatk: No such file or directory'
Vous avez une idée du problème?
-ben
Bonjour,
La bonne commande de gatk 3.8 est
gatk3 -T DepthOfCoverage
Cordialement,
1 « J'aime »