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

Missing reads with --simplified on #17

Open
maressyl opened this issue Feb 3, 2023 · 0 comments
Open

Missing reads with --simplified on #17

maressyl opened this issue Feb 3, 2023 · 0 comments

Comments

@maressyl
Copy link

maressyl commented Feb 3, 2023

Hi,

We are overall very happy with MutScan, but It seems some reads are missed when --simplified is on :

R1.fastq :

@NL500104:785:HCMC2AFX5:1:11303:22509:17003 1:N:0:ACAGTGAC+CAGTGACA
GGAAAGGGAAGGACTGGGAGAGACACAAAGACCAGAGCCAGCCTCAGGGACAAGAGATTCCAGTTTTAGGCCTTT
+
AAAAAEEEE6EEEEEEEEEEEEEEEEEEEAEEEEEEEEEEEE//EEAEEEEEEEEE/EEE/EEEE<EEEEE/<<E

R2.fastq :

@NL500104:785:HCMC2AFX5:1:11303:22509:17003 2:N:0:ACAGTGAC+CAGTGACA
TTCCCTCGTTGAACATGCTGTCAAACCAGGACACTGGCTCCAGCTTGTGTTTCTGTTCTTGACCATCGTCTGGGA
+
AAAAAEEEEEEE6EEEEEEAEAAEEEEEAEEEEEEEEEEEEEEAAEAE<EEEE/EAEEEEE/EEEEEEEEAEE/E

design.conf :

#name, left_seq_of_mutation_point, mutation_seq, right_seq_of_mutation_point, chromosome
rs8017_C, CCAGGACACTGGCTCCAGCTTGTGTTTCTG, C, TCTTGGCCATCGTCTGGGAGTGGACATGCA, chr16
rs8017_T, CCAGGACACTGGCTCCAGCTTGTGTTTCTG, T, TCTTGGCCATCGTCTGGGAGTGGACATGCA, chr16

mutscan --read1 R1.fastq --read2 R2.fastq --support 0 --html mutscan.html --mutation design.conf --simplified on


---------------
rs8017_C CCAGGACACTGGCTCCAGCTTGTGTTTCTG C TCTTGGCCATCGTCTGGGAGTGGACATGCA chr16

---------------
rs8017_T CCAGGACACTGGCTCCAGCTTGTGTTTCTG T TCTTGGCCATCGTCTGGGAGTGGACATGCA chr16

/usr/local/bin/mutscan --read1 R1.fastq --read2 R2.fastq --support 0 --html mutscan.html --mutation design.conf --simplified on
Mutscan v1.14.1, time used: 0 seconds

mutscan --read1 R1.fastq --read2 R2.fastq --support 0 --html mutscan.html --mutation design.conf


---------------
rs8017_C CCAGGACACTGGCTCCAGCTTGTGTTTCTG C TCTTGGCCATCGTCTGGGAGTGGACATGCA chr16

---------------
rs8017_T CCAGGACACTGGCTCCAGCTTGTGTTTCTG T TCTTGGCCATCGTCTGGGAGTGGACATGCA chr16
1, pos: 55, distance: 1, forward
@NL500104:785:HCMC2AFX5:1:11303:22509:17003 2:N:0:ACAGTGAC+CAGTGACA
TTCCCTCGTTGAACATGCTGTCAAA CCAGGACACTGGCTCCAGCTTGTGTTTCTG T TCTTGACCATCGTCTGGGA
+
AAAAAEEEEEEE6EEEEEEAEAAEE EEEAEEEEEEEEEEEEEEAAEAE<EEEE/E A EEEEE/EEEEEEEEAEE/E

/usr/local/bin/mutscan --read1 R1.fastq --read2 R2.fastq --support 0 --html mutscan.html --mutation design.conf
Mutscan v1.14.1, time used: 0 seconds
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

1 participant