Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BAM file to Histogram on UCSC Genome Browser

    How do I make a track of a histogram of the coverage of the bam files (Track 4 in the session provided below) I was able to load upto UCSC? Do I need to use BEDTools?

    Here is the session http://goo.gl/7o4Il

    Tracks
    1) Scale: will open on the whole of the long arm of Chr22
    2) AD32 region: Defined region of interest for this investigation
    3) Abhay-Annotations: SNP of interest from BAM file (23,466 recorded)
    4) Chr22 BAM: The bam files aligned
    5) RefSeq Genes: Chr22 genes
    6) Wiki Track: For collaborative working and also shows in purple the SeattleSeq potential mutation changes in areas with a read depth of more than 5

    The rest of the tracks are reference tracks for context and further work.
    7) Omim Genes: Just to see if there are any reported similar diseases in a particular area of interest
    8) Flagged SNP: to eliminate know clinically significant SNPs
    9) Microsatellites: Maybe useful if confirmation work.
    10) STS Markers: Defined anchors again if needed for further investigations
    11) GC%: To help with PCR primer design for further work if needed
    12) CpG Islands: Basic promoter info
    13) Evo CpG: Basic promoters across species
    14) Evolutionary conservation: To see how conserved an area is over species again to help with judging how important a mutant is.

    Any help with this would be appreciated.

    Thanks

  • #2
    Also UCSC GB has highlighted changes in the BAM Files. Is it possible to have these in a track as well. It would be good to have Seattle Seq annotations, UCSC annotations and ANNOVAR annotations as different tracks to compare and see what is in agreement with all? If there are other SNP annotations that people want to recommend that would be good to have them as tracks as well.

    And if we could also figure out how to generate a track of the gaps or once we have the histogram areas of low coverage and gaps that would be helpful.

    I have the BAMTools manual and MacOSX but not sure how to get started. Maybe there is a simple pipe for me to get used to using the program. Have limited unix skills but happy to learn.

    Once again thanks.

    Comment


    • #3
      Code:
      genomeCoverageBed -bg -ibam $file.sorted.bam -g sacCer2.fa.fai > file.bedgraph
      
      wigToBigWig file.bedgraph sacCer2.fa.fai file.bw
      There may be a way to put those two together, come to think of it.

      (bash)

      Code:
      wigToBigWig <(genomeCoverageBed -bg -ibam async.sorted.bam -g sacCer2.fa.fai) sacCer2.fa.fai file.bw
      You could also try IGV, which does the histogram for you.

      Comment


      • #4
        Thanks for this do I do this in Terminal or X11?

        Is that correct if I need to do this on another unix system let me know.

        Do I need to download any software or will this do this automatically... or do I do this from UCSC GB itself?

        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
        34 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-08-2024, 08:03 AM
        0 responses
        72 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-07-2024, 08:13 AM
        0 responses
        82 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