Question: HISAT2 alignment Error?
0
gravatar for k.kathirvel93
12 months ago by
k.kathirvel93200
India
k.kathirvel93200 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 12 months ago by k.kathirvel93200
1
gravatar for ATpoint
12 months ago by
ATpoint22k
Germany
ATpoint22k 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 12 months ago • written 12 months ago by ATpoint22k

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 12 months ago • written 12 months ago by k.kathirvel93200

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 12 months ago • written 12 months ago by ATpoint22k

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 12 months ago by k.kathirvel93200
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: 2523 users visited in the last hour