Tweet
Viral genome vs transcripts
i have a fasta file containing all the virus genome which will affect cassava. i have transcript file also in fasta. in transcript file each transcript have 100 nucleotide length. I want to know whether these transcripts are inside the virus genom. any one help me pls?
-
-
One relatively straight-forward way to attack this problem would be to make a combined fasta file for the viral genomes and cassava and then map the transcripts/transcript fragments against it. You can then get an idea how likely each fragment is to come from the host or virus. I assume that tophat or something like that would be an appropriate aligner, since presumably you have a mixture of host (i.e. spliced) and viral (I assume transcripts would be single exon, but I've only made viruses, not studied them) reads and not doing so might bias things.
I should note that I've never done what you're trying to do, but if no one replies with a better idea then this is enough to get you started. -
doubt regarding tophat
can we put the two files as input in tophat.means that genome file in fasta and also the transcript file?Originally posted by dpryan View PostComment
-
Sort of, I can read a few different meanings into what you wrote so I'll write some pseudo-commands just to ensure that my meaning is clear. Let us assume that your Cassava genome file is called "Cassava.fa" with an annotation file "Cassava.gtf", your virus genome file is called "Viruses.fa", and the transcripts whose origin you're interested in determining are in a file called "reads.fa". The work flow would then be something like the following:Originally posted by vijesh View Post
You can then see which reads originated from which chromosome/contig or whichever virus by using samtools view on the resulting accepted_hits.bam file (you'll probably want to filter out reads that map equally to both Cassava and viral genomes, for which you'll probably have to write a program).Code:cat Cassava.fa Viruses.fa > CombinedGenome.fa bowtie2-build CombinedGenome.fa Combined tophat -G Cassava.gtf Combined reads.fa
It would be good to compare the results with and without using the GTF annotation, off-hand I'm not entirely sure how or if that might bias things.Comment
-
viral genom vs transcript
i dont have the annotation file(gtf) .i have only the virus genome file (fasta) and transcript file(fasta).....Comment
-
You don't actually have to use the annotation file, I expect it might have just made things a bit more reliable (or maybe not, one would have to check).Comment
-
Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.
The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...06-02-2026, 10:05 AM -
With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.
Introduction
Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...05-22-2026, 06:42 AM
Comment