Atelier Phylogénomique Blast
From silico.biotoul.fr
(Difference between revisions)
m (→make blast database) |
m |
||
(5 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
==Liens== | ==Liens== | ||
- | *[http://silico.biotoul.fr/p/Atelier_Phylog%C3%A9nomique Atelier de Phylogénomique] | + | *[http://silico.biotoul.fr/p/Atelier_Phylog%C3%A9nomique#Pr.C3.A9liminaires Atelier de Phylogénomique Blast All-vs-All] |
+ | |||
==Préliminaires== | ==Préliminaires== | ||
<pre style="color:red;white-space: pre-wrap;white-space: -moz-pre-wrap;white-space: -pre-wrap;white-space: -o-pre-wrap"> | <pre style="color:red;white-space: pre-wrap;white-space: -moz-pre-wrap;white-space: -pre-wrap;white-space: -o-pre-wrap"> | ||
Line 27: | Line 28: | ||
MSK | MSK | ||
+ | <!-- | ||
<syntaxhighlight lang="bash"> | <syntaxhighlight lang="bash"> | ||
for i in ~/work/Prochlorococcus/peptide/*.faa; | for i in ~/work/Prochlorococcus/peptide/*.faa; | ||
Line 38: | Line 40: | ||
<pre style="color:blue;white-space: pre-wrap;white-space: -moz-pre-wrap;white-space: -pre-wrap;white-space: -o-pre-wrap"> | <pre style="color:blue;white-space: pre-wrap;white-space: -moz-pre-wrap;white-space: -pre-wrap;white-space: -o-pre-wrap"> | ||
sarray -J mkdb -o %j.out -e %j.err -t 01:00:00 --cpus-per-task=1 makeblastdb.sh | sarray -J mkdb -o %j.out -e %j.err -t 01:00:00 --cpus-per-task=1 makeblastdb.sh | ||
- | squeue -l -u | + | squeue -l -u $USER |
+ | |||
+ | ll /home/yquentin/work/Prochlorococcus/peptide/*.faa.* | ||
</pre> | </pre> | ||
+ | --> | ||
<!-- | <!-- | ||
====Intra genomes==== | ====Intra genomes==== | ||
Line 48: | Line 53: | ||
~/work/scripts/blastp_intra.pl | ~/work/scripts/blastp_intra.pl | ||
- | squeue -l -u | + | squeue -l -u $USER |
ls BlastP | wc -l | ls BlastP | wc -l | ||
Line 77: | Line 82: | ||
MSK | MSK | ||
+ | <!-- | ||
<syntaxhighlight lang="bash"> | <syntaxhighlight lang="bash"> | ||
evalue=1e-5 | evalue=1e-5 | ||
Line 96: | Line 102: | ||
<pre style="color:blue;white-space: pre-wrap;white-space: -moz-pre-wrap;white-space: -pre-wrap;white-space: -o-pre-wrap"> | <pre style="color:blue;white-space: pre-wrap;white-space: -moz-pre-wrap;white-space: -pre-wrap;white-space: -o-pre-wrap"> | ||
sarray -J mkdb -o %j.out -e %j.err -t 01:00:00 --cpus-per-task=1 blast_allall.sh | sarray -J mkdb -o %j.out -e %j.err -t 01:00:00 --cpus-per-task=1 blast_allall.sh | ||
- | squeue -l -u | + | squeue -l -u $USER |
</pre> | </pre> | ||
+ | --> | ||
<!-- | <!-- | ||
Pour les blastp inter génomes, nous allons utiliser un script similaire au précédent, mais avec une double boucle (sur -query et -db). | Pour les blastp inter génomes, nous allons utiliser un script similaire au précédent, mais avec une double boucle (sur -query et -db). | ||
Line 104: | Line 111: | ||
~/work/scripts/blastp_inter.pl | ~/work/scripts/blastp_inter.pl | ||
- | squeue -l -u | + | squeue -l -u $USER |
ls BlastP | wc -l | ls BlastP | wc -l | ||
Line 114: | Line 121: | ||
Combien de fichiers attendez-vous? | Combien de fichiers attendez-vous? | ||
</pre> | </pre> | ||
+ | |||
+ | ---- | ||
+ | *[http://silico.biotoul.fr/p/Atelier_Phylog%C3%A9nomique#Pr.C3.A9liminaires Atelier de Phylogénomique Blast All-vs-All] |
Revision as of 10:20, 26 November 2021
Contents |
Liens
Préliminaires
Question 1.5: Selon vous qu'est-ce qui guide le choix du type de séquences à utiliser dans les comparaisons (peptides ou nucléotidiques)?
Blast All-All
Nous allons utiliser NCBI_Blast+.
Nous allons copier les fichiers peptides dans un répertoire de travail:
mkdir -p ~/work/Prochlorococcus/peptide cp ~/work/Prochlorococcus/prokka/Aaa*/*.faa ~/work/Prochlorococcus/peptide/. ls -l ~/work/Prochlorococcus/peptide
make blast database
Exemple:
search_module blast module load bioinfo/ncbi-blast-2.7.1+ srun -n1 -l makeblastdb -in Aaaa.faa -dbtype prot
Vous allez procéder comme précédemment, avec un script donné à sarray, pour réaliser le makeblastdb sur tous les fichiers.
MSK
Paire de genomes
mkdir -p ~/work/Prochlorococcus/BlastP
Exemple :
module load bioinfo/ncbi-blast-2.7.1+ srun -n1 -l blastp -query ~/work/Prochlorococcus/peptide/Aaaa.faa -db ~/work/Prochlorococcus/peptide/Aaaa.faa -seg yes -dbsize 100000000 -evalue 1e-5 -outfmt 6 -num_threads 1 -out ~/work/Prochlorococcus/BlastP/Aaaa_Aaaa.tab
Question 1.6: Explicitez et justifiez les paramètres de blast utilisés dans votre script.
Boucle sur les génomes
Utilisez sarray pour réaliser les blast en toutes les paires de génomes.
MSK vérifiez que les fichiers de sorties de blast sont présents et non vides.
Question 1.7: Combien de fichiers attendez-vous?