... | ... | @@ -16,7 +16,6 @@ I added a new script in `/bioinf/transfer/marmic_NGS2022/software` called `FastA |
|
|
| U | | | |
|
|
|
| Q | | | |
|
|
|
|
|
|
|
|
|
For the assembly of the metagenomic samples we are going to use megahit. Be sure to also check out other assemblers such as IDBA-ud and SPAdes.
|
|
|
|
|
|
First, make sure it is available. It should be included in the marmic2022 conda environment.
|
... | ... | @@ -38,15 +37,15 @@ $ megahit -m 0.2 -1 G.1.fa -2 G.2.fa -o G-{your-assigned-kmer} -t 8 --min-contig |
|
|
| Student | k-mer | `#` contigs (>500bp) | N50 and L50 | Longest contig (bp) | Total length (bp) |
|
|
|
|---------|-------|----------------------|-------------|---------------------|-------------------|
|
|
|
| | | | | | |
|
|
|
| | 33 | | | | |
|
|
|
| Leonard | 33 | | | | |
|
|
|
| | 37 | | | | |
|
|
|
| | 43 | | | | |
|
|
|
| | 47 | | | | |
|
|
|
| | 53 | | | | |
|
|
|
| | 57 | | | | |
|
|
|
| | 63 | | | | |
|
|
|
| | 67 | | | | |
|
|
|
| | 73 | | | | |
|
|
|
| Louison | 43 | | | | |
|
|
|
| Ava | 47 | | | | |
|
|
|
| Ella | 53 | | | | |
|
|
|
| Kinga | 57 | | | | |
|
|
|
| Nicole | 63 | | | | |
|
|
|
| Margret | 67 | | | | |
|
|
|
| Mahum | 73 | | | | |
|
|
|
|
|
|
The assembly of samples using single k-mer sizes should take <span dir="">\~</span> 10 minutes. Once you have an assembly, you can use the stats.sh program from bbmap (/bioinf/software/bbmap) on the contigs file to get some basic information about it.
|
|
|
|
... | ... | |