Thanks!!!! :)
• 0 views
•
link
I am trying to transform Multifasta - with many sequences fasta to a csv table. But I got the following error
is there any solution to convert a fasta file to csv?
This should get you going:
from Bio import SeqIO
for re in SeqIO.parse('CHIKV1-X-gb_AB455493.fasta', 'fasta'):
print('>{}\t{}'.format(str(re.description).replace('|', '\t'), re.seq))
Prints out:
>gb:KX262887 Organism:Zika virus Strain Name:103451 Segment:null Subtype:Asian Host:Human GTTGTTGATCTGTGTGAATCA ....
>gb:KX262887 Organism:Zika virus Strain Name:103451 Segment:null Subtype:Asian Host:Human GTTGTTGATCTGTGTGAATCA ....
>gb:KX262887 Organism:Zika virus Strain Name:103451 Segment:null Subtype:Asian Host:Human GTTGTTGATCTGTGTGAATCA ....
The rest is simply a matter of possibly using comma instead of tab.
Thanks!!!! :)
Log in to answer this question.
What on earth is a multi-folder file? Given that your header is in a custom format, you're going to have to parse it using custom code.
Multifasta - with many sequences fasta. I'm going to edit sorry. Is there any solution?
You're halfway there.
re.id + redescriptionshould give you the whole header, make sure of that. Once you have it, you can split by|and then spit each element of that by:into key value pairs. Then, these key value pairs along withre.seqwill give you the columns for eachre. Write these attributes separated by,and you'll have your CSV.Thanks!!!