Questions tagged [genome]
for questions related to whole genome analyses. For questions related to individual genes, use tag gene instead. Questions specific to format (like fasta) should be tagged with a format-specific tag.
173
questions
0
votes
0
answers
27
views
What is the current state-of-the-art for classifying the ancestry groups of VCF files?
I'm wondering which tool or program is:
Most updated; doesn't use a superseded algorithm
Runs very quickly (ideally, under 30s per VCF file)
Able to be automated in some way for a large number of ...
- 86
2
votes
2
answers
65
views
Where do I get a large reference VCF?
I would like to download a large .vcf file containing many (hundreds or thousands) of samples. Ideally, I would download different population-specific .vcf files, but the ability to sort/filter by ...
- 86
4
votes
1
answer
53
views
What is the NCBI's definition of an "atypical genome"?
Using the new NCBI Datasets platform, you can browse the collection of genomes associated with one or more taxa. For example, searching Pseudomonas aeruginosa returns 19,878 genomes as of 29 March ...
- 563
2
votes
2
answers
122
views
Improving prokaryotic assembly with other contig/scaffold-level data?
I have what at first sight appears to be a high-quality MAG (~10 pieces, high completion%) that I built from a hybrid assembly (Illumina + Nanopore data) from a cyanobacterium.
Workflow:
Quality ...
- 777
2
votes
1
answer
50
views
Availability of information on genes in Gnomad VCF data
I was wondering does the gnomad genome data that is downloaded in the vcf format on variants contain information on what is the nearest gene and is the genomic location available exome/intron?
if yes ...
0
votes
1
answer
47
views
Allele frequency gnomad data
I noticed that some gnomad SNPs have allele frequency =0 and some are just not found : No genes were found in this region. what is the difference between them : for example:Y-10010520-10010520.
Is all ...
1
vote
4
answers
144
views
Grep peculiar pattern inside a file
I would like to grep everything inside the second square bracket of a file:
...
- 141
1
vote
1
answer
38
views
Why do Illumina 850k/EPIC arrays ignore CpGs which are "GC" in the forward strand?
CpGs are symmetrical, in that a CG sequence on the forward strand is hybridized to a GC --- and both dinucleotides on each opposing strand are CpGs dinucleotides which can be methylated. Conversely, ...
- 131
1
vote
1
answer
51
views
How to identify mutations in a viral genome
I have a fasta file with multiple sequences comprising 38 sequences. The length of the sequences are around 11000 bp.
How can i get changes in the genomes based in a reference genome? (aa subtitutions ...
1
vote
1
answer
88
views
Allele effects by parental-origin estimated from genotyping array data
We are dealing with relatively large family data (> 1,000 trios) with genotyping array data.
We would like to see the allele effects on children separately of mothers, fathers, and children, as in ...
4
votes
1
answer
51
views
How to promote assemblies into genomes in NCBI?
Note: I've never submitted an assembly/genome to NCBI, so excuse if my perspective is flawed.
I'm working with Drosophila subobscura. (spring fruit fly)
I see here https://www.ncbi.nlm.nih.gov/data-...
- 239
3
votes
2
answers
101
views
Searching for HLA-B in DNA results
I'm trying to find the HLA-B*15:01 variant in my DNA results, prompted from this research paper: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC8142661/#:~:text=HLA%2DB*15%3A01%20is%20strongly%...
- 131
2
votes
1
answer
52
views
What is the best way to process yeast genomes?
I have obtained several hundred raw, unassembled yeast genomes from NCBI and I am looking for advice on how to process the genomes for downstream analysis.
I have a reference genome (S288C) to use for ...
- 494
1
vote
1
answer
51
views
Lower case vs. upper case nucleotids in sequence vs. dots at the end
What is the difference between lower case and upper case nucleotides in a sequence? My other question is what are the dots at the end of the sequence? Some examples are shown below:
GGgG,GGgG,GGgG,...
2
votes
2
answers
102
views
Bacterial genomes and snpeff warnings 'WARNING_TRANSCRIPT_NO_START_CODON'
I annotated a bunch of genomes with snpeff. And although the original reads' qualities and the matches with the reference genomes were quite high, I see the above mentioned warning in almost every ...
- 153
2
votes
0
answers
39
views
Is there a way antiSMASH can be used to screen a large number of genomes just for one type of biosynthetic gene cluster?
When we use antismash command with the genome sequence file as follows:
antismash GCA_002095535.1_ASM209553v1_genomic.gbff
it screens the genome for all the ...
- 21
2
votes
1
answer
54
views
Is there a way to extract spike sequence from a Sars-CoV-2 complete genome (preferably in Python)?
I have a complete genome sequence of a Sars-CoV-2 variant, but I am only interested in the Spike sequence?
I should note that I am not a biologist (nor bioinformatics).
The complete genome looks like ...
- 121
1
vote
0
answers
36
views
System specifications for NGS data analysis
I have a 3.7 GB whole genome data of a eukaryote, for which genome assembly, gene prediction, and annotation steps have to be performed. In some time I would also need to analyze the transcriptome ...
- 11
3
votes
1
answer
46
views
How to make pan-core genome curve through command line on linux?
I´m working with a dataset of 566 genomes to analyze a pangenome. So I was working with PANWEB to create this pan core genome curve, however, there is too much sequence to work with this webserver. ...
- 185
2
votes
2
answers
42
views
How can I improve or otherwise investigate an unreliable genome tree?
Summary My genome tree doesn't agree with my gene trees and I get the feeling that my genome tree might be wrong, possibly due to long branch attraction, but I don't know how to check/fix it.
...
- 777
2
votes
0
answers
11
views
Flux balance Analysis, Result interpretation, solver Solver Heuristics
i am new to modelling of cellular networks, and I have two basic questions, regarding the result interpration.
I am using the COBRAToolbox to do a flux balance analysis of a whole genome scale model ...
- 21
0
votes
1
answer
42
views
Target genes for piRNA
Where I can find a database or tool to give me the target genes of PIWI (piRNAs) in human?
I found one but works for worm like <...
- 1,791
0
votes
0
answers
56
views
Help me to calculate Heaps Alpha value from the roary pangenome pipeline result?
I need to know whether my pan-genome is open or closed. For that, I need to calculate the ...
- 109
1
vote
0
answers
21
views
Why does gene count increase drastically after scaffolding with Hi-C data?
I have a conceptual question.
I have a diploid, outcrossing plant genome assembly of ~1.1 Gb size. The original assembly is generated from PacBio reads. After genome annotation with the ...
- 75
2
votes
1
answer
100
views
Querying The Cancer Genome Atlas (TCGA) for gene expression
The UCI gene expression cancer RNA-Seq dataset (https://archive.ics.uci.edu/ml/datasets/gene+expression+cancer+RNA-Seq) consists of 5 types of cancers (BRCA, KIRC, LUAD, COAD, PRAD) and gene ...
- 23
2
votes
2
answers
144
views
Is this a complete bacterial genome? (CheckM)
I've been working with a bacterial genome assembly. The initial stats look pretty good (3 main pieces, depth ~50x [20x - 100x] if I calculated it right). I'd expect this to be a virtually complete ...
- 777
2
votes
1
answer
339
views
How to run MD5 check for multiple fastq files in different subdirectories?
I have received Illumina sequencing reads for 100 samples. I have 8 R1.fastq.gz and 8 R2.fastq.gz files for each sample in each ...
- 75
1
vote
1
answer
492
views
tensorflow nn_model for DNA sequences: Matrix size-incompatible: In[0]: [2,1], In[1]: [784,300]
Hope anyone can help a beginner here. I'm building a proof-of concept tensorflow classifier for DNA sequences. However, the NN model does not let through train and test vectors saying the matrix size ...
- 11
0
votes
0
answers
89
views
How to manually curate a genome assembly for sequence variation or error?
I have a PacBio HiFi assembly of 1.1 Gb from a heterozygous species. I have aligned this assembly against a reference genome which is around 0.9 Gb. I can see that there are quite a few INDELs, ...
- 75
3
votes
2
answers
105
views
How to get the number of complete phage genomes available on ncbi?
I am looking to establish the total number of complete phage genomes available on NCBI.
I am not looking for any specific type, but want to understand the total diversity available.
How can I perform ...
2
votes
2
answers
73
views
Find the genomic locations of exons of certain human gene
So I know Consensus CDS has the genomic locations of exons, however, only of the exons in the CDS. Everything before or after is not included.
NIH variation viewer has those (kinda-) radio buttons ...
- 23
-1
votes
1
answer
34
views
Evolutionary speaking, what are the constraints on having new mutations?
Evolutionary speaking, some mutations lead to better fitness of an organism and its adaptation to the environment changes, but I was wondering if there are some constraints on having new mutations. ...
- 107
1
vote
2
answers
222
views
What is the most appropriate way to find the most recent common ancestor between two distantly related species
I want to specifically find the common ancestor between a lobster and a humans. I suspect it was an aquatic worm of some description. But I want to know about the nervous system of this common ...
1
vote
1
answer
102
views
Comparing homozygosity of k-mer plots
Attached are two kmer plots from two closely related species. Is that safe to say that the one on the left has higher homozygosity than the one in the right k-mer plot, due to a low to almost flat ...
0
votes
1
answer
43
views
Identifying Novel Species with ANI Taxonomy
I have sequenced a bacterial genome and after I checked the ANI scores, the scores against similar bacterias were 79 and lower. How can I confirm if this bacteria is a novel species/strain? Is there ...
- 121
0
votes
1
answer
190
views
Understand this Kmer plot from Merqury?
The attached figure is generated based on Illumina reads from multiple individuals compared to genome assembly. Looks like there are a lot of kmers are reads only (grey colored). Also, a blue peak (2x ...
1
vote
2
answers
52
views
Offline Interactive Genome Map
I use the CGView Server for constructing circular genome maps. It is really a convenient tool for me to explore circular genome maps because it's interactive. I can search a certain region and it can ...
- 121
0
votes
1
answer
80
views
Extracting Sequence Unique to a Certain Genome
I have two genomes with high similarity. But I found out that one genome has a longer sequence of about 200 kbp. I try to align the sequence with Mauve. How can I extract this unique 200 kbp?
- 121
2
votes
1
answer
32
views
Organelle genome genes
How can I get all proteins coded by mitochondrial genome or plastid genome from UniProt. I see UniProt has multiple filtering and search parameters but I could not distinguish between proteins/genes ...
- 73
0
votes
1
answer
342
views
How to Find List of Genomic Coordinates from List of Genes using UCSC Website?
Suppose I have a list of genes and I want to find their genomic coordinates. To be specific and to keep things consistent with the rest of my project, I'm using hg38 on the website genome.ucsc.edu. So ...
- 13
5
votes
1
answer
86
views
Why does the SARS-Cov2 genome has letter t [duplicate]
ATTAAAGGTT TATACCTTCC CAGGTAACAA ACCAACCAAC TTTCGAT... is part of the 5'UTR of genome of an RNA virus SARS-Cov-2.
RNA contains letters ...
- 61
0
votes
2
answers
1k
views
Create BED file from list with gene symbols
Is there a way to generate a BED file of a list of gene symbols for a given genome?
- 323
-1
votes
1
answer
84
views
Where to find rat genome for STAR alignment?
I'm trying to find the rat genome to carry out STAR alignment for several FASTQ files. Where would I find these files?
0
votes
2
answers
158
views
How to know if the DNA sequence has been assembled and why is it important to know how it was assembled?
I have downloaded my FASTA format files, that have the DNA sequences of the coding region of the genes and the DNA sequence of the complete genome, from NCBI. How can I recognize if these sequences ...
0
votes
1
answer
694
views
Converting Fastq files to Fasta files on Ubuntu
I am new to bioinformatics and programming. I tried to convert my fastq files to fasta file, but got this. What could be wrong?
...
- 1
1
vote
1
answer
66
views
How to quantifiy of specific genes from shotgun metagenome?
I have googled a "lot", couldn't find any specific answer to the question. So, I am here seeking for your guidance. My question is similar to this. I have several metagenome (n=30). But for ...
- 111
0
votes
0
answers
30
views
Search for specific sequence in group of bacterial genomes
I have a group of L. plantarum genomes. I identified the pangenome using Roary. From the pangenome, I was able to identify unique 'genes' for some of the strains. One of these (actually several were, ...
- 43
0
votes
1
answer
30
views
Comparing the genetic similarity of human genome and mice/monkey genome - where can I find a list of genes?
I am trying to start a side project where I verify the claim that human and mice/monkey share 99% of their genes. I hope to learn the basics of doing genomic analysis.
I would like guidance about ...
- 103
1
vote
2
answers
589
views
How to identify the DNA sequence of a gene in the complete genome sequence FASTA format file?
I need to identify the sequence of a gene in the complete genome sequence . I thought it was simple, instead it is not a straightforward task !
My method was the following: I downloaded the FASTA ...
2
votes
2
answers
195
views
What is the file .mums created by mummer?
I am doing genome sequence alignment using MUMmer, in particular I want to do a dotplot with mummerplot. So the passages that I did are:
1.create a file .mums with the following command line:
...