Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

loosing info when I use pairfq addinfo #14

Open
vgaout opened this issue Sep 28, 2021 · 2 comments
Open

loosing info when I use pairfq addinfo #14

vgaout opened this issue Sep 28, 2021 · 2 comments

Comments

@vgaout
Copy link

vgaout commented Sep 28, 2021

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

@sestaton
Copy link
Owner

sestaton commented Oct 1, 2021

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.

@vgaout
Copy link
Author

vgaout commented Oct 11, 2021

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


Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants