SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
HISAT2 reads aligned concordantly exactly 1 time Beyounique Bioinformatics 2 05-18-2017 03:03 PM
Bowtie2 for creating output(fastq) that are not aligned with reference zillur Bioinformatics 4 09-23-2015 12:23 AM
aligned seq > 1 times and skipping reads because of seed mismatches carolW Bioinformatics 12 05-22-2014 11:48 PM
Bowtie2 - Concordantly tags. Coryza Bioinformatics 16 02-13-2014 06:00 AM
Do I need remome the reads aligned >1 times and how? wmseq Bioinformatics 2 10-30-2013 06:00 PM

Reply
 
Thread Tools
Old 07-13-2017, 04:04 AM   #1
ge2sasag
Member
 
Location: Stockholm

Join Date: Apr 2013
Posts: 41
Default >90% aligned concordantly 0 times ChIP-seq Bowtie2

Hi,

I know this question have raised many times here and in other forums but I've tried everything other people suggested in previous posts and still can't figure it out where is the problem...... This is the output summary of Bowtie2 (I checked several times the target genome, hg19 and it's correct):

21404130 reads; of these:
21404130 (100.00%) were paired; of these:

21196512 (99.03%) aligned concordantly 0 times
104527 (0.49%) aligned concordantly exactly 1 time
103091 (0.48%) aligned concordantly >1 times
----
21196512 pairs aligned 0 times concordantly or discordantly; of these: 42393024 mates make up the pairs; of these:

906397 (2.14%) aligned 0 times
28296440 (66.75%) aligned exactly 1 time
13190187 (31.11%) aligned >1 times
97.88% overall alignment rate

I have paired-end ChIP-seq data, 50 bp reads. These are my steps (in Galaxy):

1) I groomed the fastq files to get fastqsanger (I checked if it's correct: Input FASTQ quality scores type --> Sanger & Illumina 1.8+)

2) FastQC is ok for all samples, some adapter contamination.

3) I trimmed using either TrimGalore or Trimmomatic (I get the same result using both) using forward and reverse files for each sample. I checked and I didn't missed up the forward and reverse samples (I run Bowtie2 using first file 1 and file 2 and also using first file 2 and second file 1 and I get the same result...)

4) I mapped the trimmed files to hg19 genome using Bowtie2 with option -fr and I get the result shown above. No matter what I try, I always get the same with the different samples that I have. All of them have around 32% of reads with adapters which I trimmed. I have no idea what to do with this.

Any suggestions or idea of where is the mistake? Am I missing something?

Thank you very much in advance.

Gema
ge2sasag is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:51 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO