How to change sequence name in along fasta file?

Login or Register to Ask a Question and Join Our Community

How to change sequence name in along fasta file?

Tags

beginners, fasta, unix commands

Top Forums UNIX for Dummies Questions & Answers How to change sequence name in along fasta file?

06-27-2012

Registered User

4, 0

Join Date: Jun 2012

Last Activity: 7 August 2012, 12:34 PM EDT

Posts: 4

Thanks Given: 0

Thanked 0 Times in 0 Posts

How to change sequence name in along fasta file?

Hi
I have an alignment file (.fasta) with ~80 sequences. They look like this-

>JV101.contig00066(+):25302-42404|sequence_index=0|block_index=4|species=JV101|JV101_4_0

GAGGTTAATTATCGATAACGTTTAATTAAAGTGTTTAGGTGTCATAATTT
TAAATGACGATTTCTCATTACCATACACCTAAATTATCATCAATCTGAAT
TCAGATGTTTATTATAAAAATTAGATGAAAAATATGTTAATATACAAGTA

>JV501.contig00066(+):24356-42404|sequence_index=0|block_index=4|species=JV501|JV501_4_0
AATGACGATTTAGATGAAAAATAT...

The name of the sequences are too big and I want to just keep JV101, JV501 and delete rest of the words after the dot. I am new to unix, please suggest an easy unix command to do this.

Thanks

Baika

baika

View Public Profile for baika

Find all posts by baika

06-27-2012

Registered User

3,733, 1,154

Join Date: Apr 2009

Last Activity: 3 August 2016, 11:03 AM EDT

Posts: 3,733

Thanks Given: 7

Thanked 1,154 Times in 1,124 Posts

Try:

Code:

perl -pe 's/^>[^.]+\K.*//' file

bartus11

View Public Profile for bartus11

Find all posts by bartus11

06-27-2012

Registered User

4, 0

Join Date: Jun 2012

Last Activity: 7 August 2012, 12:34 PM EDT

Posts: 4

Thanks Given: 0

Thanked 0 Times in 0 Posts

Thanks

Thank you very much for such a prompt reply. It is working

!

baika

View Public Profile for baika

Find all posts by baika

Login or Register to Ask a Question

Previous Thread | Next Thread

10 More Discussions You Might Find Interesting

1. UNIX for Beginners Questions & Answers

How to find a specific sequence pattern in a fasta file?

I have to mine the following sequence pattern from a large fasta file namely gene.fasta (contains multiple fasta sequences) along with the flanking sequences of 5 bases at starting position and ending position, AAGCZ-N16-AAGCZ Z represents A, C or G (Except T) N16 represents any of the four...

2. UNIX for Dummies Questions & Answers

Round up -FASTA file

I have the following script: awk 'FNR==NR{s+=$3;next;} { print $1 , $2, 100*$3/s }' and the following file: >P39PT-1224 Freq 900 cccctacgacggcattggtaatggctcagctgctccggatcccgcaagccatcttggatatgagggttcgtcggcctcttcagccaagg-cccccagcagaacatccagctgatcg >P39PT-784 Freq 2...

3. Shell Programming and Scripting

Count and search by sequence in multiple fasta file

Hello, I have 10 fasta files with sequenced reads information with read sizes from 15 - 35 . I have combined the reads and collapsed in to unique reads and filtered for sizes 18 - 26 bp long unique reads. Now i wanted to count each unique read appearance in all the fasta files and make a table...

4. Shell Programming and Scripting

Extract sequence from fasta file

Hi, I want to match the sequence id (sub-string of line starting with '>' and extract the information upto next '>' line ). Please help . input > fefrwefrwef X900 AGAGGGAATTGG AGGGGCCTGGAG GGTTCTCTTC > fefrwefrwef X932 AGAGGGAATTGG AGGAGGTGGAG GGTTCTCTTC > fefrwefrwef X937...

5. UNIX for Dummies Questions & Answers

Change sequence names in fasta file

I have fasta files with multiple sequences in each. I need to change the sequence name headers from: >accD:_59176-60699 ATGGAAAAGTGGAGGATTTATTCGTTTCAGAAGGAGTTCGAACGCA >atpA_(reverse_strand):_showing_revcomp_of_10525-12048 ATGGTAACCATTCAAGCCGACGAAATTAGTAATCTTATCCGGGAAC...

6. Shell Programming and Scripting

Parsing a fasta sequence with start and end coordinates

Hi.. I have a seperate chromosome sequences and i wanted to parse some regions of chromosome based on start site and end site.. how can i achieve this? For Example Chr 1 is in following format I need regions from 2 - 10 should give me AATTCCAAA and in a similar way 15- 25 should give...

7. Shell Programming and Scripting

How to change strt-up sequence of services on SLES10/11

Hi, Can anyone tell me what the following 2 lines are doing base=${0##*/} link=${base#*} I found this in a start up service script and I think it is giving the service link names which in turn will change the start up sequence of services.

8. Shell Programming and Scripting

How to change strt-up sequence of services on SLES10/11

Hi All, I want to change the start-up sequence of services on SLES10/11. I have my own start-up scripts for some services and I want them to start in a particular order(not in alphabetical order) Can anyone help me on this issue?

9. Red Hat

change the shutdown sequence in red hat linux

Hi, I use red hat linux kernel 2.6 I want to add the application shutdown in shutdown sequence . I add the K script in /etc/rc.d/ all sub directory for all running level . But the auto shutdown application is not appear when I type "shutdown -r now" .. There is no indication the application...

10. Solaris

how to change the boot sequence?

Hi all, I have solaris and xp installed... Usually solaris occupies the first boot slot. So i wanted to know if it is possible to change the boot sequence to xp first and then solaris?

Login or Register to Ask a Question