SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Correlation between sequencing depth and false positives Jane M Bioinformatics 8 07-30-2012 05:52 AM
False negative variant calling in haploids. Call variants using coverage (not stats)? Genomics101 Bioinformatics 3 06-07-2012 11:29 AM
Fusions false positives - using Breakdancer dicatoma General 0 01-23-2012 06:29 AM
Fusions false positives - using Breakdancer dicatoma Bioinformatics 0 01-23-2012 04:10 AM
remove false positives from less than 17 mer hpgala Illumina/Solexa 0 02-26-2011 12:30 PM

Reply
 
Thread Tools
Old 12-28-2013, 12:31 PM   #1
arcolombo698
Senior Member
 
Location: Los Angeles

Join Date: Nov 2013
Posts: 142
Default Method for Calling Variants, False Positives

Hello.

I have a question regarding the algorithm our lab uses for variant calling.
Essentially when calling the variants we call the variants for each chromosome individually. But this can cause errors (false positives) because if the reference map actually belongs to a different chromosome,say chromosome i, the algorithm will generate a score for the current chromosome having the variants called, say chromosome x.

However, for computational speed, splitting up the variant calling is more efficient because each node will do the process for each chromosome.

I am just looking for general advice.

The only way I can test this method is to finalize the VCF files, and use VCFtools to remove false positives.
-or I can use GATK, or Atlas2 to do a comparative analysis.
arcolombo698 is offline   Reply With Quote
Old 12-29-2013, 09:21 PM   #2
yueluo
Member
 
Location: Guangzhou China

Join Date: Aug 2013
Posts: 82
Default

Hi, just curious, did you
1) align your reads to individual chromosomes and then call variant, or 2) align reads to the whole genome and then call variant for individual chromosomes ?
yueluo is offline   Reply With Quote
Old 12-30-2013, 10:30 AM   #3
arcolombo698
Senior Member
 
Location: Los Angeles

Join Date: Nov 2013
Posts: 142
Default Calling Variants of aligned reads to the whole genome

Hello.

I aligned the reads to the entire genome, and then called the variants for each individual chromosome.

(#2)
arcolombo698 is offline   Reply With Quote
Old 12-30-2013, 11:22 AM   #4
vivek_
PhD Student
 
Location: Denmark

Join Date: Jul 2012
Posts: 164
Default

Its really not advisable to be splitting the mapping process by chromosomes. You will lose quite a bit of information and run into false positives in the latter stages, better idea to quicken this process is to split the reads into multiple files, align them to the reference and merge in the end.

Variant calling by chromosome is an entirely different issue and totally safe to divide this into multiple jobs by chromosome without worrying about false positives.
vivek_ is offline   Reply With Quote
Old 12-30-2013, 12:29 PM   #5
arcolombo698
Senior Member
 
Location: Los Angeles

Join Date: Nov 2013
Posts: 142
Default

I have a question the alignment process.

for a given sample, I normally perform quality control on the samples and afterward for a specific sample, I run the alignment on the paired trimmed samples (ie. running the alignments on the entire trimmed fastq files).

Are you suggesting that it is not ideal to run the alignment on an entire fastq file? but a partition of the fastq file? Does it matter if the project is DNA ? or RNAseq?

Thank you.
I am just confused because I normally run the alignments on a complete fastq file, and not a partition of a fastq file.

Is there a resource that you'd recommend?
arcolombo698 is offline   Reply With Quote
Old 12-30-2013, 12:45 PM   #6
vivek_
PhD Student
 
Location: Denmark

Join Date: Jul 2012
Posts: 164
Default

Quote:
Originally Posted by arcolombo698 View Post
I have a question the alignment process.

for a given sample, I normally perform quality control on the samples and afterward for a specific sample, I run the alignment on the paired trimmed samples (ie. running the alignments on the entire trimmed fastq files).

Are you suggesting that it is not ideal to run the alignment on an entire fastq file? but a partition of the fastq file? Does it matter if the project is DNA ? or RNAseq?

Thank you.
I am just confused because I normally run the alignments on a complete fastq file, and not a partition of a fastq file.

Is there a resource that you'd recommend?
No, its not an issue if you decide to run the alignment of an entire fastq file against the entire reference sequence.

I was suggesting the preferred way of doing the task in parallel to may be utilize a HPC cluster if you have access to one. You'd rather split the fastq file than split the reference file for parallelism.
vivek_ is offline   Reply With Quote
Old 12-30-2013, 01:03 PM   #7
arcolombo698
Senior Member
 
Location: Los Angeles

Join Date: Nov 2013
Posts: 142
Default

thank you. makes sense.

is it more accurate to run an alignment on a smaller fastq file input into the aligner?
arcolombo698 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 04:50 PM.


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