[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