Question: How to generate custom fasta and gtf files for RNA spike-ins
0
gravatar for 2405592M
6 months ago by
2405592M30
2405592M30 wrote:

Hi guys, I've been doing some RNA-seq on mice and I used 3x E. coli tRNA spike-ins. I'm assuming for analysis, I'll have to generate a fasta and GTF file for these spike-ins? If so, how do I do this? Also, I'm assuming this will have to be merged with the mouse fasta and GTF files. How do I do this as well? Once these steps are completed, can I move on to alignment? The spike-ins were custom made oligos and so I only have the 3 sequences. I don't know any gene loci/coordinates etc.

rna-seq spike-in • 293 views
ADD COMMENTlink modified 6 months ago by michael.ante3.3k • written 6 months ago by 2405592M30
0
gravatar for michael.ante
6 months ago by
michael.ante3.3k
Austria/Vienna
michael.ante3.3k wrote:

Assuming the tRNAs are unspliced sequences, you can add them like the ERCC92 spike in annotation: Each spike in has its own contig and each gene annotation goes from start to the end.

Something like:

>mySpikeIn1
ACGA....

and the GTF:

mySpikeIn1   spikein  exon 1  9999 0.0  +  . gene_id "mySpikeIn1"; transcript_id "mySpikeIn1.1";

Please refer to the GTF file format documentation regarding spaces/tabs and necessary fields.

Cheers,

Michael

ADD COMMENTlink written 6 months ago by michael.ante3.3k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1502 users visited in the last hour