Questions tagged [annotation]

The process of attaching contextual information to a biomolecular (DNA, RNA, or protein) sequence, such as the location of genes, promoters, or repeats in genomic DNA, or the function of a protein.

Filter by
Sorted by
Tagged with
1
vote
0answers
42 views

UTR features in gbk files

I have a gbk files for a viral genome that I have used biopython to add features of the 3' and 5' ...
3
votes
5answers
184 views

Import a tab-separated file with differing numbers of elements in each row; prokka output

I am using prokka to annotate a bacterial genome: prokka ecoli.fa Prokka is outputting a tab-separated file (called PROKKA_12142017.tsv) with differing numbers of ...
0
votes
1answer
32 views

Filtering VEP annotation file

I am filtering a VEP annotated vcf, trying to maintain just those variants classified as deleterious by SIFTand as damaging (probably or possibly included) by ...
1
vote
0answers
14 views

Finding the relevant annotations for biomarkers

I am looking to find a way to map a list of biomarkers of interest, say genes or bacterial functions, to existing annotation databases. Usually this is done in reverse, as such that based on ...
0
votes
1answer
59 views

Given a .gb file and a locus - how to get relevant annotations in Python?

Given a .gb file and a specific locus in the genome - how can I retrieve the relevant annotations in Python (i.e., annotations that include that locus)? I could ...
2
votes
0answers
20 views

Finding original papers for gene annotations

I am wondering if there's a database to look up the annotation history of an arbitrary gene/ORF, with links to the primary literature for each step of the annotation. I have found that neither NCBI ...
0
votes
1answer
33 views

Inheritance annotation for nsSNPs

I am relatively new to SNP analysis. Is there a database to find annotation about the Mendelian inheritance of SNPs? I have a small list of nsSNPs (non-synonymous SNPs) and I need to find how they are ...
2
votes
1answer
283 views

Why doesn't the mouse GRCm38/mm10 refGene genome annotation file contain non-coding transcripts?

I am new to Bioinformatics and I am exploring the refGene.txt files from the UCSC genome annotation database for several species. My question concerns the Dec. 2011 (GRCm38/mm10) assembly of the ...
1
vote
0answers
28 views

BAM file filteing to remain best isoform

I ran HiSat2, MarkDuplicate, removed reads with the lower quality score than 40 and finally only kept properly paired reads. After the BAM filtering steps, I used the Scallop results with TransDecoder....
0
votes
1answer
111 views

How to add cluster name to metadata in Seurat?

I'm working on a Seurat object and want to name the clusters according to 2 values alone (yes/no). So I want to add a new column to metadata and annotate the clusters (UMAP) with it. ...
1
vote
1answer
34 views

REMC SNP annotations for transcriptome prediction - epiXcan

I have been trying to use a pipeline for predicting gene expression for a target gene (Although you can do it for many) - the name of this method is epiXcan. Here is a link to the paper: https://www....
2
votes
0answers
39 views

What is tRNAscan-SE bit score?

I have used the command-line version of tRNAscan-SE by Lowe et al. (1997) on the reference genome sequence of C. elegans to predict tRNAs. The programme outputs a table with tRNA locations, amino acid,...
0
votes
1answer
44 views

Changing active.ident in Seurat

Im trying to change the active.ident to another column in metadata but this error keeps popping up! I recently upgraded to R version 4.0.2 from 3.6.1 The older version was working but the new one isn'...
0
votes
2answers
36 views

differential analysis of chip-seq data

I have several sets of chip-seq data. I called the peaks using Macs2. I am pretty new to the field and I will appreciate any help. I wanted to annotate the peaks and see which peaks are shared between ...
2
votes
1answer
80 views

renaming IDs in gff3 with BCBio.GFF

I wrote a script which should changes IDs in a GFF3 file. Unfortunately, the below script has two problems. It attaches the new ID to Parent which leads that ...
1
vote
1answer
71 views

Decontaminating RNA seqs for de novo transcriptome assembly and annotation of novel eukaryotes

I have raw paired-end RNA-seq reads for two novel eukaryotic species. Some background: the reads represent a copepod (arthropod) species each. The mRNA for each read set was obtained by extracting ...
0
votes
0answers
35 views

Gffcompare issue: 0 reference transcripts loaded

I am trying to use gffcompare to compare my assembled transcriptome to a reference gtf that contains information about small open reading frames (sORFs). The reference gtf was obtained by processing ...
0
votes
0answers
6 views

How does one go about developing GCMS metabolite annotation workflow?

I’m really new to untargeted metabolomics and I need to annotate a data set from a standard quadrupole GCMS. The person on the project before me tried to do the annotation manually, cleaning up the ...
0
votes
0answers
19 views

Determining a Gene's Unknown Function

I'm a beginner in this field (so pardon my rudimentary knowledge/explanation) trying to understand ways on how I could potentially identify genes of unknown function. I am trying to get as close as I ...
2
votes
2answers
129 views

RepeatMasker annotation

I collected my own denovo library using many different tools Pipeline (structure-based ) and I updated all the headers corresponding to RepeatMasker format, to obtain Repeat coverage in the genome and ...
3
votes
1answer
190 views

What does the number mean in an HGVSp annotation?

Let's take the example p.Arg452Pro that I got from an annotated VCF file that contains an HGVSp column. What does ...
1
vote
0answers
18 views

How to (efficiently) search nuccore to retrieve only annotated sequences?

