Skip to content

Commit 770db00

Browse files
committed
added pan_genome to index
1 parent 315c271 commit 770db00

2 files changed

Lines changed: 6 additions & 5 deletions

File tree

README.md

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -8,9 +8,10 @@ bioinformatics subjects:
88
* [Introduction to Unix](unix/first_steps.md)
99
* [File formats](file_formats.md)
1010
* [Quality Control](qc.md)
11-
* [Variant Calling](variant_calling.md)
1211
* [Assembly](assembly.md)
1312
* [Annotation](annotation.md)
13+
* [Variant Calling](variant_calling.md)
14+
* [Pan-genome Analysis](pan-genome.md)
1415
* [RNA-Seq](rna_seq.md)
1516
* [16s Metabarcoding Analysis](16s.md)
1617
* [Whole Metagenome Sequencing](wms.md)

pan_genome.md

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -9,7 +9,7 @@ This tutorial is inspired from [Genome annotation and Pangenome Analysis](https:
99
We'll use six *Listeria monocytogenes* genomes in this tutorial.
1010

1111
```bash
12-
wget pangenome.tar.gz
12+
wget https://github.com/HadrienG/tutorials/blob/master/data/pangenome.tar.gz
1313
tar xzf pangenome.tar.gz
1414
cd pangenome
1515
```
@@ -18,7 +18,7 @@ These genomes correspond to isolates of *L. monocytogenes* reported in
1818

1919
> Xiangyu Deng, Adam M Phillippy, Zengxin Li, Steven L Salzberg and Wei Zhang. (2010) Probing the pan-genome of Listeria monocytogenes: new insights into intraspecific niche expansion and genomic diversification. doi:10.1186/1471-2164-11-500
2020
21-
The six genomes you downaloaded were selected based on their level of completeness (finished; contigs, etc) and their genotype (type I-III):
21+
The six genomes you downloaded were selected based on their level of completeness (finished; contigs, etc) and their genotype (type I-III):
2222

2323
Genome Assembly | Genome Accession | Genotype | Sequenced by | Status
2424
--- | --- | --- | --- | ---
@@ -56,11 +56,11 @@ Roary comes with a python script that allows you to generate a few plots to grap
5656
First, we need to generate a tree file from the alignment generated by Roary:
5757

5858
```
59-
raxmlHPC -m GTRGAMMA -p 12345 -s core_gene_alignment.aln -n core.tre
59+
FastTree –nt –gtr core_gene_alignment.aln > my_tree.newick
6060
```
6161

6262
Then we can plot the Roary results:
6363

6464
```
65-
roary_plots.py RAxML_bestTree.core.tre gene_presence_absence.csv
65+
roary_plots.py my_tree.newick gene_presence_absence.csv
6666
```

0 commit comments

Comments
 (0)