How to create kallisto index that does not include ribosomal genes
Hello,
I would like to create a kallisto index that contains coding and noncoding RNA but also exclude the ribosomal genes. I know I need the cdna and ncrna fasta files from the ensembl website initially. However, I am not sure how to remove the ribosomal genes from the reference. Is there an easy way to remove the ribosomal genes?
Thank you
1 answer
SortMeRNA could be an option, too (initially designed for metagenomic data)
• 0 views
•
link
written
by
Carambakaracho
Why would that be necessary? Using isolated databases always comes at the risk of false alignments. If you quantify against a reference that includes rRNA together with all the other annotated transcripts then you can later easily remove the rRNAs from your count matrix.
Log in to answer this question.
More posts like this
-
Why am I getting such a big difference in alignment rates in kallisto?
written by bioinfo •
-
Why am I getting different results with kallisto?
written by bioinfo •
-
Do indexes made with cDNA fasta file contain only mRNAs?
written by bioinfo •
-
How do I get estimated gene counts for a reference that contains GFP using Kallisto?
written by bioinfo •
-
How to filter .fa file by gene_biotype
written by bioinfo •
-
How do I extract predicted genes from annotation from MAKER-P?
written by eennadi •
-
Are ribosomal protein coding genes depleted as well during rRNA depletion?
written by akibio •
-
Should cndas provided by Ensembl be filtered by ccds or biotype prior to running kallisto?
written by holgerbrandl •
-
what is the corresponding proteome or transcriptome (cDNA fasta in Ensembl) file in NCBI?
written by SMILE
-
How to filter Ensembl cDNA and ncRNA FASTA files by primary assembly?
written by michael.mccarthy.404 •
Did you want to remove both rRNA and ribosomal protein coding genes? Also, why did you want to remove these in the first place? It may cause some of these reads to be assigned to incorrect transcripts.
I would like to remove only the rRNA. We found out that there was an rRNA contamination from a probe we used so we would like to remove those transcripts from our data.
Then you pick out the rRNA genes after aligning to everything. Removing things you know are present from the reference could make the alignment inaccurate.