Align to custom sequences using STAR aligner
0
0
Entering edit mode
2.3 years ago
roy.granit ▴ 880

I have 3 amplicons that I have sequenced and now I wish to align them to the sequences of these three genes so I can quantify how many reads mapped to each gene.

In order to use STAR must I go through the steps of genome indexing or is there a way to align stright with the fasta file containing the 3 genes?

Thanks!

RNA-Seq STAR • 1.2k views
ADD COMMENT
1
Entering edit mode

you could make (and index) a custom file with only those 3 genes (perhaps go for the gene locus sequences, rather than CDS or such). Keep in mind their is a risk that you bias your analysis.

Personally I think STAR is fast enough to use the whole genome as reference and do some filtering after mapping.

ADD REPLY
0
Entering edit mode

I did not mention that some of these genes do not exist in the genome..

ADD REPLY
0
Entering edit mode

If you don't expect off-target amplification then why not use whole genome indexes if they are available. That can also show you off-target amplification (if there was any).

If not, you will need to index the reference sequences (which should be trivial for a small number).

BBMap allows you to align directly without creating indexes first: bbmap.sh in=your.fq.gz ref=your_multifasta.fa out=aligned.bam

ADD REPLY
0
Entering edit mode

Thanks, how reliable is bbmap? never heard of it..

ADD REPLY
1
Entering edit mode

You are missing some important knowledge in that case :-) BBMap goes toe to toe with any NGS aligner out there. I will leave this link here so you can see what the BBTools suite can do (besides alignments). A guide for BBMap the aligner is available here.

ADD REPLY

Login before adding your answer.

Traffic: 1220 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6