Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
You might try a tool such as MapSplice instead. It uses statistics based on the alignments themselves (such as distribution of reads across a junction) rather than any sequence-related measures to validate spliced alignments. It's not hard to use, but can be finicky and can consume its fair share of disk space.
-
Originally posted by gringer View PostUsing tophat for a viral genome seems a little odd -- I wasn't aware that viruses had introns.Last edited by drdna; 10-07-2013, 12:42 AM.
Leave a comment:
-
Using tophat for a viral genome seems a little odd -- I wasn't aware that viruses had introns.
Leave a comment:
-
Ah, well then those are really junk then. Predicting splicing based on a single base without a reference annotation seems like a bad idea! Thanks for the heads-up and you might consider making a bug report for tophat.
Leave a comment:
-
Tophat2 is stupid
Originally posted by dpryan View PostTo be fair, did you give tophat a reference GTF/GFF to use? If so, you're asking it to align against the transcriptome first and then convert those coordinates back to the genome. You're rather likely to get results like this by doing that (yes, it's probably better to simply soft-clip that one base, but you didn't use local-alignment and, anyway, it then matched the transcriptome).
Leave a comment:
-
To be fair, did you give tophat a reference GTF/GFF to use? If so, you're asking it to align against the transcriptome first and then convert those coordinates back to the genome. You're rather likely to get results like this by doing that (yes, it's probably better to simply soft-clip that one base, but you didn't use local-alignment and, anyway, it then matched the transcriptome).
Leave a comment:
-
Tophat2 produces thousands of invalid alignments
Users beware. Tophat2 produces thousands of CIGAR alignments that are just plain wrong. As an example, see the following .bam file line:
M01478:14:000000000-A5C8R:1:1106:17706:4022 256 CYDV_S1_L001_R1_trimmed_(paired)_contigs_50/104 4239 3 1M184N110M * CTGCTCTGCCCTATGCGATCTGTCCGATCGATCCTTCCAGACCATTGTGGAGGACGAAGATGTTGTTGATACCCCGAACGGACCGTGGCTCCCTGTGCAGGATGATGGTGT DEEEEFFFFFCFGGGGGGGGGGHHHGGGGGGHHGHHHHHHHGHHHHHHHGHHGGHGGGGGHHHHHHHHHHHHHHHGGGGGGGGGGGHGGHHGGHHHHHHHHGHHHHHHGHG AS:i:0 XN:i:0 XM:i:0 XO:i:0 XG:i:0 NM:i:0 MD:Z:111 YT:Z:UU XS:A:- NH:i:2 CC:Z:= CP:i:4423 HI:i:0
The CIGAR indicates a single match/mismatch, then 184 nucleotides that are missing from the read relative to the reference and then 110 match/mismatches. Obviously, this has to be wrong. There is no way you can predict that a single nucleotide matches a position 184 nucleotides upstream of the rest of an alignment, when there are about 40 intervening base positions that are equally valid matches. Furthermore, in this particular example, there is no mismatch. The first nucleotide precisely matches the nucleotide upstream of the 110 aligned nucleotides. In other words, Tophat is just "making stuff up"! Note: no paired end reads were used in this Tophat run.Tags: None
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:47 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Leave a comment: