(11 intermediate revisions by the same user not shown)
Line 1:
Line 1:
Creosote genome sequencing and assembly notes:
==Twig2Genome Notes==
[[Twig2Genome]]
*Sample obtained from front yard of 4951 W. McElroy Dr.
==Assembly==
*Sequences obtained from one lane of Illumina HiSeq2000
[[Creosote Assembly]]
*Fastq files delivered from UAGC
**82 files
**Headers are Sanger format (code 33)
***Description of Fastq file format with notes on specific decoding of header names generated by various technologies: http://en.wikipedia.org/wiki/FASTQ_format
**Pairend reads
***lane3_NoIndex_L003_R1_041.fastq
***lane3_NoIndex_L003_R2_041.fastq
**Need to get adapter sequences used in sequencing
***TGACCA (Not present in sequence reads)
*Check quality with fastqc: http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/
**[[Creosote First Run FastQC]]
'''Steps:'''
==Loading into CoGe==
*Merge R1 files; merge R2 files
SoapDeNovo assemly: 1,570,116 contigs (370MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12183 ('''Note:''' too many contigs to process and visualize in SynMap)
*gzip them
'''Trim sequences'''
*Get this package of Haibaos:
** git clone git://github.com/tanghaibao/jcvi.git
** SET PATH: export PYTHONPATH=/lib/python (which is the dir above jcvi)
** may need to install biopython: sudo easy_install biopython
*Run this: python -m jcvi.apps.baseclean trim R1.fastq.gz R2.fastq.gz
**Automatically trims and cleans sequences, also does the conversion to appropriate fastq format
**NOTE: This program should download trimmomatic, but may need to update the path of the timmomatic program in the program
*If the Trimmer script fails for silly reasons, you can run it from the command-line:
*Syntenic Path Assembly to peach: http://genomevolution.org/r/3xpd
==[[Syntenic Path Assembly]]==
'''Other Stuff'''
[[File:Master 12185 8400.genomic-CDS.lastz.dag.go c20 D20 g10 A2.aligncoords.gcoords ct0.w1000.ass2.cs1.csoS.nsd.png|thumb|600px|left|Syntenic path assembly with SynMap of creosote (x-axis) and peach (y-axis). Results may be regenerated at: http://genomevolution.org/r/3w95]]
'''Trimming reads'''
*Trim Paired ends with Trimmomatic: http://www.usadellab.org/cms/index.php?page=trimmomatic
*Assumes Illumina Encoding (code: 64, not code: 33)
**Need to convert for the HighSeq Reads:
** easy_install biopython
** git clone git://github.com/tanghaibao/jcvi.git
** export PYTHONPATH=/lib/python (which is the dir above jcvi)
[[File:Pseudo-assembly-creosote.png|thumb|600px|left|Pseudo-assembly of creosote (x-axis) using the peach (y-axis) genome. Syntenic comparison to the peach genome.]]
'''Converting sequences'''
[[File:Screen Shot 2014-05-14 at 8.33.27 AM.png|thumb|600px|left|Microsynteny analysis of the pseudo-assembled creosote genome to the peach genome. Orange bars are unsequenced Ns that represent contigs glued together in creosote by the [[syntenic path assembly]] method. Note the concordance of gene model coding sequences.]]
*python -m jcvi.formats.fastq convert (read help file, default conversion Sanger (code 33) to Illumina (code 64)
'''Other programs to clean sequences'''
*python -m jcvi.apps.baseclean trim fastqfile (single ended)
*python -m jcvi.apps.baseclean trim R1.fastq.gz R2.fastq.gz (paired ended)
'''keep sequences in single files (or two files for a pair of reads)'''
Syntenic path assembly with SynMap of creosote (x-axis) and peach (y-axis). Results may be regenerated at: http://genomevolution.org/r/3w95Syntenic path assembly with SynMap of creosote (y-axis) and Arabidopsis thaliana Col-0 (x-axis). Results may be regenerated at: http://genomevolution.org/r/3w96
Pseudo-Assembly
Pseudo-assembly of creosote (x-axis) using the peach (y-axis) genome. Syntenic comparison to the peach genome.Microsynteny analysis of the pseudo-assembled creosote genome to the peach genome. Orange bars are unsequenced Ns that represent contigs glued together in creosote by the syntenic path assembly method. Note the concordance of gene model coding sequences.
