SEQanswers

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 06:17 AM
SNP function prediction in non-human species spark Bioinformatics 7 01-22-2013 01:45 AM
Samtools SNP Caller versus MAQ SNP caller BertieWooster Bioinformatics 5 05-18-2012 01:41 PM
functional prediction SNP airtime Bioinformatics 1 10-24-2011 01:37 PM
Varscan for snp detection instead of gatk madsaan Bioinformatics 3 07-27-2011 04:31 AM

Reply
 
Thread Tools
Old 06-28-2012, 02:30 PM   #1
bioman1
Member
 
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
Reply

Tags
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:02 AM.


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