SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
extra high coverage regions in exome sequencing data lyz1030 Bioinformatics 1 10-22-2012 03:26 AM
Tool for analyzing low-coverage human data pravee1216 Bioinformatics 0 05-17-2012 08:11 PM
How to find coverage for RNA-seq, solid data? pinki999 Bioinformatics 2 09-28-2011 05:59 AM
Maq Ė how to filter low coverage areas? icg Bioinformatics 0 01-06-2011 01:49 AM
low 454 coverage combined with high solexa coverage strob Bioinformatics 7 10-07-2010 11:14 AM

Reply
 
Thread Tools
Old 03-29-2013, 09:29 AM   #1
willMD
Junior Member
 
Location: Maryland

Join Date: Mar 2013
Posts: 5
Default find areas of very low and extra high coverage in NGS data

Are there any tools that can report all areas of very low coverage or areas of high coverage in NGS data?
willMD is offline   Reply With Quote
Old 03-29-2013, 11:26 AM   #2
dkatzel
Junior Member
 
Location: Rockville, MD USA

Join Date: Mar 2013
Posts: 5
Default

I've written some custom code to programatically do this, but it only works on a few assembly formats (don't have time or funding to add more formats yet). But several assembler programs have ways to output the coverage levels of the contigs. It would then be up to you to parse that output and find any locations that are below or above your thresholds.
dkatzel is offline   Reply With Quote
Old 03-29-2013, 12:41 PM   #3
suryasaha
Member
 
Location: New York

Join Date: Mar 2011
Posts: 27
Default

A commonly used method is to map all your reads to contigs using bowtie2/bwa mem/bwa and calculate the coverage with samtools or bedtools.
suryasaha is offline   Reply With Quote
Old 04-01-2013, 12:39 AM   #4
jimmybee
Senior Member
 
Location: Adelaide, Australia

Join Date: Sep 2010
Posts: 119
Default

Quote:
Originally Posted by willMD View Post
Are there any tools that can report all areas of very low coverage or areas of high coverage in NGS data?
Can you give us some more information? What type of data do you have? Have you mapped reads to a reference yet?

If you do have a bam file of aligned reads to a reference sequence then you will be able to generate a plot of the per-base coverage of the reference sequence using bedtools genomeCoveragebed.
jimmybee is offline   Reply With Quote
Old 04-04-2013, 11:33 AM   #5
willMD
Junior Member
 
Location: Maryland

Join Date: Mar 2013
Posts: 5
Default

Human whole genome data aligned to hg19 in bam file. Thanks. I'll look at bedtools.
willMD 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 01:21 PM.


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