BOL: Related items

Pathway Analysis

Rahul Agarwal — Fri, 03 Oct 2014 08:51:13 -0500

Pathway Analysis is usually performed with aim to enrich the genes with their functional information and reveal the underlying biological mechanisms pursue by genes. Pathway Analysis is not only limited to what biological pathways a particular set of expressed genes follow but also to disclose the relationships between these genes. With availability of more genomics, transcriptomics and proteomics data, interactions between genes involve in multiple pathways become more clear and also relationships between the genes, their transcripts, and their gene products. However, existing tools and dbs mainly based on knowledge driven approach in which pathways will be identified by finding the correlation between the information in one of the pathway knowledge databases (KEGG,Reactome,Panther,BioCarta, Panther,GO,NCI,WikiPathways,etc) and gene expression result for a specific conditions for instance tumor, obesity , cold resistant crops/plants, etc.

Introductory Articles/ppt/sources:

http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1002375

http://bioinformatics.mdanderson.org/MicroarrayCourse/Lectures09/Pathway%20Analysis.pdf

http://gettinggeneticsdone.blogspot.de/2012/03/pathway-analysis-for-high-throughput.html

http://davetang.org/muse/tag/pathway/

https://www.biostars.org/p/42219/

http://bioinformatics.ca//files/public/Pathways_2014_Module4_v2.pdf

http://bioinformatics.ca//files/public/Pathways_2014_Module2.pdf

Impotant Database and Tools:

GeneMANIA, Cytoscape, IPA and Metacore (Commerical ), Pathway Commons, Reactome ,Panther, BioCyc, WikiPathways, Pathvisio, KEGG, NCI, Stringdb, Amigo, WebGestalt ,ConsensusPathDB ,GSEA,Blast2go

Popular R based tools:

Reactome.db, ReactomePA, ClusterProfiler, Gage, SPIA, topGO, Pathview,DOSE,GOStat

More:

http://www.bioconductor.org/help/search/index.html?q=Enrichment+analysis+

kallisto: a program for quantifying abundances of transcripts from bulk and single-cell RNA-Seq data

Jit — Mon, 07 Jan 2019 10:35:14 -0600

kallisto is a program for quantifying abundances of transcripts from bulk and single-cell RNA-Seq data, or more generally of target sequences using high-throughput sequencing reads. It is based on the novel idea of pseudoalignment for rapidly determining the compatibility of reads with targets, without the need for alignment. On benchmarks with standard RNA-Seq data, kallisto can quantify 30 million human reads in less than 3 minutes on a Mac desktop computer using only the read sequences and a transcriptome index that itself takes less than 10 minutes to build. Pseudoalignment of reads preserves the key information needed for quantification, and kallisto is therefore not only fast, but also as accurate as existing quantification tools. In fact, because the pseudoalignment procedure is robust to errors in the reads, in many benchmarks kallisto significantly outperforms existing tools. kallisto is described in detail in:

Nicolas L Bray, Harold Pimentel, Páll Melsted and Lior Pachter, Near-optimal probabilistic RNA-seq quantification, Nature Biotechnology 34, 525–527 (2016), doi:10.1038/nbt.3519

Address of the bookmark: https://pachterlab.github.io/kallisto/about

RNA-Seq Analysis: A Guide for Bioinformaticians

LEGE — Sat, 07 Dec 2024 22:22:24 -0600

RNA sequencing (RNA-Seq) has revolutionized transcriptomics, offering unprecedented insights into gene expression, splicing, and transcript diversity. For bioinformaticians, RNA-Seq analysis is a gateway to exploring the complexity of RNA biology and its implications in health and disease. This blog post provides an overview of RNA-Seq analysis, key computational steps, and tools for bioinformaticians eager to delve into this powerful technique.

What is RNA-Seq?

RNA-Seq is a next-generation sequencing (NGS) technology used to study the transcriptome—the complete set of RNA molecules in a cell. It quantifies gene expression, detects novel transcripts, and captures alternative splicing events with high sensitivity and resolution.

Workflow for RNA-Seq Analysis

RNA-Seq analysis involves several stages, each requiring computational tools and expertise.

1. Experimental Design and Data Acquisition

Before diving into analysis, bioinformaticians should consider:

Biological Replicates: Ensure statistical power to detect meaningful differences.
Sequencing Depth: Align sequencing depth to study objectives (e.g., higher depth for low-abundance transcripts).
Paired-End vs. Single-End: Paired-end sequencing provides more detailed information on transcript structure.

Once sequencing is complete, raw data is provided in FASTQ format, containing sequence reads and quality scores.

2. Quality Control and Preprocessing

Quality control (QC) ensures data integrity. Tools such as FastQC evaluate metrics like base quality, GC content, and adapter contamination.

Preprocessing Steps:

Trimming: Tools like Trimmomatic or Cutadapt remove low-quality bases and adapter sequences.
Filtering: Discard reads below a certain quality threshold or length.

3. Read Alignment

Reads are mapped to a reference genome or transcriptome to determine their origin. Alignment tools include:

HISAT2: Handles large genomes efficiently and supports spliced alignments.
STAR: High-speed aligner optimized for RNA-Seq.
Bowtie2: Suitable for short-read alignment.

Output: A SAM/BAM file containing aligned reads.

4. Transcript Assembly and Quantification

This step involves identifying transcripts and quantifying their expression levels. Tools used include:

StringTie: Assembles and quantifies transcripts from aligned reads.
Salmon/Kallisto: Perform pseudo-alignment for rapid and accurate quantification.

Expression levels are typically measured as TPM (transcripts per million) or FPKM (fragments per kilobase of transcript per million mapped reads).

5. Differential Expression Analysis

To identify genes with altered expression between conditions, bioinformaticians use tools such as:

DESeq2: Accounts for data normalization and variability.
edgeR: Handles overdispersed count data efficiently.
Limma-voom: Combines linear modeling with RNA-Seq count data.

The output includes a list of differentially expressed genes (DEGs) with statistical significance and fold-change values.

6. Functional Annotation and Pathway Analysis

Understanding the biological significance of DEGs involves:

