I can actually align to the part of transcriptome that I want (first I have to make fasta file and use it as reference file to align to). I can also make gtf file using the same file and count using that. but in both cases it returns 0 counts. so I thought maybe I should align to the genome and count on the reads that mapped to that part. but here I do not know how to do that.
counting the reads that mapped to a part of transcriptome
2 answers
Hi Sara,
Did you align the raw reads to a reference genome? if so what is the current format of your aligned reads?
Tarek
• 0 views
•
link
written
by
tarek.mohamed
Hello Sara,
once you have mapped the data against the genome of your interest, you can simply use FeatureCounts program for instance:
featureCounts -a yourGTF.gtf -t UTR5 -g transcript_id -o counts_UTR5_transcripts.txt
The -t option allows you to select a type of features to be counted. The -g option allows you to choose "how to group your data". Using transcript_id (or something equivalent in your GTF/GFF), you tell to FeatureCounts: "ouput the counts for each transcript".
Also if you prefer to use R, you can use "Rsubread" package to count your reads using featureCounts () as mention by glihm.
T
• 0 views
•
link
written
by
tarek.mohamed
Log in to answer this question.
More posts like this
-
labelling the clusters from the CD4 T cell scRNAseq data
written by Sara
-
read counts per chromosome arm
written by Sara
-
coverage of the reads over coding sequence only for one gene
written by alireza346 •
-
counting the reads that map to a special window
written by Christleemans •
-
nucleotide frequency of each gene at certain position
written by Sara
-
isolating in frame reads from a file
written by Sara
-
problem with the alignment of RNAseq data
written by Sara
-
frequency of each nucleotide at each position in sequence
written by Sara
-
counting the reads mapped only to the coding sequence
written by alirezamomeni707 •
-
alignment to a special part
written by alirezamomeni707 •
