SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat samtools snp calling JasonR Bioinformatics 5 05-13-2013 09:43 AM
problem with samtools and SNP calling lpn Bioinformatics 21 11-05-2012 01:08 PM
snp calling with bwa-samtools heiya Bioinformatics 3 08-21-2012 10:51 PM
Samtools SNP calling vidhya Bioinformatics 3 04-07-2011 07:17 AM
SAMtools and SNP calling Jan Bioinformatics 2 09-16-2010 02:01 PM

Reply
 
Thread Tools
Old 03-19-2013, 07:22 AM   #1
maize
Junior Member
 
Location: NY

Join Date: Apr 2011
Posts: 9
Default snp calling with samtools

hi all,
i am a new for snp calling with samtools in transcriptome data. recently, i am trying to call snps with samtools. I set the option as follows:
samtools mpileup -B -ugSD -f ref.fa aln_sorted.bam | bcftools view -Nbcvg - > aa.bcf
i use the default that samtools manual lists.but i think i need set the option,especially for -D(such as -D100), according to my data, but i don't know the rules or criterion clearly. i think the -D is difficult to set because the data is from RNA-seq. so i am not dry behind the ears for this.
can you help me, dear? thx~
maize is offline   Reply With Quote
Old 03-19-2013, 10:19 AM   #2
Maayanster
Member
 
Location: Vancouver, BC

Join Date: Dec 2012
Posts: 30
Default

Quote:
samtools mpileup -C50 -l location_list -d 1000000 -L 1000000 -AEfu ref.fasta my.bam > result.vcf
This is how I run mpileup on rna-seq. The large numbers for -d and -L are due to very high coverage regions is RNA seq experiments. We didn't want it to ignore high coverage regions, so we set them very high. The -l list is because I only wanted calls for certain locations. -E is extended BAQ computation... don't understand much about that but apparent;y its a good idea. -C50 is recommended for bwa alignments. -u is for getting the output as an uncompressed vcf.
Maayanster is offline   Reply With Quote
Old 03-26-2013, 06:51 PM   #3
leifive
Member
 
Location: Nanjing, China

Join Date: Mar 2013
Posts: 10
Default

Hi, all,
I got a document from , it says using pileup file from samtools and junctions.bed file generated by tophat with bedtools to filter SNPs, how do you think?
leifive 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 01:52 AM.


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