Quick start

In brief, the steps involved in getting a Trycycler consensus assembly are:

• Prepare the input files:
  • Create multiple independent assemblies of the same isolate from a long-read set. Save these as assembles/*.fasta.
  • Save the long reads used to generate those assemblies as reads.fastq.gz.
• Run Trycycler cluster to group similar contigs together:
  • trycycler cluster --assemblies assemblies/*.fasta --reads reads.fastq.gz --out_dir trycycler
• Manually inspect the clusters and decide which are valid:
  • For this example, we'll assume cluster_001, cluster_007 and cluster_008 are the good clusters which represent replicons for which we want a consensus.
  • Delete all other cluster directories (will make it easier to glob for the good clusters with *).
• Run Trycycler align on each of the clusters:
  • trycycler align --reads reads.fastq.gz --cluster_dir trycycler/cluster_001
  • trycycler align --reads reads.fastq.gz --cluster_dir trycycler/cluster_007
  • trycycler align --reads reads.fastq.gz --cluster_dir trycycler/cluster_008
  • In order for these commands to complete, it may be necessary to delete or repair some of the cluster sequences.
• Run Trycycler MSA on each of the clusters:
  • trycycler msa --cluster_dir trycycler/cluster_001
  • trycycler msa --cluster_dir trycycler/cluster_007
  • trycycler msa --cluster_dir trycycler/cluster_008
• Run Trycycler partition to divide up the reads:
  • trycycler partition --reads reads.fastq --cluster_dirs trycycler/cluster_*
• Run Trycycler consensus to make a consensus sequence for each contig cluster:
  • trycycler consensus --cluster_dir trycycler/cluster_001
  • trycycler consensus --cluster_dir trycycler/cluster_007
  • trycycler consensus --cluster_dir trycycler/cluster_008
• Combine all consensus sequences into a single FASTA:
  • cat trycycler/cluster_*/7_final_consensus.fasta > trycycler/consensus.fasta

For more information, please look at the wiki pages for each of the steps involved.