Below is an outline of the steps I followed to create the Capitella teleta Quickload:
1. Convert genome .fasta to .2bit
rsync -aP rsync:./faToTwoBit GCA_000328365.1_Capca1_genomic.fna C_teleta_Jan_2013.2bit
2. Create genome.txt
rsync -aP rsync:./twoBitInfo C_teleta_Jan_2013.2bit genome.txt
3. Get gene models from NCBI (.gff), then convert .gff to .bed
git clone git@bitbucket.org:lorainelab/genomesource.git
path+=('~/Documents/Repos/genomesource/')
export PYTHONPATH="${PYTHONPATH}:~/Documents/Repos/genomesource/"
./gff3ToBedDetail.py -g ~/Downloads/genomic.gff -b ~/Downloads/C_teleta_Jan_2013_refGene.bed
4. Check if NCBI has any information for this genome using its txid (NCBI:txid283909) (Conclusion: it didn't)
brew install tnftp
ftp ftp:
get gene2accession.gz
quit
gunzip -c gene2accession.gz | grep '^283909\t' > 283909.gene2accession.txt
5. Sort, gzip, and tabix the .bed file made in step 3
sort -k1,1 -k2,2n C_teleta_Jan_2013_refGene.bed | bgzip > C_teleta_Jan_2013_refGene.bed.gz
tabix -0 -s 1 -b 2 -e 3 C_teleta_Jan_2013_refGene.bed.gz
6. Sanity check the .bed and .2bit files - Add the .2bit file as a reference, then drag/drop the .bed file into IGB. Confirm that gene models are present, labeled correctly, and the chromosomes listed are in a logical order.
7. Create annots.xml
brew install svn
svn checkout --username=guest --password=guest https:svn mkdir C_teleta_Jan_2013
svn cp A_gambiae_Oct_2006/annots.xml C_teleta_Jan_2013/
7. Add the new genome to contents.txt and .htaccess
Tested following the instructions above. Everything looks good.
Closing ticket.