Minimal length of sequence for reference alignment.

Hi! When I am using bowtie 2 to reference alignment, console displays warnings about length of my sequences:

Warning: skipping mate #2 of read 'SRR988073.8033840 HWI-ST833:124:6:2202:2830:97445 length=101' because it was < 2 characters long Warning: skipping mate #2 of read 'SRR988073.8078016 HWI-ST833:124:6:2202:19987:137144 length=101' because length (1) <= # seed mismatches (0)

I have trimmed my data to improve quality of my sequences, but in fastq file is many short sequences now. Is there any minimal length of sequence for reference alignment? Shoud I remove short sequences? And how long sequences shoud be?

alignment next-gen sequencing

• 3,929 views • link

updated 10.2 years ago by WouterDeCoster

• written 10.2 years ago by mieszko91 •

2 answers

Typically short reads may map to multiple genome regions. Thus, it is recommended to remove the short sequences before alignment. I usually use 35 bp as the minium post-trimming length for bowtie2.

• 1 views • link

written 10.2 years ago by GZ1995

Well a length <2 characters is for sure far too short. Avoid too aggressive trimming, aligners commonly can handle low quality bases and suboptimal parts of a read.

• 1 views • link

written 10.2 years ago by WouterDeCoster

This ^. There is really no need to trim based on Q scores (which may be the case here) for re-sequencing data.

• 1 views • link

written 10.2 years ago by GenoMax

More posts like this

How to use Bowtie2 for mapping pair end reads and in what format reference genome should be used?
written 5.6 years ago by priya.saxena •

Hello Everyone I am trying to map the pair end reads using Bowtie2, I have formatted the custom reference genome using the format given in …
Warnings in the Hisat2 alignment output files
written 5.9 years ago by newbie

Before doing the alignment, I used fastqc on the fastq.gz files and observed adapter content. ![enter image description here][1] So, I removed the adapter content …
Number of reads in Fastq and Bowtie
written 7.5 years ago by yp19

Hi guys, When I check the number of reads in my fastq (paired end library), I get **4071779** paired reads. I do this using "wc …
Evaluation of HISAT2 Alignment Result
written 8.3 years ago by modarzi

Hi, I run hisat2 for one of my sample(RNA-seq fastq) but I receive lots of Warning. You can see the result of this analysis: Warning: …
Can bowtie2 mapping paired-end reads of different length?
written 9.0 years ago by whiteeagle115 •

Hello! I began to analyze some reads of a genome of a haploid organisms. The genome was made with HiSeq and during trimerization, the adapters …
Hisat2 warning: skipping mate
written 9.1 years ago by veronicaschroeder78

I am a newbie in NGS data analysis and trying to run hisat2 after successfully building a index, with the following bash script: trimmedfiles=... (list …
Bam to fastq conversion is giving error?
written 9.2 years ago by kirannbishwa01

I am trying to take the PE reads (100 bp) aligned to chr2 of one species and then convert it to PE fastq again. Then …
error in bowtie2 alignment
written 9.6 years ago by prasoonagar •

Hi I have done the trimming of my reads using cutadapt after trimming i did FASTQC and find that three modules * Per Base Sequence …
Bowtie memory warning
written 10.3 years ago by umn_bist •

I am aligning my files with bowtie 1.0.0 via our HPC. bowtie /UCSC/hg19/Sequence/BowtieIndex/genome -p 16 --chunkmbs 16 -t -m 1 -S --max ${file1:0:9}_multimap.fastq -1 ${file1} …
Error while mapping reads to assembly using bowtie
written 10.8 years ago by Biocode_user •

I am trying to map assembly back to the reads using bowtie. But with one of the libraries, i get an error Error: Warning: skipping …