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

Unanswered Questions

250 questions with no upvoted or accepted answers
9
votes
1answer
138 views

How to decide number of neighbors and resolution for Louvain clustering

I am using Louvain clustering (1,2) to cluster cells in scRNAseq data, as implemented by scanpy. One of the parameter required for this kind of clustering is the number of neighbors used to construct ...
9
votes
0answers
150 views

How to retrieve logical expressions (KO based) for reactions from KEGG?

The completeness of a module can easily be checked by looking at the Definition entry associated with the module. For example, in module M00010, it is given as: <...
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 ...
8
votes
0answers
143 views

State of the art mutation simulation software

There are many features affecting mutation probabilities, e.g. CpG mutations are 10-fold more likely than other types of mutations. Is there a model (preferably with software) which can take two ...
8
votes
0answers
191 views

Improving consensus assembly from UMI-tagged nanopore reads

The Albertsen lab has recently put out a competition/challenge for read error correction I only found out about this today, and I think the conference where high-ranking participants were going to be ...
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
0answers
103 views

Convert local alignments to spliced alignments in SAM file

I mapped RNA reads to reference genome, using LAST in split mode, and converted the MAF alignment to SAM with maf-convert. My problem is that the transcripts are not reported in a spliced manner, ...
6
votes
1answer
171 views

Pathway level analysis of single-cell gene expression

I'm looking for single-cell specific methods to construct (using gene expression data) new features that express pathway "level" or "activity", and then use these for clustering cells. One example ...
6
votes
0answers
145 views

Comparing two genome annotations

I have one-dimentional array (human genome). Also I have two annotations for it, we can think about them as different peaks (it's nucleosome and secondary structures). How can we find correlation and/...
6
votes
0answers
73 views

Where can I find summary data for how common certain mutation *types* are?

I'd like to know how common certain mutation types are in public data sets like the 1000 Genomes, ExAC, and ESP6500. Specifically, I'd like to know the distribution of stop-gains, stop-losses, ...
6
votes
0answers
57 views

How to interpret contig-alignment.psa produced by velvet

I'm using velvet to align given reads of RNA to given CDSs (i.e. coding areas and genes) of an organism, so I can generate gene-expression profiles. But after using ...
5
votes
0answers
24 views
+50

What is the minimum number of individuals per group required for ADMIXTURE analysis?

What is the minimum number of individuals required per group (population/species) for an ADMIXTURE analysis to be appropriate? Is there a way to determine this number mathematically based on the ...
5
votes
0answers
59 views

How can I read multiple different regions from a BAM file?

I’m trying to process BAM records (of a coordinate-sorted, indexed BAM file) successively in fixed-size, non-overlapping, sliding genome coordinate windows. Unfortunately I am unable to read records ...
5
votes
0answers
48 views

Run cuffcompare in strand-agnostic mode

Is there a way to run Cufflinks' cuffcompare in a strand-agnostic mode? I would like to do this because I have some RNA-seq datasets derived from an unstranded run, that should be compared to a ...

15 30 50 per page