error with bowtie2 alignment
0
0
Entering edit mode
3.8 years ago
pt.taklifi ▴ 60

Hi I am trying to align a fastq file to mitochondrial genome with bowtie 2 with the following command and store unmapped read to folder unaligned

 bowtie2 -p 1 -k 1 -D 20 -R 3 -N 1 -L 20 -i S,1,0.50 -x /home/ptaklifi/genome_folder/rCRSd/bowtie2_index/default/rCRSd --rg-id SRR10984460 -U SRR10984460_1-trimmed-pair1.fastq --un unaligned

but during alignment I get this error

Error: Read SRR10984461.33749887 33749887 length=76 has more read characters th$
terminate called after throwing an instance of 'int'
Aborted (core dumped)
(ERR): bowtie2-align exited with value 134

using the command below I looked through the read

 grep SRR10984461.33749886 -A 3 SRR10984461_1-trimmed-pair1.fastq

and here is the output :

ATTCAAACAAGTACAATTAGAAATTATAAAGGTGACATTAACAACTCATATCACAGAAATACAAAAGATCATCAGA
+
A

so first should I delete this read ? and if so, how to do that ? I appreciate your response

alignment error sequencing • 1.1k views
ADD COMMENT
1
Entering edit mode

Use cutadapt to filter out low-quality bases/reads.

ADD REPLY
1
Entering edit mode

Soemthing got corrupted during the initial trimming would be my guess. How did you do the trimming?

ADD REPLY
0
Entering edit mode
skewer -f sanger -m pe -x Nextera.fa read1.fastq read2.fastq
ADD REPLY
0
Entering edit mode

That looks about right. Can you check whether this is already in the original fastq file you downloaded? You could run repair.sh from BBMap on the original file, then trim and proceed if the original file was the problem.

ADD REPLY
0
Entering edit mode

I downloaded the fastq files from SRA . I will check if the problem is with fastq files. Thank you

ADD REPLY
0
Entering edit mode

thank you @ATpoint . as suggested , I think I should remove poor quality reads in trimming step. can I remove them with skewer ?

ADD REPLY

Login before adding your answer.

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