Question: Incorrect codon sequences reported in kissplice2reftranscriptome output?
gravatar for andres.firrincieli
2.8 years ago by
andres.firrincieli420 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,


kissplice snp k2rnt • 743 views
ADD COMMENTlink modified 2.7 years ago by vincent.lacroix130 • written 2.8 years ago by andres.firrincieli420
gravatar for vincent.lacroix
2.7 years ago by
vincent.lacroix130 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



ADD COMMENTlink written 2.7 years ago by vincent.lacroix130

Dear Vincent,

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

Best, Andrea

ADD REPLYlink written 2.7 years ago by andres.firrincieli420
gravatar for vincent.lacroix
2.8 years ago by
vincent.lacroix130 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 !


ADD COMMENTlink written 2.8 years ago by vincent.lacroix130
Please log in to add an answer.


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