Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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

  • #2
    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

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM
    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
      02-26-2024, 02:07 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-14-2024, 06:13 AM
    0 responses
    32 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-08-2024, 08:03 AM
    0 responses
    71 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-07-2024, 08:13 AM
    0 responses
    80 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-06-2024, 09:51 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X