# 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.

69 questions
Filter by
Sorted by
Tagged with
52 views

### Is there a software that allows to convert a jpeg graph into dataset?

I have a very old graph as a .jpeg picture that I would like to convert into a dataset to be able to operate with it using R/python. I have no idea if it is even possible.
• 103
46 views

### Visualization of mutations in (part of) viral genome - existing tools?

I'm trying to visualize mutations in SARS-CoV-2 sequences, somewhat like it's done on outbreak.info: but for individual sequences, probably on the base of the vcf files generated by the ARTIC ...
128 views

### How to de novo hybrid assemble with Pacbio CCS and Illumina PE reads

I would like to perform de novo genome assembly on a diploid microalgal strain. I have two datasets: PacBio CCS/HiFi reads, low coverage. Illumina PE 2x150 (standard shotgun) Does anybody have any ...
• 11
58 views

### If I have data from a full genome sequence, how can I search for a known pathogenic variant of a known gene?

I have a CS background, but am a bioinformatics neophyte I did a full genome sequence which provided me with ~100GB of files (SNP VCF, Indel VCF, BAM, Indel TBI, SNP TBI, BAM BAI, CNV VCF, CNV TBI, CV ...
• 101
168 views

### How to get a consensus sequence from a nanopore fastq files?

I am new in bioinformatic field. I would like to know a way to generate a consensus sequence from nanopore fastq files (fastq files demultiplexed). I usually generate a consensus sequence with "...
33 views

### Which GWGAS Tool is 'better', MAGMA or LDSC?

I am currently doing a project regarding gene-based analysis and gene-set analysis given a certain GWAS dataset, I only know 2 well known gene-based analysis tool which is LDSC and MAGMA. As of now, I'...
1 vote
81 views

### Finding SNP/Indel Markers of Interspecific Variation

This is a revised version of my original question, edited to include the portions already answered by @Maximilian Press. Goal: Find SNP and indel alleles which are fixed between but lack polymorphism ...
• 129
26 views

### Is is possible to predict ncRNAs from sequence and homology alone?

I'm working with a set of homologous genes (let's call it gene A) from several bacterial species. I know (from previously published research) that in gene B (a close paralogue of gene A), there is a ...
• 595
156 views

### How can I subset WGS data to the level of WES variants?

I would like to compare mutational signatures1 in patients from different studies, however some studies are based on exome seq (i.e. ~20,000 coding variants) and some are from whole genome seq (i.e. ~...
• 203
160 views

### MM/PB(GB)SA calculations for DESMOND trajectories

I did molecular dynamic simulations with DESMOND. Now, I want to obtain the interaction free energy using MM/PB(GB)SA. The free tools I found to perform such calculations are AmberTools and g_mmpbsa. ...
• 101
1 vote
40 views

### Is there a software that can predict protein-DNA binding?

Is there any software that, given a DNA sequence as input, can predict what protein/ type of protein can bind to the DNA. The sequence of interest is a potential transcriptional regulator, I want to ...
1 vote
35 views

### Gene Ontology equivalent - higher specificity?

I've recently discovered OrthoVenn2, which I'm using to compare proteomes and extract likely clade-specific genes. I'm enjoying software, and one of the things it offers is a breakdown of protein ...
• 595
226 views

### Platanus-allee phasing fail: Error(13): Error, SolveDBG exception!

I am using Platanus-allee 2.2.2 for heterozygous genome (~500mb) assembly with Illumina short reads and PacBio reads input data. I have the file contigs.fa from short reads but phasing step with ...
• 21
189 views

### Understanding some of the computational bottlenecks of Covid-19 research

I am a researcher in high-performance computing (with very little bioinformatics background), and I am trying to understand what are the current biggest computational bottlenecks of software used for ...
• 181
44 views

### Software for microbial profiling from 16S rRNA gene sequence

I have hundreds of GBs metagenomic 16S rRNA gene sequence data. I want to do microbiome composition profiling (with relative abundance) from the data. Also after that, I will do functional profiling (...
• 83
216 views

### Is it possible to do this in R?

I have a big clinical data file with 96 columns like age, gender, BMI, etc I want to see which of these clinical characteristics respond to chemotherapy. Response to chemotherapy divides patients to ...
• 1,617
199 views

### What assembler is appropriate for High-Fidelity PacBio reads

What assembler is appropriate for High-Fidelity PacBio reads? For example, canu is good for high-error PacBio reads. But what algorithm to use for HiFi reads? Would it be OK to use canu without the ...
• 2,349
1 vote
286 views

### Is there any tool for fuzzy sequence matching?

Note: this question was also asked on biostars are you aware of any tool that is able to perform error-tolerant pattern-matching search on protein FASTA files? For example, I want to know, which ...
283 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 ...
• 2,349
1 vote
34 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 ...
• 497
79 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 ...
• 2,960
105 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 ...
• 192
1 vote
117 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 ...
• 1,181
28 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 ...
• 547
487 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 ...
• 325
656 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 (...
• 2,644
308 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
122 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 ...
• 5,287
527 views

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

The dilemma: 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 (both pacbio and ONT). ...
• 3,071
1 vote
51 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 ...
100 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 ...
• 258
125 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 ...
• 4,602
1k 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. ...
• 2,644
1 vote
96 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 ...
• 2,349
1 vote
35 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 ...
• 507
87 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 ...
• 297
85 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 ...
• 297
1 vote
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 ...
• 111
2k 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,644
94 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 ...
• 69
598 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 ...
• 2,644
236 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,644
1k 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 ...
• 2,349
96 views

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 ...
263 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 ...
• 1,641
42 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: ...
• 131
698 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 : ...
895 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 ...