Question: --u option in bowtie2
0
gravatar for A
4.5 years ago by
A3.5k
A3.5k wrote:

Hi

in the rRNA filtering part i am going to map the fastq reads on the rRNA sequences and keep only the reads which has not been mapped, in cmd i should type:

 bowtie2 -x [name of the bowtie2-build indicized file containing the rRNA sequence] --un [name of the fastq file which will contain the UNMAPPED reads] -U [name of the fastq file containing the reads] -S [name of the .sam file that will contain the MAPPED and UNMAPPED reads]

I could not understand  about --un option because i don't know which i should type instead of [name of the fastq file which will contain the UNMAPPED reads]..

is the name optional??  in typed some optional names with gz Suffix but always error

 

sequencing • 1.3k views
ADD COMMENTlink modified 4.5 years ago by jivarajivaraj40 • written 4.5 years ago by A3.5k
0
gravatar for RamRS
4.5 years ago by
RamRS24k
Houston, TX
RamRS24k wrote:

You should type in the name for a file which will be created afresh to hold unmapped reads. Ensure the name you provide ends with a .fastq

ADD COMMENTlink written 4.5 years ago by RamRS24k
0
gravatar for jivarajivaraj
4.5 years ago by
jivarajivaraj40 wrote:

i typed --un unmapped.fastq

but it saying that: warning bowtie2-align run directly. it is recommended that you run the wrapper script "bowtie2" instead

i am going to have tow files: one sam file containing mapped and unmapped reads and another file containing unmapped reads by --un option but i dont know for afresh file which name i should type or totally why i get error...

please help me more if possible

ADD COMMENTlink modified 4.5 years ago • written 4.5 years ago by jivarajivaraj40
1

Please move this to a separate question. And paste your full command when you create the new question.

ADD REPLYlink written 4.5 years ago by RamRS24k
1

In addition to what Ram said, read the error message. It's telling you what you did wrong.
 

ADD REPLYlink written 4.5 years ago by Devon Ryan92k

Devon,

as i typed the error for Ram, when i used bowtie2 instead of bowtie2_align, it says not recognized as external or internal command

ADD REPLYlink written 4.5 years ago by A3.5k
1

Ensure that bowtie2 is in the same directory. You'll need to have perl installed. I actually have no idea if any of this will work on Windows, there's a reason people don't use it in Bioinformatics.
 

ADD REPLYlink written 4.5 years ago by Devon Ryan92k
1

How did you know OP was using Windows? And why does it look like "sarah" and "jivarajjivaraj" are the same person/working together?

ADD REPLYlink written 4.5 years ago by RamRS24k

hey Ram, please

i need help not being searched by bioinformatics detectives!!

you suppose she knows me well but cant help more just sharing my question

ADD REPLYlink modified 4.5 years ago • written 4.5 years ago by A3.5k

You are mistaken if you think that someone sharing your question on a separate thread helps your case. Mods will redirect the person to your question and close the new question. That is standard policy in any forum. Plus, this other user does not have the requisite reputation to add a positive effect (in fact, their usage patterns are frankly quite odd)

And if a second user posts a related question on your thread, people's attention will be diverted between the subject matters, and this is called "hijacking the thread". We do this only to ensure usage etiquette is followed so you can get quick answers.

ADD REPLYlink modified 4.5 years ago • written 4.5 years ago by RamRS24k

:(

thanks Devon

ADD REPLYlink written 4.5 years ago by A3.5k
Please log in to add an answer.

Help
Access

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