More posts like this
-
Convert pileup file to genetic map file
written by taniamahmood38 •Dear all, I have pileup file of pooled-seq samples generated using samtools mpileup feature. I want to convert it to map file to compute various …
-
Genomic intervals to transcriptomic intervals
written by EVRHI, I have a set of bed intervals from human genome and I need to convert the genomic intervals to transcriptomic intervals. Any guidance would …
-
FPKM values from normalized read counts
written by Bioin •Hi Biostars, I would like to get FPKM values from normalized read counts for that I am using this formula FPKM(Fragments Per Kilobase per Million) …
-
Splitting the Spliced reads from genome Alignment
written by EVRHi, I have a RNA-seq reads mapped onto genome using tophat2 and obtained a "bam" file which is later converted into a BED file. I …
-
Counting read orientation between ranges in a bed file
written by EVRHI, I have two bed file like follows: bed_A scaffold11296 36365 36414 scaffold11296 36471 36526 bed_B scaffold11296 36302 36334 - scaffold11296 36303 36334 + scaffold11296 …
-
Minimum and maximum Intron length of Oryza sativa
written by EVRHi, I am new to Tophat alignment. I want map my paired-end reads onto Oryza sativa Genome but I couldnt find any information for minimum …
-
Finidng the closest features near transcript boundaries in genome
written by EVRHi, I aligned both my transcriptome and reads to genome and found reads have been aligned outside the transcript boundaries. For an example, if a …
-
Counting reads on genome
written by EVRHi, I aligned my reads to de novo transcriptome using bowtie. Now the unaligned reads from the previous step was mapped to genome using tophat2. …
-
Identifying full length transcripts in a transcriptome
written by EVRHi, Is there any way to find full length and in-complete transcripts that has been assembled from a de novo transcriptome? Any guidance would be …
-
High Sequence Duplication levels in FastQC
written by EVRHi, I am new to NGS analysis. I have RNAseq samples generated from Illumina platform. When the samples loaded onto FastQC , I noticed that …
What have you tried?
I tried a tool called "Quast" but it gave only the minimum, maximum and average intron length of transcripts and overall summary of genome.
I'm not aware of a tool that does what you need (but it might exists nonetheless). I guess you need to do some scripting, and essentially I would approach this as
intron_length = transcript_length - (sum of exons and UTRs). As you can see the third field contains features such as 'exon', 'gene', 'transcript', and the last field contains information about the 'parent' transcript and gene a feature belongs to.