Rename the read ID

[jiangzy26]

I got a problem when i markduplicates using picardtools , it shows Values was put into PairInfoMap more than oce, i think it's the same reads ID problem , so i want to know how to handle this problem with using the command rename of seqtk, Thank you very much

[tseemann]

You can rename your reads as follows:

seqtk rename old.fq.gz prefix_ | gzip > new.fq.gz

This wil name all the reads prefix_1, prefix_2, .... , prefix_346341 etc