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

How to use mutscan with UMI ? #16

Open
biofer opened this issue Aug 17, 2022 · 1 comment
Open

How to use mutscan with UMI ? #16

biofer opened this issue Aug 17, 2022 · 1 comment

Comments

@biofer
Copy link

biofer commented Aug 17, 2022

when use fastp + gencore , can not use GATK SamToFastq.
gencore copy the flag of R1 reads as R2 reads (or R2 as R1)?
For example,
Before gencore

A00877:990:HL5F3DSX3:4:2111:6361:29731:UMI_CCCTGAGTGT	1187	chr1	27023101	60	2S111M	=	27023101	111	GCGGGAAAGGAGCTGCAGGACGGGGCCGAGAGCAATGGGGGTGGCGGCGGCGGCGGAGCCGGCAGCGGCGGCGGGCCCGGCGCGGAGCCGGACCTGAAGAACTCGAACGGGAA	FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF:	NM:i:0	MD:Z:111	MC:Z:2S111M	AS:i:111	XS:i:20	RG:Z:M30713-3MJ15_M30713-3MJ15_llluma_L4
A00877:990:HL5F3DSX3:4:2111:6361:29731:UMI_CCCTGAGTGT	1107	chr1	27023101	60	2S111M	=	27023101	-111	GCGGGAAAGGAGCTGCAGGACGGGGCCGAGAGCAATGGGGGTGGCGGCGGCGGCGGAGCCGGCAGCGGCGGCGGGCCCGGCGCGGAGCCGGACCTGAAGAACTCGAACGGGAA	FFFFFF,FFFFF::FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF	NM:i:0	MD:Z:111	MC:Z:2S111M	AS:i:111	XS:i:20	RG:Z:M30713-3MJ15_M30713-3MJ15_llluma_L4

After gencore

A00877:990:HL5F3DSX3:4:2111:6361:29731:UMI_CCCTGAGTGT	1187	chr1	27023101	60	2S111M	=	27023098	111	GCGGGAAAGGAGCTGCAGGACGGGGCCGAGAGCAATGGGGGTGGCGGCGGCGGCGGAGCCGGCAGCGGCGGCGGGCCCGGCGCGGAGCCGGACCTGAAGAACTCGAACGGGAA	FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF	NM:i:0	MD:Z:12G2G17T24G8G11G1G18T1G8	MC:Z:3M1D110M	AS:i:66	XS:i:22	RG:Z:M30713-3MJ15_M30713-3MJ15_llluma_L4	FR:i:3
A00877:990:HL5F3DSX3:4:2111:6361:29731:UMI_CCCTGAGTGT	1187	chr1	27023101	60	2S111M	=	27023101	111	GCGGGAAAGGAGCTGCAGGACGGGGCCGAGAGCAATGGGGGTGGCGGCGGCGGCGGAGCCGGCAGCGGCGGCGGGCCCGGCGCGGAGCCGGACCTGAAGAACTCGAACGGGAA	FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:	NM:i:0	MD:Z:111	MC:Z:2S111M	AS:i:111	XS:i:20	RG:Z:M30713-3MJ15_M30713-3MJ15_llluma_L4	FR:i:3

Am I making a mistake, or is there some other solution。
I want run fastp with the ConsensusReads.
Thank you for your great work!

@yanjunbiocsi
Copy link

same requests

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