Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [software-recommendation]

Questions asking whether a software exists for a certain purpose; questions asking which software is best among a selection for a certain purpose.

0
votes
0answers
16 views

Can I calculate maximum parsimony for a binary dataset with morphological characters?

I have a dataset with binary morphological characters. I am wondering if it is possible to calculate the maximum likelihood of a phylogenetic tree for this dataset, but it seems that all of the ...
1
vote
0answers
27 views

Colinearity analysis with respect to gene orientation

I was reading the manual of MCScanX and found that besides a all-vs-all protein blast file the tool takes a gtf file with four columns only: scaffold name, protein name, from, to. As the tool is not ...
2
votes
0answers
22 views

Is there an R package that computes homoplasy excess ratios (HER)?

I want to calculate the homoplasy excess ratio for a phylogenetic tree and a dataset with character states and am wondering if there are any R packages out there that will do this.
3
votes
1answer
63 views

How to map short sequences to long reads, recovering all multiply-mapped high-quality matches

The dilemma: I have a problem where I have around one million short sequences (21 bp to several 100s of basepairs) for which I need to identify all occurrences of in 20-30x coverage noisy long reads (...
0
votes
0answers
37 views

case/control phenotype prediction - population stratification with R

I'm struggling these last days on a case/control phenotype prediction problem based on genotype (individual x SNPs with more than 10k individuals and approx 18k SNPs encoded as 0,1,2). Thanks to PCA ...
8
votes
0answers
41 views

Best practices for deciding if two structural variants are actually the same variant?

I know that I can use bcftools isec to compare two VCF files containing single nucleotide variants (SNVs); it will generate four possible output files: one with all ...
5
votes
2answers
74 views

Deciding which samples go in which batch

I have 370 samples to sequence, we probably will end up using only 96 samples per run (due to the barcode with the primers we'll use). This means running 4 batches. To minimize the batch effect I need ...
7
votes
1answer
160 views

estimate genome size: kmer-based approach from PacBio reads

Can anyone suggest a software/method for kmer analysis using PacBio reads (RSII)? Something similar to Jellyfish, that I saw in a nice tutorial - but must be suitable for long, noisy reads. ...
1
vote
0answers
37 views

Detect gene duplication events from a whole genome in FASTA format

Are there any tools for finding duplications from a whole genome in FASTA format? Let's say the genome has gene A duplicated somewhere in the genome and gene B triplicated in tandem elsewhere in the ...
1
vote
0answers
19 views

Which secondary structure prediction algorithms don't query the PDB/DSSP for known homologous sequences?

There are several secondary structure prediction algorithms. I am interested in the algorithms which predict secondary structure without accessing the PDB to cheat by using the DSSP or STRIDE entry ...
1
vote
0answers
35 views

Do any computational phylogenetic methods enable the specification of ancestral states?

Various phylogenetic algorithms estimate ancestral states of a phylogenetic dataset. Is there a way in either maximum parsimony, distance-based methods, or Bayesian inference to indicate what the ...
0
votes
2answers
41 views

Why does `pratchet()` in `phangorn` not find the most parsimonious trees?

I called pratchet() on a dataset with binary values for the characters and it returned a parsimony value of 35. I then calculated the parsimony value for the ...
0
votes
0answers
16 views

Where can I find ChemAxon libraries for a java project?

I have to code a project using chemaxon libraries for literature review. Yet it is the first time I use these and I don't know where to find it. Here is a code example that use it but generate an ...
1
vote
0answers
17 views

XCMS isn't evaluating the proper peak due to retention time drift

I'd like to do a targeted screening analysis for trace contaminates. However, due to the nature of the sample preparations, the retention times I'm getting drift a lot from run to run and XCMS isn't ...
6
votes
3answers
154 views

visualisation of genome alignment

I was asked to check the synteny of some genes in two genome assemblies of the same species (PacBio & Illumina). I was given two scaffolds couples (so 4 scaffolds total) in the Illumina genome and ...
2
votes
0answers
92 views

Program to make a haplotype network for a specific gene

Currently my lab manually uses pop art to make the network after using PHASE to id an individual's alleles. Example: I just feel like there's got to be a better, more time effective way to do this.
7
votes
1answer
57 views

How to efficiently compute the exact percentage of non-unique k-mers in a genome for given k?

I'm looking for some software that can "efficiently" (time and memory) compute the exact percentage of non-unique k-mers in a genome for given k. I don't need the k-mers or the abundances itself, I ...
1
vote
1answer
87 views

gffread: GFaSeqGet errors on coordinate overhang

Disclaimer: I had this issue posted on Tuxedo Tools users group and shared it on Twitter, but could not get an answer to this from the developers, nor find documentation of this issue online. So, I'll ...
3
votes
2answers
85 views

Exon-exon junctions: compare experimental transcripts to reference annotation

My aim is to parse an experimental transcript set (obtained by RNAseq) to check which splice junctions are already reported in a reference annotation and which ones are new. I tried ...
3
votes
2answers
84 views

What software/tool can I use to visualise a multiple sequence alignment in the xmfa format as produced by command-line tool Mauve?

What software/tool can I use to visualise a multiple sequence alignment in the xmfa format as produced by command-line tool Mauve? I am particularly interested in seeing the single nucleotide ...
3
votes
0answers
41 views

CRISPR Sequence Finder and Database Download

I am searching for tools to pull CRISPR Spacers from Bacterial Genomes. I am aware of the CRISPRDB and the corresponding identification tool on the web server. Are there other tools for finding ...
6
votes
2answers
107 views

Estimate the length of poly-A tails from randomly-primed RNAseq data

So a poly-A tail is a long chain of adenine nucleotides that is added to a messenger RNA (mRNA) molecule during RNA processing to increase the stability of the molecule. For my project, I would like ...
2
votes
0answers
32 views

Large range parameter identification

I have a model (ODE) of an enzymatic reaction, and I'd like to estimate four of its parameters to fit some available data. To obtain the model I am just using Mass-Action Kinetics. Some more details: ...
3
votes
1answer
123 views

problem of “ordering in physical positions” phasing SNPs with Shapeit

I'm trying to phase my data (whole genome resequencing SNPs) with Shapeit (from ped file) and I get this error message : ...
2
votes
1answer
45 views

Matlab scatter diagram

I have two columns of data, x and y, how can I make a ‘thermodynamic diagram’ with matlab? The z axis is the density value of dots, like below figure. The x-axis is angle, y-axis is distance. I want ...
7
votes
0answers
62 views

Chimera Alignments

I have a structure with two subunits. I am trying to show movement of the C-terminal subunit upon ligand binding by superposition with another structure from the same strain in the apo form. I want ...
7
votes
2answers
56 views

How to measure or assign hydrophobicity score values to individual amino acids of a PDB structure?

I want measure the hydrophobicity of each amino acid within a PDB structure file. Since I have the PDB file I want to consider the 3D information, rather than sequence-only measures such as GRAVY. I ...
5
votes
2answers
131 views

Software recommendation: find DNA sequence distribution over entire transcript

I would like to create a density/histogram of the distribution of a particular DNA sequence over the entire transcript using R and/or command line tools. From here, I would like to use the coordinates ...
6
votes
2answers
521 views

RIP-seq analysis?

Given an experiment consisting of an input (baseline RNA) and IP (pulldown to find RNAs bound to certain protein of interest)... Is a DE analysis performed over the RNA-seq data from the samples (lets ...
7
votes
4answers
77 views

Introduce errors in reference transcripts according to external dataset error model

I would like to modify some reference transcripts from Ensembl (D. melanogaster) to introduce a controlled rate of random errors in the sequences. The idea would be to introduce random base ...
5
votes
4answers
2k views

How to extract / convert gff3 CDS sequences to multifasta

I would like to extract all the CDS from a batch of genomes. I have found a perl script from BioStars but this does not seem to work for me. I would preferably like to have a script/ method which will ...
9
votes
2answers
478 views

How to convert fastq to fast5

fast5 is a variant of HDF5 the native format in which raw data from Oxford Nanopore MinION are provided. You can easily extract ...
4
votes
1answer
175 views

Software recommendations - DNA composition

I'm looking to generate data based on the DNA composition of a region of my genomes (data is incomplete genomes from HiSeq runs in fasta format). I'm looking for software which will give me sliding ...
11
votes
4answers
1k views

How to convert BED to GFF3

I would like to convert a BED format to GFF3. The only useful tool that I could find via a google search seems to be Galaxy, ...
2
votes
0answers
51 views

ngs time course experiment [duplicate]

I have a question regarding the RNA-Seq analysis. We have a time course data regarding a mouse model wt and mutant treated with a drug (10uM) and the taken down at different time points: 8 time ...
5
votes
3answers
112 views

Visualization tools for 3C/Hi-C long-range interaction data?

I have identified a genomic region (CpG island) that is enriched for a characteristic of interest in the vast majority of my samples. After viewing it in the UCSC genome browser, I have confirmed ...
6
votes
2answers
623 views

How to convert GFF3 to GTF2

I would like to convert a file in gff3 format to a gtf2.2 format. The reason why I would like to do this is: I have a set of ...
6
votes
2answers
196 views

What tools can I use for a bacterial core/pan genome pipeline?

I want to perform a genome comparison on a group of isolates. I want to look into two broad groups of taxa and compare the accessory genome in each group. I have been using prokka (v1.12) and roary (...
4
votes
3answers
248 views

tools to reconcile experimental transcripts with reference annotation

Looking for tools to reconcile alignment file of experimental transcripts mapped to genome (SAM/BAM) with the reference transcriptome annotation (GTF) from Ensembl (organism: D. melanogaster). The ...
4
votes
2answers
88 views

Reordering scaffolds according to a reference without a genetic map

I am trying to reorder scaffolds of a rice species, but no genetic map is available right now. Oryza sativa Japonica is a close relative of this rice species. Mummer was used to do a whole genome ...
9
votes
7answers
3k views

Fast way to count number of reads and number of bases in a fastq file?

I am looking for a tool, preferably written in C or C++, that can quickly and efficiently count the number of reads and the number of bases in a compressed fastq file. I am currently doing this using <...
5
votes
1answer
58 views

How can I dock a protein to a nucleic acid?

I have a protein of interest and I would like to now how it interacts with RNA. I have structures of both molecules. What tool can I use?
8
votes
4answers
982 views

How do I generate a color-coded tanglegram?

I want to compare two phylogenies and colour the association lines based on some metadata I have. I have been using ape cophyloplot but I have not had any success in getting the lines to colour ...
3
votes
3answers
179 views

Analysis of differential transcript usage (DTU)

Recent breakthroughs in bioinformatics tools for quantification (e.g. Cufflinks/Kallisto/Salmon etc.) and tools which can identify differential transcript usage (DTU) (e.g. DRIMSeq, Cufflinks etc.) ...
7
votes
5answers
216 views

What motif finding software is available for multiple sequences ~10Kb?

I have around ~3,000 short sequences of approximately ~10Kb long. What are the best ways to find the motifs among all of these sequences? Is there a certain software/method recommended? There are ...
15
votes
5answers
312 views

Tools for simulating Oxford Nanopore reads

Are there any free open source software tools available for simulating Oxford Nanopore reads?