Can bowtie2 mapping paired-end reads of different length?
1
0
Entering edit mode
4.2 years ago

Hello!

I began to analyze some reads of a genome of a haploid organisms. The genome was made with HiSeq and during trimerization, the adapters were removed and this generated reads of different sizes in the paired-end files.

I used the following command:

bowtie2 -x SalI -1 121_MSO.Adtrim_1.fastq -2 121_MSO.Adtrim_2.fastq -p 20 -D 20 -R 3 -N 1 -L 20 -i S,1,0.5 -S 1_21.sam

When I try to rotate the bowtie it generates the following error:

Warning: skipping mate #2 of read 'MG00HS19:1041:CBBT1ANXX:1:1101:11867:2106 2:N:0:GAGTGG' because length (1) <= # seed mismatches (1)
Warning: skipping mate #2 of read 'MG00HS19:1041:CBBT1ANXX:1:1101:11867:2106 2:N:0:GAGTGG' because it was < 2 characters long
Warning: skipping mate #2 of read 'MG00HS19:1041:CBBT1ANXX:1:1101:14352:2179 2:N:0:GAGTGG' because length (1) <= # seed mismatches (1)
Warning: skipping mate #2 of read 'MG00HS19:1041:CBBT1ANXX:1:1101:14352:2179 2:N:0:GAGTGG' because it was < 2 characters long
Warning: skipping mate #2 of read 'MG00HS19:1041:CBBT1ANXX:1:1101:15207:2236 2:N:0:GAGTGG' because length (1) <= # seed mismatches (1)
Warning: skipping mate #2 of read 'MG00HS19:1041:CBBT1ANXX:1:1101:15207:2236 2:N

Thanks for all!

bowtie2 genome • 2.0k views
ADD COMMENT
1
Entering edit mode

It says that the reads are too short to align. When you trim the data, try to remove reads are less than 20bp long, as your -L is 20

ADD REPLY
0
Entering edit mode

Thanks! I will see that!

ADD REPLY
0
Entering edit mode
4.2 years ago
Bioaln ▴ 350

From documentation, the -L parameter dictates the lower seed length bound. This is according to warnings the problem (as noted by @geek_y)

ADD COMMENT

Login before adding your answer.

Traffic: 1969 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