Gene Ontology (GO) Analysis: Tools like DAVID or clusterProfiler categorize genes based on their biological functions.
Pathway Enrichment Analysis: Identifies pathways enriched in DEGs using tools like KEGG, Reactome, or GSEA.

7. Visualization

Visualizing results enhances interpretability. Common visualizations include:

Heatmaps: Show expression patterns across samples (e.g., pheatmap).
Volcano Plots: Highlight significant DEGs (e.g., ggplot2).
PCA/UMAP: Assess sample clustering and variability (e.g., Seurat).

Challenges in RNA-Seq Analysis

Batch Effects: Technical variability can confound biological signals. Combat this with normalization techniques or batch-correction tools like ComBat.
Low-Quality Samples: Poor-quality RNA impacts downstream analyses.
Computational Complexity: RNA-Seq generates massive datasets, requiring robust computing resources and optimized pipelines.

Key Tools and Resources

Bioconductor: A treasure trove of R packages for RNA-Seq analysis.
Galaxy: A web-based platform for running RNA-Seq workflows.
Nextflow/Snakemake: Workflow management tools to streamline analyses.

Applications of RNA-Seq

RNA-Seq is used in diverse research areas, including:

Cancer Transcriptomics: Identifying tumor-specific expression profiles.
Developmental Biology: Studying dynamic transcriptome changes.
Drug Discovery: Screening genes modulated by therapeutic compounds.

Conclusion

RNA-Seq analysis is a cornerstone of modern transcriptomics, offering bioinformaticians a versatile toolkit for unraveling gene expression and regulation. Mastering RNA-Seq workflows and tools empowers researchers to transform raw sequencing data into biological discoveries.

Whether you’re investigating disease mechanisms, exploring cellular pathways, or developing new therapeutics, RNA-Seq is a powerful ally in your bioinformatics arsenal.

GeneSpring webinar- Uncovering mechanisms of hepatotoxicity on 14 June at 8AM PST

Yeshodari — Tue, 23 May 2017 06:48:40 -0500

Uncovering Mechanisms of Hepatotoxicity for High Affinity Antisense Oligonucleotides – 3’ end RNA-seq Profiling Using GeneSpring GX

High affinity antisense oligonucleotides (ASOs) containing bicylic modifications (BNA) such as locked nucleic acid (LNA) or constrained ethyl (cEt) designed to induce target RNA cleavage have been shown to have enhanced potency along with a higher propensity to cause hepatotoxicity. In order to unravel the mechanism of this hepatotoxicity, we leveraged GeneSpring GX analysis software to analyze transcriptional profiles from the livers of mice treated with a panel of highly efficacious hepatotoxic or non-hepatotoxic LNA ASOs.

Speaker:
Sebastien A. Burel, PhD
Director, Nonclinical Development, Ionis Pharmaceuticals, California

Details:
14 June 2017, 8 AM PST

Register for this Webinar

Integrated mRNA and microRNA transcriptome analysis in Strand NGS

Strand — Tue, 04 Aug 2015 05:04:46 -0500

Using a nasopharyngeal carcinoma case study, this paper highlights the integrated transcriptome analysis capabilities of Strand NGS demonstrating the identification of miRNA – mRNA interactions in regulatory networks.
Read the application note on Integrated mRNA and microRNA transcriptome analysis in Strand NGS by Veena Hedatale and Rohit Gupta. For more information, please contact us

Calling narrow and broad peaks from ChIP-Seq data

Strand — Tue, 04 Aug 2015 05:06:03 -0500

Know about the state-of-the-art algorithms implemented in Strand NGS for detecting the binding sites of transcription factor (narrow peaks) and enriched regions of histone modification (broad peaks) from ChIP-Seq data.
Read the benchmarking study on Calling narrow and broad peaks from ChIP-Seq data in Strand NGS by Rohit Gupta and Anita Sathyanarayanan. For more information, please contact us

Webinar on Implications of Next Generation Sequencing in Molecular Diagnosis of Cancer on 28 Sep 2016

Strand — Thu, 01 Sep 2016 01:52:25 -0500

Abstract:

Genetic testing requires screening of the entire gene, which by conventional sequencing is time consuming and expensive. Next Generation Sequencing (NGS) based approaches increase the sensitivity of mutation detection, making it fast and cost-effective compared to the conventional tests performed in a reflex-testing mode. Strand NGS includes workflows with quality assessment and filter sections that do not require any manual intervention. Post-analytical workflows in Strand NGS allow users to execute sequence analysis with stringent filtering to eliminate false positive and low quality reads. This simplifies the analysis in large scale cohort settings, where every sample needs to be processed identically.

In this webinar we will discuss the implications of next generation sequencing based tests in multi-gene testing. We will also show how NGS based tests help to identify copy number variations, split read analysis and breakpoint identification. Finally, we will show a brief glimpse of Indian cohort data, where NGS based tests have shown improved mutation detection. In this webinar, we will present clinical case studies in on Hereditary Breast and Ovarian Cancer (HBOC) and Retinoblastoma patients to demonstrate how CNV analysis in Strand NGS enables researchers to detect and visualize copy number changes ranging from single exon to full gene.

Speaker:

Dr. Jaya Singh, Senior Scientist, Strand Life Sciences

Date: 28 September 2016

Session1: 2:30 PM IST

Session2: 10 PM IST

Register here: http://www.strand-ngs.com/webinar_registration

Surrogate Variable Analysis (SVA)

Jit — Thu, 30 Oct 2014 08:01:58 -0500

The sva package contains functions for removing batch effects and other unwanted variation in high-throughput experiment. Specifically, the sva package contains functions for the identifying and building surrogate variables for high-dimensional data sets. Surrogate variables are covariates constructed directly from high-dimensional data (like gene expression/RNA sequencing/methylation/brain imaging data) that can be used in subsequent analyses to adjust for unknown, unmodeled, or latent sources of noise. The sva package can be used to remove artifacts in three ways:

(1) identifying and estimating surrogate variables for unknown sources of variation in high-throughput experiments (Leek and Storey 2007 PLoS Genetics,2008 PNAS),

(2) directly removing known batch effects using ComBat (Johnson et al. 2007 Biostatistics) and

(3) removing batch effects with known control probes (Leek 2014 biorXiv).

