Welcome, Guest. Please Login or Register
UGENE Bulletin Board
  Welcome to our forum.
  HomeHelpSearchLoginRegister  
 
 
Pages: 1 2 
How to create consensus sequence from forward and reverse chromatogram (Read 16995 times)
Jul 16th, 2019 at 12:20am

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Hello,

I am a German mycologist and new to this forum.

My problem: I want to create the consensus sequence from a forward and a reverse chromatogram.
To do this I start with Tools -> Sanger data analysis -> Reads de novo assembly (with CAP3)...
But when I add the two chromatograms with "Add... > Open... > Run" then the following error message is created: "[19:11:19] 'Load document: 'Russula exalbicans B9-18461_1F-ITS4.cap.ac..."

Can you please help me?
Forward chromatogram *ITS1F.ab1 and reverse chromogram *.ITS4.ab1 are attached.

Bernie
« Last Edit: Jul 16th, 2019 at 2:54am by Bernie »  
IP Logged
 
Reply #1 - Jul 16th, 2019 at 1:11pm

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
Hello,

I'm unable to open the attached files. Could you please try to share them again or, alternatively, share the files in another way (e.g. using Google Drive)?

 
IP Logged
 
Reply #2 - Jul 16th, 2019 at 8:47pm

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Hello,

when I try to send a dropbox link then I get an error message

Sorry, you are not allowed to post messages containing active links to websites or images before you have posted 5 normal messages. 

Bernie
 
IP Logged
 
Reply #3 - Jul 16th, 2019 at 10:55pm

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Hello,
I have sent you the chromatogram files via researchgate.net. Smiley
Bernie
 
IP Logged
 
Reply #4 - Jul 17th, 2019 at 2:49pm

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
Hello Bernie.

It appears that that there are too many mismatches between the reads. I think, this is the reason why CAP3 was not able to produce the result.

I had to make a workaround to align the reads:
* Open both *.ab1 files in UGENE.
* In the Project View at the left part of the UGENE window select "[s]" objects of the *.ab1 files. In the context menu select "Export sequences as alignment". Thus, the reads are opened in the Alignment Editor.
* Select the second sequence in the Alignment Editor and convert it to it's reverse-complement (select "Actions > Edit > Replace selected rows with reverse-complement").
* On the "Pairwise Alignment" tab of the options panel in the Alignment Editor select the both sequences, optionally configure other parameters and click "Align".

See the result alignment in the attachment.

Note that you may select different options for calculating consensus of the reads. And export the consensus sequence, if required.
 
IP Logged
 
Reply #5 - Jul 17th, 2019 at 11:21pm

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Hello Olga,

thank you for the hints!
Is it possible to show both chromatograms aligned one above the other and to edit them manually to optimize the consensus sequence?  Smiley

Bernie
 
IP Logged
 
Reply #6 - Jul 18th, 2019 at 11:30am

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
Bernie, you're welcome!

Quote:
Is it possible to show both chromatograms aligned one above the other and to edit them manually to optimize the consensus sequence?

In the current UGENE version you will have to export the consensus sequence and then map the reads to it as reference ("Tools > Sanger data analysis > Map reads to reference"), tweak the input parameters, if required.

The result will be opened in the Sanger Reads Editor (https://www.youtube.com/watch?v=lDovNM1oZEw).
 
IP Logged
 
Reply #7 - Jul 18th, 2019 at 12:04pm

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
When I try to do this it results in an error (please see attachment)
Bernie
 
IP Logged
 
Reply #8 - Jul 18th, 2019 at 2:40pm

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
Which consensus sequence did you use? I've tried to export the consensus, calculated with the "Levitsky" algorithm and the reads were successfully mapped to it. See in the attached archive the consensus and the *.ugenedb file that can be opened with the Sanger Reads Editor.
 

Archive.zip (47 KB | 266 )
IP Logged
 
Reply #9 - Jul 21st, 2019 at 4:12am

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Hello Olga,

now it works well!! But it is a long and "dangerous" way to reach the Sanger reads mapping!  I have chosen the Levintsky Consensus type.

I would like to ask you two more questions:

a) When exporting consensus there is a small box "Keep gaps". I left it open. Is that correct?


b) When at the end editing the chromatograms to optimize the consensus sequence: Is it possible to search a string whithin a chromatogram? For instance the string "GGATCATTA" or "TTGACCTCAAATC".

Thank you in advance!

Bernie
 
IP Logged
 
Reply #10 - Jul 22nd, 2019 at 3:26pm

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
Quote:
now it works well!!

Great!

Quote:
But it is a long and "dangerous" way to reach the Sanger reads mapping!

In future we plan to open the Sanger reads de novo assembly result directly in the Sanger Reads Editor.
Maybe we should also choose some other algorithm for de novo assembly, so that reads with low similarity (like in this case) are also taken into account.

Quote:
a) When exporting consensus there is a small box "Keep gaps". I left it open. Is that correct?

It doesn't matter for this scenario. Gaps ("-") may be inserted into the reference on the mapping step, if required (i.e. if the reads contain some insertion in comparison with the reference).

Quote:
b) When at the end editing the chromatograms to optimize the consensus sequence: Is it possible to search a string whithin a chromatogram? For instance the string "GGATCATTA" or "TTGACCTCAAATC".

Sorry, but this is one more issue we need to improve. You can do the search in the Sequence View, for example, but for now you can't do it in the Sanger Reads Editor.
 
IP Logged
 
Reply #11 - Jul 22nd, 2019 at 5:48pm

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Thank you once more!
« Last Edit: Jul 23rd, 2019 at 12:56am by Bernie »  
IP Logged
 
Reply #12 - Jul 22nd, 2019 at 5:50pm

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
You're welcome!
 
IP Logged
 
Reply #13 - Sep 26th, 2019 at 5:21pm

Bernie   Offline
YaBB Newbies

Gender: male
Posts: 9
*
 
Hello,

what do you think, when will the next version of UGENE be on line?

Bernie
 
IP Logged
 
Reply #14 - Sep 26th, 2019 at 5:27pm

Olga Golosova   Offline
YaBB Administrator

Posts: 338
*****
 
The next version 33 will come out very soon, but it will not contain changes for Sanger reads de novo assembly.
 
IP Logged
 
Pages: 1 2