I have direct RNA data mapped to the gencode transcriptome with minimap2. Finding the 'true' transcript of origin for a read is nontrivial as there are many secondary alignments with very close alignment scores to the primary. After visualising I can see some alignments are to transcripts which start further 3 prime than my alignment. However, due to the mechanism of direct RNA sequencing, the three prime ends of reads are the true end site.
I want to discard alignments to transcripts that have a 3' start site over 100nt prior to my read start site.
I've thought about simply extracting TES from the gencode gtf but these are genomic coordinates and I need to use the transcriptome mapping. Another way I've been thinking is if the query end site is over 100nt of my read end site, to discard the alignment. But I am not sure how to do this, any ideas? Thanks.
•
written
by
josiegleeson1 •
0 answers
No answers yet.
Log in to answer this question.
-
lncRNA quantification using featureCounts
written by kuntalasb •
-
Selection of start codons in Ensembl
written by i.sudbery
-
mapToTranscripts() with + and - strand regions
written by newscient •
-
how bowtie deals with reads cowering simultaneously 5` and 3` ends of DNA fragment
written by k.panov •
-
Prediction of TSS from tuxedo pipeline
written by bioinfo_ga •
-
problem determining read directions
written by blooria •
-
Bowtie2 or RSEM (with Bowtie2) alignment of PE reads with reference transcriptome?
written by deanpett •
-
SAM format paired-end alignment handling
written by darxsys
-
Get all the translation start site for human genome
written by Ming Tommy Tang
-
What preliminary steps should I do before checking the quality of sequence for RNA-seq
written by M K
Did you end up solving this issue? I am facing it now... direct RNA sequencing is tough!
the problem looks simple but I would need a example bam with a few reads to test.
are u the dude from jvarkit? cuz if so, nice work bro