FU Logo
  • Startseite
  • Kontakt
  • Impressum
  • Home
  • Listenauswahl
  • Anleitungen

[Seqan-dev] Fwd: Re: bug ?!

<-- thread -->
<-- date -->
  • From: Laurent MANCHON <lmanchon@univ-montp2.fr>
  • To: SeqAn Development <seqan-dev@lists.fu-berlin.de>
  • Date: Tue, 17 Jan 2012 23:01:53 +0100
  • Organization: SKT
  • Reply-to: lmanchon@um2.fr, SeqAn Development <seqan-dev@lists.fu-berlin.de>
  • Subject: [Seqan-dev] Fwd: Re: bug ?!

--Hi,

is anyone able to fix the bug in the FAR software developed by Matthias Dodt FAR?

thank you --



Hi Laurent!

Sorry for the delay in reply. Yes you are right, there is a bug in the
software. Unfortunately i cannot fix this at the moment since i have
moved out of sience several month ago. I contacted my group but
currently there seems to be nobody taking care of this. I think it
will be fixed in the future, but i can't tell you when. I will get
back to you if i have further information.

thanks,

mat

2012/1/6 laurent manchon<lmanchon@users.sourceforge.net>:
 --hello,

 i have some problem to clean my sequences using this adaptor string:
 AATCTCGTATGCCGTCTTCTGCTTGC

 this is my input file:
sequence_100104
 ATAGTCAACGGCTGAGGATTTGGCATCTCGTA
sequence_100105
 ATAGTCAAGATAGAGCTACTGGAGAATCTCGA
sequence_100106
 ATAGTCAAGATAGAGCTACTGGAGAATCTCGT
sequence_100107
 ATAGTCAAGATAGAGCTACTGGGGAATCTCGT


 i use this command to clean:

 far --source reads.fa --target reads_no_adapters.fa --adapters myAdapters.fa --format
 fasta --cut-off 4 --min-overlap 2 --min-readlength 16 --trim-end right

 and this is the output:

sequence_100104
 ATAGTCAACGGCTGAGGATTTGGCATCTCGTA
sequence_100105
 ATAGTCAAGATAGAGCTACTGGAGAATCTCGA
sequence_100106
 ATAGTCAAGATAGAGCTACTGGAG
sequence_100107
 ATAGTCAAGATAGAGCTACTGGGG


 So, sequence_100106 and sequence_100107 are well cleaned but not
 sequence_100104, the good results should be:

sequence_100104
 ATAGTCAACGGCTGAGGATTTGGC
sequence_100105
 ATAGTCAAGATAGAGCTACTGGAGAATCTCGA
sequence_100106
 ATAGTCAAGATAGAGCTACTGGAG
sequence_100107
 ATAGTCAAGATAGAGCTACTGGGG

 substring \"ATCTCGTA\" is a part of adaptor sequence !

 Somebody has a solution to resolve this problem ?

 thank you.







 --
 This message was sent to your SourceForge.net email alias via the web mail form.  You may reply to this message directly, or via https://sourceforge.net/sendmessage.php?touser=2041456
 To update your email alias preferences, please visit https://sourceforge.net/account





<-- thread -->
<-- date -->
  • seqan-dev - January 2012 - Archives indexes sorted by:
    [ thread ] [ subject ] [ author ] [ date ]
  • Complete archive of the seqan-dev mailing list
  • More info on this list...

Hilfe

  • FAQ
  • Dienstbeschreibung
  • ZEDAT Beratung
  • postmaster@lists.fu-berlin.de

Service-Navigation

  • Startseite
  • Listenauswahl

Einrichtung Mailingliste

  • ZEDAT-Portal
  • Mailinglisten Portal