[BioRuby] Fwd: mask on Ilumina-format

遠藤大二 daijiendoh at gmail.com
Tue Mar 29 09:51:25 UTC 2011

Thank you Nishiyama-san and Bonnal-san

1 I use bioruby on Ubuntu10.10 with synnaptic-mediated installation.
The version is 1.40-2. So I cannot use mask. And I will install
bioruby manually.
2 To setup fastq_illumina -format I had not been know. The scripts below.

q = Bio::Fastq.new("@test\natgcgataatgc\n+test\nCCCCDDDDBXXX\n")
 q.format = :fastq_illumina

As it is very difficult for me how to make a new object with certain
format. Mr Bonnal's suggestion is very helpful for me.

Thanks again

---------- 転送メッセージ ----------
From: Raoul Bonnal <bonnal at ingm.org>
日付: 2011年3月29日16:51
件名: Re: [BioRuby] mask on Ilumina-format
To: 遠藤大二 <daijiendoh at gmail.com>
Cc: bioruby at lists.open-bio.org

>From a Goto-san's mail:
 q = Bio::Fastq.new("@test\natgcgataatgc\n+test\nCCCCDDDDBXXX\n")
 q.format = :fastq_illumina
 puts q.mask(your_quality_threshold, your_masking_char)
keep in mind mask masks the bases < your_quality_threshold
by default the format is :fastq_sanger if you call q.quality_scores (
without specify q.format)
On 29/mar/2011, at 01.45, 遠藤大二 wrote:

Dear All
I start to use illumina-formated data.

On the start, data was parsed as below.


Parsing script*****

require 'bio'
quality_threshold = 60
Bio::FlatFile.open('text.txt').each do |entry|
      p entry.format
      p entry.entry_id
      p entry.seq
      hq_seq = entry.mask(quality_threshold, mask_char='n')


p entry.format   returns nil
entry.mask    returns    undefined method `mask' for
#<Bio::Fastq:0x7f1eb5d9dc28> (NoMethodError)

How I set format to FASTQ_ILLUMINA and use 'mask' method.
'http://bioruby.open-bio.org/rdoc/'  is so difficult for me and need
help to use it.

With best wishes,
Daiji Endoh
BioRuby Project - http://www.bioruby.org/
BioRuby mailing list
BioRuby at lists.open-bio.org

linkedin: http://it.linkedin.com/in/raoulbonnal
twitter: http://twitter.com/ilpuccio
skype: ilpuccio
irc.freenode.net: Helius
github: https://github.com/helios

More information about the BioRuby mailing list