Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
alignment of bisulfite treated reads fadista Bioinformatics 40 04-03-2012 12:32 PM
MethylCoder: software for bisulfite treated reads brentp Bioinformatics 15 07-07-2011 12:17 PM
Illumina mapping with bwa bpetersen Bioinformatics 5 03-14-2011 11:04 PM
Question about Illumina reads, SNPs and mapping assemblies allcreation Bioinformatics 8 12-09-2010 11:53 AM
increase in 3' %T after filtering BS-Treated reads brentp Bioinformatics 2 09-29-2010 06:10 AM

Thread Tools
Old 12-06-2011, 04:02 AM   #1
Location: Pittsburgh

Join Date: Aug 2011
Posts: 72
Default help mapping BS treated illumina reads using rmapbs

I recently ran some WGBS treated DNA on a hiseq2000 and I am trying to align it using the rmapbs software that I saw people using in the forums. Warning:I am a benchwork scientist with not much programming experience, but I am trying to learn. My question is about specifying the reference chromosome directories and the BS reads that I sequenced.

1. For specifying the chromosome directory (-c), can I point it to a folder containing all my chromosome.fa files?

2. To specify what reads I want it to analyze, can I also point it to a folder containing fastq files or does the path have to end in a specific .fastq file. Just as an FYI, illumina puts reads in separate files (I guess for data management purposes) automatically, so I have like 266 .fastq gzipped files. This might also be a problem.

Basically I am getting an error and I think these might be the source but I am not sure. I will paste the code I was running below. I probably did something stupid. I have read the manual but haven't figured out what the problem is.

/usr/local/bin/rmap/rmapbs -o /111123_SN874_0071_AD0F4CACXX/22647_aligned -c /UCSC/hg19/Sequence/Chromosomes/ /111123_SN874_0071_AD0F4CACXX/unaligned/Project_kinome_WGBS/Sample_22647_BS/read_1/

Any advice would be greatly appreciated!
shawpa is offline   Reply With Quote

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 12:22 AM.

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