20 September 2013
Aphaenogaster, mason, and genome assemblyNext post Previous post
Gotelli Lab Meeting
First official Burlington UseR Meetup!
Last week’s script failed because mason user quota limited to 10GB. Changed directory to save intermediate files to data capacitor
/N/dc/scratch/tg-johnsg/results and ran fine.
Suggestion from khmer list to use RSEM for comparative expression analysis. Tutorial here. Note that it’s bundled with Trinity!
Ant hunting. Cold morning (39F overnight) and four of us found no colonies in nearly two hours of searching East Woods. Andrew and I went back out to Centennial Woods in afternoon and still only found one small colony. Not surprisingly, the cold weather seems to have pushed the ants underground. May try next week if weather warms back up…but may just be too late in the year to collect colonies.
Installed sabre to demultiplex reads.
sudo make failed with multiple errors
demulti_paired.c:(.text+0x78): undefined reference togzread’` and I found the solution on stackoverflow
Change this line:
$(CC) $(CFLAGS) $(LDFLAGS) $(OPT) $? -o $(PROGRAM_NAME) to $(CC) $(CFLAGS) $(OPT) $? -o $(PROGRAM_NAME) $(LDFLAGS) the -lz option has to be at the end.
and installed successfully!
Who: JSG, SHC, Bruce Kingham, Erin Bernberg, Olga Shevchenko
About: Genome assembly for Aphaenogaster
de novo genome assembly should definitely include Pac Bio
Pac Bio - single molecule sequencer - large 20-30kb library - avg. read lengths 5-6kpb - 60-80k reads per smart cell - use for scaffold - error rate is error rate of polymerase, can be high - lower error for smaller pieces, error - very good combined w
Hybrid assembly - PacBio and Illumina - requires knowledge of programs
Illumina - avoid mate-pair libraries can be problematic. Library prep lends itself to generating inverted or incorrectly-paired reads
What we need to provide for sequencing:
DGE - not commonly done. Cost of Illumina so cheap that not really necessary.
Down the road, for sequencing other genomes.
Tentative plan: 1) Start with Illumina sequencing (use Genomic-tip for DNA extraction to get idea of yield. DO NOT VORTEX. Be nice to DNA. Use wide-bore pipette tips after cell lysis) 2) Wait for males to develop in colony 3) When have enough males for 10ug DNA for PacBio, do that. May see major upgrade in sequencing chemistry by next spring so timing may work out. 4) BioNanoGenomics (for optical mapping) could be useful in structural variants or repeats.
Problem with oases_pipeline script because does not expand regular expression " -fastq -shortPaired A22-**-interleaved.fq " to include all A22 interleaved files, though this does work when running velveth. Must be an issue with the python wrapper. Specificed each file individually and running!
ssh authentication for for github so can
push commits without requiring login.
Generate ssh key and verify if doesn’t exist already
Change git remote origin using this command:
git remote set-url origin email@example.com:USERNAME/REPO.git
Should be able to push commits without password entry!
This work is licensed under a Creative Commons Attribution 4.0 International License.