Also include the read ID in the fastq header
authorDiane Trout <diane@caltech.edu>
Tue, 1 Dec 2009 02:03:04 +0000 (02:03 +0000)
committerDiane Trout <diane@caltech.edu>
Tue, 1 Dec 2009 02:03:04 +0000 (02:03 +0000)
scripts/qseq2fastq

index 56589ee0c7d00322f03455aa5eb3de2df64c19f4..42628980d4f8b78809791ebd9bff1df92a3280c6 100644 (file)
@@ -56,13 +56,14 @@ def qseq2fastq(destination, qseqs, trim=None):
             
             pass_qc = record[10]
 
-            destination.write('@%s_%s:%s:%s:%s:%s pf=%s%s' % ( \
+            destination.write('@%s_%s:%s:%s:%s:%s/%s pf=%s%s' % ( \
                 machine_name,
                 run_number,
                 lane_number,
                 tile,
                 x,
                 y,
+                read,
                 pass_qc,
                 os.linesep))
             destination.write(sequence[trim])