Question: bowtie2 alignment to rRNA, ERCC spike ins, Repetitive Elements and index creation
gravatar for Melonhead
11 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 • 246 views
ADD COMMENTlink modified 9 months ago by Biostar ♦♦ 20 • written 11 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 11 months ago • written 11 months ago by ATpoint32k

Apologies, will do next time Thank you

ADD REPLYlink written 11 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: 1290 users visited in the last hour