Questions tagged [genomics]
Genomics is an interdisciplinary field of science focusing on the structure, function, evolution, mapping, and editing of genomes.
80
questions
2
votes
0answers
11 views
When to use Ariba vs Roary or other genomics tools to call drug-resistance variants?
I have a large dataset of assemblies of bacterial S. aureus genomes (~1000, same species). I am looking at different tools to estimate the presence of drug-resistance factors such as blaZ or ampC. ...
3
votes
1answer
42 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 ...
1
vote
1answer
66 views
0
votes
0answers
26 views
How to find Genotype, from list of mutations:lineage list?
I'm very new to python and I'm getting along nicely (I'd like to believe); however, I must be missing something here. I'm looking to read each and every file and compare each line against one another.
...
0
votes
0answers
31 views
How to get % similarity between strains and mutation files
I’m very new to python, and having some difficulty getting hang of some more complicated things
I have multiple files which look like so:
hCoV-19/Singapore/4/2020|EPI_ISL_410535|2020-02-03
hCoV-19/...
1
vote
1answer
34 views
How to identify an unknow species from a fasta containing genetic sequences
I am a Biochemist that is unfamiliar with bioinformatic tools and new to academia as a whole.
I am currently using ILLUMINA PE data, which I trimmed (Trimmomatic), corrected (Rcorrector) and assembled ...
0
votes
1answer
29 views
Pointing Cromwell to local EBS mounted storage on AWS Batch
I'm trying to run the WDL tasks from the GTEx RNA-Seq pipeline with Cromwell using AWS Batch as a backend.
I store the STAR alignment index on an Elastic Block Store since my Cromwell server instance ...
2
votes
1answer
39 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, ...
0
votes
0answers
15 views
Filtering out nonarthropod sequences using Blobtools2
I am using Blobtools2 and am working my way through the sample dataset that they provide before using my own. The filtering tutorial I am following is --> https://blobtoolkit.genomehubs.org/...
0
votes
1answer
33 views
How to annotate gene length to a list of gene symbols using UCSC data?
I have a list of HGNC gene symbols, I am looking to get the gene length of each gene. Although I also describe these genes with lots of UCSC datasets as features, so I am wondering if there is a ...
0
votes
0answers
11 views
Swapping to effect increasing allele in case/control studies
I want to find the effect-increasing allele for a given set of SNPs. For quantitive traits, I understand the effect sign is changed when the A1 and A2 are swapped. How does this work for odds ratios?
...
1
vote
0answers
34 views
“How the clustal omega can be reverse engineered, to trace ancestral inversion mutations via the guide tree?”
I apologise for very basic question but I am very new to biology and have very specific knowledge of this field, I am currently working on Bioinformatics in my machine learning project which is why I ...
1
vote
1answer
61 views
How to compact variant data to their genes without bias?
I have a dataset of genes I am trying to collect data on from public databases, to use as features in machine learning. I am trying to take some features from UCSC genome browser (e.g. number of CpG ...
0
votes
2answers
56 views
How identifiable are human omics data and how to mitigate their identifying features?
Say a database were to store human omics datasets. The human subjects are known and the sample size is rather small in size initially (n=500). The database contains genomics, transcriptomics, ...
1
vote
0answers
79 views
Kraken2 > OTU format > Phyloseq
A collaborator has passed me over Kraken2 outputs *.report and *.kraken, from a metatranscriptomic sequencing experiment conducted on the minION.
I would like to use this data in phyloseq, however, I ...
0
votes
0answers
10 views
What is the origin of the 'source Gene ID' references given in the 'gene_presence_absence.csv' output of Roary?
I am learning to use Roary for preparing a pan genome for some lactobacillus strains. In the 'gene_presence_absence.csv' output of Roary (which i view in excel), a 'source Gene ID' is given for each ...
0
votes
0answers
28 views
How do I delete ligands from a list of 50,000
I need to delete ligands from a list of over 50,000 ligands in a .sdf file format. Is there any easy and quick way of deleting the ligands from the list.
Thank you,
Dr. Hussain.
1
vote
0answers
46 views
How to aggregate multiple MAF files together?
I have cancer genomic data (tumor/normal whole exome sequencing) from 50
patients that received the same type of treatment, half of whom responded. These come in the form of 50 .maf files, along with ...
0
votes
1answer
19 views
Is the number of contigs and scaffolds an artifact of genome assembly and rather arbitrary, or is it canonical for each species?
E.g. bowhead whale genome has X contigs and Y scaffolds, but couldn't these easily be different numbers?
1
vote
2answers
61 views
How to combine two Genome-wide Association Study (GWAS)? [closed]
I did a GWAS analysis in the past for antibiotic resistance of E. Coli and the results were interesting (matching the literature). I did a new GWAS analysis for some new samples, but the results are ...
1
vote
1answer
69 views
Problem with classification model of genomic data: every machine learing model predicts wrongly almost always the same subset of dataset
First of all, I'd like to apologize for any spelling or grammar mistakes.
I'm having a problem using R for a classification problem. My dataset contains ~300.000 genomic data, and the features are ...
0
votes
0answers
27 views
Phasing partially phased genomic data
I am trying to create an integrated callset with coding regions from phased exomes and noncoding regions from a genotyping array. I would like to merge the data together in the correct phase. Since ...
0
votes
0answers
24 views
Coding object structure properties into a sequence
Hopefully, I found the right place to ask. Please, note that I'm not a specialist in the current field.
Is there an algorithm to code information about object structural properties into a sequence (...
1
vote
1answer
39 views
How to compute Shanon entropy
I want to analyze some viral sequences. I want to use shanon entropy as a measure to analyze the sequences. However I do not understand the concept of shannon entropy. Would anyone explain please how ...
2
votes
0answers
25 views
What Cellbase resources are 0-indexed, and which are 1-indexed?
The UCSC Genome Browser Team is clear about that the genome browser GUI indexing which is 1-based, closed interval. Genomic sequence retrieval via the UCSC REST API is 1-based, closed interval. All ...
0
votes
0answers
59 views
Mapping a disease integrome on the whole genome
I'm currently studying a disease which has complex traits (non mendelian disease). Very few data are available on humans. But data collected from cell/animal models have been used to identify human ...
0
votes
0answers
24 views
How to decide the parameters and parameter values when running ChIP-seq Peak Calling for H3K27ac data in MACS?
I have alignment data in SAM files from ChIP-seq analysis for H3K27ac mark. I need to call peaks using MACS2. But I'm having issues with selecting parameters and values of the parameters. I will be ...
0
votes
1answer
25 views
What are the largest fully publicly available molecular QTL datasets?
I'm familiar with the sQTL and eQTL GTEx data that can be downloaded from the
GTEx Portal, but I'm interested in gathering at least two other types of molecular QTL datasets:
Larger eQTL and/or sQTL ...
1
vote
1answer
37 views
What does these terms mean?
I am working on cancer genomics failed to follow some terms even by googling
What is the difference of focal genome amplifications versus convergent amplification of given region of genome for ...
0
votes
2answers
77 views
CRISPR screening and systems biology
I am not sure how to tackle this one, so I’ll explain the general idea I have in mind.
Given results of a CRISPR knockout screening experiment (like in this paper) with two types of assays: control ...
2
votes
1answer
52 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 ...
2
votes
1answer
68 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
vote
3answers
119 views
Python module for fetching NCBI id for a list of species
I have a list of scientific names of species. Is there a python module that can fetch NCBI taxonomy IDs?
1
vote
0answers
57 views
How to apply RDP, Greengenes and other special taxonomies in Krona?
I used Kraken 2 to classify my 16S metagenomic data using both RDP and Greengenes database. As these are special databases the taxonomic ids assigned do not match with the NCBI taxonomic ids. I found ...
1
vote
2answers
1k views
What is the name of this type of figure?
This figure comes from this wiki page.
I googled "gene location figure" and "gene location plot", none gets similar results. I also tried search by image, none of the results is close.
so, what is ...
1
vote
1answer
56 views
Error Viewing file on HiGlass
I am implementing hi-glass python on Jupyter, and upon running the test data code
...
3
votes
2answers
444 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 minimap on the FASTQ files, resuling in ...
0
votes
1answer
120 views
creating population AF from set of individuals in 1000 genomes population?
What's the recommended way to create population-level Allele Frequency variant entries from a set of individuals (vcfs?) in the 1000 genomes (or similar) project?
For example, given the IBS ...
0
votes
1answer
85 views
what is the most complete vcf file of population allele frequencies that can be built/downloaded from public datasets?
What is the most complete vcf file of population allele frequencies that can be built/downloaded from public datasets nowadays?
About 5 years or so ago, it used to be the latest release of the CSHL ...
4
votes
1answer
85 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 ...
1
vote
2answers
93 views
Transcription factor targets for Human and Mouse genes
I want to get transcription factor targets for my list of genes (For Human and Mouse)
I have found ENCODE Transcription Factor Targets database (The link is given below). But I am looking for some ...
1
vote
1answer
51 views
Highly variable elements?
Which type of genomic elements in the UCSC browser can be used as an example of the highly labile = highly variable = low conservation regions?
4
votes
1answer
75 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
votes
1answer
56 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 ...
4
votes
1answer
63 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
votes
2answers
241 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
votes
2answers
118 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 ...
6
votes
1answer
318 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 ...
1
vote
1answer
79 views
Are there free cloud computing platforms for biology projects? [duplicate]
I want to implement the analysis found in the paper RNA-Seq of Tumor-Educated Platelets Enables Blood-Based Pan-Cancer, Multiclass, and Molecular Pathway Cancer Diagnostics The project id is 281708. ...
5
votes
1answer
137 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:
<...
