Skip to main content

Questions tagged [sequence-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.

37 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
6 votes
0 answers
65 views

Why does repeatmasker annotate transposons of length 1?

I am working on Mac OSX. I have downloaded RepeatMasker from www.repeatmasker.org. I have downloaded RepBase from www.girinst.org. I have download the WS266 version of the C. elegans genome from ...
Biomagician's user avatar
  • 2,459
4 votes
0 answers
36 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 ...
Life_Searching_Steps's user avatar
3 votes
0 answers
258 views

How to remove redundancy from a gtf file?

I have an annotation file. I would like to remove redundancy, as shown in the example (in the real file, I have a lot of these redundant cases). I would like to consider only one of the following ...
Marco's user avatar
  • 161
3 votes
1 answer
239 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. ...
lizaveta's user avatar
  • 203
2 votes
0 answers
17 views

Confusing naming in MANE GFF file

I'm working with the MANE v1.2 annotation (with RefSeq keys) and am puzzled about the ID naming convention being used. Specifically, there are several transcripts (18 by my count) that have hyphenated ...
merv's user avatar
  • 651
2 votes
0 answers
251 views

PROKKA setupdb / cmpress error

I am trying to run --setupdb in Prokka, but I'm getting this error that I can't seem to find the answer for. When I run prokka --setupdb I get ...
Carol Luz's user avatar
2 votes
1 answer
75 views

Sanger sequencing annotation error

I am a student in a Cancer lab. Working with sanger is new to me. While analyzing a report we found an insertion that has not been reported in any databases so far, we were working on checking if the ...
user avatar
2 votes
0 answers
55 views

In VEP annotation, how is the codon field interpreted?

After annotating with VEP a VCF file, we obtain different fields. One of them is called Codons which represents the affected codon in the transcript of the gene. ...
user324810's user avatar
  • 1,115
2 votes
0 answers
25 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 ...
biohacker's user avatar
  • 131
2 votes
0 answers
199 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: ...
Biomagician's user avatar
  • 2,459
2 votes
0 answers
78 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,...
Biomagician's user avatar
  • 2,459
1 vote
0 answers
25 views

How to compare CDS, CNE, miRNA, UCNE across new related genomes?

I have 5 sets of raw Illumina whole genome shotgun data from 5 different species, and 3 reference genomes that are assembled/annotated, I want to produce a similar table to Table S2 of this paper: ...
JohnDoe23's user avatar
1 vote
0 answers
29 views

David tool online not converting most genes

I'm using DAVID tool online, in order to perform functional annotation chart. I extracted the most relevant genes from a dataset (GSE73304), but have the problem that DAVID is not able to understand ...
Riesz98's user avatar
  • 41
1 vote
0 answers
159 views

Prokka + table2asn (NCBI submission)

I have a couple of MAGs of interest that I am trying to submit to the NCBI (with annotations). The process of submitting the annotations is rather challenging to me. I have the output of Prokka, which ...
Laura's user avatar
  • 909
1 vote
0 answers
19 views

How Can You Reformat an OTU Table from qiime2 pathway to Include Sequence IDs?

I am trying to run Tax4Fun in MicrobiomeAnalyst, but am running into trouble with my table formatting. From the qiime2 pathway, I have an OTU, metadata and taxonomic table formatted like so: In ...
shannonevelyne's user avatar
1 vote
0 answers
793 views

conversion of a gbff file to a gff file

I want to get the annotations in gff format for a assembled genome. But for some species there is no available annotation in gff format in NCBI but only in gbff format. And I don’t know how to convert ...
Kishor Kumar Sarker's user avatar
1 vote
0 answers
137 views

snpEff annotation and variant in a gene

I'm using snpEff to annotate genetic variants. However, the annotation often returns more than one annotation (see below for an example): ...
Anti's user avatar
  • 153
1 vote
0 answers
25 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 ...
Anik Dutta's user avatar
1 vote
0 answers
23 views

How to get a list of all human morphogens?

Gene Ontology provides many genes annotated as taking part in various morphogenesis processes but I want to get a list of all morphogen coding genes specifically. Uniprot does not have "...
Sashko Lykhenko's user avatar
1 vote
0 answers
33 views

How to train annotations tools

I would like to train Augustus, SNAP and GlimmerHMM. I found protein sequences in GenBank and orthodb.org. Furthermore, I found HMM files on busco-data.ezlab.org. ...
user3523406's user avatar
1 vote
0 answers
138 views

SnpEff not correctly annotating multiple-nucleotide polymorphisms

I have some SARS-CoV-2 sequencing data that I'm trying to annotate with SnpEff, however SnpEff doesn't appear to be recognizing multiple adjacent SNPs within the same codon and correctly calling them ...
Set's user avatar
  • 241
1 vote
0 answers
494 views

Prokka error message

I am trying to run Prokka but the following error message appears: Could not run command: makeblastdb -hash_index -dbtype prot -in /cluster/software/prokka/1.13.7-gompi-2019a/db/kingdom/Archaea/sprot ...
ahmedbar's user avatar
1 vote
0 answers
88 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' ...
Lamma's user avatar
  • 175
1 vote
0 answers
17 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 ...
alexanderjansma's user avatar
1 vote
0 answers
38 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....
user977828's user avatar
1 vote
0 answers
157 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 ...
nattzy's user avatar
  • 31
1 vote
0 answers
155 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 ...
Oren Milman's user avatar
1 vote
0 answers
76 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 ...
Oren Milman's user avatar
1 vote
0 answers
17 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,
user977828's user avatar
1 vote
0 answers
20 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 ...
Some student's user avatar
1 vote
0 answers
218 views

Downloading full search results from BlastKOALA?

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 ...
Scientist's user avatar
  • 111
1 vote
0 answers
102 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 ...
Biomagician's user avatar
  • 2,459
0 votes
0 answers
18 views

The current status of human genome annotation

I was wondering if you can suggest something on the current status of the human genome annotation efforts. I know about ENSEMBL/HAVANA/GENCODE, and I've read the relevant papers; however, I was ...
Alexander Predeus's user avatar
0 votes
0 answers
42 views

What tool(s) can I use to use a .faa file to feed into JSMol?

JSMol, the molecular modelling software, requires PDB file format to work and render the graphics. Currently, I am writing an annotation software - backend in Perl - that takes in a .faa file (fasta ...
Harish Prabhakar's user avatar
0 votes
0 answers
93 views

An error with annovar

I want to install annovar but in one particular part I always get this error ...
Zizogolu's user avatar
  • 2,160
0 votes
0 answers
32 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, ...
Mark's user avatar
  • 143
0 votes
0 answers
50 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 ...
Andrew Hollins's user avatar