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
          Exploring the Dynamics of the Tumor Microenvironment
          by seqadmin




          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
          07-08-2024, 03:19 PM
        • seqadmin
          Exploring Human Diversity Through Large-Scale Omics
          by seqadmin


          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
          06-25-2024, 06:43 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 07-16-2024, 05:49 AM
        0 responses
        26 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-15-2024, 06:53 AM
        0 responses
        32 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-10-2024, 07:30 AM
        0 responses
        40 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-03-2024, 09:45 AM
        0 responses
        205 views
        0 likes
        Last Post seqadmin  
        Working...
        X