Features in GenBank files · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Features in GenBank files

Hi,

I'm trying to retrieve GenBank files for genes with older assembly version (GRCh37) using eutils Bash

For example: efetch -db nucleotide -id NC_000012.11 -seq_start 53700940 -seq_stop 53715954 -format gb and I noticed in theFEATURE section in the .gb file it does not have features such as mRNA, CDS region. It only contains source, does anyone know how can I retrieve mRNA, CDS altogether with FEATURES?

Thanks very much!

gene mrna genbank

• 1,556 views • link

updated 6.3 years ago by Sej Modha

• written 6.3 years ago by speycast •

1 answer

You can use -format gbwithparts for that.

efetch -db nucleotide -id NC_000012.11 -seq_start 53700940 -seq_stop 53715954 -format gbwithparts

• 0 views • link

written 6.3 years ago by Sej Modha

Thanks so much Sej!!

• 0 views • link

written 6.3 years ago by speycast •

Log in to answer this question.

More posts like this

Best way for finding annotated counterparts of unknown transcript after BLAST
written 5.4 years ago by mwittep •

Hi all! I am building a local pipeline in order to identify unknown transcripts. One part of this pipeline is identifying if the unknown sequences …
Is there a way to do a batch extraction of DNA Coding Sequences from List of Protein Accessions (NC…
written 5.8 years ago by biobaker •

Hello! I have been banging my head against a wall for the past 2 days trying to get a fasta file of DNA sequences from …
Retrieve all annotated features in a GenBank record
written 5.9 years ago by al-ash

Is edirect supporting retrieval of all annotated features from a record? Example: I can retrieve all features annotated as "gene" in GenBank record of a …
Retrieve reference sequence GenBank for mRNA transcript of interest
written 6.2 years ago by speycast •

Hi, I want to retrieve reference sequence GenBank files for genes with transcript of interest using E-utils For example: GFM1 gene has two mRNA transcripts …
Retrieve exon ranges from E-utils
written 6.3 years ago by speycast •

Hi, I'm trying to retrieve exon ranges for a gene from GenBank format `efetch -id NG_007086.2 -format gb -mode xml -db nucleotide | xtract -pattern …
Code for splitting long string into Genbank record gives error
written 7.5 years ago by toth.joe •

I am trying to read a long amino acid string into a Biopython Genbank record object so that a Genbank file can be written. Here …
Strand not working in biopython efetch
written 7.9 years ago by tim.ivanov.92 •

I noticed i can't get efetch to retrieve sequence on a reverse strand: def fasta_from_chords(x): import sys sys.path.append('/uge_mnt/home/tim_ivanov/pythonlibs/biopython-1.71') from Bio import Entrez, SeqIO Entrez.email = …
Efetch from nuccore db using URL / Perl fails with Bad Gateway error for GenBank, but succeeds for …
written 8.6 years ago by senthil.ramamurthy •

Hello all, My previously working perl script to fetch a given region suddenly has failed. I use the Bio::DB::EUtilities library as shown below. use strict; …
Retrivie a large set of genomic sequence with custom range each
written 9.8 years ago by guedes.aureliano •

Hello, I have a list with genome ID's and genomic region of each gene, and ID of protein seq of each gene. INSDC CP002471.1 675883 …
How to pass multiple seq_start and seq_stop statements to NCBI efetch
written 9.9 years ago by massa.kassa.sc3na

Hi, </br> I wasn't able to find anywhere how to pass several seq_start and seq_stop optional arguments to list of queries for NCBI efetch. <br> …