SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
setting parameters for Novoalign 2-6 bp insertions deletions 12jrowley2 RNA Sequencing 2 11-01-2013 09:48 AM
Illumina insertions deletions andreitudor Bioinformatics 2 02-14-2012 07:20 AM
PubMed: Strategy for Robust Detection of Insertions, Deletions, and Point Mutations i Newsbot! Literature Watch 0 03-01-2011 03:30 AM
Coverage and deletions Hena Bioinformatics 0 01-05-2011 11:05 AM
How to find insertions in the genome? dina Bioinformatics 0 10-24-2009 11:55 AM

Reply
 
Thread Tools
Old 05-27-2010, 02:42 AM   #1
joa_ds
Member
 
Location: belgium

Join Date: Dec 2008
Posts: 52
Default SOAP2 and deletions/insertions

I was told that SOAP should be able to detect deletions and insertions up to 3bp.

So i made a reference file and a file with a normal read and some reads with 1,2,3 bp deleted and inserted. Just to test how the output format etc would look like and to get optimal parameter settings.

i have used different -e settings etc, but it appears that not a single read with deletion/insertions gets aligned...

All i get is

Begin Alignment ...
8 ok 0.00 sec
Total Reads: 8
Alignment: 1 (12.50%)

-> Only the normal read gets aligned...

Is there a secret parameter i should set somewhere?
joa_ds is offline   Reply With Quote
Old 05-27-2010, 02:53 AM   #2
francesco.vezzi
Member
 
Location: Udine (Italy)

Join Date: Jan 2009
Posts: 50
Default

Hi
I can't help you so much but I can report another strange behaviour of SOAP2. Despite the fact that you vary the -v option that is the total amount of errors (mismatches) allowed in one read the output that you obtain is always the same. I Try to change several parameters like the seed length without understand how to specify an exact number of errors.

Despite the fact that SOAP2 is one of the most used aligners it has some troubles: there is no a clear documentation, the source code is not available and hence we actually don't know what kind of heuristics are applied in order to speed up the query time...

Francesco
francesco.vezzi is offline   Reply With Quote
Old 05-27-2010, 02:56 AM   #3
joa_ds
Member
 
Location: belgium

Join Date: Dec 2008
Posts: 52
Default

Well, I can tell you the following.

I have simulated loads of data. I have generated datasets with more than a billion reads where error profiles/PCR errors/sequencing errors were simulated and acted as a background for real SNPs that needed to be discovered. The real SNPs always come out on top. So far so good, but with deletions/insertions, I am completely confused.

Indeed, the manual/documentation is not clear. Did you have any success contacting the authors?
joa_ds is offline   Reply With Quote
Old 05-27-2010, 03:00 AM   #4
francesco.vezzi
Member
 
Location: Udine (Italy)

Join Date: Jan 2009
Posts: 50
Default

Quote:
Originally Posted by joa_ds View Post
Well, I can tell you the following.
Indeed, the manual/documentation is not clear. Did you have any success contacting the authors?
Yes I contact them more than once.... but they never reply...
It is really strange that the most used aligner (or at least one of the most used) is so obscure.....
francesco.vezzi is offline   Reply With Quote
Old 05-27-2010, 03:28 AM   #5
joa_ds
Member
 
Location: belgium

Join Date: Dec 2008
Posts: 52
Default

I just checked out their site and found the following thing.

SOAPindel: SOAPindel is developed to find the insertion and deletion specially for re-sequence technology : Coming soon.

I am just wondering if indel detection is already included in SOAP2? Has anyone ever seen an indel being detected by using SOAP2?
joa_ds is offline   Reply With Quote
Old 09-20-2010, 01:44 AM   #6
dingxiaofan1
Member
 
Location: Hong kong

Join Date: Jul 2010
Posts: 17
Default Bwa-maq

How about bwa-maq (http://bio-bwa.sourceforge.net/bwa.shtml), Pindel(seems not available now)
dingxiaofan1 is offline   Reply With Quote
Old 09-20-2010, 05:50 AM   #7
drio
Senior Member
 
Location: 4117'49"N / 24'42"E

Join Date: Oct 2008
Posts: 323
Default

Quote:
Originally Posted by joa_ds View Post
Well, I can tell you the following.

I have simulated loads of data. I have generated datasets with more than a billion reads where error profiles/PCR errors/sequencing errors were simulated and acted as a background for real SNPs that needed to be discovered. The real SNPs always come out on top. So far so good, but with deletions/insertions, I am completely confused.

Indeed, the manual/documentation is not clear. Did you have any success contacting the authors?
In your simulations do you introduce indels or it is only single base changes?
How does that affects your results?
__________________
-drd
drio is offline   Reply With Quote
Old 09-21-2010, 04:59 AM   #8
KaiYe
Senior Member
 
Location: amsterdam

Join Date: Jun 2009
Posts: 133
Default

Quote:
Originally Posted by dingxiaofan1 View Post
How about bwa-maq (http://bio-bwa.sourceforge.net/bwa.shtml), Pindel(seems not available now)
Pindel is still being actively updated.
http://www.ebi.ac.uk/~kye/pindel/release/
New version of Pindel is being tested with the following additional functions:
1. Allow sequence errors/SNPs in the same reads containing INDELs/SVs
2. non-template sequence in deletions
3. inversions
4. tandem duplications
5. breakpoints of large insertions
KaiYe 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 05:31 PM.


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