I suggest you have a look to export2sam.pl script.
In this file you could find a section "quality value conversion table" explaining quality format start point and end point.
So you could know which quality format is used in result sam file, i.e. phred+33
Hope this help.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Base qualities during qseq2fastq and export2sam conversion
I started off with Illumina qseq format - converted to fastq.
I used the sol2std conversion from Maq and used the fastq to align against reference using BWA.
Used the original qseq files to align against reference using ELAND.
I got the output of BWA in SAM format and ELAND in export (which I converted to SAM using export2sam.pl)
Now, the purpose is to detect SNPs and indels using samtools which expects (or interprets) the base qualities in a particular phred scale (Phred+33 correct?).
Could someone help me sort the out the base quality scales in these file formats and what is the correct way to format the input SAM files before running pileup?
Same question for GATK unified caller - what kind of base qualities should I ensure is in the input SAM files?
If I use a standard converter such as sol2std or export2sam how can I be sure the base qualities are in the correct scale?
Thanks in advance.Last edited by hrajasim; 05-12-2011, 11:18 AM.Tags: None
Latest Articles
Collapse
-
by seqadmin
Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.
3D Genomics
While spatial biology often involves studying proteins and RNAs in their...-
Channel: Articles
01-01-2025, 07:30 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 07:35 AM
|
0 responses
20 views
0 likes
|
Last Post
by seqadmin
Yesterday, 07:35 AM
|
||
Started by seqadmin, 01-23-2025, 09:43 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
01-23-2025, 09:43 AM
|
||
Started by seqadmin, 01-23-2025, 08:36 AM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
01-23-2025, 08:36 AM
|
||
Started by seqadmin, 01-17-2025, 09:38 AM
|
0 responses
37 views
0 likes
|
Last Post
by seqadmin
01-17-2025, 09:38 AM
|
Leave a comment: