SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
SNP calling on low coverage NGS data JdeBruin Bioinformatics 4 04-22-2015 01:01 AM
find areas of very low and extra high coverage in NGS data willMD Bioinformatics 4 04-04-2013 10:33 AM
Tool for analyzing low-coverage human data pravee1216 Bioinformatics 0 05-17-2012 07:11 PM
low 454 coverage combined with high solexa coverage strob Bioinformatics 7 10-07-2010 10:14 AM

Reply
 
Thread Tools
Old 01-02-2018, 10:44 PM   #1
raman91
Junior Member
 
Location: Singapore

Join Date: May 2016
Posts: 6
Default No coverage or very low coverage in the Complete Genomics data

Hi,

Has anyone worked on Complete Genomics data? Most of the exons in the data have no coverage or very low coverage (5-10 reads only) even though the reads are sequenced at 40X coverage. Can anyone explain why is it so?

Thanks

Regards
raman91 is offline   Reply With Quote
Old 01-03-2018, 02:56 AM   #2
Gorgon_VZ
Member
 
Location: Germany

Join Date: Oct 2016
Posts: 10
Default

Hi! Is it really the exons only, that show reduced coverage and the rest of the genome shows 40x? Have you checked the mean of all exons or for one particular gene only? What kind of genome and data are you working with? In general I would say it's kind of normal to see differences in coverage over the genome, but I am not aware of a exon specific bias. However, if you look on a particular gene or even a gene-family, the coverage might be reduced because of homology regions and multimapped reads with low mapping quality that are filtered out. Maybe it is worth to have a look on the percentage of reads you are able to map against your reference. A high number of discarded sequences could be a hind of such an effect.
Gorgon_VZ is offline   Reply With Quote
Old 01-09-2018, 12:20 AM   #3
raman91
Junior Member
 
Location: Singapore

Join Date: May 2016
Posts: 6
Default

Iím sorry, I didnít explain the issue properly. Iím working on CG WGS data with avg read depth of 40x. The QC metrics looks good for all parameters and the alignment rate is 97.43%. However, I had used cgaTools to convert tsv files provided by CG to BAM. When I visualize these BAM files on IGV, I see minimal coverage at all exons and major parts of introns for all genes (See example image attached). The trend is normally scant coverage hills at junction of intron and exon. I believe something went wrong in my conversion step. Can anyone please suggest a solution for this?
Or is there any other visualization tool specific for CG data that I should be using?Ē

The IGV snapshots is attached herewith.

Thanks in advance
Attached Files
File Type: ppt IGV_IMAGE.ppt (243.0 KB, 5 views)

Last edited by raman91; 01-09-2018 at 12:48 AM.
raman91 is offline   Reply With Quote
Old 01-09-2018, 09:16 AM   #4
Gorgon_VZ
Member
 
Location: Germany

Join Date: Oct 2016
Posts: 10
Default

Hi! It could be an issue due to differing genome builds. Ich guess igv is using hg19 and reads are mapped to hg38.
Gorgon_VZ is offline   Reply With Quote
Old 01-09-2018, 09:27 AM   #5
Gorgon_VZ
Member
 
Location: Germany

Join Date: Oct 2016
Posts: 10
Default

You should be able to Check this by zoom in to nucleotide Level in igv. If the reads Do Not Match the reference this would be a Hind. I am Not a 100% Sure but i believe igv only uses the coordinates and cigar of the bam and does Not care about matching nucleotides. So maybe it is just a slippage of coordinates between hg19 and 38.
Gorgon_VZ is offline   Reply With Quote
Old 01-10-2018, 10:42 PM   #6
raman91
Junior Member
 
Location: Singapore

Join Date: May 2016
Posts: 6
Default

Thanks for your reply. I am pretty sure that the reads are mapped to hg19 build only. I checked the IGV browser too and the read sequence match to the reference sequence.
raman91 is offline   Reply With Quote
Reply

Tags
complete genomics, coverage, sequence alignment

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 04:31 AM.


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