SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Confusion about single end vs paired end read output on Illumina HiSeq acidcoated Illumina/Solexa 1 03-19-2015 08:56 AM
Identify differentially expressed genes from single-end Vs. paired end samples DAnand Bioinformatics 2 11-04-2014 11:35 AM
Given BAM/SAM file, how to see if it's single-end or paired-end sequencing? xxatbio Bioinformatics 2 08-11-2013 02:51 AM
Forcing paired-end data mapped as single-end in SAM puggie Bioinformatics 1 03-16-2013 10:50 AM
Paired-end Bam from single-end aligned sam ramouz87 Bioinformatics 4 08-17-2011 12:55 PM

Reply
 
Thread Tools
Old 03-01-2018, 12:48 AM   #1
jennorocks
Junior Member
 
Location: UK

Join Date: Oct 2017
Posts: 3
Default Paired end SAM to single end

I am performing PAR-CLIP which is typically done using single end however we ran paired because it fit in with our in house set up.

I've mapped the reads (as a pair) with bowtie and generated a SAM file that has the pairs alternating.

However, it seems the program I am using to peak call (PARalyzer) doesn't seem to work with the paired file. I get barely any clusters whereas if I just map the forward reads and use those I get 1000s.

Since these are short reads anyway with almost complete overlap the obvious approach is to just work with R1 but is there any way to combine them so I dont lose the mapping info from having R1 and R2

Hope that makes sense
jennorocks is offline   Reply With Quote
Old 03-01-2018, 12:53 AM   #2
jennorocks
Junior Member
 
Location: UK

Join Date: Oct 2017
Posts: 3
Default

Just for an example, a pair of reads might looks like this in my SAM

NTTCTGTAAGCCTGTGAAGGTGTGCTGTGAGGCAT NGCCTCACAGCACACCTTCACAGGCTTACAGAAC

So between them theres good overlap but they both add 1 to 2 bases so only using 1 of them loses some information
jennorocks is offline   Reply With Quote
Reply

Tags
clip, paired end read, par-clip, sam, sam file

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 08:59 PM.


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