Pairfq
Pairfq copied to clipboard
loosing info when I use pairfq addinfo
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
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.
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