Magic-BLAST incorporates within the NCBI BLAST code framework ideas developed in the NCBI Magic pipeline, in particular hit extensions by local walk and jump (http://www.ncbi.nlm.nih.gov/pubmed/26109056), and recursive clipping of mismatches near the edges of the reads, which avoids accumulating artefactual mismatches near splice sites and is needed to distinguish short indels from substitutions near the edges.

Address of the bookmark: https://ncbi.github.io/magicblast/

Ref https://www.lncrnablog.com/lncpipe-a-nextflow-based-pipeline-for-identification-and-analysis-of-long-non-coding-rnas-from-rna-seq-data/

Address of the bookmark: https://github.com/likelet/LncPipe

ORLANDO, USA, Oct 17, 2017/ PRNewswire/

Strand NGS now supports large scale RNA- and small-RNA-Seq and Unique Molecular Identifiers (UMIs) for DNA-, RNA-, and small-RNA-Seq.

Strand Life Sciences announced the latest version release of its bioinformatics flagship product, Strand NGS, at the Annual Meeting of the American Society of Human Genetics today. Two major themes in Strand NGS v3.1 address recent challenges in next generation sequencing (NGS).

The first theme is large-scale RNA-Seq data analysis. Current cross-cohort RNA- and small-RNA-Seq studies span tens of replicates and batches across hundreds of samples, sometimes conducted across several different institutions. For such studies, Strand NGS v3.1 includes confounding variable analysis to eliminate technical effects, including batch effects; the t-SNE plot; profile and heat-map plots of gene-body coverage; and several other notable visual enhancements.

The second new feature is support for Unique Molecular Identifiers, or UMIs, for DNA-, RNA- and small-RNA-Seq. UMI support in Strand NGS is end-to-end, spanning alignment to variant calling in DNA-Seq, and alignment to quantification in RNA- and small-RNA-Seq. The Bioo Scientific, Qiagen, and Rubicon UMI protocols are natively supported, and an intuitive interface allows the specification of custom UMI protocols.

“For liquid biopsies and low-grade FFPE samples, UMI support in DNA-Seq enables the detection of somatic variants at low concentrations. In RNA-Seq, large-scale and UMI support can be used in single-cell-based studies that reveal tumor-cell heterogeneity, even at low concentrations”, says Dr. Vamsi Veeramachaneni, Chief Scientific Officer, Strand Life Sciences.

“At Strand, we are continuously working towards improving the accuracy and efficiency of NGS data analysis. Customers can look forward to Strand NGS becoming available on the cloud in the near future”, says Dr. Ramesh Hariharan, Chief Executive Officer, Strand Life Sciences.

Visit Strand Life Sciences at ASHG booth #1017 to know more about Strand NGS v3.1 and other products and service offerings from Strand Life Sciences. Click here to access detailed agenda and v3.1 release notes.

About Strand Life Sciences

Strand Life Sciences is a premier life science informatics innovation company. Founded in 2000, Strand is a leader in technology innovations for healthcare using genomics. By enhancing sequence-based diagnostics and clinical genomic data interpretation using a strong foundation of computational, scientific, and medical expertise, Strand is bringing individualized medicine to the world. To know more, visit www.strandls.com

Address of the bookmark: http://www.sthda.com/english/wiki/rna-seq-differential-expression-work-flow-using-deseq2

What Are Kallisto and Salmon?

At their core, both Kallisto and Salmon are tools for quantifying transcript abundance from RNA-Seq reads. They bypass traditional alignment-based methods, replacing them with pseudoalignment or quasi-mapping, which drastically speeds up the process.

• Kallisto was developed by Lior Pachter’s lab and introduced the concept of pseudoalignment using a de Bruijn graph.
• Salmon, developed by Rob Patro’s group, builds on this idea with quasi-mapping and offers additional features like advanced bias correction.

Head-to-Head Comparison

1. Algorithm

• Kallisto uses pseudoalignment, focusing on matching k-mers from reads to a transcriptome index.
• Salmon uses quasi-mapping, which adds more flexibility and can also work with aligned reads (BAM files).

2. Input and Flexibility

• Kallisto works with raw FASTQ reads and requires a custom transcriptome index.
• Salmon accepts FASTQ or pre-aligned BAM files, giving you more workflow options.

3. Bias Correction

One of Salmon’s major advantages is its sophisticated bias correction system. It corrects for:

• Sequence-specific bias
• Positional bias
• GC-content bias

Kallisto offers basic sequence bias correction but lacks the comprehensive models found in Salmon.

4. Speed and Resources

• Kallisto is blazing fast and slightly more memory-efficient.
• Salmon is still very fast, but the added features can come at a small computational cost.

5. Output and Downstream Analysis

• Both tools provide transcript-level quantifications and support bootstrapping for variance estimation.
• Salmon can also summarize counts at the gene level if provided with a mapping file (--geneMap).
• Kallisto integrates seamlessly with Sleuth for differential expression analysis.
• Salmon works well with tximport, DESeq2, edgeR, and other Bioconductor tools.

Choosing the Right Tool

Example Command Lines

Kallisto (paired-end):

kallisto quant -i transcriptome.idx -o output -b 100 sample_R1.fastq sample_R2.fastq

Salmon (paired-end, bias correction):

salmon quant -i salmon_index -l A -1 sample_R1.fastq -2 sample_R2.fastq \
  -p 8 --validateMappings --seqBias --gcBias -o output

Conclusion

Both Kallisto and Salmon are exceptional tools that have transformed RNA-Seq analysis. Your choice largely depends on your priorities—whether it's speed, accuracy, flexibility, or compatibility with downstream tools.

For many users, Salmon offers a more complete and flexible solution, especially when bias correction and gene-level outputs are essential. However, Kallisto remains a favorite for quick, accurate quantification, especially when paired with the Sleuth pipeline.

Qualifications

PhD or Post-Graduated Master degree is required. Successful candidates will have some expertise in data analysis of NGS data by using methods running in R and UNIX environment. Familiarity with genome databases and browsers is required.

Application

Candidates should send a CV and a brief personal statement focusing on their skills and interests related to the research project.

Contacts

Start date: 1° April 2014
Salary on grant: 25,000 euros per year.
Contact Person (Referent): Mario Capasso
Ref. Email: mario.capasso@unina.it and achille.iolascon@unina.it
Tel: +39 081 3737889

The closing date for application will be 26 June 2015.

Address of the bookmark: http://omega.omicsbio.org/

So far miniasm is in early development stage. It has only been tested on a dozen of PacBio and Oxford Nanopore (ONT) bacterial data sets. Including the mapping step, it takes about 3 minutes to assemble a bacterial genome. Under the default setting, miniasm assembles 9 out of 12 PacBio datasets and 3 out of 4 ONT datasets into a single contig. The 12 PacBio data sets are PacBio E. coli sample, ERS473430, ERS544009, ERS554120, ERS605484, ERS617393, ERS646601, ERS659581, ERS670327, ERS685285, ERS743109 and a deprecated PacBio E. coli data set. ONT data are acquired from the Loman Lab.

For a C. elegans PacBio data set (only 40X are used, not the whole dataset), miniasm finishes the assembly, including reads overlapping, in ~10 minutes with 16 CPUs. The total assembly size is 105Mb; the N50 is 1.94Mb. In comparison, the HGAP3produces a 104Mb assembly with N50 1.61Mb. This dotter plot gives a global view of the miniasm assembly (on the X axis) and the HGAP3 assembly (on Y). They are broadly comparable. Of course, the HGAP3 consensus sequences are much more accurate. In addition, on the whole data set (assembled in ~30 min), the miniasm N50 is reduced to 1.79Mb. Miniasm still needs improvements.

Miniasm confirms that at least for high-coverage bacterial genomes, it is possible to generate long contigs from raw PacBio or ONT reads without error correction. It also shows that minimap can be used as a read overlapper, even though it is probably not as sensitive as the more sophisticated overlapers such as MHAP and DALIGNER. Coupled with long-read error correctors and consensus tools, miniasm may also be useful to produce high-quality assemblies.

Minimap and miniasm are ultrafast tools for (i) mapping and (ii) assembly. Designed for long, noisy reads, they do not have a correction or consensus step, and therefore the resulting assemblies are contiguous (i.e. long) but very noisy (i.e. full of errors)

We start with an all against all comparison:

minimap -Sw5 -L100 -m0 -t8 reads.fq reads.fq | gzip -1 > reads.paf.gz

Then we can assemble

miniasm -f reads.fq reads.paf.gz > reads.gfa

Convert GFA to FASTA:

awk '/^S/{print ">"$2"\n"$3}' reads.gfa | fold > reads.fa

And then count how many contigs:

grep ">" reads.fa | wc -l

# Download sample PacBio from the PBcR website
wget -O- http://www.cbcb.umd.edu/software/PBcR/data/selfSampleData.tar.gz | tar zxf -
ln -s selfSampleData/pacbio_filtered.fastq reads.fq
# Install minimap and miniasm (requiring gcc and zlib)
git clone https://github.com/lh3/minimap && (cd minimap && make)
git clone https://github.com/lh3/miniasm && (cd miniasm && make)
# Overlap
minimap/minimap -Sw5 -L100 -m0 -t8 reads.fq reads.fq | gzip -1 > reads.paf.gz
# Layout
miniasm/miniasm -f reads.fq reads.paf.gz > reads.gfa

Address of the bookmark: https://github.com/lh3/miniasm

Address of the bookmark: https://github.com/marbl/MashMap