Build canine dna.primary_assembly.fa for kb ref
1
0
Entering edit mode
23 months ago
Candice • 0

I am following the instruction to build the index using kb.

The command should be

kb ref -i index.idx -g t2g.txt -f1 cdna.fa -f2 intron.fa -c1 cdna_t2c.txt -c2 /
intron_t2c.txt --workflow lamanno -n 8 /
Mus_musculus.GRCm38.dna.primary_assembly.fa.gz /
Mus_musculus.GRCm38.98.gtf.gz

In Ensembl dog genome, there is no such file as dna.primary_assembly.fa.gz. They're separate files (by chromosome) as shown here.

How could I create a dna.primary_assembly.fa.gz for canine genome that includes all chromosomes?

Thanks!

bustools Ensembl kallisto • 548 views
ADD COMMENT
1
Entering edit mode

You should be able to cat the files together.

ADD REPLY
1
Entering edit mode
23 months ago
Candice • 0

I've found the answer. It turns out the Readme in Ensemble says "If the primary assembly file is not present, that indicates that there are no haplotype/patch regions, and the 'toplevel' file is equivalent."

Therefore, I can use the 'Canis_lupus_familiaris.ROS_Cfam_1.0.dna.topleve.fasta.gz' here

ADD COMMENT

Login before adding your answer.

Traffic: 2549 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6