This is a test version of Biostars. For the public version, visit https://www.biostars.org.
making a multi-fasta file out of 9 separate fasta file

Hi

I want to do a python task which involves extracting the sequence from a multi-fasta file using a file of IDs for Mycosphaerella graminicola genome. I do not have data(multi-fasta file) for this task. Therefore, I downloaded from ENA using the following link. https://www.ebi.ac.uk/ena/data/search?query=Mycosphaerella+graminicola. Now I have 9 separate fasta files containing Mycosphaerella chr 3,chr 4,chr 5, chr 6, chr 10, chr 12, chr 13, chr 15, chr 19,chr 21 whole genome short gun sequencing.

My question 1)how to make a single multi-fasta file using these separate 9 fasta files. 2) Is the multi-fasta file available on web straightway?

Any Help would be highly appreciated! Thanks

sequence
• 8,292 views • link
written by rupashree.choudhury • rupashree.choudhury

2 answers

1) how to make a single multi-fasta file using these separate 9 fasta files. 

cat *.fasta > Mgraminicola.fasta

But why don't you use GCA_000219625.1, it seems the assembly is at chromosome level, with 21 chromosomes.

• 0 views • link
written by h.mon h.mon

May be better to do cat chr*.fasta > Mgraminicola.fasta.

• 0 views • link
written by GenoMax GenoMax

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted. Upvote|Bookmark|Accept

• 0 views • link
written by WouterDeCoster WouterDeCoster

May be this also works:

$ awk "1" chr*.fa > out.fa
• 0 views • link
written by cpad0112 cpad0112

Log in to answer this question.

More posts like this