Today I attempted to take my BLAST output on Hummingbird and work through this blast-2-slim notebook of Steven’s. I had gone through this before using Cgigas. I got stuck today at the part where you
GO-GOslim.sorted. Will ask tomorrow. Additionally, today I worked on using
transrate to look at my C. bairdi assembly using Laura’s notebook as a guide.
BLAST to GO slim
The notebook is in the Downloads on Hummingbird in
grace-Cbairdi-transcriptome/notebooks/20181105-blast-Cbairdi_swiss-prot.ipynb. I couldn’t figure out how to
push it successfully to my GitHub repo notebook…
rsync-ed the notebook to Owl from Hummingbird: /scaphapoda/grace/Crab-project/20181025-blast-Cbairdi_swiss-prot.ipynb.
Anyway, I got stuck where you have to
GO-GOslim.sorted. I checked the two files in TextWrangler on Hummingbird, and they both looked good. Additionally, I looked at the files using
!tail. The one thing that I noticed, which I’m not sure is typical, is that in the
_blast-GO-unfolded.sorted file, there were a lot more rows of the “TRINITY_DN#…” column, than there where of the “GO….” column. Maybe that is the problem?
Will ask Steven or Sam tomorrow morning.
This Thursday is GSS. I’m hoping to have some basic assembly data to share on my poster, and I’ll have some stuff after I get the above issues sorted out.