snp calling with samtools

[leifive]

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?

[Maayanster]

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.

[maize]

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~