Skip to content

Missing reads with --simplified on #17

@maressyl

Description

@maressyl

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

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions