Question: Incorrect codon sequences reported in kissplice2reftranscriptome output?
0
gravatar for andres.firrincieli
3 months ago by
andres.firrincieli0 wrote:

Dear Kissplice users,

after running kissplice2reftranscriptome

kissplice2reftranscriptome -b assemblychec095n10_okay_okalt_newnames.fasta.transdecoder.bed results_mDr1_EcoC_mDr2_EcoC_mDr3_EcoC_WW1_EcoC_WW2_EcoC_WW3_EcoC_k25_coherents_type_0a.fa -t out_blat_SNP_on_transcripts_80.psl -s kissDE_output_C_pvalue05

and reading the mainOutput[timestamp].tsv file, somenthing is not clear for me. As reported in the UserGuide, in the columns Codon_1 and Codon_2 should be reported
the SNP(s) from the upper and lower paths respectively. However, the codons reported in both columns (Codon_1 & Codon_2) are the same. This happen for all the SNPs; is this normal?

Thanks in advance for your time,

Best

kissplice snp k2rnt • 180 views
ADD COMMENTlink modified 11 weeks ago by vincent.lacroix40 • written 3 months ago by andres.firrincieli0
1
gravatar for vincent.lacroix
11 weeks ago by
vincent.lacroix40 wrote:

Dear Andres, Thanks for sending the files. It seems that the problem comes from the fact that you used a different kmer value in KisSplice and KisSplice2RefTranscriptome. We are not catching this error correctly in kissplice2reftranscriptome and it was indeed not obvious to know what the problem was. Thanks for reporting this bug. In the next release, we will clarify the manual, and add an error message.

Using the current version, the solution for you is simply to add the option -K 25 in kissplice2reftranscriptome

Best,

Vincent

ADD COMMENTlink written 11 weeks ago by vincent.lacroix40

Dear Vincent,

thanks a lot for your help. I will try the single -K solution.

Best, Andrea

ADD REPLYlink written 11 weeks ago by andres.firrincieli0
0
gravatar for vincent.lacroix
12 weeks ago by
vincent.lacroix40 wrote:

Dear Andres,

Columns Codon 1 and Codon 2 should be distinct, otherwise there would be no SNP. Could you send me a sample of your data so that we try to reproduce your bug ? You can send directly at vincent dot lacroix at univ-lyon1 dot fr

Thanks !

Vincent

ADD COMMENTlink written 12 weeks ago by vincent.lacroix40
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: 761 users visited in the last hour