Question: How to align .fastq files to an Addgene vector sequence?
gravatar for Kristin Muench
2.2 years ago by
United States
Kristin Muench480 wrote:


I would like to find out if any reads in my .fastq file were transcribed from a vector sequence that is not in the human reference genome.

The sequence is available here:

I thought I might try aligning the .fastq files to the sequence using TopHat, as if the sequence were a human reference genome, and seeing if any alignments pop up.

However, I'm not sure how to go about doing this.

Should I make the sequence above into a .gtf file somehow?

How do I make the corresponding annotation (.gff) file?

Is there an easier way to go about doing this? E.g., isolating every sequence in the .fastq file and using grep to search for it in the vector sequence?

Thank you!

EDIT: I am aligning Illumina RNA-Seq .fastq files. Also, I'd appreciate any resources folks have regarding how to modify or add chromosomes to a reference file!

rna-seq alignment • 883 views
ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by Kristin Muench480

My standard answer when someone mentions Tophat:
You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. (If you can't get access to that publication, let me know and I'll -cough- help you.) There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon.

ADD REPLYlink written 2.2 years ago by WouterDeCoster43k

Eek, thank you for the heads-up! I'll take a look at those other methods.

ADD REPLYlink written 2.2 years ago by Kristin Muench480
gravatar for WouterDeCoster
2.2 years ago by
WouterDeCoster43k wrote:

I believe the most correct would be to add an additional chromosome, your vector, to the human genome and index that for alignment using e.g. STAR or HISAT2 (assuming you have Illumina RNA-seq data, which you did not specify).

Aligning only to the vector (without the human genome) or vice versa might lead to false positive/negative alignments.

ADD COMMENTlink written 2.2 years ago by WouterDeCoster43k

Great point! And yes, I am using Illumina RNA-Seq data.

ADD REPLYlink written 2.2 years ago by Kristin Muench480
gravatar for genomax
2.2 years ago by
United States
genomax80k wrote:

You could use from BBMap suite to bin the reads to quickly find out which ones are from human genome and which are from vector. You can decide how to handle those reads that multi-map (both within and across the genomes).

ADD COMMENTlink written 2.2 years ago by genomax80k

Oh, interesting! Thank you!

ADD REPLYlink written 2.2 years ago by Kristin Muench480
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1813 users visited in the last hour