Questions tagged [software-usage]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1 vote
0 answers
101 views

How to modify dot plot in MUMmer 3 for bacteria comparative genomics?

I was trying to make a dotplot to visualize genome-genome sequence alignment by using the MUMmer 3 software (it is typically used to compare whole genome sequences of bacteria). I runned the same ...
user avatar
2 votes
2 answers
103 views

What is the file .mums created by mummer?

I am doing genome sequence alignment using MUMmer, in particular I want to do a dotplot with mummerplot. So the passages that I did are: 1.create a file .mums with the following command line: ...
user avatar
1 vote
0 answers
253 views

Prokka error message

I am trying to run Prokka but the following error message appears: Could not run command: makeblastdb -hash_index -dbtype prot -in /cluster/software/prokka/1.13.7-gompi-2019a/db/kingdom/Archaea/sprot ...
user avatar
2 votes
2 answers
228 views

Genome assembly of SRR12196449 with SPAdes

I am trying to assemble the run SRR12196449 with SPAdes. The description of their project is: This project expected to standardize a method for amplification and ...
user avatar
0 votes
0 answers
30 views

Blast templates not found in PSI-TM Coffee

I'm trying to run psicoffe mode on Ubuntu 18.04, but it returns the error: 17633 -- WARNING: Impossible to find BLAST Templates Check that your blast server is properly installed [See documentation][...
user avatar
  • 1
0 votes
2 answers
449 views

Permanent error in installing R package

I want to install NNLM ...
user avatar
  • 1,617
1 vote
1 answer
174 views

What is the perferred method of optimization or energy minimization of small molecules downloaded from PubChem?

I will be docking a small library of molecules that I have downloaded from PubChem using AutoDock. I thought why not minimize their energy before docking using ChemDraw. I am thinking; what is the ...
user avatar
  • 25
1 vote
0 answers
473 views

Create a genotype matrix from VCF file using bcftools

I have a vcf file that contains SNPs from chromosome 2 for 2000 samples. I want to create a genotype matrix from this vcf file. The file is in zipped format as well so that I can use ...
user avatar
  • 331
1 vote
1 answer
46 views

A box plot of qualitative variables

I have 2 groups of patients; responders to chemotherapy and non-responders I have calculated cancer cell fraction (CCF) of a set of driver genes (by variant allele frequency) for each group and I ...
user avatar
  • 1,617
-2 votes
1 answer
104 views

Running a script in python [closed]

I am trying to run a script in python but I getting this error ...
user avatar
  • 1,617
1 vote
1 answer
194 views

Interpreting quantitive outputs from maximum likelihood phylogenetic trees

I ran a calculation in RAxML to determine the majority consensus phylogeny of a maximum likelihood bootstrap (How to show bootstrap values on a phylogenetic tree constructed with RAxML), and I got ...
user avatar
  • 103
1 vote
2 answers
5k views

How to show bootstrap values on a phylogenetic tree constructed with RAxML

I apologize for this very basic question. With @Michael G.'s help, I have been learning how to use RAxML to build a phylogenetic tree (How to define an outgroup to build a robust amino acid tree) and ...
user avatar
  • 103
0 votes
1 answer
30 views

AutoDock 4 -soft docking

How could I turn on/off the smooth potential in AutoDock4? I need to do two docking. The first one is soft docking with smooth potential and the second is re-docking without the smooth potential.
user avatar
  • 1
2 votes
0 answers
29 views

How to interprete PlasmidSeeker results concerning plasmids?

I've been interested in plasmid detection in bacterial genomes. I tried plasmidfinder and this software found only one plasmid in each of my two bacterial genomes. I was dubious about these results, ...
user avatar
2 votes
1 answer
229 views

BUSCO does not find buscos in reference sequence of C. elegans (technical problem)

I am working on a GNU/Linux machine. I have installed BUSCO using conda. I have downloaded the nematode dataset nematoda_odb9.tar.gzhere. I downloaded the reference ...
user avatar
0 votes
1 answer
220 views

Plot a circos plot to show the consistency between 2 samples

I have called SNP and INDEL in two matched samples by strelka and extract this information from .vcf file and I have these ...
user avatar
  • 1,617
14 votes
1 answer
3k views

What is the difference between a Bioinformatics pipeline and workflow?

I want to understand the difference between pipeline systems and workflow engines. After reading A Review of Scalable Bioinformatics Pipelines I had a good overview of current bioinformatics ...
user avatar
  • 437
0 votes
1 answer
374 views

How do I tell featureCounts to ignore fragments that overlap with more than one feature?

I am using featureCounts from the subread package to count the number of fragments per gene in an RNA-seq experiment. I feed it a BAM file with paired-end reads aligned to a reference genome. In the ...
user avatar
1 vote
0 answers
92 views

BPGA - functional analysis

I have been using BPGA to perform pan-genome analysis on my set of isolates. Now, I need the functional analysis as well. However I am having problem with the function, as the program says I am ...
user avatar
  • 11
4 votes
1 answer
280 views

Can blat use more than one core/CPU to speed up the alignment?

I am using BLAT to align two versions of the genome of C. elegans. I can see in the Activity Monitor of my Mac Book Pro High Sierra that blat is using 100% of a CPU....
user avatar
4 votes
1 answer
92 views

Codon usage analysis for whole genomes

I am new to bioinformatics. So if these questions seem you to a bit childish please forgive me. I have two queries. I am intending to perform a codon usage analysis followed by correspondence ...
user avatar
  • 87
1 vote
0 answers
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 ...
user avatar
  • 507
1 vote
0 answers
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 ...
user avatar
1 vote
0 answers
95 views

Mugsy error: Can't find species II dna at output/software/mugsy_x86-64-v1r2.3/mugsy line 501

I have installed mugsy in order to create a multiple genome alignment and a phylogenetic tree of several species of nematodes. The following command successfully pulls out the help: ...
user avatar