Question: HISAT2 alignment Error?
0
gravatar for k.kathirvel93
2.2 years ago by
k.kathirvel93260
India
k.kathirvel93260 wrote:

Hi EveryOne,

I am working in Paired-wnd RNA-seq data set, When i am trying to align my reads against reference HG19 genome by using HISAT2, am getting this error:

Warning: skipping mate #2 of read 'ST-E00142:733:HTKKMCCXY:1:1207:16163:5229 2:N:0:ACTTGA' because it was < 2 characters long

This is happening only after trimming. Thanks in advance.

ADD COMMENTlink written 2.2 years ago by k.kathirvel93260
1
gravatar for ATpoint
2.2 years ago by
ATpoint42k
Germany
ATpoint42k wrote:

This is not an error but a warning. Apparently your trimmer cuts of most of the read because of an adapter detection or low-quality of the base pairs. It would be smart to discard reads shorter than a certain cutoff, e.g. 20bp categorically as they have little change to be properly aligned (which most trimmers do by default), simply because of the size. The shorter the size, the more likely they will be multimappers. A 2bp read cannot be aligned anyway, as one might imagine quiet well. Please post the full command of the trimming.

ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by ATpoint42k

Thanks for the reply ATpoint, Here is my Adapter removal command:

cutadapt3 -a AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC  -A AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTAGATCTCGGTGGTCGCCGTATCATT -o Sample_1.fastq -p Sample_2.fastq Sample1_1fastq Sample1_2.fastq

I am not getting any warning, while aligning with raw unfilterd data, with filtered data only am getting warning.

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by k.kathirvel93260

Cutadapt has no default minimum read length. I never noticed that because I always use -m 18. Can also be 25 or so. In any case, it should prevent you from these warnings.

ADD REPLYlink modified 2.2 years ago • written 2.2 years ago by ATpoint42k

Should I do something for this warning OR i can go further ?. Because after this, i am getting more than 95 % alignment rate. Thanks

ADD REPLYlink written 2.2 years ago by k.kathirvel93260
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: 1135 users visited in the last hour