paired-end files with reads on different strand
1
0
Entering edit mode
8.8 years ago

I want to align reads in 2 different files using paired-end alignment with the help of bowtie. when I try to align it with bowtie1 it fails and a little portion of it just align to the reference genome. I aligned each file separately as a single read and when checked 2 SAM output files I saw that the flag values for entries in first file and second file are 0 for some entries and 16 for others,which I think means that some reads for each file were on the + strand and others on - strand also when flag value for entry number N is 0 for first file it's corresponding flag value for N'th entry in second file is 16 . I want to know what should to do if I want to align it bowtie1? Is it true if I reverse all entries in first file and make all reads of first file on the same strand and do the same thing for second strand and then align it?

The other question is that bowtie2 can map these two files, how can I be sure that mapping is true?

alignment sequencing bowtie paired-end • 2.2k views
ADD COMMENT
3
Entering edit mode
8.8 years ago

Don't reverse your reads, paired end reads are supposed to be on opposite strands. Bowtie is smart enough to know this. The first thing I'd try is changing the settings to tell it to expect a very large potential insert size between the reads; I think bowtie often refuses to report alignments if it thinks the reads are too far apart, and the default distance is very small.

ADD COMMENT
0
Entering edit mode

Thanks, The problem solved :)

ADD REPLY
2
Entering edit mode

If your problem got solved, please accept the answer by clicking on the green tick, which only you could see below the answer.

ADD REPLY

Login before adding your answer.

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