[Biopython] reading two fastq files at the same time

xyz mitlox at op.pl
Mon Sep 6 09:24:14 EDT 2010


Hi,
How is it possible to read two fastq files at the same time in 
BioPython? I have the following BioRuby example:

require 'bio'

begin
   fq1 = Bio::FlatFile.open(Bio::Fastq, ARGV[2])
   fq2 = Bio::FlatFile.open(Bio::Fastq, ARGV[3])

   while (entry1 = fq1.next_entry) and (entry2 = fq2.next_entry)

     fastq_A1 = entry1.entry_id
     fastq_A2 = entry1.seq

     fastq_B1 = entry2.entry_id
     fastq_B2 = entry2.seq
   end

rescue => err
   raise "Exception: #{err}"
end

Thank you in advance.


More information about the Biopython mailing list