Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
To find the difference in base pair position of two fasta file rajeshkmaurya08 Bioinformatics 1 02-09-2017 09:33 PM
A strange size difference of fastq file qinhao13 Bioinformatics 1 01-06-2015 02:56 AM
Is it possible to do DE analysis on single tumor-normal pair with EdgeR? ymc RNA Sequencing 3 07-28-2013 01:52 AM
Normal size of .sam output file from Tophat GiladZil RNA Sequencing 0 07-04-2011 09:02 AM

Thread Tools
Old 03-08-2018, 02:38 PM   #1
Junior Member
Location: US

Join Date: Mar 2018
Posts: 1
Cool WXS: High amount of duplications! Enormous file size difference in normal-tumor pair!

I wanna preface this by saying I'm relatively new to NGS analysis.

I recently received raw data from WXS (paired end 100bp reads with 100X coverage (12Gb data) with Agilent SureSelect All Human Exon V5 kit). I noticed something's really off from the getgo.

The file size between the normal and tumor pair are enormous. Read 1 and 2 fastq of the normal sample are about 7GB each, and the reads 1 and 2 from the tumor sample is about 40GB each. I've worked with exome data before, and they were usually close in size.

Anyway, I assume everything is okay, and put the raw data through the "pipeline." You know the usual, alignment, sorting, duplicate marking, indel align, mate-info fixing, Base-Recalibration, etc. Almost 45% of the reads failed the "DuplicateReadFilter"(GATK-baserecalibrator
), and another 45% failed the mappingQualityzero filter(GATK-baserecalibrator
)! I have to filter out >90% of my sequence reads!

In my previous runs, I've filtered out at most ~10%.

This made me run the FastQC on the reads, I should have done in the beginning. In the sequence duplication level section, FastQC reports that 25% of the seqs will remain if deduplicated! I see double peaks in the "per sequence GC content"!

All this is new to me. I'm used to seeing yellow bars in the sequence quality section, I don't see it here for some reason.

Would someone show me the ropes?
FallingBabyEagle is offline   Reply With Quote

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 03:38 AM.

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