What does Illumina raw data look like?

[alig]

Hi,

Re : There's a really good tool in the maq package (latest release) called fq_all2std.pl

I tried to use

fq2fa Convert various FASTQ-like format to FASTA

to convert my illumina seq data from fastq to fasta as I want the quality in fasta format to run Mosaik's gigBayes program.

But the Maq perl script fq_all2std.pl fq2fa <in.txt> command

just seemed to print the results to the screen & not place them in a fasta file.

Am I doing something really silly here?

Only I've got a 1.8 Gb illumina seq text file so this process takes a while & I need it in a file, not printed to the screen

thanks alig

[hannat]

I have similar data, seq.txt and prb just like you,
seq.txt
........................................................................
6 1 914 893 GCTACTGCCGTGACCTCATTTCTCTTA
6 1 898 905 GAAAAAGAGAAAGTTTAGGAGATCGAT
.....................................................................................
prob.txt
.....................................................................................
-30 -30 30 -30 -30 30 -30 -30 -30 -30 -30 30 30 -30 -30
-30 -30 30 -30 -30 -30 -30 -30 30 -30 -30 30 -30 -30 3
0 -30 -30 -30 30 -30 -30 -30 -30 30 -30 -30 -30 -30 30
-30 -30 30 -30 30 -30 -30 -30 -30 30 -30 -30 -30 30 -30
-30 -30 -30 -30...
.............................................................................

but when i run
fq_all2std.pl seqprb2std seq.txt prb.txt
The output is like following,
...........................................
@6:1:914:893
GCTACTGCCGTGACCTCATTTCTCTTA
+
???????????????????????????
..................................................

And i had lots of the warnings, similar things like this

but there is other problems, i got lots of warning message like this:
Argument "\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0..." isn't numeric in numeric gt (>) at /usr/local/bin/fq_all2std.pl line 152, <$fhq> line 6609.
....................................................................................................................................................................................................... line ......


i wonder if this kind of warning is happening to others too, if so, what do you think the problem is?
now i am checking my prb.txt, i guess there is some lines which was not accpeted.

[kwebb]

Great tool!

Thanks for the info!

[zee]

There's a really good tool in the maq package (latest release) called fq_all2std.pl

See below:


Usage: fq_all2std.pl <command> <in.txt>

Command: scarf2std Convert SCARF format to the standard/Sanger FASTQ
fqint2std Convert FASTQ-int format to the standard/Sanger FASTQ
sol2std Convert Solexa/Illumina FASTQ to the standard FASTQ
fa2std Convert FASTA to the standard FASTQ
seqprb2std Convert .seq and .prb files to the standard FASTQ
fq2fa Convert various FASTQ-like format to FASTA
export2sol Convert Solexa export format to Solexa FASTQ
export2std Convert Solexa export format to Sanger FASTQ
csfa2std Convert AB SOLiD read format to Sanger FASTQ
instruction Explanation to different format
example Show examples of various formats

Note: Read/quality sequences MUST be presented in one line.

[zee]

There's a really good tool in the maq package (latest release) called fq_all2std.pl

See below:


Usage: fq_all2std.pl <command> <in.txt>

Command: scarf2std Convert SCARF format to the standard/Sanger FASTQ
fqint2std Convert FASTQ-int format to the standard/Sanger FASTQ
sol2std Convert Solexa/Illumina FASTQ to the standard FASTQ
fa2std Convert FASTA to the standard FASTQ
seqprb2std Convert .seq and .prb files to the standard FASTQ
fq2fa Convert various FASTQ-like format to FASTA
export2sol Convert Solexa export format to Solexa FASTQ
export2std Convert Solexa export format to Sanger FASTQ
csfa2std Convert AB SOLiD read format to Sanger FASTQ
instruction Explanation to different format
example Show examples of various formats

Note: Read/quality sequences MUST be presented in one line.

[kwebb]

Update

I've managed to convert my data using the solexa2fasta.pl script. However the tool included with Maq, sol2sanger, does not work with my data. Can someone please explain?

Thank you!

[kwebb]

Hi

I'm trying to work through some of the various assembler programs before actually collecting my own Illumina data. I've found some test datasets here:

http://sharcgs.molgen.mpg.de/download.shtml

but I'm not sure if the file formats are the same as raw data from the Genome Analzyer.

The files are s_4_seq.txt and s_4_prb.txt and the first few lines look like this:
s_4_seq.txt
4 1 56 910 AACTTACAATTGAAAATATAAACTCAT
4 1 64 716 AAGATGATTATATGTCTTCCTTTTCGA
4 1 890 894 TCAAACCAATCAGACCTATGTTTCATA

s_4_prb.txt
40 -40 -40 -40 40 -40 -40 -40 -40 40 -40 -40 -40 -4
0 -40 40 -40 -40 -40 40 40 -40 -40 -40 -40 40 -40
-40 40 -40 -40 -40 40 -40 -40 -40 -40 -40 -40 40

So my questions are
1. Is this the raw data format from the machine?
2. How do I get these files into fastq format? The maq converter and sanger perl scripts previously mentioned do not seem to work.

Thank you!