Hello,
Assume you have a Fasta file from Genbank with one or more entries. Each header shows the accession number first, and the description second.
Say I want to annotate/extract/store a sub sequence of one or more entries and not lose the association with the main sequence. Is there a standard way of doing this?
My first intuition is to create a fasta file that has headers as:
>{Accession} {Description} {subsequence description}
Is there a better way?
Thanks!
1 answer
No there is no standard way other than the requirement that fasta header line start with a >. You can store anything you want on the header line. Be aware that many programs will drop things after first white space, so if you want the header to be carried forward in entirety then make sure you keep the things connected as a single string.
> {Accession}_{Description}_{subsequence description}
Log in to answer this question.
-
convert list to fasta format
written by pinheirofabiano •
-
How to change headers of fasta files with the help of another file (nucl_wgs.accession2taxid.gz) co…
written by Sebastian •
-
Batch Annotation of Bacterial Gene Clusters
written by kayrouz.1 •
-
What is unique identifier for organism in NCBI gbff file?
written by Naren
-
Dowload all complete genomes grom GenBank .gbk full format
written by mmart12 •
-
Extracting subsequence from FASTA using Samtools faidx not working
written by hmg •
-
Mauve Headers / Add DNA sequence to genbank file
written by Cricket •
-
Best annotation GTFs for transcript detection vs novel transcript discovery
written by trakhtenberg
-
Refseq Genbank To Fasta Format Failing With Contig Fields
written by Käpt'N Silico •
-
How Can I Save Bioperl Sequence Nested Features In Genbank Or Embl Format?
written by Ryan Thompson
