SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Snp discovery without a reference lletourn De novo discovery 38 01-12-2016 05:15 AM
reference-free SNP discovery Marius De novo discovery 5 03-30-2011 11:23 AM
compare CLC SNP's with dbSNP boetsie Bioinformatics 1 02-08-2011 01:25 AM
HELP: some suggestions for SNP discovery in 454? linikujp Bioinformatics 1 04-07-2010 12:39 AM
Nonsynonymous SNP (nsSNP) discovery tools? jpeaco02 Bioinformatics 2 11-08-2009 01:13 PM

Reply
 
Thread Tools
Old 04-11-2011, 05:27 AM   #1
extari
Junior Member
 
Location: Bilbao (Spain)

Join Date: Mar 2011
Posts: 9
Default CLC and SNP discovery

Hi there,

I am using CLC 4.5 software for doing a de novo assembly and SNP discovery.
I have already perform my de novo assembly, and I have mapped my reads to contigs with the main objective of SNP discovery.
This software (CLC) has the option "SNP Detection", I have done it and I have discovered SNPs (after filtering).
Does anyone of all know how can I get my consensus sequences with the ambiguities in one fasta file, or txt format?

I mean, imagine this sequence with this SNP (*), and these ambiguities (-)
5` ... A G G G A T A G C C C G A T A G G G T A C A G G T T A ... 3`
- * - -
When I try to save my consensus sequence as a txt file, I can only save as one of these:
a) if in the assembly selected VOTE (for ambiguities):
5` ... A G G G A T A G C C C G A T A G G G T A C A G G T T A ... 3`

b) if in the assembly selected Unknown (for ambiguities):
5` ... A G G G N T A G C C N G A T N G G G T A N A G G T T A ... 3`

c) if in the assembly selected Ambiguity (for ambiguities):
5` ... A G G G A R A G C C C R A T M G G G T A Y A G G T T A ... 3`

But the one I want to save is this one: 5` ... A G G G A T A G C C C R A T A G G G T A C A G G T T A ... 3` (with the only annotated SNP).

Could you help me, please???


Thank you!
extari is offline   Reply With Quote
Old 04-11-2011, 03:46 PM   #2
flipwell
Member
 
Location: Queensland

Join Date: Feb 2011
Posts: 14
Default

I haven't tried it, but during SNP detection you can select to annotate the consensus with SNPs or not - if you choose vote for ambiguities during alignment then annotate with SNPs during detection I'm guessing it should give you the result you are after
flipwell is offline   Reply With Quote
Old 04-13-2011, 05:45 AM   #3
extari
Junior Member
 
Location: Bilbao (Spain)

Join Date: Mar 2011
Posts: 9
Thumbs down

Hi ya,

I have tried to do what you have told me, but that one is not the solution, look what happens when I do that (see the attachment).
The SNP is annotated, but I cannt have the fasta file as ambiguities (R) or as [A/G], what is the thing I would like, there must be a way to do that, but, how?
Attached Images
File Type: png Screenshot.png (7.6 KB, 29 views)
extari is offline   Reply With Quote
Old 04-13-2011, 07:55 AM   #4
NextGenSeq
Senior Member
 
Location: USA

Join Date: Apr 2009
Posts: 482
Default

That's a terrible alignment. You should increase the gap penalty.
NextGenSeq is offline   Reply With Quote
Old 04-13-2011, 10:01 AM   #5
extari
Junior Member
 
Location: Bilbao (Spain)

Join Date: Mar 2011
Posts: 9
Default

no, don't worry about the aligment, that was a bad example.
the question is how can i have the fasta file with the polymorphic site?
extari is offline   Reply With Quote
Old 04-13-2011, 10:20 PM   #6
szilva
Member
 
Location: Hungary

Join Date: Aug 2009
Posts: 16
Default

If you point to the sequence name on the left side of the panel and right click on the name, you should get a menu like "Open Copy of Sequence" . Once you have got a copy, you can export that into FASTA.
szilva is offline   Reply With Quote
Old 04-14-2011, 12:24 AM   #7
extari
Junior Member
 
Location: Bilbao (Spain)

Join Date: Mar 2011
Posts: 9
Default

Yes, that is true, but the problem persists, in this way I cannot have the polymorphism annotation in the sequence (as "R" or "[A/G]") in the FASTA file.
I don't know how to do this, I am desperate!!!
extari is offline   Reply With Quote
Old 04-14-2011, 01:10 AM   #8
sklages
Senior Member
 
Location: Berlin, DE

Join Date: May 2008
Posts: 628
Default

what does CLC recommend? AFAIK their support is quite helpful.
sklages is offline   Reply With Quote
Old 04-14-2011, 01:13 AM   #9
extari
Junior Member
 
Location: Bilbao (Spain)

Join Date: Mar 2011
Posts: 9
Default

In the user guide there is nothing about this issue.
I have sent them an e-mail, I am waiting for response. If I can solve this issue, I will tell you, but I think there is no way to do this...
extari is offline   Reply With Quote
Old 04-15-2011, 01:32 AM   #10
extari
Junior Member
 
Location: Bilbao (Spain)

Join Date: Mar 2011
Posts: 9
Unhappy CLC support has answered

Hello all,

CLC support has answered to my question in a e-mail. They have told me this: "I'm afraid there is no way to export fasta sequence with ambiguities or IUPAC letters from the CLC software. However, if you have chosen to annotate your Consensus, and/or your Reference sequence, you could then export to GFF or Genbank format. The SNP variants that have been called will appear in the annotations. At that point though, you would have to script something to generate files of the format that suits you needs."

So, there is no way to do this...

If anyone knows about these type of script, please, let me know.

Thank you!
extari 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:08 PM.


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