Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Is it possible to convert a SNP.txt to a bed file or get a SNP.bed from samtools? Ling Bioinformatics 7 04-02-2015 07:17 AM
SNP function prediction in non-human species spark Bioinformatics 7 01-22-2013 02:45 AM
Samtools SNP Caller versus MAQ SNP caller BertieWooster Bioinformatics 5 05-18-2012 02:41 PM
functional prediction SNP airtime Bioinformatics 1 10-24-2011 02:37 PM
Varscan for snp detection instead of gatk madsaan Bioinformatics 3 07-27-2011 05:31 AM

Thread Tools
Old 06-28-2012, 03:30 PM   #1
Location: US

Join Date: May 2012
Posts: 80
Default samtools,varscan-SNP prediction

Hi all,
I am new to NGS analysis. I have illumina RNA-seq paird end reads (1.pe1.fastq,2.pe2.fastq,..4.pe1.fastq,4.pe2.fastq) and assembled transcriptome reference sequence(reference.fa). I have aligned the illumina reads (1.pe1.fastq,1.pe2.fastq..4.pe1.fastq,4.pe2.fastq) to the reference transcriptome (reference.fa) using bowtie. Then created mpileup file using samtools. I have predicted SNPs using variant caller varscan. I have found SNPs predicted for each contigs, also with homozygosity & heterozygosity predicted.

1.How can I get distribution of SNPs in each chromosome? (because varscan reports for each contig)
2.From varscan prediction, I get heterozygosity predicted more (about 95%) than heterozygous (5%). I would like know whether my prediction is biologically relevant?
bioman1 is offline   Reply With Quote

mpileup, samtools, varscan

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 10:21 PM.

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