Removing batch effects and using surrogate variables in differential expression analysis have been shown to reduce dependence, stabilize error rate estimates, and improve reproducibility, see (Leek and Storey 2007 PLoS Genetics, 2008 PNAS or Leek et al. 2011 Nat. Reviews Genetics).

More at http://www.bioconductor.org/packages/release/bioc/html/sva.html

Installing Trinity !

Rahul Nayak — Tue, 03 Jul 2018 04:31:15 -0500

Download this https://github.com/trinityrnaseq/trinityrnaseq/releases/tag/Trinity-v2.6.6

➜ cd trinityrnaseq-Trinity-v2.6.6

➜ trinityrnaseq-Trinity-v2.6.6 git:(master) ✗ make
Using gnu compiler for Inchworm and Chrysalis
cd Inchworm && (test -e configure || autoreconf) \
&& sh ./configure --prefix=`pwd` && make install
checking for a BSD-compatible install... /usr/bin/install -c
checking whether build environment is sane... yes
checking for a thread-safe mkdir -p... /bin/mkdir -p
checking for gawk... gawk
checking whether make sets $(MAKE)... yes
checking for g++... g++
checking whether the C++ compiler works... yes
checking for C++ compiler default output file name... a.out
checking for suffix of executables...
checking whether we are cross compiling... no
checking for suffix of object files... o
checking whether we are using the GNU C++ compiler... yes
checking whether g++ accepts -g... yes
checking for style of include used by make... GNU
checking dependency style of g++... gcc3
checking for library containing cos... none required
configure: creating ./config.status
config.status: creating Makefile
config.status: creating src/Makefile
config.status: creating config.h
config.status: config.h is unchanged
config.status: executing depfiles commands
make[1]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm'
Making install in src
make[2]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/src'
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT Fasta_entry.o -MD -MP -MF .deps/Fasta_entry.Tpo -c -o Fasta_entry.o Fasta_entry.cpp
mv -f .deps/Fasta_entry.Tpo .deps/Fasta_entry.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT IRKE_run.o -MD -MP -MF .deps/IRKE_run.Tpo -c -o IRKE_run.o IRKE_run.cpp
mv -f .deps/IRKE_run.Tpo .deps/IRKE_run.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT sequenceUtil.o -MD -MP -MF .deps/sequenceUtil.Tpo -c -o sequenceUtil.o sequenceUtil.cpp
mv -f .deps/sequenceUtil.Tpo .deps/sequenceUtil.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT IRKE.o -MD -MP -MF .deps/IRKE.Tpo -c -o IRKE.o IRKE.cpp
mv -f .deps/IRKE.Tpo .deps/IRKE.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT KmerCounter.o -MD -MP -MF .deps/KmerCounter.Tpo -c -o KmerCounter.o KmerCounter.cpp
mv -f .deps/KmerCounter.Tpo .deps/KmerCounter.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT string_util.o -MD -MP -MF .deps/string_util.Tpo -c -o string_util.o string_util.cpp
mv -f .deps/string_util.Tpo .deps/string_util.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT Fasta_reader.o -MD -MP -MF .deps/Fasta_reader.Tpo -c -o Fasta_reader.o Fasta_reader.cpp
mv -f .deps/Fasta_reader.Tpo .deps/Fasta_reader.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT stacktrace.o -MD -MP -MF .deps/stacktrace.Tpo -c -o stacktrace.o stacktrace.cpp
mv -f .deps/stacktrace.Tpo .deps/stacktrace.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT argProcessor.o -MD -MP -MF .deps/argProcessor.Tpo -c -o argProcessor.o argProcessor.cpp
mv -f .deps/argProcessor.Tpo .deps/argProcessor.Po
g++ -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -o inchworm Fasta_entry.o IRKE_run.o sequenceUtil.o IRKE.o KmerCounter.o string_util.o Fasta_reader.o stacktrace.o argProcessor.o
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT cigar_tweaker.o -MD -MP -MF .deps/cigar_tweaker.Tpo -c -o cigar_tweaker.o cigar_tweaker.cpp
mv -f .deps/cigar_tweaker.Tpo .deps/cigar_tweaker.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT SAM_entry.o -MD -MP -MF .deps/SAM_entry.Tpo -c -o SAM_entry.o SAM_entry.cpp
mv -f .deps/SAM_entry.Tpo .deps/SAM_entry.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT SAM_reader.o -MD -MP -MF .deps/SAM_reader.Tpo -c -o SAM_reader.o SAM_reader.cpp
mv -f .deps/SAM_reader.Tpo .deps/SAM_reader.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT Cigar.o -MD -MP -MF .deps/Cigar.Tpo -c -o Cigar.o Cigar.cpp
mv -f .deps/Cigar.Tpo .deps/Cigar.Po
g++ -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -o cigar_tweaker cigar_tweaker.o SAM_entry.o sequenceUtil.o SAM_reader.o Cigar.o string_util.o stacktrace.o Fasta_reader.o Fasta_entry.o
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT pull_reads_with_kmers.o -MD -MP -MF .deps/pull_reads_with_kmers.Tpo -c -o pull_reads_with_kmers.o pull_reads_with_kmers.cpp
mv -f .deps/pull_reads_with_kmers.Tpo .deps/pull_reads_with_kmers.Po
g++ -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -o pull_reads_with_kmers pull_reads_with_kmers.o Fasta_reader.o Fasta_entry.o sequenceUtil.o KmerCounter.o string_util.o stacktrace.o argProcessor.o
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT FastaToDeBruijn.o -MD -MP -MF .deps/FastaToDeBruijn.Tpo -c -o FastaToDeBruijn.o FastaToDeBruijn.cpp
FastaToDeBruijn.cpp: In function ‘void createGraphPerRecord(std::vector >, int, bool, ArgProcessor)’:
FastaToDeBruijn.cpp:199:39: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (seq_region.size() < kmer_length) { continue; } // can be e
^
mv -f .deps/FastaToDeBruijn.Tpo .deps/FastaToDeBruijn.Po
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT DeBruijnGraph.o -MD -MP -MF .deps/DeBruijnGraph.Tpo -c -o DeBruijnGraph.o DeBruijnGraph.cpp
mv -f .deps/DeBruijnGraph.Tpo .deps/DeBruijnGraph.Po
g++ -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -o FastaToDeBruijn FastaToDeBruijn.o argProcessor.o Fasta_reader.o Fasta_entry.o sequenceUtil.o string_util.o stacktrace.o DeBruijnGraph.o
g++ -DHAVE_CONFIG_H -I. -I.. -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -MT fastaToKmerCoverageStats.o -MD -MP -MF .deps/fastaToKmerCoverageStats.Tpo -c -o fastaToKmerCoverageStats.o fastaToKmerCoverageStats.cpp
fastaToKmerCoverageStats.cpp: In function ‘void populate_kmer_counter_from_reads(KmerCounter&, std::__cxx11::string&)’:
fastaToKmerCoverageStats.cpp:226:18: warning: unused variable ‘kmer_length’ [-Wunused-variable]
unsigned int kmer_length = kcounter.get_kmer_length();
^
mv -f .deps/fastaToKmerCoverageStats.Tpo .deps/fastaToKmerCoverageStats.Po
g++ -std=c++0x -pedantic -fopenmp -Wall -Wextra -Wno-deprecated -m64 -g -O2 -o fastaToKmerCoverageStats fastaToKmerCoverageStats.o argProcessor.o Fasta_reader.o Fasta_entry.o sequenceUtil.o string_util.o stacktrace.o KmerCounter.o
make[3]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/src'
test -z "/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/bin" || /bin/mkdir -p "/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/bin"
/usr/bin/install -c inchworm cigar_tweaker pull_reads_with_kmers FastaToDeBruijn fastaToKmerCoverageStats '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/bin'
make[3]: Nothing to be done for 'install-data-am'.
make[3]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/src'
make[2]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm/src'
make[2]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm'
make[3]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm'
make[3]: Nothing to be done for 'install-exec-am'.
make[3]: Nothing to be done for 'install-data-am'.
make[3]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm'
make[2]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm'
make[1]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Inchworm'
cd Chrysalis && make UNSUPPORTED=yes
make[1]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Chrysalis'
./MakeDepend.cc: In member function ‘bool makefile_builder::dependency::operator<(const makefile_builder::dependency&) const’:
./MakeDepend.cc:93:47: warning: suggest parentheses around ‘&&’ within ‘||’ [-Wparentheses]
dependent_ == other.dependent_ &&
^
./MakeDepend.cc: In member function ‘bool makefile_builder::subdir_and_target::operator<(const makefile_builder::subdir_and_target&) const’:
./MakeDepend.cc:130:41: warning: suggest parentheses around ‘&&’ within ‘||’ [-Wparentheses]
subdir_ == other.subdir_ &&
^
./MakeDepend.cc: In function ‘int main(int, char**)’:
./MakeDepend.cc:223:20: warning: deprecated conversion from string constant to ‘char*’ [-Wwrite-strings]
char * options = "Ad:g:u:fh";
^
./MakeDepend.cc: In member function ‘void makefile_builder::DumpDependencies(const string&)’:
./MakeDepend.cc:491:46: warning: suggest parentheses around ‘&&’ within ‘||’ [-Wparentheses]
if ( ! this->FileIsSource( dependent ) && ! this->FileIsHeader( dependent
^
./MakeDepend.cc:493:44: warning: suggest parentheses around ‘&&’ within ‘||’ [-Wparentheses]
this->FileIsHeader( dependent ) && this->FileIsSource( provider ) )
^
./MakeDepend.cc: In member function ‘bool makefile_builder::FileIsHeader(const string&) const’:
./MakeDepend.cc:1339:32: warning: suggest parentheses around ‘&&’ within ‘||’ [-Wparentheses]
return ( filename.size() > 2 &&
^
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:04 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./aligns/KmerAlignCore.cc -o obj/aligns/KmerAlignCore.o
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:06 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/AACodons.cc -o obj/analysis/AACodons.o
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:09 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/BreakTransByPairs.cc -o obj/analysis/BreakTransByPairs.o
./analysis/BreakTransByPairs.cc: In function ‘int main(int, char**)’:
./analysis/BreakTransByPairs.cc:105:14: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (i=0; i > DeBruijnGraph::get_candidate_weldmers(kmer_int_type_t, int)’:
./analysis/DeBruijnGraph.cc:473:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i = 0; i < left_extensions.size(); i++) {
^
./analysis/DeBruijnGraph.cc:476:27: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int j = 0; j < right_extensions.size(); j++) {
^
./analysis/DeBruijnGraph.cc: In member function ‘void DeBruijnGraph::recursively_construct_kmer_extensions(kmer_int_type_t, std::vector&, std::vector >&, char, std::map&, int)’:
./analysis/DeBruijnGraph.cc:513:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i = 0; i < adjacent_kmers.size(); i++) {
^
./analysis/DeBruijnGraph.cc:527:41: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (kmer_extension_chars.size() == flank_extension_length) {
^
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:30 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/sequenceUtil.cc -o obj/analysis/sequenceUtil.o
./analysis/sequenceUtil.cc: In function ‘bool contains_non_gatc(std::__cxx11::string)’:
./analysis/sequenceUtil.cc:33:20: warning: array subscript has type ‘char’ [-Wchar-subscripts]
if (_base_to_int[c] > 3)
^
./analysis/sequenceUtil.cc: In function ‘kmer_int_type_t kmer_to_intval(std::__cxx11::string)’:
./analysis/sequenceUtil.cc:264:26: warning: array subscript has type ‘char’ [-Wchar-subscripts]
int val = _base_to_int[c];
^
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:31 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/stacktrace.cc -o obj/analysis/stacktrace.o
/bin/rm -f lib_BubbleUpClustering_temp.a
ar -qc lib_BubbleUpClustering_temp.a obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/BubbleUpClustering.o obj/analysis/DNAVector.o obj/analysis/DeBruijnGraph.o obj/analysis/KmerTable.o obj/analysis/NonRedKmerTable.o obj/analysis/sequenceUtil.o obj/analysis/stacktrace.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./BubbleUpClustering
g++ -ggdb3 -fopenmp -o ./BubbleUpClustering obj/analysis/BubbleUpClustering.o -L. -lm -pthread -l_BubbleUpClustering_temp
/bin/rm lib_BubbleUpClustering_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:32 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/Chrysalis.cc -o obj/analysis/Chrysalis.o
./analysis/Chrysalis.cc: In function ‘int main(int, char**)’:
./analysis/Chrysalis.cc:384:28: warning: ignoring return value of ‘int system(const char*)’, declared with attribute warn_unused_result [-Wunused-result]
system(command.c_str());
^
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:35 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/TranscriptomeGraph.cc -o obj/analysis/TranscriptomeGraph.o
./analysis/TranscriptomeGraph.cc: In function ‘int TranscriptomeGraph(vecDNAVector&, FILE*, int, bool)’:
./analysis/TranscriptomeGraph.cc:668:20: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (i=0; i<=d.isize()-k; i++) {
^
/bin/rm -f lib_Chrysalis_temp.a
ar -qc lib_Chrysalis_temp.a obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/Chrysalis.o obj/analysis/DNAVector.o obj/analysis/TranscriptomeGraph.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./Chrysalis
g++ -ggdb3 -fopenmp -o ./Chrysalis obj/analysis/Chrysalis.o -L. -lm -pthread -l_Chrysalis_temp
/bin/rm lib_Chrysalis_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:39 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/CreateIwormFastaBundle.cc -o obj/analysis/CreateIwormFastaBundle.o
/bin/rm -f lib_CreateIwormFastaBundle_temp.a
ar -qc lib_CreateIwormFastaBundle_temp.a obj/analysis/AACodons.o obj/analysis/CreateIwormFastaBundle.o obj/analysis/DNAVector.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./CreateIwormFastaBundle
g++ -ggdb3 -fopenmp -o ./CreateIwormFastaBundle obj/analysis/CreateIwormFastaBundle.o -L. -lm -pthread -l_CreateIwormFastaBundle_temp
/bin/rm lib_CreateIwormFastaBundle_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:42 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/GraphFromFasta.cc -o obj/analysis/GraphFromFasta.o
./analysis/GraphFromFasta.cc: In function ‘bool SimpleHalves(const DNAVector&)’:
./analysis/GraphFromFasta.cc:255:15: warning: suggest parentheses around ‘&&’ within ‘||’ [-Wparentheses]
&&
^
./analysis/GraphFromFasta.cc: In function ‘void report_iworm_graph(std::map&, std::map&)’:
./analysis/GraphFromFasta.cc:631:31: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int j = 0; j < adjacent_nodes.size(); j++) {
^
./analysis/GraphFromFasta.cc: In function ‘int main(int, char**)’:
./analysis/GraphFromFasta.cc:1164:30: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (i=0; i bubble_up_cluster_growth(std::map&, std::map&)’:
./GraphFromFasta_MPI.cc:706:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i = 0; i < pool_vec.size(); i++) {
^
./GraphFromFasta_MPI.cc:727:27: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i = 0; i < pool_vec.size(); i++) {
^
./GraphFromFasta_MPI.cc:834:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int i = 0; i < pool_vec.size(); i++) {
^
./GraphFromFasta_MPI.cc: In function ‘void add_unclustered_iworm_contigs(svec&, vecDNAVector&)’:
./GraphFromFasta_MPI.cc:1144:30: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (size_t i = 0; i < p.size(); i++) {
^
./GraphFromFasta_MPI.cc: In function ‘int main(int, char**)’:
./GraphFromFasta_MPI.cc:1453:60: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (iworm_counter % 1000 == 0 || iworm_counter == dna.size()-1) {
^
./GraphFromFasta_MPI.cc:1458:24: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (j=0; j<=d.isize()-k; j++) {
^
./GraphFromFasta_MPI.cc:1479:27: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (c == i)
^
./GraphFromFasta_MPI.cc:1503:27: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (c == i)
^
./GraphFromFasta_MPI.cc:2235:24: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (j=0; j<=d.isize()-k; j++) {
^
./GraphFromFasta_MPI.cc:2264:27: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (c == i) {
^
./GraphFromFasta_MPI.cc:2359:27: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
if (c == i) {
^
./GraphFromFasta_MPI.cc:2475:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (i=0; i ^
./GraphFromFasta_MPI.cc:2483:30: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (size_t j = 0; j < p.size(); j++) {
^
./GraphFromFasta_MPI.cc:2496:30: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (size_t j = 0; j < p.size(); j++) {
^
/bin/rm -f lib_GraphFromFasta_MPI_temp.a
ar -qc lib_GraphFromFasta_MPI_temp.a obj/GraphFromFasta_MPI.o obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/DeBruijnGraph.o obj/analysis/KmerTable.o obj/analysis/NonRedKmerTable.o obj/analysis/sequenceUtil.o obj/analysis/stacktrace.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./GraphFromFasta_MPI
g++ -ggdb3 -fopenmp -o ./GraphFromFasta_MPI obj/GraphFromFasta_MPI.o -L. -lm -pthread -l_GraphFromFasta_MPI_temp
/bin/rm lib_GraphFromFasta_MPI_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:53 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/IsoformAugment.cc -o obj/analysis/IsoformAugment.o
./analysis/IsoformAugment.cc: In function ‘int main(int, char**)’:
./analysis/IsoformAugment.cc:87:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (j=0; j ^
./analysis/IsoformAugment.cc:122:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (j=0; j ^
/bin/rm -f lib_IsoformAugment_temp.a
ar -qc lib_IsoformAugment_temp.a obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/IsoformAugment.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./IsoformAugment
g++ -ggdb3 -fopenmp -o ./IsoformAugment obj/analysis/IsoformAugment.o -L. -lm -pthread -l_IsoformAugment_temp
/bin/rm lib_IsoformAugment_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:27:56 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/JoinTransByPairs.cc -o obj/analysis/JoinTransByPairs.o
./analysis/JoinTransByPairs.cc: In function ‘int main(int, char**)’:
./analysis/JoinTransByPairs.cc:293:14: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (i=0; i ^
./analysis/JoinTransByPairs.cc:341:14: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (i=0; i ^
/bin/rm -f lib_JoinTransByPairs_temp.a
ar -qc lib_JoinTransByPairs_temp.a obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/JoinTransByPairs.o obj/analysis/KmerTable.o obj/analysis/NonRedKmerTable.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./JoinTransByPairs
g++ -ggdb3 -fopenmp -o ./JoinTransByPairs obj/analysis/JoinTransByPairs.o -L. -lm -pthread -l_JoinTransByPairs_temp
/bin/rm lib_JoinTransByPairs_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:28:00 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/QuantifyGraph.cc -o obj/analysis/QuantifyGraph.o
/bin/rm -f lib_QuantifyGraph_temp.a
ar -qc lib_QuantifyGraph_temp.a obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/KmerTable.o obj/analysis/NonRedKmerTable.o obj/analysis/QuantifyGraph.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./QuantifyGraph
g++ -ggdb3 -fopenmp -o ./QuantifyGraph obj/analysis/QuantifyGraph.o -L. -lm -pthread -l_QuantifyGraph_temp
/bin/rm lib_QuantifyGraph_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:28:03 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/ReadsToTranscripts.cc -o obj/analysis/ReadsToTranscripts.o
/bin/rm -f lib_ReadsToTranscripts_temp.a
ar -qc lib_ReadsToTranscripts_temp.a obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/NonRedKmerTable.o obj/analysis/ReadsToTranscripts.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./ReadsToTranscripts
g++ -ggdb3 -fopenmp -o ./ReadsToTranscripts obj/analysis/ReadsToTranscripts.o -L. -lm -pthread -l_ReadsToTranscripts_temp
/bin/rm lib_ReadsToTranscripts_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:28:09 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./ReadsToTranscripts_MPI.cc -o obj/ReadsToTranscripts_MPI.o
/bin/rm -f lib_ReadsToTranscripts_MPI_temp.a
ar -qc lib_ReadsToTranscripts_MPI_temp.a obj/ReadsToTranscripts_MPI.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/NonRedKmerTable.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./ReadsToTranscripts_MPI
g++ -ggdb3 -fopenmp -o ./ReadsToTranscripts_MPI obj/ReadsToTranscripts_MPI.o -L. -lm -pthread -l_ReadsToTranscripts_MPI_temp
/bin/rm lib_ReadsToTranscripts_MPI_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:28:14 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./ReadsToTranscripts_MPI_chang.cc -o obj/ReadsToTranscripts_MPI_chang.o
/bin/rm -f lib_ReadsToTranscripts_MPI_chang_temp.a
ar -qc lib_ReadsToTranscripts_MPI_chang_temp.a obj/ReadsToTranscripts_MPI_chang.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/NonRedKmerTable.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./ReadsToTranscripts_MPI_chang
g++ -ggdb3 -fopenmp -o ./ReadsToTranscripts_MPI_chang obj/ReadsToTranscripts_MPI_chang.o -L. -lm -pthread -l_ReadsToTranscripts_MPI_chang_temp
/bin/rm lib_ReadsToTranscripts_MPI_chang_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:28:18 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/RunButterfly.cc -o obj/analysis/RunButterfly.o
./analysis/RunButterfly.cc: In member function ‘void JobControl::ClearFiles()’:
./analysis/RunButterfly.cc:107:25: warning: ignoring return value of ‘int system(const char*)’, declared with attribute warn_unused_result [-Wunused-result]
system(rm.c_str());
^
./analysis/RunButterfly.cc: In member function ‘void JobControl::Submit()’:
./analysis/RunButterfly.cc:77:27: warning: ignoring return value of ‘int system(const char*)’, declared with attribute warn_unused_result [-Wunused-result]
system(theCmmd.c_str());
^
./analysis/RunButterfly.cc:85:27: warning: ignoring return value of ‘int system(const char*)’, declared with attribute warn_unused_result [-Wunused-result]
system(theCmmd.c_str());
^
/bin/rm -f lib_RunButterfly_temp.a
ar -qc lib_RunButterfly_temp.a obj/analysis/RunButterfly.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./RunButterfly
g++ -ggdb3 -fopenmp -o ./RunButterfly obj/analysis/RunButterfly.o -L. -lm -pthread -l_RunButterfly_temp
/bin/rm lib_RunButterfly_temp.a
g++ -W -Wall -Wno-unused -Wno-deprecated -ansi -pedantic -Wno-long-long -fno-nonansi-builtins -Wctor-dtor-privacy -Wsign-promo -Woverloaded-virtual -Wendif-labels -O3 -ggdb3 -DMAKE_DATE='"Die Jul 3 11:28:23 CEST 2018 "' -DMAKE_OS_RELEASE='"4.4.0-128-generic"' -DMAKE_RELEASE='"3.0"' -DNEW_MAKEFILE -imacros system/BigFileDefines.h -pthread -ftemplate-depth-30 -fno-strict-aliasing -mieee-fp -iquote . -ISpines -fopenmp -c ./analysis/TranscriptomeFromVaryK.cc -o obj/analysis/TranscriptomeFromVaryK.o
./analysis/TranscriptomeFromVaryK.cc: In function ‘int main(int, char**)’:
./analysis/TranscriptomeFromVaryK.cc:216:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (j=0; j<=d.lsize()-k; j++) {
^
./analysis/TranscriptomeFromVaryK.cc:218:22: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
for (int x=j; x ^
/bin/rm -f lib_TranscriptomeFromVaryK_temp.a
ar -qc lib_TranscriptomeFromVaryK_temp.a obj/aligns/KmerAlignCore.o obj/analysis/AACodons.o obj/analysis/DNAVector.o obj/analysis/KmerTable.o obj/analysis/NonRedKmerTable.o obj/analysis/TranscriptomeFromVaryK.o obj/base/ErrorHandling.o obj/base/FileParser.o obj/base/StringUtil.o obj/util/mutil.o
./checkLock ./TranscriptomeFromVaryK
g++ -ggdb3 -fopenmp -o ./TranscriptomeFromVaryK obj/analysis/TranscriptomeFromVaryK.o -L. -lm -pthread -l_TranscriptomeFromVaryK_temp
/bin/rm lib_TranscriptomeFromVaryK_temp.a
make[1]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/Chrysalis'
cd trinity-plugins && make trinity_essentials
make[1]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins'
tar xvf seqtk-trinity.v0.0.2.tar.gz
seqtk-trinity-0.0.2/
seqtk-trinity-0.0.2/.gitignore
seqtk-trinity-0.0.2/LICENSE
seqtk-trinity-0.0.2/Makefile
seqtk-trinity-0.0.2/README.md
seqtk-trinity-0.0.2/khash.h
seqtk-trinity-0.0.2/kseq.h
seqtk-trinity-0.0.2/seqtk.c
seqtk-trinity-0.0.2/testing/
seqtk-trinity-0.0.2/testing/newformat_1.fq
seqtk-trinity-0.0.2/testing/newformat_2.fq
seqtk-trinity-0.0.2/testing/oldformat_1.fq
seqtk-trinity-0.0.2/testing/oldformat_2.fq
cd seqtk-trinity-0.0.2 && make
make[2]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/seqtk-trinity-0.0.2'
gcc -g -Wall -O2 -Wno-unused-function seqtk.c -o seqtk-trinity -lz -lm
seqtk.c: In function ‘stk_comp’:
seqtk.c:444:16: warning: variable ‘lc’ set but not used [-Wunused-but-set-variable]
int la, lb, lc, na, nb, nc, cnt[11];
^
make[2]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/seqtk-trinity-0.0.2'
mv seqtk-trinity-0.0.2/seqtk-trinity ./BIN/.
tar -zxvf "ParaFly-0.1.0".tar.gz && \
cd "ParaFly-0.1.0" && sh ./configure --prefix=`pwd` && make install && \
cp bin/ParaFly ../BIN/
ParaFly-0.1.0/
ParaFly-0.1.0/LICENSE
ParaFly-0.1.0/Makefile
ParaFly-0.1.0/Makefile.am
ParaFly-0.1.0/Makefile.in
ParaFly-0.1.0/README.md
ParaFly-0.1.0/aclocal.m4
ParaFly-0.1.0/bin/
ParaFly-0.1.0/bin/ParaFly
ParaFly-0.1.0/config.h
ParaFly-0.1.0/config.h.in
ParaFly-0.1.0/config.log
ParaFly-0.1.0/config.status
ParaFly-0.1.0/configure
ParaFly-0.1.0/configure.ac
ParaFly-0.1.0/depcomp
ParaFly-0.1.0/install-sh
ParaFly-0.1.0/missing
ParaFly-0.1.0/src/
ParaFly-0.1.0/src/.deps/
ParaFly-0.1.0/src/.deps/ParaFly.Po
ParaFly-0.1.0/src/.deps/argProcessor.Po
ParaFly-0.1.0/src/Makefile
ParaFly-0.1.0/src/Makefile.am
ParaFly-0.1.0/src/Makefile.in
ParaFly-0.1.0/src/ParaFly.cpp
ParaFly-0.1.0/src/argProcessor.cpp
ParaFly-0.1.0/src/argProcessor.hpp
ParaFly-0.1.0/stamp-h1
ParaFly-0.1.0/test/
ParaFly-0.1.0/test/Makefile
ParaFly-0.1.0/test/cmds.txt
ParaFly-0.1.0/test/rand_exit.pl
checking for a BSD-compatible install... /usr/bin/install -c
checking whether build environment is sane... yes
checking for a thread-safe mkdir -p... /bin/mkdir -p
checking for gawk... gawk
checking whether make sets $(MAKE)... yes
checking for g++... g++
checking for C++ compiler default output file name... a.out
checking whether the C++ compiler works... yes
checking whether we are cross compiling... no
checking for suffix of executables...
checking for suffix of object files... o
checking whether we are using the GNU C++ compiler... yes
checking whether g++ accepts -g... yes
checking for style of include used by make... GNU
checking dependency style of g++... gcc3
checking for library containing cos... none required
configure: creating ./config.status
config.status: creating Makefile
config.status: creating src/Makefile
config.status: creating config.h
config.status: config.h is unchanged
config.status: executing depfiles commands
make[2]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0'
Making install in src
make[3]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/src'
g++ -DHAVE_CONFIG_H -I. -I.. -pedantic -fopenmp -Wall -Wextra -Wno-long-long -Wno-deprecated -m64 -g -O2 -MT ParaFly.o -MD -MP -MF .deps/ParaFly.Tpo -c -o ParaFly.o ParaFly.cpp
mv -f .deps/ParaFly.Tpo .deps/ParaFly.Po
g++ -DHAVE_CONFIG_H -I. -I.. -pedantic -fopenmp -Wall -Wextra -Wno-long-long -Wno-deprecated -m64 -g -O2 -MT argProcessor.o -MD -MP -MF .deps/argProcessor.Tpo -c -o argProcessor.o argProcessor.cpp
mv -f .deps/argProcessor.Tpo .deps/argProcessor.Po
g++ -pedantic -fopenmp -Wall -Wextra -Wno-long-long -Wno-deprecated -m64 -g -O2 -o ParaFly ParaFly.o argProcessor.o
make[4]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/src'
test -z "/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/bin" || /bin/mkdir -p "/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/bin"
/usr/bin/install -c ParaFly '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/bin'
make[4]: Nothing to be done for 'install-data-am'.
make[4]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/src'
make[3]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0/src'
make[3]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0'
make[4]: Entering directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0'
make[4]: Nothing to be done for 'install-exec-am'.
make[4]: Nothing to be done for 'install-data-am'.
make[4]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0'
make[3]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0'
make[2]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins/ParaFly-0.1.0'
ln -sf Trimmomatic-0.36 Trimmomatic
make[1]: Leaving directory '/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity-plugins'
sh ./util/support_scripts/trinity_install_tests.sh
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Performing Unit Tests of Build

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Inchworm: has been Installed Properly
Chrysalis: has been Installed Properly
QuantifyGraph: has been Installed Properly
GraphFromFasta: has been Installed Properly
ReadsToTranscripts: has been Installed Properly
parafly: has been Installed Properly
➜ trinityrnaseq-Trinity-v2.6.6 git:(master) ✗
➜ trinityrnaseq-Trinity-v2.6.6 git:(master) ✗ ./Trinity

###############################################################################
#

______ ____ ____ ____ ____ ______ __ __
| || \ | || \ | || || | |
| || D ) | | | _ | | | | || | |
|_| |_|| / | | | | | | | |_| |_|| ~ |
| | | \ | | | | | | | | | |___, |
| | | . \ | | | | | | | | | | |
|__| |__|\_||____||__|__||____| |__| |____/

#
#
# Required:
#
# --seqType :type of reads: ('fa' or 'fq')
#
# --max_memory :suggested max memory to use by Trinity where limiting can be enabled. (jellyfish, sorting, etc)
# provided in Gb of RAM, ie. '--max_memory 10G'
#
# If paired reads:
# --left :left reads, one or more file names (separated by commas, no spaces)
# --right :right reads, one or more file names (separated by commas, no spaces)
#
# Or, if unpaired reads:
# --single :single reads, one or more file names, comma-delimited (note, if single file contains pairs, can use flag: --run_as_paired )
#
# Or,
# --samples_file tab-delimited text file indicating biological replicate relationships.
# ex.
# cond_A cond_A_rep1 A_rep1_left.fq A_rep1_right.fq
# cond_A cond_A_rep2 A_rep2_left.fq A_rep2_right.fq
# cond_B cond_B_rep1 B_rep1_left.fq B_rep1_right.fq
# cond_B cond_B_rep2 B_rep2_left.fq B_rep2_right.fq
#
# # if single-end instead of paired-end, then leave the 4th column above empty.
#
####################################
## Misc: #########################
#
# --SS_lib_type :Strand-specific RNA-Seq read orientation.
# if paired: RF or FR,
# if single: F or R. (dUTP method = RF)
# See web documentation.
#
# --CPU :number of CPUs to use, default: 2
# --min_contig_length :minimum assembled contig length to report
# (def=200)
#
# --long_reads :fasta file containing error-corrected or circular consensus (CCS) pac bio reads
# (** note: experimental parameter **, this functionality continues to be under development)
#
# --genome_guided_bam :genome guided mode, provide path to coordinate-sorted bam file.
# (see genome-guided param section under --show_full_usage_info)
#
# --jaccard_clip :option, set if you have paired reads and
# you expect high gene density with UTR
# overlap (use FASTQ input file format
# for reads).
# (note: jaccard_clip is an expensive
# operation, so avoid using it unless
# necessary due to finding excessive fusion
# transcripts w/o it.)
#
# --trimmomatic :run Trimmomatic to quality trim reads
# see '--quality_trimming_params' under full usage info for tailored settings.
#
#
# --no_normalize_reads :Do *not* run in silico normalization of reads. Defaults to max. read coverage of 50.
# see '--normalize_max_read_cov' under full usage info for tailored settings.
# (note, as of Sept 21, 2016, normalization is on by default)
#
# --no_distributed_trinity_exec :do not run Trinity phase 2 (assembly of partitioned reads), and stop after generating command list.
#
#
# --output :name of directory for output (will be
# created if it doesn't already exist)
# default( your current working directory: "/home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/trinity_out_dir"
# note: must include 'trinity' in the name as a safety precaution! )
#
# --workdir :where Trinity phase-2 assembly computation takes place (defaults to --output setting).
# (can set this to a node-local drive or RAM disk)
#
# --full_cleanup :only retain the Trinity fasta file, rename as ${output_dir}.Trinity.fasta
#
# --cite :show the Trinity literature citation
#
# --verbose :provide additional job status info during the run.
#
# --version :reports Trinity version (Trinity-v2.6.6) and exits.
#
# --show_full_usage_info :show the many many more options available for running Trinity (expert usage).
#
#
###############################################################################
#
# *Note, a typical Trinity command might be:
#
# Trinity --seqType fq --max_memory 50G --left reads_1.fq --right reads_2.fq --CPU 6
#
#
# and for Genome-guided Trinity:
#
# Trinity --genome_guided_bam rnaseq_alignments.csorted.bam --max_memory 50G
# --genome_guided_max_intron 10000 --CPU 6
#
# see: /home/urbe/Tools/trinityrnaseq-Trinity-v2.6.6/sample_data/test_Trinity_Assembly/
# for sample data and 'runMe.sh' for example Trinity execution
#
# For more details, visit: http://trinityrnaseq.github.io
#
###############################################################################

Rcorrector: efficient and accurate error correction for Illumina RNA-seq reads

BioStar — Tue, 04 Feb 2020 23:23:16 -0600

Rcorrector has an accuracy higher than or comparable to existing methods, including the only other method (SEECER) designed for RNA-seq reads, and is more time and memory efficient. With a 5 GB memory footprint for 100 million reads, it can be run on virtually any desktop or server. The software is available free of charge under the GNU General Public License from https://github.com/mourisl/Rcorrector/.

Usage: perl run_rcorrector.pl [OPTIONS]
OPTIONS:
	Required
	-s seq_files: comma separated files for single-end data sets
	-1 seq_files_left: comma separated files for the first mate in the paried-end data sets
	-2 seq_files_right: comma separated files for the second mate in the paired-end data sets
	-i seq_files_interleaved: comma sperated files for interleaved paired-end data sets
	Optional
	-k INT: kmer_length (<=32, default: 23)
	-od STRING: output_file_directory (default: ./)
	-t INT: number of threads to use (default: 1)
	-trim : allow trimming (default: false)
	-maxcorK INT: the maximum number of correction within k-bp window (default: 4)
	-wk FLOAT: the proportion of kmers that are used to estimate weak kmer count threshold, lower for more divergent genome (default: 0.95)
	-ek INT: expected number of kmers; does not affect the correctness of program but affects the memory usage (default: 100000000)
	-stdout: output the corrected reads to stdout (default: not used)
	-verbose: output some correction information to stdout (default: not used)
	-stage INT: start from which stage (default: 0)
		0-start from begining(storing kmers in bloom filter) ;
		1-start from count kmers showed up in bloom filter;
		2-start from dumping kmer counts into a jf_dump file;
		3-start from error correction.

Address of the bookmark: https://github.com/mourisl/Rcorrector/