Questions tagged [rna-alignment]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1
vote
1answer
55 views

What is the meaning of split read?

I want to use rna seq data to later perform functional tests on fusion genes. so before that I need to filter the "best results" (of rnaseq) for deciding which candidates I actually want to ...
2
votes
1answer
23 views

Understanding ViennaRNA RNAdistance scoring table

I'm trying to compare the output of 2 different algorithms of RNA structure prediction (my implementation of Nussinov vs RNA-mfold algorithm) using the RNAdistance algorithm that is part of ViennaRNA ...
0
votes
1answer
39 views

What is the best way to map/align my reads on a given genome?

I am frequently using a ballgown package for my rnaseq analysis, but recently I have had a new task to have my reads mapped on two different genomes to understand the level of alignment between the ...
0
votes
0answers
39 views

Best way to align to find inserted sequence

We have some RNA from knock-in mice, there are two different sequences we're looking for. We have aligned to the mouse genome using STAR but the sequence isn't there which isn't too surprising What is ...
0
votes
2answers
149 views

Does rRNA depletion protocol give higher number of mapped reads in Intronic regions?

Recently, I have downloaded a publicly available dataset, which are 350 tumor samples. I see the following information from the published paper. They used Ribo Zero Gold and rRNA was depleted. Strand ...
2
votes
2answers
37 views

Is there a computational tool or possibility to identify mRNA isoforms from the count matrix of a bulk RNA sequencing dataset?

I have the counts matrix of an RNA sequencing dataset of fibroblasts and I wish to identify isoforms of a particular gene of interest in it. Can anyone please hint me on a bioinformatics method to ...
2
votes
1answer
521 views

What pitfalls exist with running FastQC on a bam file?

I was analysing a bad sequencing run of some RNA data using FastQC, I supplied it RNA-STAR-aligned bam files (and the human reference). The output of MultiQC had many more "unique reads" in ...
3
votes
1answer
262 views

cDNA and alignment mapping

I am confused with RNA seq alignment. My understanding is that after Mature mRNA is isolated from the cell, it is then fragmented and using reverse transcriptase enzyme a cDNA copy is created which is ...
9
votes
1answer
3k views

What are "split reads" and "intron clusters?"

I'm working through the example data set for LeafCutter and the documentation mentions "split reads": This will cluster together the introns fond [sic] in the junc files listed in test_juncfiles....
1
vote
0answers
97 views

How to correctly parallelise RSeQC scripts with GNU parallel?

I have a .bam, as ouput from STAR aligner, from which I need to extract some info using RSeQC while using all the computational resources available to increase ...
4
votes
1answer
173 views

Merging ssu-align alignments

ssu-align is an HMM-based MSA tool for ribosome's small subunit (SSU) sequences. It's bundled with three SSU models: archaeal, bacterial eukaryotic Approach Given a set of archaeal and ...
2
votes
1answer
1k views

Disk space error while aligning reads using STAR

Hi am trying to align RNA-seq reads using STAR through following commands, ...
4
votes
1answer
2k views

Can I run STAR without an annotation file?

I wish to use Rascaf to scaffold a fragmented draft genome. For this, I need to provide a BAM file of aligned RNA-seq reads and the draft genome. So, I indexed the draft genome with STAR like this: <...
3
votes
2answers
1k views

Error creating indices using STAR

I am trying to index wheat genome using STAR through following command ...
0
votes
1answer
345 views

Why is my STAR reference genome indexing aborting on my GNU/Linux server but not on my Mac OS X laptop?

I am running the following command to index my genome: ...
6
votes
2answers
3k views

5' and 3' bias in Rna-seq data

I'm working with rna-seq samples. I see 5' bias and also 3' bias in the per-base sequence content plot. From this link I see that the bias at the start of the sequences appears to be the result of ...
5
votes
3answers
189 views

Show presence of known mutation in RNA-seq data

We have RNA-seq fastq data from control (WT) patients and a patient with a point mutation at a known location in one gene. I'd like to retrieve the reads aligning to that gene and show the presence of ...
6
votes
1answer
1k views

STAR-long parameters for aligning RNA ONT reads to genome

Are there any suggested parameters to align ONT reads to the reference genome using STAR-long? For now, I used the parameters suggested here, but I noticed a weird behaviour. I have RNA reads (D. ...
8
votes
1answer
218 views

Convert local alignments to spliced alignments in SAM file

I mapped RNA reads to reference genome, using LAST in split mode, and converted the MAF alignment to SAM with maf-convert. My problem is that the transcripts are not reported in a spliced manner, ...
6
votes
3answers
118 views

How can I compute gene expression for a set of RNA reads?

I'm trying to compute a gene expression profile for an organism. I have gene nucleotide sequences of the mentioned organism stored in a fasta file and a set of paired reads stored in two separate ...
7
votes
1answer
91 views

How to interpret contig-alignment.psa produced by velvet

I'm using velvet to align given reads of RNA to given CDSs (i.e. coding areas and genes) of an organism, so I can generate gene-expression profiles. But after using ...
11
votes
2answers
473 views

What to use to edit RNA alignments?

I have many alignments from Rfam Database, and I would like to edit them. I saw that many tools are used for Protein sequence alignments, but there is something specific to edit RNA alignments ? e....