Pairfq icon indicating copy to clipboard operation
Pairfq copied to clipboard

loosing info when I use pairfq addinfo

Open vgaout opened this issue 3 years ago • 2 comments

When I apply pairfq addinfo to trim-galore-cleaned-reads (8.8GB) I obtain outputs with a weight of 4.4GB. What could be the reason of this?

perl pairfq addinfo -i sample1.fq -o sample1_pairinfo.fq -p 1

perl pairfq addinfo -i sample2.fq -o sample2_pairinfo.fq -p 2

After that, if I use pairinfo makepairs, the program does not find any paired read in the samples:

perl pairfq makepairs -f sample1_pairinfo.fq -r sample2_pairinfo.fq -fp for_paired.fq -rp rev_paired.fq -fs for_unpaired.fq -rs rev_unpaired.fq

vgaout avatar Sep 28 '21 11:09 vgaout

It is hard to understand from the commands alone. I will need to see the format. Please show a sample of the first file.

head -4 sample1.fq

Also, the same with the output:

head -4 sample1_pairinfo.fq

That should be informative enough to understand the problem.

sestaton avatar Oct 01 '21 19:10 sestaton

When I use that command in the sample to which I applied pairfq it appears like:

@SRR3583357.3.1/1 NAACTACTGTTAGTTTGTTAAGAAGTGATTTTTTTGGGGCGAGTTGAGTGGTGTGTAGTGGGGTGTTATTTGTTAGTTTATATAGTGATGTAGGTAGTGTG + #1=DFDFFHHHFHHHIJIIJIJJJJHIEIJJJJJJIIIIJJGHBDFBD?CDCCABDDDEDBDD8;;BADDEDDCCCCDEDDEEEADDDEDDCDDCDDDDCA

while the original was:

@SRR3583357.3.1 SN7001427:168:C7T7WACXX:2:2201:1499:2047 length=101
NAACTACTGTTAGTTTGTTAAGAAGTGATTTTTTTGGGGCGAGTTGAGTGGTGTGTAGTGGGGTGTTATTTGTTAGTTTATATAGTGATGTAGGTAGTGTG
+SRR3583357.3.1 SN7001427:168:C7T7WACXX:2:2201:1499:2047 length=101


vgaout avatar Oct 11 '21 14:10 vgaout