Questions tagged [genbank]

NCBI sequence annotation flat file.

Filter by
Sorted by
Tagged with
3 votes
1 answer
65 views

What does this accession NCBI code mean: 6MWN_B?

According to this article, accession codes should consist from a combination of uppercase letters following a combination of digits. If this is a RefSeq, it can have a prefix as a combination of ...
user avatar
  • 355
2 votes
2 answers
76 views

How to properly annotate sequence variants and errors in a GenBank file format and how to keep track of successive versions of a GenBank file

EDIT By thinking over the issue, I realized it was mixing two linked but separate questions: Q1. How to handle the situation when there are divergences between a GenBank file and an actual or ...
user avatar
  • 121
2 votes
3 answers
87 views

How to programatically download SARS-CoV-2 fasta from NCBI/Genbank via API?

I would like to download a number of SARS-CoV-2 fastas from NCBI/Genbank using a web endpoint. I looked on their website but it's confusing. I have the accession numbers, but can't figure out how to ...
user avatar
1 vote
1 answer
127 views

Entrez (Biopython) esearch and efetch not returning sequence as expected

I'm trying to use Entrez (through Biopython) to download the sequence of a TMV replicase gene. I have the following code: ...
user avatar
1 vote
0 answers
45 views

Convert genbank file into tbl file [closed]

I want to submit some phage sequences to NCBI. So I want to transform a genbank file into a tbl file. Can someone help me please? Thanks
user avatar
  • 341
5 votes
1 answer
148 views

KeyError when getting features from a genbank file with biopython with some accessions but not others

I'm very new to python but i've been using it to extract the sequence of a gene from a genbank file. The issue is is that sometimes i'll get the output i want (prints the sequence to a file) and ...
user avatar
  • 53
0 votes
0 answers
26 views

How can I retreive specific qualifiers from genbank files in R?

I'm trying to extract the host field from the features table using the biofiles package. From what I understand, the function filter (from the same package) could be of help but I can't extract the ...
user avatar
  • 101
0 votes
2 answers
234 views

Splitting a GenBank file into smaller files

I have a genbank file that I have been trying to run in HMMER using the following command: ...
user avatar
2 votes
1 answer
64 views

How to programmatically classify a protein according to its genbank feature

Say I found an interesting protein in a genbank file, e.g.: ...
user avatar
1 vote
1 answer
42 views

Convert sequence in MS doc to fasta or genbank

How can I convert a sequence provided in a Microsoft doc file into a fasta or genbank format?
user avatar
  • 2,574
3 votes
0 answers
29 views

Template files for bioattributes and meta-data [closed]

I am in process of submiting my FASTQ files to the SRA database. I am a little confused about how to fill the biosample attributes and SRA metadata. I would really appreciate if you anyone could share ...
user avatar
  • 109
1 vote
0 answers
74 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' ...
user avatar
  • 143
1 vote
0 answers
128 views

GCF VS GCA or combination of both for pangenomic studies

I have just started to learn bioinformatics and pangenomics. So if this question seems to you pretty basic then I apologize in advance. As we know, NCBI for the genome database, there are two kinds of ...
user avatar
  • 87
1 vote
1 answer
51 views

Why does the Genbank phylogenetic tree of 4000 sars-cov-2 sequences only display 200 nodes?

I manually selected all 4000 sequences on the 20 pages of the SARS-CoV-2 sequences on genbank here. When I click on Build phylogenetic tree, it only has 200 nodes. ...
user avatar
  • 113
1 vote
2 answers
174 views

Parsing a GenBank file with multiple gene entries

Basically a GenBank file consists of gene entries (announced by 'gene') followed by its corresponding 'CDS' entry (only one per gene) like the two shown here below. I would like to extract part of the ...
user avatar
  • 135
3 votes
2 answers
4k views

Existing tool for converting gff3 to genbank (gbk)

I want to convert my gff3 annotation files to genbank format for use in Mauve. I found the seqret tool here https://www.ebi.ac.uk/Tools/sfc/emboss_seqret/ which can perform this task, but my files (...
user avatar
  • 33
5 votes
1 answer
240 views

which NCBI tool is optimized to identify a species from a DNA fragment?

“The Iceman” was a man who lived 5300 years ago and whose body was recovered from an Alpine glacier in 1991. Some fungal material was recovered from his clothing and sequenced. Ice man : found as a ...
user avatar
1 vote
1 answer
57 views

Biopython Genbank writer not splitting long lines

I am parsing a csv file of annotated sequences and using Biopython to generate Genbank files for each. I want to add annotations of the sequence features. My output file shows features listed without ...
user avatar
  • 21
2 votes
1 answer
73 views

Parsing specific features from Genbank by label?

I'm trying to parse a Genbank file to find a specific feature. I can pull it out if I know the feature type (e.g. repeat_region). E.g. if I'm looking for this feature: ...
user avatar
4 votes
2 answers
336 views

Entrez.efetch returns incomplete genbank records

I am using the biopython Entrez.efetch command to retrieve all features (CDS, mRNA, ...) of some genomes. In this case (NC_014649, Acanthamoeba polyphaga mimivirus), it works as expected: ...
user avatar
  • 595
1 vote
0 answers
18 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 ...
user avatar
7 votes
3 answers
1k views

Importing GFF file with Biopython

Is there a way to import a GFF file for an organism with Biopython in the same way you can for a Genbank file? For example, ...
user avatar
7 votes
1 answer
321 views

ANY (technical) reason behind submitting sequences to GenBank versus ENA Sequence

The DNA sequence sections of the three INSDC databases (i.e., DDBJ, ENA Sequence and GenBank) are synchronized periodically and strive to keep their stored data as ubiquitously accessible as possible. ...
user avatar
1 vote
1 answer
78 views

Genomics: Downloading genomic data from the NCBI Genbank database

I'm taking a genomics course at school, and we have a final project due in a few weeks. Now, I can write bash scripts to sort and analyze the data, etc. However, I'm not quite sure how to actually ...
user avatar
7 votes
1 answer
81 views

Source for whole genome comparisons at NCBI Genomes

The NCBI Genomes database has these dendrograms for (presumably) whole genome comparisons for certain species, e.g. Pseudonomas aeruginosa or Escherichia coli. How were these comparisons done? ...
user avatar
2 votes
0 answers
105 views

How to extract Genbank sequence and identify ORFs from the sequence in R?

I used ape::read.GenBank and then used the function ORFik::findORFs to find ORFs in the sequence of interest. But the former ...
user avatar
  • 428
1 vote
1 answer
638 views

Extract certain fields using from GenBank file using Bash script

I need to extract certain lines from a GenBank file on Linux. Input file: ...
user avatar
3 votes
3 answers
113 views

How to get the product of a CDS

I need the name of the protein in /product="protein_name" using bash commands. Beware, there is a lot of whitespace between lines. ...
user avatar
3 votes
1 answer
2k views

How to extract the protein fasta file from a genbank file?

Let's say I have a genbank file, e.g. the one downloaded from here, that contains entries as ...
user avatar
  • 743
-1 votes
1 answer
105 views

How do I get the gene order from genbank file using gene names?

I have a list of gene names which I know are present in a genbank file. I want to order these genes based on the order of their presence in the genbank. This should also account for whether they are ...
user avatar
  • 925