Questions tagged [genomics]
Genomics is an interdisciplinary field of science focusing on the structure, function, evolution, mapping, and editing of genomes.
138
questions
20
votes
1
answer
465
views
The state, limitations and comparisons of large variant stores
Background: We're increasingly needing some way of storing lots of variant data associated with lots of subjects: think clinical trials and hospital patients, looking for disease-causing or relevant ...
- 788
13
votes
1
answer
186
views
Is it possible to use SNP heterozygosity as a proxy for Indel heterozygosity?
I have estimated genome-wide heterozygosity levels using maximum likelihood and classical substitution model implemented in package atlas. These estimates are way more robust than classical SNP ...
- 5,517
7
votes
1
answer
783
views
ethnicity check either from bam or vcf files
What tool could I use to check the ethnicity of a human bam or vcf file? I would like to use the results as a QC check to know whether a given sample or set of samples match the ethnicity information ...
- 2,324
7
votes
1
answer
401
views
ANY (technical) reason behind submitting sequences to GenBank versus ENA Sequence
The DNA sequence sections of the three INSDC databases (i.e., DDBJ, ENA Sequence and GenBank) are synchronized periodically and strive to keep their stored data as ubiquitously accessible as possible. ...
6
votes
3
answers
1k
views
Somatic tumor only variant calling?
I'm evaluating possibilites for somatic tumor variant calling without paired-normal samples. I'm aware of the consequences without a normal sample.
All the popular variant calling tools such as ...
- 2,739
6
votes
1
answer
493
views
6
votes
2
answers
395
views
Meaning of the FORMAT fields of the VCF file coming from GIAB project
After reading the GIAB paper in https://www.biorxiv.org/content/early/2018/05/25/281006 and its Figure 1, I am still having trouble understanding the data inside the GIAB VCF file for HG001 (...
- 161
6
votes
1
answer
739
views
Proof of Breakpoint Reversal Sorting Approximation Algorithm
I'm kind of new to bioinformatics and trying to self-study. I'm reading a bioinformatics book: An Introduction to Bioinformatics Algorithms and ran into some problems about understanding the proof of ...
- 61
5
votes
4
answers
2k
views
Range overlap python error with genomic regions
I have two files
s3.txt :
1 10 20
1 5 20
2 20 30
2 25 30
1 10 50
2 20 60
1 14 17
s4.txt:
1 10 20
2 20 30
I am trying to match col0 ...
5
votes
3
answers
283
views
How to adapt the fgseaL function to perform rapidGSEA computation of gene ranks across 9 different phenotype labels?
I wish to adapt the r language function fgseaL, https://github.com/ctlab/fgsea , to perform rapidGSEA, https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1244-x , computation of ...
- 191
5
votes
1
answer
180
views
How to find all variable-length seqs with an exact 5' and 3' match in a FASTA file
Context
I am interested in finding all of the promotors specific to a particular sigma factor. I have identified the -35 and -10 sites from the literature, bold denotes -10, -35, binding sites:
<...
4
votes
5
answers
3k
views
using python to write bioinformatics pipelines tutorial
I was wondering if there is a tutorial or a small code snippet to understand how to write bioinformatics pipeline using python, for example
use a aligner (say hisat)
get the output and process it ...
4
votes
3
answers
2k
views
What do the FASTQ file names mean here?
I would like to run an analysis on Nanopore data. I'm trying to download sample data from https://github.com/nanopore-wgs-consortium/chm13. I would like to run minimap2 on the FASTQ files, resulting ...
- 2,739
4
votes
1
answer
652
views
Using BLAT command line tool to blat split sequences
I have the nucleotide sequence:
AATTGAGGCACATTTTTTTTTAGACAGTCTTGCTCTGTTGCCCAGGCTGGAGTGCAGTGGTGTGATCATAGCTCACTGCAGCCTCGACCTCCTGGGCTCAACAAAGCACACAGTGGGCGGATCCCCACCAG
When I blat this on UCSC Genome ...
- 241
4
votes
1
answer
115
views
Generating DNA sequences with constraints
I would like some advice on potential strategies to address the following problem.
I want to write a program that will generate DNA sequences that are optimized on two constraints based on an input ...
- 595
4
votes
1
answer
92
views
Does UCSC definition for telomere making sense?
UCSC provides a database of telomere regions for each chromosome:
ftp://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/gap.txt.gz
It reported the telomere section in chr21 as:
...
- 2,739
4
votes
1
answer
71
views
Are there any Genome in a Bottle-like resources for non-humans (especially for invertebrates)?
Genome in a Bottle is an excellent resource that provides many types of DNA and RNA sequencing reads for a single individual/cell line to test genome assembly and analysis tools. For example there are ...
- 3,141
4
votes
1
answer
55
views
Identifying somatic mutations in cell lines
I would like to identify the somatic mutations present in a cell line and characterise the genes that are potentially affected by those mutations. For example, are there oncogenes mutated in a ...
- 160
3
votes
2
answers
275
views
Selecting 65000 SNPs where AF is close to 0.5 in all or most populations
I am evaluating the tool somalier (https://github.com/brentp/somalier) and I need to create a list of about 65,000 SNPs where the allele frequency (AF) is as close to 0.5 as possible across the most ...
- 2,324
3
votes
1
answer
1k
views
What are the coordinates of the extended major histocompatibility complex (xMHC) in GRCh38?
Is there an accepted boundary of the extended major histocompatibility complex (xMHC) in human? Specifically, I am interested in the boundary coordinates for Genome Reference Consortium Human Build 38 ...
3
votes
1
answer
132
views
Installing Sibelia for Ragout on Mac OSX
I am trying to use Ragout: https://github.com/fenderglass/Ragout
to fill the gaps in my de novo genome assembly. You can access the article freely here: https://www.ncbi.nlm.nih.gov/pubmed/24931998
...
- 2,459
3
votes
1
answer
188
views
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 ...
- 41
3
votes
1
answer
9k
views
What are the differences between GWAS and QTL mapping? [closed]
What are the differences between the two methods?
What advantages does one have over the other, and what are their limitations relative to one another?
- 69
3
votes
1
answer
2k
views
Searching for start and stop codons for protein sequencing of contigs
I need to convert contigs into their respective protein sequences given a reference genome (i.e. I need to take a substring, whose position is already known on the string, and I need to locate the ...
3
votes
1
answer
101
views
Difference between copy number neutral reads and active reads
In this paper, the authors talk about copy number neutral reads (as reads that equally represent parental chromosomes) and active reads (as reads from only one parent chromosome):
We reasoned that ...
- 391
3
votes
1
answer
72
views
predict the foldability of single-stranded DNA molecules
I have a list of regions of the human genome and I want to predict if single-stranded molecules in a buffer would tend to fold and create pin structures by sequence self-complementarity. What's the ...
- 2,324
3
votes
1
answer
73
views
How to query the Human Microbiome Project (HMP) to find all subjects with both 16s and WGS workups?
I am looking for a query to run on the HMP database that will return all subjects who have had BOTH 16s and whole genomes sequence (WGS) workups. I am currently using this query...
...
- 265
2
votes
2
answers
129
views
Genomics Statistics Problem
I have used a python script to identify target sequences in a DNA sequence file. There are two classes of sequence: coding and non-coding.
I have identified 728 sequences of interest. 597 of these ...
2
votes
2
answers
174
views
BWA-mem and sambamba read group line error
this question has been asked [and answered] on Stack Overflow
This is a two-part question:
help interpreting an error;
help with coding.
I'm trying to run bwa-mem ...
2
votes
2
answers
71
views
Genomic coordinates of predicted binding pocket
I would like to know the genomic coordinates of a proteins ligand-binding pocket. For example, on this PDB page you can see that there is xray crystallography data for how CYP2C9 interacts with two ...
- 841
2
votes
1
answer
202
views
What does "motif 2 bp or longer" mean?
I'm reading a preprint on a new somatic variant discovery tool, and the first sentence of the Discussion section has me confused.
Across the four datasets analyzed in this study, we discovered that ...
- 5,060
2
votes
1
answer
64
views
At what point is a gene different between species that we call it a different gene?
I am new to genetics and I know humans share many genes with mice for instance but that there are slight differences in conserved nucleotide sequences. Is there a community consensus around at what ...
- 376
2
votes
1
answer
87
views
Are there genomics databases for neurological and psychiatric disorders?
Are there databases that relate genetic disorders with neurological disorders?
I am thinking of something like COSMICS for cancer? AFAIK cancer usually develops due to a combination of 5-10 ...
- 1,427
2
votes
1
answer
49
views
How are the "effect_allele" and "other_allele" bases chosen in a PGS file from the PGS Catalog?
How are the "effect_allele" and "other_allele" bases chosen in a PGS file from the PGS Catalog?
For example, for PGS002723 (https://www.pgscatalog.org/score/PGS002723/), if I set ...
- 125
2
votes
1
answer
39
views
Inconsistent replicate numbers in RNA-seq
I have 3 groups for the RNA-seq analysis (Control, treatment A and treatment B). There are 2 replicates for control and treatment A and 3 replicates for treatment B (lost 2 replicates due to a mistake ...
- 101
2
votes
1
answer
38
views
Using very closely related strains to increase coverage for short read, de novo assembly
Do you think it's possible to combine short read Illumina libraries (WGS) from multiple closely related eukaryotic microbial strains (e.g. libraries from a re-sequencing study, >99% ITS1 sequence) ...
- 31
2
votes
1
answer
655
views
how to retain reads with low mapping quality (MAPQ) scores when using samtools view -q
I have been using the -q option of samtools view to filter out reads whose mapping quality (MAPQ) scores are below a given ...
2
votes
1
answer
282
views
Transfer annotations from one genome assembly to another
I've got an annotated draft genome assembly made using short and long read strategies. I've also done optical mapping to stitch some of the contigs together, with the goal being chromsome-length ...
- 21
2
votes
1
answer
2k
views
How to download data from SRA in Linux systems via the command line?
My workflow for downloading data from SRA has been the following:
Access SRA Run Selector.
Enter the accession number for the project of interest.
Download "Accession List" for the "...
2
votes
1
answer
134
views
COVID-19 GWAS: what is a cross-replicating association?
In a genomic study of patients infected with SARS-CoV-2, the authors detected cross-replicating associations with rs11385942 at locus 3p21.31 , the association signal spanned the genes SLC6A20, ...
- 153
2
votes
1
answer
106
views
How to extract hypervariable region of mitochondrial dna from its fasta file?
Recently I've done some analysis on human mitochondrial DNA and now I want to run the same analysis on only hypervariable regions of mitochondrial DNA to see whether the changes that I have observed ...
- 145
2
votes
1
answer
35
views
tensorQTL interaction issue
I hope that the question has its place on this board.
I am trying to run the mapping of cis eQTL using an interaction term with tensorQTL. I am running into the following index error when running the ...
- 31
2
votes
1
answer
74
views
Can I discard missing alleles?
I am converting a biallelic VCF into a SQL table and one of my tables will be something like:
...
- 123
2
votes
1
answer
507
views
Extracting base sequences from ABI/AB1 sanger sequencing chromatogram
I am trying to understand the sanger sequencing ABI/AB1 file format better, and extract base calls from given signal intensities over time.
As I understand, reading in a raw AB1/ABI file into python, ...
- 21
2
votes
1
answer
345
views
Difference between bioinformatics and cheminformatics [closed]
Anyone can answer, and should be in simple way
- 29
2
votes
1
answer
62
views
Statistical or performance test for genomic intervals between 2 samples
Is there an appropriate statistical or performance test to perform when I want to compare genomic intervals between 2 samples. Sample A is the the benchmark of comparing the intervals and sample B is ...
- 141
2
votes
1
answer
62
views
Repetition-free Longest Common Subsequence as similarity measure
I'm starting a thesis in Computer Science with topic the RFLCS (Repetition-free Longest Common Subsequence). I know very little about Bioinformatics problems and Genomics. I would like to understand ...
- 121
2
votes
1
answer
559
views
How can I improve ESS values for ucldMean in beast?
Using starBeast I ran a 50 locus tree, with 34 individuals, using evolutionary models indicated by AICc in jModel test. I ran this tree with a MCMC chain length of 500 million, storing every 10,000 ...
- 163
2
votes
0
answers
35
views
How to perform liftover from 38 to 37 in R?
I have some gwas summary statistics in GRCh38 that I want to lift to GRCh37. I am trying to liftover in R using this code:
...
- 45
2
votes
0
answers
20
views
Nucleotide alignments from mmseqs "tblastn"
I am searching a protein sequences against a nucleotide database ("tblastn") using mmseqs and need nucleotide sequence alignments as an output.
I start by the following (pg 23 of https://...
- 121