I have to admit, I haven't tried this since I switched to CRISPick.
Hey everyone, is anyone using the CLI version of CHOPCHOP for sgRNA design? I am trying to get it to work, but I have some questions about it.
I need a .gene_table file that I can extract from the UCSC. However, the tracks
RefSeq GenesorEnsemble Genesare not available. forRefSeq GenesNCBI RefSeqseems to be the correct track. However, I can't find an properEnsemble Genestrack here. Any ideas where to find this?Is there a way to directly use
HGNC Symbolsas IDs for the query?Do all reference files (
.2bit,.ebwt,.gene_table) need to be in different folders?
1 answer
This is described in README.
":gene_table" is simple genePred table with named columns e.g. https://chopchop.cbu.uib.no/genomes/isoforms/hg38.gene_table You could transform gff>f into genePred by using tools like http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/gtfToGenePred
To use the symbols you have to have them in the proper column in the genePred table, you could process the table to replace the ensemble gene names, or something like this might work
gtfToGenePred -genePredExt -geneNameAsName2 gencode.v29.annotation.gtf hg38.genePred- Paths to the reference files are specified in the config.json
{
"PATH": {
"PRIMER3": "./primer3_core",
"BOWTIE": "bowtie/bowtie",
"TWOBITTOFA": "./twoBitToFa",
"TWOBIT_INDEX_DIR": "/your/full/path/to/2bit_folder",
"BOWTIE_INDEX_DIR": "/your/full/path/to/ebwt_folder",
"ISOFORMS_INDEX_DIR": "/your/full/path/to/ebwt_transcriptome_folder_and_2bit_of_genome",
"ISOFORMS_MT_DIR": "/your/full/path/to/vienna_MT_folder",
"GENE_TABLE_INDEX_DIR": "/your/full/path/to/genePred_folder"
},
"THREADS": 1
}
Log in to answer this question.