Tweet
Converting .sam to .vcf
I have some files in the .sam format which I obtained after I ran Stampy on some Cancer reads. I now have to convert the .sam files to .vcf files to do some analysis and am clueless on how to do the conversion. Any help or hints would be much appreciated.
-
-
Use samtools with mpileup and use bcftools command:
samtools and bcftools are part of the samtools collection.
PLEASE DO examine the options to each command and understand them.
samtools mpileup -R -d 1000 -l rangehg18.chr -uf hg18.fa test.bam | bcftools view -g - > test.vcf
additional note: you will want to convert your sam into a sorted bam ("test.bam" in the example above)Last edited by Richard Finney; 06-19-2012, 07:34 AM.
-
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...05-06-2024, 07:48 AM -
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...04-22-2024, 07:01 AM
Comment