Tool: DiscoSnp becomes DiscoSnp++
4
gravatar for pierre.peterlongo
2.5 years ago by
France
pierre.peterlongo600 wrote:

DiscoSnp major update: DiscoSnp becomes DiscoSnp++

Major modifications were made on DiscoSnp these last few weeks. Improvements are the following:

  • Thanks to the recoding using the GATB library:
    • An even quicker execution speed, and the parallelization of the kissnp2 module.
    • An improved progression message
    • A unique file for storing the graph. This file (.h5) may be used in any GATB tool.
  • DiscoSnp is no more limited to isolated SNP detection:
    • Up to P (parameter) close SNPs may be found within a unique bubble
    • Insertions and deletions of length lower or equal to D (parameter) are also detected.

Regarding these modifications, and as DiscoSnp is not limited only to SNP prediction anymore, we decided to change also its name.
Thus DiscoSnp becomes DiscoSnp++. The first release version is DiscoSNP++-2.0.0.

Original DiscoSnp paper: Reference-free detection of isolated SNPs

Official DiscoSnp++ page: http://colibread.inria.fr/software/discosnp/

Comments and feedbacks are welcome.

indels gatb tool snp discosnp • 1.3k views
ADD COMMENTlink modified 17 months ago by acros10 • written 2.5 years ago by pierre.peterlongo600
2
gravatar for pierre.peterlongo
2.0 years ago by
France
pierre.peterlongo600 wrote:

Hello Paolo,

Thanks for your message.

We don't have much  hindsight with these kind of datasets. We are currently investigating the application of discoSnp on RAD seq data (take care of heterogeneous coverage), but we cannot assess yet precisely the results quality.

We have never conducted any experiments nor had any feedback on applying discoSnp on GBS (I'd be curious to know results).

Pierre

ADD COMMENTlink written 2.0 years ago by pierre.peterlongo600
1
gravatar for ruggeripaolo82
2.0 years ago by
United States
ruggeripaolo8210 wrote:

Hi Pierre,

I'm new of SNPs discovery pipelines and I'm in the need to use a free-reference tool to obtain SNPs from a non-model species. I wonder if DiscoSNP can be able to handle GBS or RAD data to obtain a set of new SNPs loci.

Many thanks! 

Paolo

ADD COMMENTlink written 2.0 years ago by ruggeripaolo8210
0
gravatar for ruggeripaolo82
24 months ago by
United States
ruggeripaolo8210 wrote:

Hi Pierre,

Thanks a lot to have provided me a reply about my previous inquiry! Maybe, I could try to run DiscoSNP with my current dataset but I read the User Guide and I'm not sure it will work well......but if there will be something of interesting I'll let you know! In any case, I hope to find new updates for RAD tags on DiscoSNP, soon, it would be quite adequate for my current data.

Thanks a lot again!!

Cheers!

Paolo

ADD COMMENTlink written 24 months ago by ruggeripaolo8210
0
gravatar for acros
17 months ago by
acros10
acros10 wrote:

Hi Pierre,

DISCOSNP++ seems to be exactly what I am looking for to call SNPs in microsatellite flanking regions for a coral species. However it doesn't seem to detect any SNPs even though I can visually seem them. I am wondering if the sequences are too short? I am inputing a fasta file with sequences at a single loci to test: Here is a short example fof_test.txt: data_test/align_P1.fasta

>S10_C10_1_a1

ATGTGAGAAAGCTAATCCATTTTTAAT

>S10_C12_1_a1

ATGTGAGAAAGCTAATCCATTTTTAAT

>S9_C8_1_a1

ATGTGAGAAAGCTAATCCATTTTTAAT

>S9_C9_1_a2

ATGTGAGAAAGCTAATCCATTTTTAAT

>S10_C18_1_a1

ATAATCCATTTTTAAT

>S10_C19_1_a1

ATAATCCATTTTTAAT

>S7_C13_1_a1

ACGTGAGAAAGCTAATCCATTCTAAT

>S7_C14_1_a2

ACGTGAGAAAGCTAATCCATTCTAAT

>S11_C21_1_a3

ATGTGAGAAAGCTAATCCATTTTTAAC

However I don't get a result file and the results on the terminal show empty bubbles. Any idea of what I could be doing wrong?

Thank you, Annick

 

 

ADD COMMENTlink written 17 months ago by acros10

Hi Annick, 

 

Thanks for your message and your patience :) 

disco finds variants surrounded by at least k nucleotides left and right. By default, k=31 explaining why disco does not find results on your example.

 

Pierre

ADD REPLYlink written 17 months ago by pierre.peterlongo600
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 645 users visited in the last hour