Newest Questions
6,549 questions
1
vote
0
answers
32
views
STACKS process_radtags for demultiplexing RAD-seq data returning 57.5% "barcode not found" drops, retaining only 41.2% of reads
I'm having a very specific issue with demultiplexing my single-digest RAD-seq data using STACKS process_radtags. I have 389 individuals across 5 different libraries,...
- 11
1
vote
0
answers
5
views
BioID SAINT downstream analysis
I have SAINT results which contain 2 types of baits CCER2-C and CCER2-N. I need to determine proteins which interact with CCER2-C vs CCER2-N. I have come across heat maps and dot plots by Prohits-viz ...
- 153
0
votes
0
answers
8
views
ChimeraX AlphaFold3 prediction pseudobond naming
I had a question about how to label the residues between which pseudobonds occur in ChimeraX 1.9.
I generated pseudobonds by using the PAE output from AlphaFold3 and the command
...
- 1
0
votes
0
answers
28
views
Choosing between Intel and AMD for bioinformatics and molecular dynamics workflows on Linux
I’m in the process of building a workstation for scientific computing, primarily focused on bioinformatics and molecular dynamics simulations. My daily tasks include:
• Genome assembly and ...
1
vote
0
answers
11
views
Anomalies while docking and visualising using Python on Colab
I need some help with what I am trying to make. This program/tool aims to perform the entire molecular docking and visualisation process on a Jupyter Notebook in Colab, without having to manually do ...
- 11
1
vote
0
answers
10
views
GRaNIE package: the output contains the wrong rat genome assembly
I initialize the gene regulatory network with the appropriate genome assembly and ensembl version for my data:
...
- 11
0
votes
0
answers
11
views
Association (Protein interaction) in Intact database results what can really mean?
I'm looking for every possible interaction between proteins. Many results were association (not physical just "association"). The interaction detection method were anti tag ...
0
votes
0
answers
20
views
Warning in bcftools when subsetting a VCF by variant ID
I have successfully subset a VCF by variant ID after phasing and imputing using Beagle 5.5, however, I am getting the warning "[W::vcf_parse_info] INFO 'END' is not defined in the header, ...
- 1
3
votes
1
answer
15
views
Getting organism information from the ENA API
I have been generating Snakemake sample tables from a GEO identifier (GSExxxxxx) directly by making requests to the NCBI and ENA APIs. It works, but with some caveats (explained in this open question: ...
- 767
0
votes
1
answer
19
views
Combining SCTransform, FindIntegrationAnchors, and Harmony
I'm trying to integrate 9 scRNA datasets from different publications. When using SCTranform combined with IntegrateData using anchors, the resulting UMAP is fairly noisy with some clear batch effect. ...
1
vote
1
answer
12
views
Installation of biogrid2 plugin in cytoscape
I'm trying to install the Biogrid2 plugin in the latest version of cytoscape (v3.10.3), but I get this error:
...
- 11
3
votes
1
answer
31
views
Construction of 1000 genome phase 3 reference panel
Currently I'm learning to merge reference panel with WGS data to attain a more population-specific panel for imputation of SNP array data.
But I face 1 issue, it is not clear how to construct ...
- 31
1
vote
0
answers
7
views
Finding matched RNA-seq and Ribo-seq datasets for Nicotiana benthamiana under the same condition
I am working on translation efficiency analysis in Nicotiana benthamiana. To do this properly, I need paired RNA-seq and Ribo-seq datasets collected under the same biological condition (same tissue, ...
-1
votes
0
answers
19
views
How to respond the comment on the sample size in qPCR
In my bioinformatics article, I discuss the identification of MEG3 and MAPK3 as potential therapeutic targets for osteoarthritis prognosis. In the experimental stage, Tissue samples from three case ...
- 1
1
vote
0
answers
14
views
Unable to install cytoscape on the new version of sequoia (15.4.1)
I get this error message when trying to open the installer:
You can’t open the application “Cytoscape Installer” because this application is not supported on this Mac.
I am pretty sure this is ...
- 11
2
votes
1
answer
19
views
What tool to use for visualizing combined TF binding sites and H3K27me3 enrichment across multiple cell lines using deepTools?
I'm analyzing ChIP-seq data for the transcription factor of interest (TF) and the H3K27me3 histone modification across three different cell lines: H9, HSMM, and iPSC. For each cell line, I have:
TF ...
1
vote
0
answers
36
views
CRISPR Screen Alignments Mystery
I am analyzing a CRISPR screen. Illumina MiSeq, 5 million 75PE reads, and a PhiX spike in of 10%. The experiment is a depletion CRISPR library over time of a cell line.
It is a subset of 700 guide ...
1
vote
2
answers
18
views
Will Cytoscape allow me to draw a network piece by piece and test claimed interactions?
As a scientific editor, I sometimes need to envisage proposed gene–enzyme–metabolite networks while editing papers. Being a visual thinker, I have to plot the claimed interactions and effects (as ...
- 11
2
votes
0
answers
25
views
What is the difference between JASPAR IDs MA0668.3 and MA1993.2 for Neurod2?
I was looking into the JASPAR database for transcription factor binding motifs related to Neurod2, and I found two ids:
MA0668.3
MA1993.2
Both IDs correspond to the same transcription factor (...
- 21
4
votes
0
answers
43
views
Understanding NCBI eutils API results
I have made a Python function that can turn an NCBI GEO identifier (GSExxxxxx) into an ENA bioproject accession using the NCBI API. I need this to get sample tables which I then feed into Snakemake. ...
- 767
0
votes
0
answers
20
views
Cytoscape stringapp takes too much time getting protein queries and cloudflare's servers time out
I'm trying to load a list of 19153 genes in Cytoscape with the stringapp protein query, but since the list is gigantic, it takes a lot of time to gather all the protein names and cloudflare's server ...
- 1
1
vote
0
answers
12
views
Automated Molecular Docking Script for Library of Receptor Proteins and Singular Ligand
My experiment is to see if a lipid molecule binds to the surface proteins of a class of cells, to see if there is potentially an interaction.
I am relatively new to molecular docking, but am curious ...
- 11
0
votes
0
answers
9
views
Inquiry about using EMBOSS MATCHER with Python
I am new to Bioinformatics and want to know if this is possible. I have a .txt file with a series of 50 gene sequences and I want to know if I can find the identity of them iteratively through a ...
1
vote
0
answers
45
views
Add statistical test to this ggplot boxplot in R
Good morning,
I was wondering if anyone could assist me with the following: I am struggling to add a statistical test to this boxplot like a bar (preferably with Wilcox test). I tried with adding ...
- 65
1
vote
1
answer
126
views
How to subset a large BAM file using BED file?
I have few NIST BAM files which require to be filtered for reads that are overlapping regions with the BED file. I have tried the following command using samtools:
...
2
votes
1
answer
23
views
Weighted mean allele frequency (x-chromosome) calculation
This is a well-known formula that I found in a printed book:
$AF = \frac{(AF_{female} \hspace{4mm} \times 2 \times N_{females} \hspace{4mm})~+~(AF_{male} \hspace{4mm} \times N_{males} \hspace{4mm})}{2 ...
- 131
2
votes
0
answers
12
views
BEAST Low Posterior ESS and Posterior Support Values for MCC nodes
I've been working with BEAST for the past couple weeks to try to make a comprehensive phylogeny of 16S for a genus of frogs. There are two groups of interest to me: the mainland clade and the ...
2
votes
1
answer
8
views
PSIMI 25 style has disappeared from Style menu in Cytoscape [Cytoscape]
I have used the PSIMI 25 style, which was a built-in style, to be able to easily distinguish human proteins/nodes from non-human ones in Cytoscape (I'm using 3.10.3) networks. This seems to have ...
3
votes
1
answer
37
views
analysis of a few genes in RNA-Seq datasets
I'm analyzing the gene expression of three genes in several RNA-Seq datasets (each dataset investigates a different disease - I'm evaluating the difference of RNA counts in healthy vs disease).
I'm ...
- 119
2
votes
0
answers
16
views
How to match output alleles of modkit and sniffles2 outputs in the wf human variation pipeline?
I am using the wf-human-variation pipeline to obtain (1) DNA methylation data and (2) structural variation data. According to their documentation, the methylation results are labelled according to ...
- 21