This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tool: AliNe (Alignment in Nextflow)

AliNe is a pipeline written in Nextflow that focuses on the alignment process. It aligns reads to a reference genome using the tools of your choice among 16 different alignment tools (more to come).

AliNe supports short reads (both paired-end and single-end) as well as long reads generated by PacBio and Oxford Nanopore Technologies (ONT). It is designed to minimize user inputs and avoid common parameter mistakes ( e.g. scoring system, strandedness, orientation).

To try it you need Nextflow + Docker or Singularity.
Usage:

nextflow run Juke34/AliNe -r v1.1.0 -profile docker \
  --reads https://github.com/Juke34/AliNe/raw/refs/heads/main/test/illumina/yeast_R1.fastq.gz \
  --genome https://raw.githubusercontent.com/Juke34/AliNe/refs/heads/main/test/yeast.fa \
  --read_type short_single \
  --aligner bbmap,bowtie,bowtie2,bwaaln,bwamem,bwasw,graphmap2,hisat2,kallisto,minimap2,ngmlr,nucmer,star,subread,sublong \
  --trimming_fastp \
  --fastqc \
  --samtools_stats

It has been published in JOSS: publication here.
Full documentation here.

Feedbacks and contributions are Welcome !

nextflow alignment

I standardize the output for all aligner in sorted bam file. It is why it is specific to “read” alignment. I don’t think we can output bam with last. It exists already the pairgenomealign nf-core pipeline for whole genome alignment.

Indeed, LAST outputs to MAF. However it also provides a maf-convert command that does SAM, whose output can be piped in samtools view for BAM or CRAM output. LAST has been used to align Nanopore RNA reads to genomes, for instance.

Then yes it is feasible. Thx for the information

Lastis now present in AliNe v1.3.0
Before that bwa-mem2 and salmon have been added in AliNe v1.2.0

1 answer

Release AliNe v1.6.0 !

  • cram output available
  • parameter to filter unmapped reads
  • Dragmap aligner added
  • tsv comparaison table made automatically at then end when statistics activated

One thing that would make the documentation even more useful is to generate a bunch of stats from data across different organisms and platforms. And published that with the examples/docs.

In essence, the documentation could serve as a reference for choosing an aligner.

How long will it take to run aligner A on data obtained from organism B with coverage C and sequencing platform X. How good is that alignment? Are there any notable differences between aligners? Can I run that test myself - most of the time hunting down the data is the tedious part.

Currently, we are always guessing and usually go with past experience...

Log in to answer this question.