This seems that should be done for each chromosome individually, isn't it ?!
I will try to do so, thanks a lot ^_^
Hi all :)
Where can I find a credible and reliable source (website) for downloading genomes on chromosome level only ?!
I don't want contigs or scaffolds, or even cDNAs !!
Thanks in advance
2 answers
Download any genome you want from any source, then subset that file based on the chromosome. Most are distributed as FASTA sequences. You could select among the chromosomes easily (and within a loop!) by using Biostrings in R, amongst a ton of other programming options.
library(Biostrings)
a <- readDNAStringSet("your_genome.fa", format="fasta")
#show the data structure
a
#the chromosome names are stored in
names(a)
#and you could write to a new file with:
b <- a[whatever chromosome you want]
writeXStringSet(b, file="your_chromosome.fa", format="fasta")
You got it. Wrap it in a loop across indices/chromosome names.
See the answer to this question:
@Erik Wright ...
I took a look at the answers in that thread, and yours seems not applicable to me since I'm not familiar with R language.
Is there any other options to download the full genomes in chromosomes level ?!
I would be very grateful if you can help me out in this .
Log in to answer this question.
Hi all :) I'm really sorry for what I'm about to do, but I have to cross-reference the same question on some other specialized forums since I have waited so long for the appropriate answer and no one seems willing to help , unfortunately :/ Greetings !!
Both my answer and Erik's answer are valid solutions for doing what you are trying to do. My guess is nobody suggested other options because ours are sufficient. You asked on a bioinformatics group and were given bioinformatics-based solutions. If you wanted an application that would download these for you, you should have explicitly said so.
Thanks for the clarification, I truly appreciate your help and concern but I wrote earlier that I'm not familiar with R language and I don't want applications either ... just some kinda straightforward solution with the minimal confusion, I'm still a beginner in this field and need guidance !!