Comments (2)
Hello Trandamere,
I am also new to RNA-SEQ analysis. I also followed this same protocol for my dataset. And i am having some problem for identifying novel transcripts and genes from the below gffcompare ".stats" file
i was running the following command $ gffcompare –r chrX_data/genes/chrX.gtf –G –o merged stringtie_merged.gtf. and i got the result file as
= Summary for dataset: stringtie_merged.gtf
Query mRNAs : 253181 in 70635 loci (216728 multi-exon transcripts)
(23264 multi-transcript loci, ~3.6 transcripts per locus)
Reference mRNAs : 216257 in 60158 loci (189357 multi-exon)
Super-loci w/ reference transcripts: 51791
-----------------| Sensitivity | Precision |
Base level: 100.0 | 93.0 |
Exon level: 99.9 | 93.6 |
Intron level: 99.4 | 94.0 |
Intron chain level: 99.7 | 87.1 | Transcript level: 99.8 | 85.2 | Locus level: 100.0 | 84.4 |
Matching intron chains: 188838
Matching transcripts: 215729
Matching loci: 60158
Missed exons: 0/623537 ( 0.0%)
Novel exons: 23741/674273 ( 3.5%)
Missed introns: 2160/383827 ( 0.6%)
Novel introns: 4747/405880 ( 1.2%)
Missed loci: 0/60158 ( 0.0%)
Novel loci: 10239/70635 ( 14.5%)
Total union super-loci across all input datasets: 70632 253181 out of 253181 consensus transcripts written in merged.annotated.gtf (0 discarded as redundant).
from this result where i can find the novel genes and transcripts?
from gffcompare.
Closing this as it is generally addressed in the answer to #24 and the online documentation.
from gffcompare.
Related Issues (20)
- .tracking file does not always show the 'closest' reference transcript for SETs HOT 4
- replace '.' class code in .tracking and combined files with a concatenation of participant codes
- How is 'tss_id' assigned?
- Exon-level sensitivity is not 100%, but there are zero missed exons.
- disable discard duplicate
- gffcompare -d option not having effect HOT 1
- Number of samples reported exceeds the maximum HOT 3
- keep reference transcripts
- Where can I get gffcompare v0.12.5? HOT 1
- Comparing two different mouse strain gtf files
- Definition of "Locus Level" HOT 2
- Availability of new version on Conda and biocontainers
- gff file and the XLOC numbers are not ordered correctly
- Extract transposon HOT 1
- Missing tmap and refmap files in query gff directory + missing statistics in stats file
- How to keep the CDS info for each transcript? HOT 2
- Transcript classification codes ? HOT 1
- class_code 'u' in all transcripts HOT 1
- Merging different gtf files
- How is class code X is identified in stranded data?
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from gffcompare.