Most Popular

1500 questions
3
votes
1 answer

Finding URL for UCSC data hub

I am trying to build a UCSC track hub (a private one). I have all my hub files on a university server in a simple folder called /hub. I am using Ubuntu 16.04 and the server is also Ubuntu (16.04.3 LTS (GNU/Linux 4.4.0-98-generic x86_64). I have…
mf94
  • 203
  • 1
  • 4
3
votes
1 answer

Why is cell count stated relative to animal size?

Why do some single-cell RNA-seq papers give the number of cells sampled as a fraction of the size of the organ, animal, or embryo, as opposed to just saying how many cells they sampled? For instance, Wagner et al.'s paper "Single-cell mapping of…
eric_kernfeld
  • 380
  • 1
  • 11
3
votes
1 answer

Why DoHeatmap Does not show all genes in genes.use?

I am heatmaping a list of genes by DoHeatmap function in Seurat R package. I am sure I have 212 genes but heat map shows only a few of my genes > DoHeatmap( + object = seurat, + genes.use = genes, + slim.col.label = TRUE, + …
Zizogolu
  • 2,148
  • 11
  • 44
3
votes
1 answer

Do people today with a disease caused by single gene mutation that is inheritable - all share a common ancestor?

First question here so sorry if it is about both genetics and evolution. So, if I get a disease which can be clearly identified on the gene level and all my descendants inherit it, after 500 years do all people with the disease share the common…
Perf Lar
  • 33
  • 2
3
votes
2 answers

How to map reads shorter than 32bp with minimap2?

I mapped fasta sequences to a reference. The fasta sequences range in length from 17bp up to several hundred bp. I used minimap2 with the following command: minimap2 -ax sr ref.fa seqs.fa, et cetera. Only the sequences greater than 32bp len mapped…
conchoecia
  • 3,141
  • 2
  • 16
  • 40
3
votes
1 answer

how to filter out cell after doing QC

I was wondering that the threshold that we set after doing QC in Seurat. there is an example # GenePlot is typically used to visualize gene-gene relationships, but can be used for anything # calculated by the object, i.e. columns in…
hua
  • 441
  • 1
  • 7
  • 14
3
votes
1 answer

Glocal\semi-local\Hybride Globale-Local alignment with Python

I was looking for a simple way to do a glocal alignment. The case I have is I have a small sequence which should be find in a bigger one, thus typically a glocal alignment. Also I can not Install anything on the target machine (Python 3.6 with…
Ckln
  • 43
  • 4
3
votes
0 answers

RNA-Seq type and and optimal fusion detection

There are several popular types of RNA-Seq library prep which are frequently used: total RNA (with and without ribosomal depletion), mRNA-Seq/poly-A, and targeted mRNA-Seq/RNA exome. What would the impact of choice of these various types of RNA-Seq…
Matt Bashton
  • 1,069
  • 6
  • 16
3
votes
1 answer

Path to career in bioinformatics?

I've been completing various courses on bioinformatics through edX and Coursera, and I've made a strong footing in the material. I'm not sure what to pursue next in bioinformatics that would bring me closer to research in the field or employment.…
Rock910
  • 41
  • 1
3
votes
2 answers

Hemoglobin subunits genes in scRNA-seq

In one scRNA-seq sample I encountered the genes: Hbb-bs, Hba-a1 and Hba-a2. These genes appear on top of the list of the highest expressed genes but the 75% percentile of cells have very low counts. Therefore the expression is highly…
gc5
  • 1,783
  • 18
  • 32
3
votes
1 answer

How to assess the quality of assembled .fasta genome files?

I have assembled 3 .fasta files from contigs infastq format of 3 different Homo sapiens. I would like to see if the assembled genomes which are in .fasta format are in good quality or not. In addition, I would be happy to see heritages of genes…
0x90
  • 1,437
  • 9
  • 18
3
votes
1 answer

The algorithm behind makeblastdb?

What is the algorithm/concept behind‚makeblastdb? I know makeblastdb produces a .phr, .pin, and a .psq file, that represents an index or a data structure that is better optimised for searching a given query. What is exactly stored in these…
Paul
  • 327
  • 2
  • 8
3
votes
1 answer

Change sample names while using cummeRbund

I want to analyze RNA seq data using R package cummeRbund. Using this package I have a problem changing sample names. When I run this command in cummeRbund samples(cuff) sample_index sample_name sample_name parameter value 1 1 …
3
votes
3 answers

Annotation with Prokka or RAST

I was experimenting Prokka and RAST annotation tools. So, I took a well-annotated swinepox virus genome from genebank (NCBI Reference Sequence: NC_003389.1). I ran those sequences on Prokka and RAST Seed server at the same time. I can see that only…
L R Joshi
  • 719
  • 3
  • 11
3
votes
3 answers

How to visualize genome track of gene in specific cell-lines?

I'm trying to make a plot showing genome tracks of specific genes in specific cell-lines of RNA-seq and Chip-seq data. It should look something like this I have recently seen this encode, but in this I didn't find any MDAMB231 cell-lines. May I…
stack_learner
  • 1,262
  • 14
  • 26