The Stack Overflow podcast is back! Listen to an interview with our new CEO.

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.

Filter by
Sorted by
Tagged with
0
votes
0answers
15 views

Is there a way to extract reads assigned to a particular taxonomic node from SRA data?

I've noticed that SRR submissions have an 'Analysis' section, where one can see what organisms are enriched in the sample; there is even a great graphic display that uses Krona to highlight what ...
2
votes
2answers
49 views

How can I transform a mapped BAM file into an unmapped BAM file?

In order to use MergeBamAlignment (Picard), I need an unmapped BAM file and a mapped BAM file. I have two mapped BAM files: one with reads mapped to the reference I want but without metadata such as ...
1
vote
0answers
16 views

Is there a computational method/method/way to predict if a cell population in a tissue immigrated or is enriched locally?

I have single cell data which I have analysed for differential expression. In my experiment, I subjected two groups of mice (a control and a treatment group) to treatment that will lead to the ...
3
votes
1answer
31 views

Detecting broad peaks in sRNA-seq data

What kind of tool would be appropriate do detect "broad peaks" in small RNA-seq sequencing data? MACS2 appears to be developed for ChIP-seq data, but I see that there is a ...
3
votes
1answer
71 views

programmatic secondary structure prediction for >36-mer DNA oligonucleotides

I'm writing a tool to automate Sanger sequencing primer design for a production lab that uses a universal-tail chemistry Sanger sequencing to verify NGS results. Essentially, the template DNA is ...
1
vote
0answers
48 views

How to do protein quantification with Peptideshaker output (TMT labeled, MS2)?

I have TMT tagged MS2 mass spec run processed with Peptideshaker. I saved the project in a peptideshaker format cpsx. Now, I would like to do the quantification ...
0
votes
0answers
16 views

Number of recombination events from GWAS data per person

Is there a way to calculate number of recombination events from GWAS data per person? Autosomal chromosomes. I looked online to see if any tool/package is available but none. I could phase data ...
3
votes
0answers
15 views

Iupred definition of long/short form disorder

I need to predict disorder of proteins and here there is a description of Iupred in its two variants, long and short-disorder. This is a manuscript defining the method. I couldn't find a precise ...
4
votes
2answers
162 views

Plotting coverage of annotation over collection of region

I'm trying to plot "meta" coverage of annotation: i.e. features (eg. gene class) over certain regions. It is similar to read coverage plots over gene body, except my input is two bed files (both in ...
5
votes
2answers
221 views

full visualisation of draft genomes alignment

I have two draft genomes (aprox. 500Mb each), one with ~10'000 scaffolds (genome1.fasta) and one with ~1'000 contigs (...
0
votes
1answer
33 views

Exporting nt's differences from MEGA alignment

I have an alignment of some full genome sequences of DENV1 that we have recently sequenced from the serum of patients. I have already cleaned the sequences, built the alignments and done the SNPs ...
1
vote
0answers
51 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
39 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.
4
votes
1answer
207 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 (...
1
vote
0answers
45 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
1answer
73 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
87 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 ...
8
votes
1answer
602 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
67 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
26 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 ...
4
votes
2answers
81 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
50 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 ...
1
vote
0answers
20 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 ...
7
votes
3answers
562 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 ...
6
votes
1answer
63 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
196 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 ...
4
votes
2answers
125 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 ...
2
votes
2answers
279 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 ...
5
votes
1answer
69 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
163 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 ...
3
votes
0answers
37 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
251 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 : ...
3
votes
1answer
92 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 ...
9
votes
0answers
103 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 ...
8
votes
2answers
71 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
153 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
982 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 ...
8
votes
4answers
92 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
3k 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 ...
11
votes
2answers
973 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
196 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
2k 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
60 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
154 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 ...
7
votes
2answers
1k 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
248 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 (...
5
votes
3answers
286 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
219 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 ...
10
votes
7answers
5k 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
68 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?