SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Uploading BAM Track to UCSC Genome Browser jjw14 Bioinformatics 4 02-01-2019 07:15 AM
Error message when uploading a file into UCSC genome browser petrosul General 2 02-14-2012 10:08 AM
not able to load a bed file into UCSC genome browser ataaillah Introductions 0 01-18-2012 07:46 AM
Upload Bam file to custom track UCSC Genome Browser gabrielw Bioinformatics 4 06-15-2011 11:26 AM
RNA-Seq: ENCODE whole-genome data in the UCSC genome browser (2011 update). Newsbot! Literature Watch 1 11-24-2010 01:08 PM

Reply
 
Thread Tools
Old 10-14-2011, 03:06 AM   #1
qnc
Member
 
Location: London

Join Date: Jun 2010
Posts: 51
Default 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
qnc is offline   Reply With Quote
Old 10-14-2011, 03:10 AM   #2
qnc
Member
 
Location: London

Join Date: Jun 2010
Posts: 51
Default

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.
qnc is offline   Reply With Quote
Old 10-14-2011, 05:54 AM   #3
mgogol
Senior Member
 
Location: Kansas City

Join Date: Mar 2008
Posts: 197
Default

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.
mgogol is offline   Reply With Quote
Old 10-14-2011, 06:11 AM   #4
qnc
Member
 
Location: London

Join Date: Jun 2010
Posts: 51
Default

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?
qnc is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:48 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO