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.

Filter by
Sorted by
Tagged with
2 votes
0 answers
9 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 ...
user avatar
0 votes
1 answer
27 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 <...
user avatar
  • 1,513
0 votes
0 answers
14 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 ...
user avatar
  • 109
1 vote
0 answers
17 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 ...
user avatar
2 votes
1 answer
61 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 ...
user avatar
  • 23
2 votes
2 answers
41 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 ...
user avatar
  • 539
2 votes
1 answer
78 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 ...
user avatar
1 vote
1 answer
90 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 ...
user avatar
  • 11
0 votes
0 answers
29 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, ...
user avatar
3 votes
2 answers
82 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 ...
user avatar
2 votes
2 answers
45 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 ...
user avatar
-1 votes
1 answer
32 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. ...
user avatar
  • 107
0 votes
2 answers
79 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 ...
user avatar
0 votes
0 answers
17 views

Why does Nucmer Genome alignment produce faulty Dot plot?

I am using Mummer v.3.23 to make a dot plot between two genomes of the same species. One genome was assembled using PacBio and one from Dovetail. None of them is ...
user avatar
1 vote
1 answer
53 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 ...
user avatar
0 votes
1 answer
23 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 ...
user avatar
  • 121
0 votes
1 answer
54 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 ...
user avatar
1 vote
2 answers
44 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 ...
user avatar
  • 121
0 votes
1 answer
40 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?
user avatar
  • 121
0 votes
0 answers
18 views

Is there a newer/better book than Genome-Scale Algorithm Design available? It's from 2015

I'm considering buying the book "Genome-Scale Algorithm Design: Biological Sequence Analysis in the Era of High-Throughput Sequencing" by Veli Mäkinnen and others. However, as this book was ...
user avatar
  • 129
2 votes
1 answer
31 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 ...
user avatar
  • 73
0 votes
1 answer
77 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 ...
user avatar
  • 13
6 votes
1 answer
83 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 ...
user avatar
-1 votes
2 answers
175 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?
user avatar
-3 votes
1 answer
38 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?
user avatar
0 votes
1 answer
46 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 ...
user avatar
0 votes
1 answer
319 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? ...
user avatar
1 vote
0 answers
35 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 ...
user avatar
0 votes
0 answers
25 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, ...
user avatar
  • 33
0 votes
1 answer
26 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 ...
user avatar
1 vote
2 answers
121 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 ...
user avatar
2 votes
2 answers
93 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: ...
user avatar
0 votes
1 answer
374 views

STAR aligner multiple fastq files

I’m using STAR to align fastq files from SMART-seq2. I have raw data folder containing sub-folders with samples names the sub-folders each contain fastq file. How can I make a bash command in order ...
user avatar
4 votes
0 answers
27 views

Do RepeatModeler results contain functional domains?

The repeat families predicted by RepeatModeler contain known transposable elements (TEs) and unknown ones. How do we know whether some of these may actually: within a functional domain of a gene or ...
user avatar
0 votes
3 answers
45 views

Genome QC + Assembly Pipeline semantics

I’m trying to create a pipeline for genome assembly. How best can I “redirect/pipe” from existing fasta files (or files in general) to other steps of the pipeline? I was thinking of going from the SRA ...
user avatar
2 votes
1 answer
76 views

What tools can I use to look up my alleles, genotypes, or phenotypes from my sequenced DNA (WGS)?

(I am not a Bioinformatics expert, please forgive and educate me if I've used any wrong terms or assumptions here) I bought a "Whole Genome Sequence" kit, which gave me the following VCF ...
user avatar
  • 121
3 votes
2 answers
72 views

Gold standard benchmark

This page is claimed to contain a gold standard benchmark for viral genome assembly. https://github.com/cbg-ethz/5-virus-mix The claim is here: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5411778/ ...
user avatar
0 votes
0 answers
84 views

How to filter a genome assembly consistsing of a large number of contigs?

I did some de novo genome assemblies with Illumina PE data using SPAdes, whereas most of them consisting of a large number of contigs(>1000). I have several questions below. Do we need to filter ...
user avatar
  • 59
1 vote
0 answers
73 views

How to extract certain genes including non-coding regions?

I would like to look at the non-coding regions of some genes. I have around 1000 full genome assemblies and I was able to extract the nucleotide sequences for certain proteins with ...
user avatar
  • 1,181
1 vote
1 answer
30 views

Mapping protein refseq to Gene ID

I have a protein refseq (eg, NP_000029). How can I get the corresponding gene ID (ag, APC) from NCBI using an R package?
user avatar
1 vote
1 answer
95 views

how to train a gene dataset with a nearest shrunken centroid classifier?

I have a data file named "geneexp.csv". the data contains information about gene expression of three different cell types (CD4 and CD8, CD19) I want to classify cells by performing the ...
user avatar
1 vote
2 answers
217 views

Is it possible to filter contaminated reads for raw PacBio sequences (not HiFi reads) before assembly?

De novo genome assembly for non-model organisms face the issue of bacterial contamination. For assembled contigs with mostly bacterial-like sequences (based on BLAST search), the entire contig can be ...
user avatar
2 votes
1 answer
105 views

Merging several fq.gz files or R1 and R2 classes into a single one

I have the following files: ...
user avatar
2 votes
2 answers
150 views

How can I subset WGS data to the level of WES variants?

I would like to compare mutational signatures1 in patients from different studies, however some studies are based on exome seq (i.e. ~20,000 coding variants) and some are from whole genome seq (i.e. ~...
user avatar
0 votes
1 answer
53 views

Is it possible to convert BAM file from one genome assembly to the other?

I Have multiple BAM files that are referenced to UCSC genome assembly GRCh37/hg19 that are read in different time frames. Now, I am planning a different studies that require assembling all the data ...
user avatar
0 votes
1 answer
478 views

kraken2 to OTU table creation

I did run kraken2 and I get two kinds of output How do I generate OTU table from these outputs for phyloseq usage?
user avatar
  • 11
1 vote
2 answers
52 views

Why does Coronavirus genome contain poly A region at the 3' end in GenBank? [duplicate]

If you look at the Coronavirus genome below, it contains poly-A region at the end. I know that the viral RNA contains a poly-A tail. But, the genome in GenBank is stored in cDNA format. That means the ...
user avatar
  • 709
1 vote
1 answer
51 views

Get raw genome sequence

I'm currently working on testing some matching algorithms for strings. I would like to do some tests on raw genomic data as I expect different results from random strings given the lower entropy. How ...
user avatar
  • 149
1 vote
1 answer
53 views

How do I download the mitochondrial haplogroup datasets for human genetics online?

I seem to have landed on the mitomap.org site, but I don't know what to make of it or what do with it / how to get the genomes onto my computer. It sounds like the genomes are stored in GenBank, but ...
user avatar
  • 133
1 vote
1 answer
425 views

How to fix RepeatMasker.lib.nsq missing error when running RepeatModelor?

I have run RepeatModelor and I don't see the classified output file. I got this error. ...
user avatar