Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • bnfoguy
    Member
    • May 2011
    • 17

    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.
  • Richard Finney
    Senior Member
    • Feb 2009
    • 701

    #2
    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.

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 06-09-2026, 11:58 AM
    0 responses
    17 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-05-2026, 10:09 AM
    0 responses
    27 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-04-2026, 08:59 AM
    0 responses
    38 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 06-02-2026, 12:03 PM
    0 responses
    61 views
    0 reactions
    Last Post SEQadmin2  
    Working...