Hello,
I want to make a script that will extract a sequence from a multifasta file. I have made a file with 19 genomes from NCBI and I want to be able to extract a genome that I need from it. I have played a little with Bio::SeqIO and learned how to create a object that has the hole file in it and then right the whole information in another file.
What i do not know how to do is how to select only the fasta sequence from the genome that I want and write it onto a file.
I think I have to create a third object that would contain only that genome, but I still don't know how to read that part from my sequence object...
Thanks in advance.
I want to make a script that will extract a sequence from a multifasta file. I have made a file with 19 genomes from NCBI and I want to be able to extract a genome that I need from it. I have played a little with Bio::SeqIO and learned how to create a object that has the hole file in it and then right the whole information in another file.
What i do not know how to do is how to select only the fasta sequence from the genome that I want and write it onto a file.
I think I have to create a third object that would contain only that genome, but I still don't know how to read that part from my sequence object...
Thanks in advance.
Comment