Mapping Phred+64 encoded reads with Bowtie2
0
0
Entering edit mode
3.0 years ago
Denis ▴ 210

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?

genome alignment software error • 673 views
ADD COMMENT
1
Entering edit mode

You should add --phred64 which indicates that qualities are Phred+64.

ADD REPLY
0
Entering edit mode

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 resulted SAM file i usually got)?

ADD REPLY
1
Entering edit mode

If I had to guess, the mapping results might be skewed as bowtie2 might (notice I said might) think that you have just really good higher quality sequences than you actually have. Not sure if bowtie2 can automatically determine the encoding of the reads or not.

ADD REPLY

Login before adding your answer.

Traffic: 865 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6