Header Leaderboard Ad

Collapse

What type of phred quality scores does samtools 1.3 expect

Collapse

Announcement

Collapse

SEQanswers June Challenge Has Begun!

The competition has begun! We're giving away a $50 Amazon gift card to the member who answers the most questions on our site during the month. We want to encourage our community members to share their knowledge and help each other out by answering questions related to sequencing technologies, genomics, and bioinformatics. The competition is open to all members of the site, and the winner will be announced at the beginning of July. Best of luck!

For a list of the official rules, visit (https://www.seqanswers.com/forum/sit...wledge-and-win)
See more
See less
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • What type of phred quality scores does samtools 1.3 expect

    I am working with NG sequence data that was generated in 2012 using the Illumina HiSeq 1000. I have had some trouble getting samtools to accurately identify heterozygotes and phase bam files generated from these reads. I am wondering if said problems are a result of samtools misreading the phred scores in the fastq/bam files. Does anyone know the default method of phred-score encoding that samtools 1.3 expects?

    e.g.
    Sanger/Illumina 1.8 (ASCII 33 to 126)
    Solexa/Illumina 1.0 (ASCII 59 to 126)
    Illumina 1.3 (ASCII 64 to 126)

  • #2
    All SAM/BAM/CRAM files expect Sanger/Illumina 1.8 format (aka, "phred + 33"). If you aligned phred+64 fastq files and didn't tell your aligner to fix the phred scores then you're going to have problems. I happen to have written a smaller conversion program if you happen to have done this (https://github.com/dpryan79/Answers/...iostars_133825).

    Comment

    Latest Articles

    Collapse

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Today, 07:14 AM
    0 responses
    4 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 01:08 PM
    0 responses
    6 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-01-2023, 08:56 PM
    0 responses
    120 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 06-01-2023, 07:33 AM
    0 responses
    257 views
    0 likes
    Last Post seqadmin  
    Working...
    X