I wish to efficiently search nuccore to retrieve only annotated sequences. By "annotated" I mean that the nuccore entry contains annotation information. An equivalent definition (if I understand ...
1
vote
0answers
19 views

Is there a built-in parser for the textual feature table that Entrez.efetch returns?

Entrez.efetch helpfully gives me the feature table of a nuccore entry. If I understand correctly, efetch can only return the ...
0
votes
0answers
12 views

miRNA Homologues (Orthologues) between Mouse and Human

Where may I find a list of mature miRNA names that are orthologues (homologues) between human and mouse? It has to be the mature miRNA, and not just the precursor/hairpin/gene names - i.e. I need the -...
3
votes
4answers
173 views

Purpose of ### (3 consecutive pound signs / hashtags / octothorps) in GFF3

I downloaded the annotation of the C. elegans genome in GFF3 format from Ensembl. I typed the following command, hoping to get the header of the file (lines starting with ...
2
votes
3answers
222 views

Modifing a GFF3 file and writting to a new file

I modified the GFF3 file with the below code: ...
1
vote
1answer
69 views

Getting sequences of one transcript per gene out of annotation and genome

I got a genome (data/genome/genome.fasta) and braker-based genome annotation (data/genome/annotation.gff3), now I would like to ...
1
vote
1answer
66 views

Genome annotation of a Bacillus strain

I finished the sequencing of the genome of a biocontrol activity strain Bacillus amyloliquefaciens by Illumina hiseq. I have already gone through some papers but there were no details of processing in ...
0
votes
0answers
9 views

Batch annotation with RAST

I have been looking into ways to batch annotate my bacterial genomes via RAST, however I cannot access the tutorials on blog.theseed.org as my browser prevents connection for security reasons. Does ...
1
vote
0answers
15 views

SwissProt Proteins into my Jbrowser

How is it possible to get the SwissProt Proteins into my Jbrowser in the same way like shown below? Thank you in advance,
0
votes
1answer
69 views

Ensembl protein identifiers from different assemblies

In order to obtain orthologs for human genes, I am working with InParanoid version 7.0. This version of the database uses Ensembl Protein Identifiers (ENSP) from ENSEMBL version v54, which is based on ...
7
votes
1answer
814 views

Where to download JASPAR TFBS motif bed file?

I am interested in determining if any transcription factor binding site motifs are enriched in some BED files from a DNA methylation experiment. I am looking for a database that has BED Files ...
1
vote
2answers
273 views

How to retrieve the subcellular location info from uniprot?

I have a list of proteins with their gene names/symbols. I would like to identify which of them are extracellular proteins or plasma membrane proteins. I can get this information in UniprotKB by ...
3
votes
2answers
140 views

Generate SNP/indel annotation in Arabidopsis

I have a genome-wide list of germline SNPs and short indels for Arabidopsis thaliana, which I generated with Varscan. Regardless of the tool used to generate them, I would like to annotate them, i.e. ...
2
votes
1answer
55 views
4
votes
2answers
116 views

How to predict stop codons in Illumina reads?

I have Illumina MiSeq paired-end reads from 150bp amplicons mapped to my reference genome (> 1000X coverage). These reads have indels that may or may not induce frameshifts. If the indel induces a ...
2
votes
1answer
54 views

Adding an attribute to GFF3 file

I failed to add Note=Gene description to mRNA attribute with the below code: ...
1
vote
1answer
53 views

Are there databases to annotate non-coding mutations?

I would like to ask if there are some known databases that allows to annotate non-coding mutations in a given region for the human genome? Preferably, if the database is downloadable so I am able to ...
2
votes
4answers
119 views

Merging two .txt files

I have two .txt files and let's call it 'File A' and 'File B.Here are examples of what file A and file B looks like: FILE A ...
5
votes
3answers
3k views

Retrieve detailed gene descriptions

Given a list of gene IDs, how do you retrieve the gene description, summary and other detailed information in R?
40
votes
5answers
8k views

Feature annotation: RefSeq vs Ensembl vs Gencode, what's the difference?

What are the actual differences between different annotation databases? My lab, for reasons still unknown to me, prefers Ensembl annotations (we're working with transcript/exon expression estimation)...
2
votes
1answer
71 views

hgsql not found when running a script despite being accessible on command-line

I am working on macOS. I am following the steps described here to lift over an annotation from the reference genome of C. elegans to another assembly of C. elegans. I have run this code: ...
2
votes
1answer
24 views

In SBOL, when should I annotate a DNA sequence vs. making sub-components?

When representing a DNA sequence in SBOL as a ComponentDefinition, you can mark things like promoters and coding sequences in two different ways: either as a ...
1
vote
1answer
73 views

ANNOVAR, 1000g database calling

I have recently downloaded Annovar, and succesfully filtered my files with a few databases, downloaded from Annovar. But when I try to do the exact same thing with the 1000 genomes project database, I ...
8
votes
2answers
278 views

Comparing two genome annotations

I have one-dimentional array (human genome). Also I have two annotations for it, we can think about them as different peaks (it's nucleosome and secondary structures). How can we find correlation and/...
0
votes
1answer
96 views

What is the `table` parameter in the UCSU ldHgGene tool?

I want to lift over the annotation from E. coli release 29 to E. coli release 42. I have followed the steps for an annotation liftover mentioned here. I got the chain file for the final step and ran: ...
1
vote
1answer
12 views

SnpEff gene annotation with @ such as PCDHB@

I'm using snpEff to annotate some VCFs and getting some spurious annotations that don't have transcript ids. Some of these annotations have a gene name followed by ampersand such as PCDHB@ and I'm not ...
2
votes
1answer
47 views

No refseq transcripts in 2nd half of Y chromosome for HG38

Obviously there is a gap in my understanding of how sex chromosomes are annotated. I've been working on some CNV calls and noticed that most of the 2nd half of ChrY is missing annotations and no reads ...
0
votes
2answers
42 views

Evidence of emergence of “genuine” novel protein interactions by aa mutation [closed]

Anyone knows of an example of a protein that, without coming from a recent duplication event, underwent a mutation(s) that caused it to have a novel interaction with a new ligand, substrate, other ...
2
votes
3answers
962 views

Annotation with Prokka or RAST.

I was experimenting Prokka and RAST annotation tools. So, I took a well-annotated swinepox virus genome from genebank (NCBI Reference Sequence: NC_003389.1). I ran those sequences on Prokka and RAST ...