Question: bowtie2 alignment to rRNA, ERCC spike ins, Repetitive Elements and index creation
gravatar for Melonhead
4 months ago by
Melonhead10 wrote:

I have downloaded .fastqfiles from SRA and now trying to filter out reads aligning to ERCC spike ins, rRNA and Repetitive Elements (from RepeatMasker).

I would like to do this step by step, repeating the command on the unaligned results from the previous step using something like this:

bowtie2 -p 2 - N1 --un-conc <unaligned fastq> --al-conc <aligned fastq> \
    -1 <dowloaded mate1 fastq. -2 <downloaded mate 2 fastq> \
    -S <output alignment sam>

As I got from the manual, first I need to build .bt2 indices for each alignment (namely, for ERCC, rRNA, RepetitiveElements) How can I create indices for the aforementioned command and which data should I download for it?

Would be happy for any help!

rna-seq alignment • 147 views
ADD COMMENTlink modified 10 weeks ago by Biostar ♦♦ 20 • written 4 months ago by Melonhead10

Hi artem.baranowsky, welcome.

Please use the formatting bar to highlight code. I made the changes for you this time.

enter image description here

ADD REPLYlink modified 4 months ago • written 4 months ago by ATpoint22k

Apologies, will do next time Thank you

ADD REPLYlink written 4 months ago by Melonhead10
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: 2174 users visited in the last hour