SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   General (http://seqanswers.com/forums/forumdisplay.php?f=16)
-   -   Chromosome sizes don't match (http://seqanswers.com/forums/showthread.php?t=65661)

sjb343 01-19-2016 07:35 PM

Chromosome sizes don't match
 
Hi,

I'm having a heck of a time finding the chromosome sizes for mm9 that match the mm9 Bowtie pre-built indexes on their website. I have ChIP-seq data aligned using Bowtie2 to mm9 and used Macs2 to call peaks generating bedgraph files. Now I'm trying to convert these to BigWig files using bedgraphtobigwig but the UCSC chromosome sizes are not matching up. I can't find any threads with this problem but I'm not sure what I'm doing wrong. Since I can't seem to find anything, my questions is, is there a way to extract the chromosome sizes from the available Bowtie mm9 pre-built indexes?

Thanks,
sj

dariober 01-20-2016 12:49 AM

You can get chromosome sizes from UCSC with something like this (credit: bedtools docs):

Code:

mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -e \
        "select chrom, size from mm9.chromInfo"  > mm9.genome

Then you could check if these sizes match the sizes used by bowtie2 to build the index (they should). I think the utility bowtie2-inspect with -s option should give you the chromosome lengths.

Finally, you could extract chrom sizes from the header of a bam file with e.g:

Code:

samtools view -H aln.bam \
    | grep -P "@SQ\tSN:" \
    | sed 's/@SQ\tSN://' \
    | sed 's/\tLN:/\t/' > genome.txt



All times are GMT -8. The time now is 09:11 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.