SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
parameter set of bwa-samtools for snp calling heiya De novo discovery 10 04-25-2012 06:27 PM
Set ssaha2 parameter for mismatch xpchen1011 Bioinformatics 0 03-12-2012 09:02 PM
How can i set the "mismatch" parameter in bwa genelab Bioinformatics 14 01-12-2011 02:08 PM
where can I get ref seq hg18 culmen Bioinformatics 2 07-20-2010 12:28 PM
what can be done with micro rna seq from SOLiD of an organism with no ref? KevinLam Bioinformatics 0 06-02-2010 07:45 PM

Reply
 
Thread Tools
Old 09-09-2014, 10:51 AM   #1
david_d
Junior Member
 
Location: Utah

Join Date: Sep 2014
Posts: 1
Default cnv-seq - what do I set the --ref parameter to?

I am trying to detect cnv's in a bam file of an entire human genome using cnv-seq. I was able to create the "hits" file from the bam file using samtools view.

My problem is that I do not know what to put in for the --ref parameter. From what I understand I need to have a hits file from the human reference genome. I have the human reference genome in fasta file format but that will not help me to make a "hits" file. What do I do?

Thank you in advance for your help!

David
david_d is offline   Reply With Quote
Old 07-01-2015, 08:38 AM   #2
vbajic
Junior Member
 
Location: Europe

Join Date: Jul 2015
Posts: 2
Default

Quote:
Originally Posted by david_d View Post
My problem is that I do not know what to put in for the --ref parameter. From what I understand I need to have a hits file from the human reference genome. I have the human reference genome in fasta file format but that will not help me to make a "hits" file. What do I do?

Thank you in advance for your help!
David
I was looking for the explanation of parameters in CNV-seq and I came across your post.
I guess it's already way too late for the answer... however it might be useful to someone.

--ref should be another hits file. For that you need to have another genome sequenced and same as for your --test you need to map the reads to the same reference genome and extract hits file. CNV-seq will compare coverage for --ref and --test hits files. You can download raw sequence reads from for example NCBI site

Ideally for your purpose would be to have raw sequence reads from same human reference as the one you have in fasta file and onto which you mapped your reads.

Cheers!
vbajic is offline   Reply With Quote
Reply

Tags
cnvseq, cnvseq reference

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 05:37 AM.


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