SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
Strand bias in variant calls wangzkai Bioinformatics 0 06-27-2016 04:19 PM
Manual Quality Control for Variant Calls onconaut General 1 01-22-2015 08:51 AM
Variant calls with a low fraction of alt reads Jeremy37 Bioinformatics 9 04-17-2012 06:18 PM
convert CASAVA variant calls to VCF? krish Bioinformatics 0 12-01-2011 08:44 PM
Complete Genomics Variant Calls quicksand21 Bioinformatics 2 10-11-2011 07:21 AM

Reply
 
Thread Tools
Old 02-27-2017, 02:14 PM   #1
larangela
Junior Member
 
Location: Orange County, CA

Join Date: Feb 2017
Posts: 3
Default Help with Repeat Variant Calls

Hello!

We conducted library prep and sequencing with Illumina. During data analysis, bioinformatics claims that the variants for all 40 samples were exactly the same. The %PF was >80, % aligned was about 3.64 (read 1) and 3.54 (read 4), with a density of 235. Illumina tech support stated that the run quality is passing but what would cause the same variants to be called for every sample in the run?

Last edited by larangela; 02-27-2017 at 02:27 PM.
larangela is offline   Reply With Quote
Old 02-27-2017, 03:49 PM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,695
Default

Perhaps you can give more details, like the sequencing platform, experiment type, whether the samples are from the same individual, how multiplexing was performed, and so forth? Also, why you expect only 3.64% of the reads to align, which is incredibly low.
Brian Bushnell is offline   Reply With Quote
Old 02-27-2017, 03:56 PM   #3
larangela
Junior Member
 
Location: Orange County, CA

Join Date: Feb 2017
Posts: 3
Default

The TruSeq custom amplicon low input run (from PCR) was sequenced on a MiSeq for various individuals. We used Illumina's library prep kit as well. The align % is very low which was the only portion that seemed strange, per the Illumina Scientist. I'm not quite sure what could have caused this issue but it's occurred in a couple of runs.
larangela is offline   Reply With Quote
Old 02-27-2017, 06:50 PM   #4
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,695
Default

I would not use the data unless you can determine why the mapping rate was so low and ensure it would not cause incorrect results. Most likely, the samples were highly contaminated, which would also explain why all gave the same variant calls. I suggest you BLAST some of the 96% of nonhuman reads to public databases like nt and find out what they hit, and report back. Also, I'm glad it aroused your suspicion that all samples indicated identical variants; that does not usually occur outside of highly inbred villages, unless the variants are close together. It would be helpful if you could describe the methodology for variant-calling, and sample selection. Are the samples all from the same race, or same village? And how big is the amplicon region?
Brian Bushnell is offline   Reply With Quote
Old 02-27-2017, 08:15 PM   #5
nucacidhunter
Senior Member
 
Location: Iran

Join Date: Jan 2013
Posts: 1,055
Default

Quote:
Originally Posted by larangela View Post
Hello!

The %PF was >80, % aligned was about 3.64 (read 1) and 3.54 (read 4), with a density of 235.
I wonder if those figures are from secondary analysis. Those to me looks like the sequencing metrics reported by Illumina SAV. %PF would be % of clusters passing filter and 3.64 and 3.54 %reads aligned to PhiX spike in.

PS. Analysis can be done using the TruSeq Amplicon BaseSpace App or MiSeq Reporter TruSeq Amplicon Workflow (on the MiSeq).

Last edited by nucacidhunter; 02-27-2017 at 11:09 PM.
nucacidhunter is offline   Reply With Quote
Old 02-28-2017, 07:52 AM   #6
larangela
Junior Member
 
Location: Orange County, CA

Join Date: Feb 2017
Posts: 3
Default

Thank you for your responses! Those figures are from a secondary analysis derived from SAV. Our primary analysis is done by a qualified third party.
larangela is offline   Reply With Quote
Old 02-28-2017, 11:25 AM   #7
nucacidhunter
Senior Member
 
Location: Iran

Join Date: Jan 2013
Posts: 1,055
Default

I wonder if you have run data through MiSeq Reporter as well and checked design manifest for accuracy of target loci.
nucacidhunter 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:50 PM.


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