search

DMU-lilab / pTrimmer

Used to trim off the primer sequence from mutiplex amplicon sequencing
GNU General Public License v3.0
21 stars 5 forks source link

Add primer name to fastq read names #26

Open OscarAspelin opened 5 months ago

OscarAspelin commented 5 months ago

Hi, very nice tool and nice addition with the --info flag. I'd like to ask you if it would be possible to also output the primer name (e.g,. the aux info) to each classified read name instead of the index on the "+" line?

XLZH commented 5 months ago

Thanks for your suggestion. However, adding information in the read name might affect the processing by some software. Therefore, it's better to place the additional information in the comment (start with '+' line). If you need to obtain the primer name corresponding to each trimmed read, it is easily to extract it from the amplicon file using the primer index given by pTrimmer (e.g., 'P=342' indicate the 342-th primer in the amplicon file).