Creosote: Difference between revisions

From CoGepedia
Jump to navigation Jump to search
← Older edit
No edit summary
No edit summary
 
(13 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
***lane3_NoIndex_L003_R1_041.fastq
***lane3_NoIndex_L003_R2_041.fastq
**Need to understand if these are paired-end reads
**Need to get adapter sequences used in sequencing
**Description of Fastq file format with notes on specific decoding of header names generated by various technologies: http://en.wikipedia.org/wiki/FASTQ_format
*Check quality with fastqc: http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/
**[[Creosote First Run FastQC]]


*Sequences cleaned using trimReads by Haibao Tang: https://github.com/tanghaibao/trimReads/tree/
==Loading into CoGe==
**NOte: Only use on single reads
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)
**Ran with supplied adapter sequence file:
>Adapter 4
TGACCA
>Adapter 4 rc
TGGTCA
**Command-line run:
Running /home/elyons/bin/trimReads  -Q 33 -f /home/elyons/projects/genome/data/creosote/Sample_lane3/adapter/adapter.faa ./lane3_NoIndex_L003_R2_015.fastq
**Output of trimReads:


SoapDeNovo assembly of contigs >= 2000nt: 6,976 contigs (20MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12185


ABySS assembly (bpsize=64, 2kb minimum contig size): 122,972 contigs (392MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12275
*Syntenic Path Assembly to peach: http://genomevolution.org/r/3xmq


''New Notes for processing"
Velvet assembly 515,190 contigs (241MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12245


*Trim Paired ends with Trimmomatic: http://www.usadellab.org/cms/index.php?page=trimmomatic
CLC4 assembly 685,475 contigs (508MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12244
*Assumes Illumina Encoding (code: 64, not code: 33)
*Syntenic Path Assembly to peach: http://genomevolution.org/r/3xpd
**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)
** python -m jcvi.formats.fastq  (Install missing packages)


'''Steps:'''
==[[Syntenic Path Assembly]]==
*Merge R1 files; merge R2 files
*gzip them
*Run this: python -m jcvi.apps.baseclean trim R1.fastq.gz R2.fastq.gz
**NOTE: This program should download trimmomatic, but may need to update the path of the timmomatic program in the program
*Note:  Bao recommends CLC for genome assembly.  Runs faster, less memory, less sensitive to bad data.  Compute intensive.
*If the Trimmer script fails for silly reasons, you can run it from the command-line:
java -Xmx4g -cp Trimmomatic-0.13/trimmomatic-0.13.jar org.usadellab.trimmomatic.TrimmomaticPE lane3_NoIndex_L003_R1_001.b64.fastq.gz lane3_NoIndex_L003_R2_001.b64.fastq.gz lane3_NoIndex_L003_R1_001.pairs.fastq.gz lane3_NoIndex_L003_R1_001.frags.fastq.gz lane3_NoIndex_L003_R2_001.pairs.fastq.gz lane3_NoIndex_L003_R2_001.frags.fastq.gz ILLUMINACLIP:adapters.fasta:2:40:15 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
*Running SOAPdenovo:
SOAPdenovo31mer all -s ../../soap.config.eric -o SoapAssem -K 25 -p 16 -R -d -D -F
**Note: if SOAP crashes, try another XXmer binary (e.g. 63mer)


'''Running Velvet'''
[[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]]
**Need to interleave reads:
~/src/velvet_1.1.04/shuffleSequences_fastq.pl lane3_NoIndex_L003_R1_001.pairs.fastq lane3_NoIndex_L003_R2_001.pairs.fastq merged_pairs.fastq
**set threading of velvet with env var
export OMP_NUM_THREADS=32
**running velveth
OMP_NUM_THREADS=32 velveth VelvetAssem 31 -shortPaired -fastq.gz merged_pairs.fastq.gz -short -fastq.gz lane3_NoIndex_L003_R1_001.frags.fastq.gz -short -fastq.gz lane3_NoIndex_L003_R2_001.frags.fastq.gz
OMP_NUM_THREADS=32 velvetg VelvetAssem -scaffolding yes -exp_cov auto -cov_cutoff auto -min_contig_lgth 200 -ins_length 150


[[File:Master 11022 12185.CDS-genomic.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 (y-axis) and Arabidopsis thaliana Col-0 (x-axis).  Results may be regenerated at: http://genomevolution.org/r/3w96]]


==Pseudo-Assembly==
[[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.]]


 
[[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.]]
 
'''Other Stuff'''
*python -m jcvi.formats.fastq convert  (read help file, default converstion
*python -m jcvi.apps.baseclean trim fastqfile (single ended)
*python -m jcvi.apps.baseclean trim R1.fastq.gz R2.fastq.gz (paired ended)
 
*Cat all the R1s together
*Cat all the R2s together

Latest revision as of 13:35, 14 May 2014

Twig2Genome Notes

Twig2Genome

Assembly

Creosote Assembly

Loading into CoGe

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)

SoapDeNovo assembly of contigs >= 2000nt: 6,976 contigs (20MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12185

ABySS assembly (bpsize=64, 2kb minimum contig size): 122,972 contigs (392MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12275

Velvet assembly 515,190 contigs (241MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12245

CLC4 assembly 685,475 contigs (508MB): http://genomevolution.org/CoGe/OrganismView.pl?dsgid=12244

Syntenic Path Assembly

Syntenic path assembly with SynMap of creosote (x-axis) and peach (y-axis). Results may be regenerated at: http://genomevolution.org/r/3w95
Syntenic 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.
Retrieved from "https://genomevolution.org/wiki/index.php?title=Creosote&oldid=6452"