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

Hi I have sequences like this

> abc.ff
ccgtac

> qwe.oo
aattcg

> dee.ff
aatt

> rte.oo
ttcc

and I want to be like this

> ff 
ccgtacaatt

> oo
aattcgttcc

can you tell me how to do it thsanks

genome sequence alignment r
• 1,150 views • link
updated by GenoMax GenoMax • written by a.moner • a.moner

In R: mutate a ID column from name column, group_by ID column and paste sequence column.

• 0 views • link
written by shenwei356 shenwei356

In shell: Step 1) split sequences by FASTA header using seqkit

seqkit split --id-regexp "\.(.+?)$" -i seqs.fa

Step 2) combining sequences in a file to one sequence. for one splited file:

$ echo ">ff" > combined.fa
$ cat seqs.id_ff.fa | seqkit fx2tab | cut -f 2 |   paste -s -d "" >> combined.fa
• 0 views • link
written by shenwei356 shenwei356

Hello a.moner!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.

Cheers!

PS: Not really appropriate to delete old question without following up and just repeating the question in a new thread.
• 0 views • link
written by WouterDeCoster WouterDeCoster

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Want to reverse the fasta sequence
    written by harry • harry
    >ENST00000014935 TAGGTTTTTAATTACCTTTATTTTATTTTGCCAAACATACCTGGGAATACCTTTTATTTTTTTTTTACCTTGGGGTGATGGTTCCAAACCATAAATGTGATTATAGTTAACACATGACCCTTCTAGCGTCCCAGCCAGTGTTTTTCCTGACCTCTCTTCTTTGGAGAGGAGGATGGAAGGGAGGGGTCCGGCATGCTGCTGGCATTTTGCTGTGTCCTGCAGCCCCTTTCCGGGACACCTGGGTTCACACAGCTTTTTAGCTTACATAACTGGTGCAGATTTTCTGTGTGGAGATGTTGCCTTGACCAGCCTTGGCTGGACTTTACCAGGCATGCAGAAGCCTGTACCAACACAGACTACAGCACCCAGGAGGTGCGAGTGTGGCTGCTCAGCGGTTATAACAGGCCTGACTGCATTGTTCACCGGATTATAATGAGCCAAAATGTTTCCCGGTGTTTGCTGGTTTCAGGGAAGGAGTTTGATATAGCAGATTAACCACCCTCCTTGTAGCTATTGGGGCTTAATGGTTTCCTGGTGATTCTTACCAATCCACAATAAACATGGCCCATTGGCATATCTGCTGCACAAGTGTCCTATCTCACCAATCTGGGTTTTTGTTCTCAGTAACTTTCCTTCTTGTCATACAACATCTTCATTCCTCTTTCTGAACCCTCCCTTCCCCTACCCCAACCCAGAGCCCACTTTGTCTCCACTCCTGATACTACACTACCTGGCAGGTGGCATGAGTGCAGGGCCCCTGGCTTCCTCTCCTAATCTAGGCACAAGCCCAACCAAAGAACAAGAGCCAAATCAAACAAGGCAGGCAGGGGTGGACTACAGTCACAGGGCAACTATAGTTGAAGCCCCCCAGCCCCAGGGCTGGATGGACGGGGGAGGCTGGGGTTTAAGTCCCAAAAGGCAGCAGGCCCTGGGGGGGTAGGGGGACGCTCAGGCAGCAGGGCACAGCTGAGGGGACAGGAGTGATAGCAGCAACAGAACAGTGAGGCTGAGAGGCTGGACGCTGTGCGCCTGGCTCAGCTTCAGCTCCACCTCCACGGGGTAGTGGTCACTGATGTTGAGGGC >ENST00000014944 TCCTCCTCGGTGAGCTGGAAGCTCGTGGGGAAGTCAAAGGCAGCCGCAGTGTGCAGCAGACTCCGGCAGCGCTCCCCGTGCAGCACGACGCGGTCATAGGTGCAGTGGGTGCTGGCCCGCACTGTGGTGTCCTCCCCATCGGCAATCACCCAGTGGAAGCCTGGCTCAGTCCGCAGCTCCAGCTTGTCCAGGCGCTTTTTGGTCAGTGAAGCGCAGTCAGCATTGAAGTCCCCAAGCAGGATCACGTC Likewise, I have a huge fasta sequence and I want to reverse the sequence of those fasta sequences like below:- >ENST00000014935 …
  • Script for making exon file
    written by harry • harry
    Hi everyone, I have an exon.fasta file which looks like --- >ENSE00002234944 TTAACTTGCCGTCAGCCTTTTCTTTGACCTCTTCTTTCTGTTCATGTGTATTTGCTGTCTCTTAGCCCAGACTTCCCGTGTCCTTTCCACCGGGCCTTTGAGAGGTCACAGGGTCTTGATGCTGTGGTCTTCATCTGCAGGTGTCTGACTTCCAGCAACTGC TGGCCTGTGCCAGGGTGCAAGCTGAGCACTGGAGTGGAGTTTTCCTGTGGAGAGGAGCCATGCCTAGAGTGGGATGGGCCATTGTTCATCTTCTGGCCCCTGTTGTCTGCATGTAACTTAATACCACAACCAGGCATAGGGGAAAGATTGGAGGAAAGATGA GTGAGAGCATCAACTTCTCTCACAACCTAGGCCA >ENSE00001746346 CAGGGCCATCAGGCACCAAAGGGATTCTGCCAGCATAGTGCTCCTGGACCAGTGATACACCCGGCACCCTGTCCTGGACACGCTGTTGGCCTGGATCTGAGCCCTGGTGGAGGTCAAAGCCACCTTTGGTTCTGCCATTGCTGCTGTGTGGAA >ENSE00002312635 CAGGGCCATCAGGCACCAAAGGGATTCTGCCAGCATAGTGCTCCTGGACCAGTGATACACCCGGCACCCTGTCCTGGACACGCTGTTGGCCTGGATCTGAGCCCTGGTGGAGGTCAAAGCCACCTTTGGTTCTGCCATTGCTGCTGTGTGGAAGTTCACTCC TGCCTTTTCCTTTCCCTAGAGCCTCCACCACCCCGAGATCACATTTCTCACTGCCTTTTGTCTGCCCAGTTTCACCAGAAGTAGGCCTCTTCCTGACAGGCAGCTGCACCACTGCCTGGCGCTGTGCCCTTCCTTTGCTCTGCCCGCTGGAGACGGTGTTTG TCATGGGCCTGGTCTGCAGGGATCCTGCTACAAAGGTGAAACCCAGGAGAGTGTGGAGTCCAGAGTGTTGCCAGGACCCAGGCACAGGCATTAGTGCCCGTTGGAGAAAACAGGGGAATCCCGAAGAAATGGTGGGTCCTGGCCATCCGTGAGATCTTCCCA GGGCAGCTCCCCTCTGTGGAATCCAATCTGTCTTCCATCCTGCGTGGCCGAGGGCCAGGCTTCTCACTGGGCCTCTGCAGGAGGCTGCCATTTGTCCTGCCCACCTTCTTAGAAGCGAGACGGAGCAGACCCATCTGCTACTGCCCTTTCTATAATAACTAA AGTTAGCTGCCCTGGACTATTCACCCCCTAGTCTCAATTTAAGAAGATCCCCATGGCCACAGGGCCCCTGCCTGGGGGCTTGTCACCTCCCCCACCTTCTTCCTGAGTCATTCCTGCAGCCTTGCTCCCTAACCTGCCCCACAGCCTTGCCTGGATTTCTAT CTCCCTGGCTTGGTGCCAGTTCCTCCAAGTCGATGGCACCTCCCTCCCTCTCAACCACTTGAGCAAACTCCAAGACATCTTCTACCCCAACACCAGCAATTGTGCCAAGGGCCATTAGGCTCTCAGCATGACTATTTTTAGAGACCCCGTGTCTGTCACTGA AACCTTTTTTGTGGGAGACTATTCCTCCCATCTGCAACAGCTGCCCCTGCTGACTGCCCTTCTCTCCTCCCTCTCATCCCAGAGAAACAGGTCAGCTGGGAGCTTCTGCCCCCACTGCCTAGGGACCAACAGGGGCAGGAGGCAGTCACTGACCCCGAGACG …
  • Phages sequencing anaysis
    written by anasofiamoreira94 anasofiamoreira94
    Hello, Im sorry for this question, but i have trimmed fastq sequences of my phages. And now I would like to know the similarity between …
  • Number of genomes sequences in NCBI
    written by misterie misterie
    Hi, Could you tell me how can I check number of organisms for which full genome sequences are available within NCBI? Also I would like …
  • Number of genomes sequences in NCBI
    written by misterie misterie
    Hi, Could you tell me how can I check number of organisms for which full genome sequences are available within NCBI? Also I would like …
  • How To Turn A Two Column Data File Containing Pairs Into A series of clusters
    written by stefano.campanaro • stefano.campanaro
    Hi All, I intersected to data sets and I have a tabulated file with two columns like: AA BB BB CC CC DD BB AA …
  • Remove ambiguous bases from alignment sequence
    written by a.moner • a.moner
    Hi, could you please help me with this I have very big alignment file and there are many ambiguous bases, 1- I want to replace …
  • parse data in excel or using perl
    written by sharmatina189059 • sharmatina189059
    Hi can anyone tell me how to parse a file which looks like this. file1-- 1 2 3 4 5 file 2 1 1234 1 …
  • combine multiple fasta sequences
    written by a.moner • a.moner
    Hi I have sequences like this >abc.ff ccgtac >qwe.oo aattcg >dee.ff aatt >rte.oo ttcc and I want to be like this >ff ccgtacaatt >oo aattcgttcc …
  • How Could I Find The Software For "Divide And Conquer" Approach For Sequence Assembly?
    written by Farid De La Ossa • Farid De La Ossa
    <p>Hi everybody! Could any of you tell me how could I find software in internet like "Divide and Conquer" in order to put together oligonucleotide …