Entering edit mode
5.5 years ago
Denis
▴
310
Hi,
I faced with strange behavior of Bowtie2
. I mapped reads with Phred+64
encoding by command:
bowtie2 --very-sensitive -x index -U R1.fastq -S Alignment.sam
I'm wondering why i didn't get any error message in that case? Is Bowtie2
able to determine reads encoding format automatically?
You should add
--phred64
which indicates thatqualities are Phred+64
.Thanks! But why i didn't get error by using this command line? I use this command a long time and only recently realized that the reads have
Phred+64 encoding
. How it may influence on the mapping results (everything looked well in resultedSAM
file i usually got)?If I had to guess, the mapping results might be skewed as
bowtie2
might (notice I saidmight
) think that you have just really good higher quality sequences than you actually have. Not sure ifbowtie2
can automatically determine the encoding of the reads or not.