SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
ERROR: [bcf_sync] incorrect number of fields.. vinodhsri Bioinformatics 21 08-11-2017 01:42 PM
CIGAR string from BWA-SW output incorrect ? robs Bioinformatics 13 01-13-2012 04:07 AM
Cufflinks calls incorrect read length sclindsay Bioinformatics 1 10-16-2011 01:22 PM
BWA generating incorrect CIGAR string? foxyg Bioinformatics 6 09-16-2011 11:22 AM
SAMTools showing incorrect depth for a given position rdeborja Bioinformatics 1 05-07-2011 03:25 PM

Reply
 
Thread Tools
Old 10-12-2011, 04:27 PM   #1
Phillip Morin
Junior Member
 
Location: San Diego, CA

Join Date: May 2011
Posts: 4
Question incorrect fastq calls by vcf2fq

We are doing SNP discovery from Illumina reads assembled to a set of ~80 loci, from ~25 individuals. After generating the BAM and VCF files with Samtools, we use vcf2fq to generate the Fastq files from the VCF files. In general this works well, but I have found that when there is only a single read at a site, if it differs from the reference sequence, it is still called as the reference sequence rather than the alternate base (or N). Does anyone know what parameter would be changed to generate the correct Fastq sequence even when there is only one read, or to change it to an N when it differs from the reference with only one read?
I'm not the person who put all of these scripts together, so I may have explained our analysis pipeline incorrectly, but the important part is getting the appropriate function out of vcf2fq.
Phillip Morin 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 08:31 AM.


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