Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Compare VCF Files lyz1030 Bioinformatics 8 11-13-2013 03:59 AM
how to compare tophat output files with and without "_random" sequences EA01 Illumina/Solexa 2 06-21-2013 01:05 AM
compare Cuffdiff output to microarray data? wellswu RNA Sequencing 0 02-28-2013 05:55 AM
looking for a software to compare genotypes between VCF files? caswater Bioinformatics 5 06-14-2012 03:22 PM
VarScan output to VCF OlgaMikh Bioinformatics 4 06-08-2012 12:46 AM

Thread Tools
Old 10-15-2013, 03:59 PM   #1
Senior Member
Location: .

Join Date: Mar 2011
Posts: 157
Default vcf-compare output

I am looking at some VCF files; I like the output from vcf-compare, eg:

#SN Summary Numbers. Use `grep ^SN | cut -f 2-` to extract this part.
SN Number of REF matches: 297982
SN Number of ALT matches: 262298
SN Number of REF mismatches: 4131
SN Number of ALT mismatches: 35684
SN Number of samples in GT comparison: 0

but now want to get the actual lines from my VCF that correspond to what is of interest (eg number of REF mismatches). I am intermediate at Perl and had a look at the script but it is WAY over my head! Can anyone give me any ideas on how to cajole the REF mismatches from my VCFs?
bruce01 is offline   Reply With Quote

perl, variant analysis, vcf file, vcftools

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 06:07 AM.

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