Hi, this solution works for me while changing the sequence labels does not. Can anyone direct me to the documentation that shows the difference between applying this “–decontaminate-pairs=lenient” option compared to the standard output?
Hi all,
Thank you for reaching out to the bioBakery Lab forum.
Form Kneaddata>=v0.11.0, there is an option for --input1 and --input2 so if the read identifiers are the same in both pairs, no modification to the input files are necessary as Kneaddata will automatically take care of it.
Regards,
Sagun
Hello,
I am using kneaddata v0.12.0 with the options --input1 and --input2 specified, and have encountered the same issue others have described. My sequence identifiers are @A00674:488:HGYKGDSX5:1:1101:3965:1000 1:N:0:ACGATCAG+GAGACGAT for R1, and @A00674:488:HGYKGDSX5:1:1101:3965:1000 2:N:0:ACGATCAG+GAGACGAT for R2.
In other words, the read identifiers appear to be the same in both pairs.
I used default parameters for bowtie2.
2 Likes