This is a test version of Biostars. For the public version, visit https://www.biostars.org.
UCSC Fasta files by chromosome for rn6

Hi,

Ive noticed UCSC downloads page has fasta files by chromosome for hg19 and mm9. But that isn't available for rn6. Is there some place from where i can download fasta files by chromosome for rn6 genome assembly ?

Apoorva

ucsc fasta rn6

3 answers

The UCSC utility twoBitToFa can get the sequence for an individual chromosome at a time via the -seq option, which you can run in a loop for each chromosome:

$ wget  http://hgdownload.cse.ucsc.edu/goldenPath/rn6/bigZips/rn6.chrom.sizes
$ for chr in $(cut -f1 rn6.chrom.sizes); do twoBitToFa -seq="${chr}" http://hgdownload.soe.ucsc.edu/goldenPath/rn6/bigZips/rn6.2bit $chr.fa; done
$ # verify with the faSize util:
$ head -1 rn6.chrom.sizes
chr1    282763074
$ faSize chr1.fa
282763074 bases (14711797 N's 268051277 real 153502700 upper 114548577 lower) in 1 sequences in 1 files
%40.51 masked total, %42.73 masked real

You can download twoBitToFa and faSize from the following the directory appropriate for your operating system here: http://hgdownload.soe.ucsc.edu/admin/exe/

If you have further questions about UCSC data or tools feel free to send your question to one of the below mailing lists:

  • General questions: genome@soe.ucsc.edu
  • Questions involving private data: genome-www@soe.ucsc.edu
  • Questions involving mirror sites: genome-mirror@ose.ucsc.edu

ChrisL from the UCSC Genome Browser

@ChrisL: Do you know why *.chromFa.tar.gz (assembly sequence in one file per chromosome.) are not made available for all genomes at UCSC? They are for human genome.

We make these individual chromosome files when the chromosome count is under 100. Most assemblies these days have over 100 chromosomes and so we don't run them up anymore.

FYI, NCBI Rno6 chromosome wise sequences are available at: ftp://ftp.ncbi.nlm.nih.gov/genomes/R_norvegicus/ and ftp://ftp.ncbi.nlm.nih.gov/genomes/R_norvegicus/Assembled_chromosomes/seq/

Thank you very much

Apoorva

Log in to answer this question.