Modify Fasta headers

Can I modify the following headers

>gi|392981410|ref|NC_018080|pseudocap|273_483_2027_+
>gi|392981410|ref|NC_018080|pseudocap|273_2056_3159_+

so on

>protein 1 [Pseudomonas aeruginosa 19BR]
>protien 2 [Pseudomonas aeruginosa 19BR]

so on?

sequence

• 2,170 views • link

updated 4.0 years ago by Ram

• written 10.7 years ago by utkarsh.sood •

1 answer

awk 'BEGIN{counts=1} { if ($0 ~ /^>/)  { print ">protein "counts" [Pseudomonas aeruginosa 19BR]"; counts++ } else { print }}' <input_file>

• 0 views • link

updated 4.0 years ago by Ram

• written 10.7 years ago by GouthamAtla

More posts like this

select longest isoform from RefSeq proteome fasta
written 7.1 years ago by wd •

Hi I would like to extract the longest isoform from a RefSeq proteome fasta. I looked at suggested solutions on this and other fora but …
How to calculate the fraction of bases that are different from the reference base using samtools mp…
written 7.6 years ago by M.O.L.S •

**Questions:** - How are sites identified to find the fraction of bases that are different from the reference base at each position? It seems from …
Using HTseq-count for Pseudomonas Aeruginosa data
written 8.0 years ago by jspainhour3 •

I am trying to use HTseq to generate count data from a Pseudomonas Aeruginosa rna-seq sample. I end up with all of my output going …
HTSeq-counts can not map any reads to gtf file
written 8.2 years ago by nazaninhoseinkhan

Dear all, I need to get the raw read counts of a *pseudomonas aeruginosa UCBPP-PA14 * RNASeq data set. After mapping with "Bowtie for Illumina", …
StringTie cannot be run as a problem in naming convention
written 8.5 years ago by nazaninhoseinkhan

Hi all, I know this question has been asked previously, however I still cannot solve my problem. I have obtained both gtf and corresponding fasta …
SOAP aligner removes samfile headers
written 8.6 years ago by fiona.newberry

Hi, I am using samtools to convert the sam file outputs from SOAP aligner to bam file. However, the samfiles do not have any sort …
Fetching NCBI gene symbols from NCBI protein ids or GI identifiers
written 8.8 years ago by lakhujanivijay

I have the following information from blastx annotations of bacterial genes predicted by prodigal: Sequence name Sequence desc. Sequence length Hit desc. Hit ACC gene_1_contig_1 …
blast results: grepping and replacing
written 8.8 years ago by kinggang •

Hi, I am stuck with my blast result (outfmt 6). In first column i got the contig id;like second column also i got subject id. …
Populate data frame column X with substring from column Y using R
written 9.0 years ago by ddowlin

Hi all, I have a data frame in R with a column of gene identifiers (taken from fasta headers). These vary and include both Ensembl …
Modify Fasta header
written 10.7 years ago by utkarsh.sood •

Hello I have a fasta file with sequence headers written as ``` >0|quiver|1..2075|- >0|quiver|2210..3058|- >0|quiver|3112..4169|- ``` and so on till around 1000 sequences. I want …