Does bowtie need to have reads in the standard sanger format or can it accept the default file created from the 1.4 illumina pipeline in which the quals are not standard sanger?
Cheers
L
Cheers
L
You are currently viewing the SEQanswers forums as a guest, which limits your access. Click here to register now, and join the discussion
man ascii
#!/usr/bin/env perl ################################################ # prints a table with phred, ASCII, phred+33, phred+64, p ################################################ use strict; use warnings; my @phreds = (0..62); my $step = 2; printf "%6s %6s %6s %6s %10s\n", 'phred', 'ASCII', 'Ill33', 'Ill64', 'p'; for (my $i = 0; $i < @phreds; $i+=$step ){ my $phred = $phreds[$i]; printf "%6d %6d %6s %6s %10f\n", $phred, $phred+64, chr($phred+33), chr($phred+64), phred2p($phred); } sub phred2p{ return 10 ** (-(shift) / 10.0 ); }
--phred64-quals input quals are Phred+64 (same as --solexa1.3-quals)
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 05-02-2024, 08:06 AM
|
0 responses
16 views
0 likes
|
Last Post
by seqadmin
05-02-2024, 08:06 AM
|
||
Started by seqadmin, 04-30-2024, 12:17 PM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
04-30-2024, 12:17 PM
|
||
Started by seqadmin, 04-29-2024, 10:49 AM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
04-29-2024, 10:49 AM
|
||
Started by seqadmin, 04-25-2024, 11:49 AM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
04-25-2024, 11:49 AM
|
Comment