SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   How to convert .txt file to .bed or .gff, How can we use chip seq data in R software (http://seqanswers.com/forums/showthread.php?t=1096)

forevermark4 01-23-2009 06:11 PM

How to convert .txt file to .bed or .gff, How can we use chip seq data in R software
 
Hi Everyone,

Myself Yogesh kumar and I am new in illumina solexa work for chip seq analysis (cisgenome software).. I need only help to conert this dataset in format which supported by UCSC browser or cisgenome software..Other analysis part I can do.. I need your help I attached my datasets links

http://www.well.ox.ac.uk/htseq/qXZxn...s_1_export.txt

http://www.well.ox.ac.uk/htseq/qXZxn...s_3_export.txt

About these datasets some information

The “s_#” suffix stands for the lane number used in the flow cell.
The “export” format is defined by multiple, tab separated, columns.
For the definition of each column see below.


EXPORT file definitions

Export files are generated by the GERALD step* of the Illumina pipeline. The program called ELAND aligns each read from each lane to the reference genome (in our case the Mus Musculus genome, UCSC release MM8).
When a match is found, then the relative position on the genome is reported on the “s_#_export.txt” file (only chromosome, start position and strand F or R).
If no match is found a “NM” (no match) is reported.
There is a line for each read, whether it aligns or not, and multiple lines for the same read if it aligns in multiple positions.
You can either parse yourself these files (R-Bioconductor, Perl script …) or use public available software (eg CisGenome).
The file contains information on the physical position of the reads in the flow cell, the nucleotide sequence of the read itself, a string of the quality call for each nucleotide in the read (a code developed by Illumina), various flags and the genomic position (if found).
For many purposes not all these information are needed.

*the illumina pipeline consists in several steps, starting from cluster recognition, passing through basecalling and ending with the alignment of the bases to the reference genome.
Not all fields are relevant to a single-read analysis.
1. Machine (Parsed from Run Folder name)
2. Run Number (Parsed from Run Folder name)
3. Lane
4. Tile
5. X Coordinate of cluster
6. Y Coordinate of cluster
7. Index string (Bland for a non-indexed run)
8. Read number (1 or 2 for paired-read analysis, blank for a single-read analysis)
9. Read
10. Quality string—In symbolic ASCII format (ASCII character code = quality value + 64) by default (Set QUALITY_FORMAT --numeric in theGERALD config file for numeric values)
11. Match chromosome—Name of chromosome match OR code indicating why no
match resulted
12. Match Contig—Gives the contig name if there is a match and the match
chromosome is split into contigs (Blank if no match found)
13. Match Position—Always with respect to forward strand, numbering starts at 1 (Blank if no match found)
14. Match Strand—“F” for forward, “R” for reverse (Blank if no match found)
15. Match Descriptor—Concise description of alignment (Blank if no match found)
• A numeral denotes a run of matching bases
• A letter denotes substitution of a nucleotide:
For a 35 base read, “35” denotes an exact match and “32C2” denotes substitution
of a “C” at the 33rd position
16. Single-Read Alignment Score—Alignment score of a single-read match, or for a paired read, alignment score of a read if it were treated as a single read (Blank if no match found)
17. Paired-Read Alignment Score—Alignment score of a paired read and its partner, taken as a pair (Blank for single-read analysis)
18. Partner Chromosome—Name of the chromosome if the read is paired and its partner aligns to another chromosome (Blank for single-read analysis)
19. Partner Contig—Not blank if read is paired and its partner aligns to another
chromosome and that partner is split into contigs (Blank for single-read analysis)
20. Partner Offset—If a partner of a paired read aligns to the same chromosome and contig, this number, added to the Match Position, gives the alignment position of the partner (Blank for single-read analysis)
21. Partner Strand—To which strand did the partner of the paired read align? “F” for forward, “R” for reverse (Blank if no match found, blank for single-read analysis)
22. Filtering—Did the read pass quality filtering? “Y” for yes, “N” for no

Are you have any idea How can we convert these files to WIG, BED and GFF for the UCSC. Any one format is sufficient. for me . otherwise how can we convert .txt file to .BED file.. I am planning to use cisgenome (two sample analysis) software

to look at data mapped on their genomic original contest.

It 'll be great favour for me

Thanks
Yogesh Kumar

ECO 01-23-2009 09:52 PM

Hi Yogesh,

Your URL's are badly formed.

forevermark4 01-23-2009 11:06 PM

Hi Everyone,

Myself Yogesh kumar and I am new in illumina solexa work for chip seq analysis (cisgenome software).. I need only help to conert this dataset in format which supported by UCSC browser or cisgenome software..Other analysis part I can do.. I need your help I attached my datasets links

http://www.well.ox.ac.uk/htseq/qXZxnawHREhej1c2lFut

http://www.well.ox.ac.uk/htseq/qXZxn...s_1_export.txt

http://www.well.ox.ac.uk/htseq/qXZxn...s_3_export.txt

About these datasets some information

The “s_#” suffix stands for the lane number used in the flow cell.
The “export” format is defined by multiple, tab separated, columns.
For the definition of each column see below.


EXPORT file definitions

Export files are generated by the GERALD step* of the Illumina pipeline. The program called ELAND aligns each read from each lane to the reference genome (in our case the Mus Musculus genome, UCSC release MM8).
When a match is found, then the relative position on the genome is reported on the “s_#_export.txt” file (only chromosome, start position and strand F or R).
If no match is found a “NM” (no match) is reported.
There is a line for each read, whether it aligns or not, and multiple lines for the same read if it aligns in multiple positions.
You can either parse yourself these files (R-Bioconductor, Perl script …) or use public available software (eg CisGenome).
The file contains information on the physical position of the reads in the flow cell, the nucleotide sequence of the read itself, a string of the quality call for each nucleotide in the read (a code developed by Illumina), various flags and the genomic position (if found).
For many purposes not all these information are needed.

*the illumina pipeline consists in several steps, starting from cluster recognition, passing through basecalling and ending with the alignment of the bases to the reference genome.
Not all fields are relevant to a single-read analysis.
1. Machine (Parsed from Run Folder name)
2. Run Number (Parsed from Run Folder name)
3. Lane
4. Tile
5. X Coordinate of cluster
6. Y Coordinate of cluster
7. Index string (Bland for a non-indexed run)
8. Read number (1 or 2 for paired-read analysis, blank for a single-read analysis)
9. Read
10. Quality string—In symbolic ASCII format (ASCII character code = quality value + 64) by default (Set QUALITY_FORMAT --numeric in theGERALD config file for numeric values)
11. Match chromosome—Name of chromosome match OR code indicating why no
match resulted
12. Match Contig—Gives the contig name if there is a match and the match
chromosome is split into contigs (Blank if no match found)
13. Match Position—Always with respect to forward strand, numbering starts at 1 (Blank if no match found)
14. Match Strand—“F” for forward, “R” for reverse (Blank if no match found)
15. Match Descriptor—Concise description of alignment (Blank if no match found)
• A numeral denotes a run of matching bases
• A letter denotes substitution of a nucleotide:
For a 35 base read, “35” denotes an exact match and “32C2” denotes substitution
of a “C” at the 33rd position
16. Single-Read Alignment Score—Alignment score of a single-read match, or for a paired read, alignment score of a read if it were treated as a single read (Blank if no match found)
17. Paired-Read Alignment Score—Alignment score of a paired read and its partner, taken as a pair (Blank for single-read analysis)
18. Partner Chromosome—Name of the chromosome if the read is paired and its partner aligns to another chromosome (Blank for single-read analysis)
19. Partner Contig—Not blank if read is paired and its partner aligns to another
chromosome and that partner is split into contigs (Blank for single-read analysis)
20. Partner Offset—If a partner of a paired read aligns to the same chromosome and contig, this number, added to the Match Position, gives the alignment position of the partner (Blank for single-read analysis)
21. Partner Strand—To which strand did the partner of the paired read align? “F” for forward, “R” for reverse (Blank if no match found, blank for single-read analysis)
22. Filtering—Did the read pass quality filtering? “Y” for yes, “N” for no

Are you have any idea How can we convert these files to WIG, BED and GFF for the UCSC. Any one format is sufficient. for me . otherwise how can we convert .txt file to .BED file.. I am planning to use cisgenome (two sample analysis) software

to look at data mapped on their genomic original contest.

It 'll be great favour for me

Thanks
Yogesh Kumar

ECO 01-24-2009 07:29 AM

OK Yogesh, you don't need to repost the same message. Just click "edit" and fix the URLs.

graveley 01-24-2009 10:23 AM

Dear Yogesh,

We do this by writing a perl script that reads in the alignment information and writes a new file in the appropriate format. I would send you what we use, but we do not use export.txt files. We are currently doing alignments with Bowtie and then converting the output to .gff and .wig files.

Brent

forevermark4 01-26-2009 01:02 AM

Hi Brent,

Thanks.. If you dont mind can you send me that perl script .. So I can try here to txt fle and to convert in .gff or .wif format .. Perl script source code which you are using to convert .txt file to .gff or .wig format or alignment script because I already know perl script how to convert .fasta to .embl or other formats

Agent47 01-27-2009 02:53 PM

Hi Brent,
I am stuck in almost same position as Yogesh.
I am using MAQ for the alignment of SOLEXA data but i am not able to convert it into .WIG and .GFF format, if you can provide me some directions for this it would be a great help
Thanks!

Arpit

tabascoj 01-30-2009 02:31 PM

Bowtie to .wig
 
Brent,
I would really appreciate any perl suggestions for getting the Bowtie alignment into a WIG file (which I intend to use in Gbrowse). I have no problem with the perl conversion of tabulated documents (e.g. Bowtie-->GFF), but I need help getting a pileup and getting the values into the WIG file.

Thanks very much.
Joe

apfejes 01-30-2009 02:40 PM

I feel silly promoting my own software, but Maq to wig and eland to wig are both handled well by FindPeaks.

http://vancouvershortr.wiki.sourcefo...indPeaksManual

You may not need the ChIP-Seq features, but you can certainly just use it for a quick conversion. (There are converters in the package for creating bed files as well)

As for bowtie, you can always have it produce a .map file and then do the same conversion.

Good luck.

jperin 02-05-2009 06:22 AM

this may be the wrong place to ask, but I've just tried findPeaks for creating our wig files and it works great. The only problem is that the wig file appears to be offset to the very beginning of the chromosome. Our reference sequence is only a small piece of chromosome 10, in this case. It appears that at some stage in performing the maq alignment, a tag isn't set properly and causes the wig file to insert a label of "hg18_dna" instead of the "chr10", and then the start position for the first base pair in the first header starts at 1, instead of 17M something... where it 'should' be.

The fasta reference file has the correct tag in it, with the right reference, but at some stage this doesn't get passed to findPeaks and the offset is not correctly inserted. I can see the beautiful wig image in a browser, but it displays at the beginning of chr10. I had to manually change the hg18_dna to chr10 for the first to work, but changing the offset isn't as simple since the multiple headers each have their own offset position and it would be hard to calculate, also assuming there's probably a simple way to fix this??

Thanks for any advice.
Juan

mudshark 02-06-2009 05:07 AM

hi

i tested several published ChipSeq applications. therefore i would like to mention the spp package for R (http://compbio.med.harvard.edu/Supplements/ChIP-seq/) which was the only piece of software that directly produced kind of convincing output on 2 sample comparison data (i also tested PeakSeq, MACS, cisgenome).

apfejes 02-06-2009 09:41 AM

I'm adding several control modes to FindPeaks this week. I hope you'll revisit that list at some point (=

Anthony

mudshark 02-07-2009 01:54 AM

hi Anthony,
afaik FindPeaks does not (yet?) support 2-sample analysis, i.e. IP vs. Input. is that correct?
T.

apfejes 02-07-2009 04:50 AM

That's correct - The code is currently in development, but should be ready shortly. (I hate saying that about software, but a lot of the code has already been written.)

Cheers,
Anthony

alperyilmaz 02-11-2009 10:00 AM

R package from BioConductor
 
There's another R package to be released, which will be available thru BioConductor. It's mentioned in a recent workshop. The workshop material can be viewed here.

Layla 06-26-2009 01:17 AM

maq to GFF2 format
 
Hi,

I know this is quite an old thread but useful for me at this moment.
Searching for a method to convert a maq alignment (to human ref genome) output from me-dip seq to GFF2 format for batman but no luck so far. Is there a script already available for this or is there a conversion tool to do this?

GFF2 example:
col1 col2 col3 col4 col5 col6 col7 col8 col9
Chr1 Homo_build36 Reference 1 45 . + 1 Sequence1

How does one define column 6 (score) and column 7 (phase) from the maq alignment file?

Help appreciated

L

kmcarr 06-26-2009 04:12 AM

Layla,

For the score I would use the mapping quality (column 7 of the MAQ .aln file). Phase is meaningless in this context, it only applies to coding sequence features. In your case just put a '.' in column 7.

joseph 06-28-2009 08:27 AM

Quote:

Originally Posted by graveley (Post 3128)
Dear Yogesh,

We do this by writing a perl script that reads in the alignment information and writes a new file in the appropriate format. I would send you what we use, but we do not use export.txt files. We are currently doing alignments with Bowtie and then converting the output to .gff and .wig files.

Brent

Hi Brent,
I would really appreciate it if you send me your perl scripts for conversion of bowtie output to .gff and .wig files.
my e-mail: jdsandjd@yahoo.com
Thanks
Joseph

Ka123$ 07-16-2009 04:48 PM

Bowtie to Bed files
 
Hi Brent
could you also help me with the perl script to convert bowtie to BED? I am not a perl person. Also can find peak convert the bowtie-align files to a cisGenome compatible file?
THanks
Ka

Ka123$ 07-16-2009 04:50 PM

HI Brent,
My email is com_bio@yahoo.com.Thanks
What kind of format are (GERALD) aligned output files?Aligned files are mapped files I suppose......just two terminologies.....?


All times are GMT -8. The time now is 06:10 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.