Review mRNA primer pairs positions on genome assembly
Dear forum, I would like to review a set of ~ 50 human RNA primer pairs (in fasta) on a genome browser. Their size range from 18 to 25nt. Some of them overlap splicing junctions. I have tried two approaches:
- blast on refseq RNA: fine but need a conversion to grch37 assembly (allowing split of the query over both side of splice junction)
- HiSat2 alignment: many are not found, I don't get why. (parameter L was set to 14, lower than min query length by 4nt).
I'm sure there is a simple solution as my objective is I guess not very original, but I don't find out. Do you have any guidance for generating a bed file fitting my needs? Thanks, Manu
• 1,826 views
•
link
written
by
Manu Prestat
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Merge overlapping paired end reads from BAM file.
written by chrisk •
-
Inserting Junctions Into Genome after Indexing with STAR
written by hkarakurt
-
How to Filter out False positive junctions detected from RNA-seq data?
written by julbtn •
-
junction_annotation.py: How many 'novel' splice junctions/splice events are resonably expected from…
written by RNAseqer
-
Standard options for mapping with STAR in context of RNA-seq data analysis regarding splicing
written by caggtaagtat
-
MapSplice Error: filter alignment failed
written by liguanghaoli •
-
STAR using splice junctions from SJ.out to find flanking exons (GRCh37)
written by stefanos.bamopoulos •
-
hisat2-build failed with 'is not reverse-deterministic, so reverse-determinize'2
written by ranzhou1005 •
-
Analyzing types of Junctions from BAM and GTF files
written by luisccleto •
-
Sjdboverhang Option In Star
written by Martombo
