SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
extremely large vcf file generated by GATK slowsmile Bioinformatics 2 12-09-2012 08:23 PM
How to get statistics data from sam/bam file generated by trinity packages taoxiang180 Bioinformatics 1 07-05-2012 12:46 AM
Problem with vcf file generated amruta.bn Bioinformatics 0 01-05-2012 03:22 AM
BAMseek large file viewer - now with VCF and FASTQ support BAMseek Bioinformatics 1 07-25-2011 06:57 PM
wig file generated by MACS doesn`t fit input sam tujchl Bioinformatics 2 06-16-2011 12:23 AM

Reply
 
Thread Tools
Old 12-18-2012, 08:43 PM   #1
janaahan13
Junior Member
 
Location: CA,USA

Join Date: Sep 2012
Posts: 6
Default extremely large vcf file generated using sam tools

Hi,

This is my first time to Samtools. I tried to create .vcf files for my bam files of my paired-end RNA seq samples. I used

Code:
samtools mpileup -D -B -Q 20 -u -f /WholeGenomeFasta/genome.fa Sample.bam | bcftools view -gc - > Sample_positions.vcf
But for ~4GB of bam file it is creating(still running) ~ 75GB size of vcf file.

Is this normal or am I making a mistake some where?

Thanks in advance.

Last edited by janaahan13; 12-18-2012 at 08:46 PM.
janaahan13 is offline   Reply With Quote
Reply

Tags
rna seq, sam tools, vcf

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:13 AM.


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