[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.

Data*****
@ILLUMINA-554790:34:6:107:18082:1570#0/1
AATCGTTGTAGATTCATTACTGCTCCTGAGCGGTTTACAGGAGTAAGTGTAGACGTTGAGAAATGATATTTGTGAA
+ILLUMINA-554790:34:6:107:18082:1570#0/1
efgggfffffffffffggaggggggggg[ggg[gggggagfdfafffaff_aadc^dfYaffcffKaca_^aXcBB
@ILLUMINA-554790:34:6:107:18268:1571#0/1
AGCAAGCAATCGCCAGCAGCAGAGTGCCAAGACGCCCCAGCGCAAACCACAGTCTGCCCGCTTTGCTGTTGCTGGT
+ILLUMINA-554790:34:6:107:18268:1571#0/1
cccSc]\LYSaYYbKcc[ccccccacY_c^aRaZSca\\OYaaaaJR`X\]]FHX^aQaHR\PWSaaccaaaILca
@ILLUMINA-554790:34:6:107:18774:1566#0/1
CCAGGCAAAGAGATCGCGCATAGCGGTAAATTTATTCGACAACAGCAGCCAGATAGCAAAATCACGCAGCAGCCAG
+ILLUMINA-554790:34:6:107:18774:1566#0/1
hhhhhhhhhhhghggghfhhhehhhhffhfhhhhhhhdahghchfaggafhcah[ffhefhcehffahhaffdh_a

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')
end

*******************

But,
p entry.format   returns nil
and
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
http://lists.open-bio.org/mailman/listinfo/bioruby

--
Ra
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