# 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.

83 questions
Filter by
Sorted by
Tagged with
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,
43 views

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 ...
34 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 ...
28 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 ...
43 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 ...
44 views

### where to find pathogenicity metadata for E.coli genomes

I have a list of about 900 E.coli genome ids (Genbank plus NCTC ids), e.g.: ...
121 views

### Modifing a GFF3 file and writting to a new file

I modified the GFF3 file with the below code: ...
40 views

### Adding an attribute to GFF3 file

I failed to add Note=Gene description to mRNA attribute with the below code: ...
50 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 ...
109 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 ...
143 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 ...
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 ...
10 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 ...
56 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 ...
64 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: ...
46 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 ...
40 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 ...
40 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 ...
120 views

### A GTF annotation for GRCH375d version

I have been given a bunch of .bam files, that come from the alignment on GRCh37d5 of human genome by Star tool. I want to extract raw read counts from these bam files by FEATURECOUNTS tool so I need a ...
13 views

### The best tool to count the number of insertion sequences in a set of bacterial genome

I have a set of 414 GenBank genomes and I need to count a number of insertion sequences (ISs, transposases) in each of them. Is there any traditional way to perform such an analysis? I came up with ...
192 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 ...
106 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 ...
59 views

### scRNASeq expression matrix with decimal values

I am trying to replicate some results of a scRNASeq experiment and, when I looked at the data provided by the author, I noticed that some of the counts in the expression matrix are represented as ...
368 views

### Identify non-coding regions from a genome annotation

I have this GTF file and I use the command below on a Linux machine to extract the coding regions of the genome: ...
148 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 ...
245 views

### Derive a GTF containing protein coding genes from a GTF file with Exons and CDS

Why I need a compatible file I’m trying to run velocyto with the R package to analyse RNA velocity (cell trajectories) with single cell RNASeq data. I have performed single cell analysis from 10x ...
64 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: ...
63 views

### Accessing UCSC genome via ssh results in a validation error

I am working on macOS High Sierra. I am following the steps described here to lift an annotation over from one version of a genome to another. I am now just using the example genomes provided in the ...
39 views

148 views

### Mapping RefSeq-ID or NCBI-ID (Entrez-ID) to EC number

I annotated my bacterial genomes using the new NCBI Prokaryotic Genome Annotation Pipeline and now, I want to annotate EC-numbers. In the master annotation file (.gff), I get the RefSeq-ID of every ...
315 views

### Transcript Coordinate Ranges to Genomic Coordinates

I have 2 GFF3 files: Features using transcript IDs as the landmarks. i.e. "CDS" feature types using coordinates from transcript space. Features using chromosome IDs as the landmarks. i.e. "exon" ...
86 views

I have just started using BlastKOALA KEGG which has been useful in annotating (aminoacid) sequences. This is their website: https://www.kegg.jp/blastkoala/ When you get results, there are links for ...
133 views

### Get genomic coordinates using GenomicFeatures by HGNC gene names

I want to get coordinates of human genes from my list (consisting of hgnc genes id) using GenomicFeatures and TxDb.Hsapiens.UCSC.hg19.knownGene R packages from Bioconductor. ...
37 views

### How to generate Genbank format out of paired end fastq data?

I've unmapped cleaned paired end sequence data in fastq format of a bacterial genome. I want to get a sequence data in Genbank format in the end. What are the exact steps that I have to follow in ...
137 views

### RATT works on example bacterial sequence and other bacterial genome but not on C. elegans genome and annotation

I am trying to use a tool called RATT described in this paper and available here. I ran the example using the following command: ...
24 views

### Variants in Exome that affect miRNAs related to Schizophrenia

I have 3 exomes results (VCF files), I would like to study variants (SNP, INDEL) that can affect the miRNAs, which are related to a complex disease (Schizophrenia). What is the best way to assess the ...
619 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 ...
79 views

### Detect gene duplication events from a whole genome in FASTA format

Are there any tools for finding duplications from a whole genome in FASTA format? Let's say the genome has gene A duplicated somewhere in the genome and gene B triplicated in tandem elsewhere in the ...
48 views

### TRAL does not find “Phobos result file”

I want to use TRAL to annotate tandem repeats in the reference genome of Caenorhabditis elegans. For this, I need to install some external software, such as Phobos. I've downloaded Phobos and I am ...
57 views

### RNA seq fasta file annotation from alignment to reference matches

I've got a fasta file with some RNA seq data and another csv file with the output from plast where I've aligned it to a reference using plastn. I'm struggling with figuring out a command to append my ...
29 views

### What is tRNAscan-SE bit score?

I have used the command-line version of tRNAscan-SE on the reference genome sequence of C. elegans to predict tRNAs. The programme outputs a table with tRNA locations, amino acid, anti codon, intro ...
845 views

### tRNAscan-SE error: FATAL: Unable to find /usr/local/bin/cmsearch executable

I have downloaded tRNAscan-SE from here. After decompressing and untaring the file, I installed it using: ./configure make make install When I type ...
28 views

### Batch detection of CRISP proteins in fasta file

Probably a naive question. I am inexperienced. I am interested in identifying potential CRISP (Cysteine-rich secretory proteins) in a certain tissue transcriptome (ca. 20k sequences in fasta). I ...
232 views

### Expand granges object different amounts upstream vs. downstream

I am attempting to get gene regions and their immediate neighborhoods using Bioconductor GenomicRanges packages. It is very easy for me to extend the regions around genes if I want to extend it the ...
93 views

### GFF file too big to load into R

I want to load a GFF (annotation) file into an R data frame in order to extract some information from it, e.g. the location of transposons in the genome. I used the following code: ...
437 views

### update dbSNP ID

I have used dbSNP build 138 to tag 'rs' ID in my VCF file. But later I came to know that now dbSNP build 150 has been released. I want to update the rs IDs in my VCF file. I tried but I found a ...
390 views

### get gene lines from gtf file

I would like to retrieve gene lines from a GTF file for which I only have exons & transcripts lines (output from Cufflinks) and alternative splicing possible. I need gene lines for compatibility ...
198 views

### Keep Format and Individual fields when annotating VCF with VEP

I'm currently updating my Variant Calling Pipeline by switching the VCF annotating software from Annovar to VEP for a variety of regions, not least how easy it is to annotate with HGVS notation and ...