I've been making genome assemblies for some Oxford Nanopore Technology (ONT) sequencing data using the Dragonflye package by Robert A. Petit III.
It was super easy to run!
Here's how I ran it on the Sanger compute farm:
First I found the version of Dragonflye on the farm:
% module avail -t | grep -i dragon
Then I loaded it:
% module load dragonflye/1.2.1
Then I assembled sequence reads for a Vibrio cholerae isolate into an assembly using Dragonflye:
% dragonflye --reads SRR31024125_1.fastq.gz --outdir SRR31024125_1.fastq_dragonflye --gsize 4000000
where SRR31024125_1.fastq.gz was my input fastq file of ONT reads,
SRR31024125_1.fastq_dragonflye was the name that I wanted to give to the output directory,
--gsize 4000000 specifies that the Vibrio cholerae genome is about 4.0 Mbase.
The output file was called SRR31024125_1.fastq_dragonflye/contigs.fa.
It took about 20 minutes to make the assembly. The input file of ONT reads was about 93 Megabytes (SRR31024125_1.fastq.gz).
No comments:
Post a Comment