View Single Post
Old 01-23-2018, 12:15 AM   #20
Gopo
Member
 
Location: Louisiana

Join Date: Nov 2013
Posts: 23
Default

I am not sure this is a bug or not, but when I try to use reformat.sh (version 37.76) to add fake qualities of Q30 to a PacBio Sequel fastq file (produced with
Code:
bamtools convert -format fastq -in sequel.subreads.bam -out sequel.subreads.fastq
), which has default quality of "!", I don't get quality of ">" in the output rather "#".

Code:
/opt/bbmap/reformat.sh qin=33 qout=33 qfake=30 in=sequel.subreads.fastq out=sequel.subreads.fqual.fastq
Gopo is offline   Reply With Quote