This is a test version of Biostars. For the public version, visit https://www.biostars.org.
-N and -L in bowtie2

Sorry friends,

I used the below syntax but the result seems strange for me, what is the reason please

bowtie2 -x [name] -N 0 -L 15 -U [file_trimmed.fastq -S file.sam

[izadi@lbox161 ~]$ bash
[izadi@lbox161 ~]$ cd /usr/data/nfs6/izadi/project
[izadi@lbox161 project]$ export BT2=/usr/people/home/izadi/fereshteh/bowtie2-2.2.5
[izadi@lbox161 project]$ echo $BT2
/usr/people/home/izadi/fereshteh/bowtie2-2.2.5
[izadi@lbox161 project]$ $BT2/bowtie2 -N 0 -L 15 -x rRNA --un SRR1211041_trimmed_unmapped.fastq -U SRR1211041_trimmed_unmapped.fastq -S mapped_and_unmapped.sam
0 reads
0.00% overall alignment rate
bowtie2 mismatch seed
• 2,435 views • link
updated by Ram Ram • written by zizigolu zizigolu

1 answer

-N 0 is the default, so you needn't specify it unless you're changing things. Often as -L is decreased you need to increase -D as well, so try the defaults and see if they produce any alignments and then tweak from there.

Also, you really don't want to give the same file to --un and -U. You'll end up over writing the file you're trying to align!

• 0 views • link
updated by Ram Ram • written by Devon Ryan Devon Ryan

Thank you Devon

Actually my adviser asked me to do like this. I'm too new yet and I can't get the reason of these arguments.

• 0 views • link
updated by Ram Ram • written by zizigolu zizigolu

Log in to answer this question.

More posts like this
  • Error: Encountered internal Bowtie 2 exception (#1) (ERR): bowtie2-align exited with value 1
    written by salamandra salamandra
    I'm trying to run bowtie2. These are the commands: UNMAP=/pathtofile/unmappedreads.txt BWTID=/pathtofile/btw2index/hg19 READS=/pathtofile/merged.fastq OUTPUT=/pathtofile/merged.SAM bowtie2 -p 4 —un $UNMAP —no-unal -x $BWTID -U $READS -S $OUTPUT …
  • alignment rate in bowtie2
    written by zizigolu zizigolu
    Hi, I used bowtie2 to map GSE69414 data sets but rate of alignment is 0 or max 0.02 How is that possible? Am I doing …
  • error when running trimmomatic
    written by zizigolu zizigolu
    Hi friends, I was going to adapter trimming by trimmomatic then I did so but ``` [izadi@lbox161 ~]$ cd /usr/data/nfs6/izadi/Trimmomatic-0.33/ [izadi@lbox161 Trimmomatic-0.33]$ java -jar trimmomatic-0.30.jar …
  • What is the meaning of 'Error running 'prep_reads'?
    written by zizigolu zizigolu
    Hi friends, I was running the tophat2 but I got error. Do you know the reason please? ``` [izadi@lbox161 bowtie2-2.2.5]$ $TOP/tophat2 -p 8 -G genes.gtf …
  • error in bowtie2 alignment
    written by zizigolu zizigolu
    Hey friends, For months I'm working with bowtie2 but today I am getting error My reference genome fasta is in bowtie2 folder, I indexed the …
  • how I can add a path to my directory?
    written by zizigolu zizigolu
    Hello, I used these commands to manage the folder and starting analysis ``` [izadi@lbox161 tophat-2.1.0.Linux_x86_64]$ cd /usr/data/nfs6/izadi/test [izadi@lbox161 test]$ mkdir $HOME/bin1 [izadi@lbox161 test]$ ln -s …
  • no happen when runing tophat in fedora
    written by zizigolu zizigolu
    friends, I downloaded [TopHat Linux x86_64 binary][1] and unzipped with `tar -xzvf tophat-2.1.0.Linux_x86_64.tar.gz` [1]: http://ccb.jhu.edu/software/tophat/downloads/tophat-2.1.0.Linux_x86_64.tar.gz I have already installed bowtie2 and samtools I did like …
  • failed to open BAM file
    written by zizigolu zizigolu
    Friends, ``` [izadi@lbox161 ~]$ bash [izadi@lbox161 ~]$ cd /usr/data/nfs6/izadi/gij [izadi@lbox161 gij]$ pwd /usr/data/nfs6/izadi/gij [izadi@lbox161 gij]$ ls S288C_rRNA.fsa my.sam SRR1944914.fastq my1.sam SRR1944914.sra my1_highquality.bam SRR1944914_trimmed.fastq my1_highquality.sam SRR1944914_trimmed_unmapped.fastq …
  • export function in linux
    written by zizigolu zizigolu
    Hello, I have a file containing the reads and reference genome. How I can tell bowtie2 where are the which it need now? ``` [izadi@lbox161 …
  • unexpected bowtie2 unpaired alignment behavior
    written by mark.rose • mark.rose
    Hello Seemingly according to the manual, bowtie2 will attempt to align paired-end reads first as pairs but then as unpaired reads if necessary. If this …