SEQanswers

Go Back   SEQanswers > Introductions



Similar Threads
Thread Thread Starter Forum Replies Last Post
bwa multiple alignment scami Bioinformatics 0 08-28-2014 10:27 PM
SNP calling from multiple genotypes sgcsd Bioinformatics 4 11-21-2013 01:52 PM
BWA alignment of multiple read files, the order of processing matters? gene_x Bioinformatics 8 07-12-2013 12:20 PM
Order of genotypes likelihoods for multiple alleles in vcf Karen Chait Berman Bioinformatics 0 07-18-2012 11:38 PM
Calling genotypes on multiple trios a11msp Bioinformatics 0 04-17-2012 01:50 AM

Reply
 
Thread Tools
Old 12-14-2016, 01:06 PM   #1
shis
Member
 
Location: USA

Join Date: Apr 2014
Posts: 16
Default BWA alignment for multiple genotypes

Hello expert members,
I would like to get your advice regarding BWA alignment options. I found in the bwa manual that it uses only two read files for paired-end and one read file for single-end data at a time. However, I have 3 genotypes (each has two fastq files) to identify and compare SNPs among them. According to the BWA guideline, I used two (read1.fastq and read2.fastq) files of one genotype at a time to make alignment against the reference genome using the following command line:

bwa mem ref.fasta BH_1.fastq_filtered BH_2.fastq_filtered > aln-BH.sam
bwa mem ref.fasta Cyp_1.fastq_filtered Cyp_2.fastq_filtered > aln-Cyp.sam
bwa mem ref.fasta SR_1.fastq_filtered SR_2.fastq_filtered > aln-SR.sam

Now I have 3 *.sam files for the 3 genotypes for SNPs variant calling. My question - is it possible to use all 6 fastq files from 3 genotypes in the single command to get a single *.sam file like below:

bwa mem ref.fasta BH_1.fastq_filtered BH_2.fastq_filtered Cyp_1.fastq_filtered Cyp_2.fastq_filtered SR_1.fastq_filtered SR_1.fastq_filtered > aln-bwa.sam

Thank you so much for taking your valuable time and I would appreciate your advice.

Thanks
Shofi
shis is offline   Reply With Quote
Old 01-12-2017, 12:58 PM   #2
shis
Member
 
Location: USA

Join Date: Apr 2014
Posts: 16
Default How to compare vcf file containing snps between two cultivars

Hi all,
I have three vcf files containing SNPs of three different cultivars. The SNPs were called using GATK's UnifiedGenotyper program. After stringent filtering with GATK's hard filters, I have approx. approx. 3.5 millions SNPs in the three cultivars. Now I would like to compare between two cultivars based on the vcf files to identify common and specific SNPs in vcf1/vcf2, vcf1/vcf3 and vcf2/vcf3, and also like to generate venn diagram.

Can anyone suggest me - which tool I can use to identify common and specific SNPs between two cultivars? I do appreciate your help.

Last edited by shis; 01-21-2017 at 08:18 AM.
shis is offline   Reply With Quote
Old 01-20-2017, 10:58 AM   #3
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,552
Default

Coincidentally, I just released a program for this purpose, in the BBMap package. You can use it like this:

Code:
comparevcf.sh in=sample1.vcf,sample2.vcf out=difference.vcf subtract
It supports union, intersection, and subtraction. Note that subtraction is order-dependent while the other operations are not.
Brian Bushnell 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 02:41 AM.


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