SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Genome coverage bioinf newbie Bioinformatics 5 01-30-2013 03:34 AM
PubMed: What can you do with 0.1x genome coverage? A case study based on a genome sur Newsbot! Literature Watch 1 04-11-2012 01:18 PM
How to obtain coverage % of the genome? liu_xt005 Bioinformatics 12 10-06-2011 08:14 AM
"nucleotide coverage" to genome feature coverage sheremey Bioinformatics 3 11-02-2010 11:24 AM
genome (%) coverage - what people usually get? dukevn Bioinformatics 15 08-15-2010 09:51 AM

Reply
 
Thread Tools
Old 12-10-2009, 06:06 PM   #1
jjk
Member
 
Location: Brisbane

Join Date: Nov 2009
Posts: 21
Default Genome coverage

I used bowtie to create sam files that are loaded in a tablet program for chloroplast assembly by tags. I guess there is an overall coverage of about 80%. But is there a way to let the program calculate the precise coverage by the tags in comparison to the chloroplast sequence?

Thx
jjk is offline   Reply With Quote
Old 12-10-2009, 11:18 PM   #2
dawe
Senior Member
 
Location: 45°30'25.22"N / 9°15'53.00"E

Join Date: Apr 2009
Posts: 258
Default

You may check BEDTools. There's a coverageBed utility. Of course, you must convert sam to bed first!
d
dawe is offline   Reply With Quote
Old 01-12-2011, 03:28 AM   #3
andrehorta
Member
 
Location: Brazil - Belo Horizonte - UFMG

Join Date: Jan 2011
Posts: 14
Default

Hi!

Are you using Linux or Windows. I can help you!

André
andrehorta is offline   Reply With Quote
Old 04-13-2011, 01:12 AM   #4
zhlu1986
Junior Member
 
Location: Beijing

Join Date: Aug 2010
Posts: 4
Default how to use coverageBed

How to use coverageBed to caculate the coverage of reads to reference?
I have changed sam to bam and then to bed.
I've got only one sam files after aln reads to reference.
Why does the coveragebed programme need to bed files a and b like this:
coverageBed [OPTIONS] -a <bed/gff/vcf> -b <bed/gff/vcf>
If this method does not work.
Is there any other software or script to caculate assembled contigs' coverage to reference?

ps.
the assemble programme is not the popular ones
zhlu1986 is offline   Reply With Quote
Old 04-13-2011, 04:56 AM   #5
chariko
Member
 
Location: Spain

Join Date: Jun 2010
Posts: 56
Default

Maybe this thread may help you.

http://seqanswers.com/forums/showthread.php?t=7679
chariko is offline   Reply With Quote
Old 04-13-2011, 04:10 PM   #6
zhlu1986
Junior Member
 
Location: Beijing

Join Date: Aug 2010
Posts: 4
Default

Quote:
Originally Posted by chariko View Post
Thank you .chariko
But I want this kind of results:

I want to caculate the rate, for example:
the length of reference is 10000bp
the length of reads mapped reagion is 2000bp in all
and the coverage rate is 20%
zhlu1986 is offline   Reply With Quote
Old 04-13-2011, 04:32 PM   #7
Michael.James.Clark
Senior Member
 
Location: Palo Alto

Join Date: Apr 2009
Posts: 213
Default

Read the documentation for BEDtools, it gives an example of how to do this exactly. Very nice PDF on their website!

It doesn't require the first input to be a bed file, either. It can be a BAM file. See the -abam option.
__________________
Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
Projects: U87MG whole genome sequence [Website] [Paper]
Michael.James.Clark is offline   Reply With Quote
Reply

Tags
bowtie, coverage, genome, sam

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 06:02 PM.


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