[Bioperl-l] Reading and writing fastq files

shalabh sharma shalabh.sharma7 at gmail.com
Sun Feb 6 20:42:02 UTC 2011


Hi,
    i am trying to read and write fastq files.
I can read them and can change format (like fastq files to fasta) but when i
try to write them back as 'fastq' format i am getting a warning:

--------------------- WARNING ---------------------
MSG: You can't write FASTQ without supplying a Bio::Seq::Quality object!

and i am not getting getting any output.
I am using bioperl 1.61

My part of code looks like this:

#!/usr/bin/perl -w
use Bio::SeqIO;
$in = Bio::SeqIO->new(-file => "$ARGV[0]", -format => 'fastq');
- - - - -- -
 - - - -- - -

while(my $seq = $in->next_seq){
         -- - - - -  -
         - - - - - - -
         $outr->write_seq($in);
}

fastq file looks like this:

@1477:2:1:1143:901/1
NTCGGTACAGCGACAAACAGACGATATCACCGGCTAAGCTCGATGGTGGTTACGGATGCGAAACAACGTGGTAGCTCAGGTAAGGATTTAAGGCCTTCTATTACTTTGGTTAATGAAGGCCGTGAACCAATTTGTGTGCCTGGACTCAATA
+
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
@1477:2:1:1143:901/2
ATTAACCACCGCACCTGCAGGCATTACATAATGCACCGCGATATTGGTTCCAGCCACCCAAATTGGNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNTNNNNNNNNNNNNNNNNNGCTTCCCATCGNTAACCACCATC
+
\cf_ff[fdcfebcad^e\YadcYdceWe\^deaed_Y[c\dd[ce^bfbdbRLWa]R^BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
@1477:2:1:1166:923/1
NTACTCCAGCGGAAAATGCTACGCTTCGATCATTGCTAATATCAAATAACGTTTTTTGCTCAACCGATGAGCTTTCCAGTCGGTAAGGAAGCGGTTCATTAGCCTGAGCGAGCGGGTCAAAAACGATATCTTCGCGAGCTTCATACTTAAC
+
BKOHONNNMN_______________BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB

I would really appreciate if some one can help me out.

Thanks
Shalabh



More information about the Bioperl-l mailing list