De-novo Genome Assembly

Last updated: Feb, 2022

Short-read assembly

Codes
  • SPAdes

    spades.py --pe1-1 <sample.R1.fastq> --pe1-2 <sample.R2.fastq> -o <output_dir>

  • Velvet (Velvet Optimiser)

    VelvetOptimiser.pl -f "-shortPaired -fastq -separate sample.R1.fastq sample.R2.fastq" \
                       -s <starting_hash_value> \
                       -e <end_hash_value> \
                       -d <output_dir>

  • SOAPdenovo

    soapdenovo all -s <config_file> -K <K-mer_size> -R -o <graph_prefix>


  • References
  • SPAdes
  • Velvet
  • SOAPdenovo
  • Long-read assembly

    Codes
  • Canu

    canu -p <assembly-prefix> \
         -d <assembly-directory> \
         genomeSize=<genome_size> \
         -nanopore <sample_guppy.fastq>

  • Necat

    necat.pl config <sample_config.txt>
    necat.pl correct <sample_config.txt>
    necat.pl assemble <sample_config.txt>
    necat.pl bridge <sample_config.txt>


  • References
  • Canu
  • Necat
  • Assembly Evaluation

    Codes
  • QUAST: assembly statistics

    quast.py <scaffolds.fasta> -o <output>

  • BUSCO

    busco -i <SEQUENCE_FILE> -l <LINEAGE> -o <OUTPUT_NAME> -m <genome>


  • References
  • QUAST
  • BUSCO