Understand why MICADo do not call any alteration on the sample with 1000reads, but does so for the with 750reads

cbib / MICADo

Looking for mutations in PacBio cancer data: an alignment-free method

Other

8 stars 2 forks source link

Closed massyah closed 5 years ago

massyah commented 8 years ago

These were generated by the read sampler
They contain the exact same alteration
For the sample with 1000 reads, MICADo do not test for the ref vs alt TAACAGTTCCTGCATGGGCGGCATGAACCGGAGGCCC vs alt TAACAGTTCCTGCATGGGTGGCATGAACCGGAGGCCC

Note that the same apply for several other synthetic samples. Run logs are attached:

massyah commented 8 years ago

Expected mutation is Ref: TAACAGTTCCTGCATGGG C GGCATGAACCGGAGGCCC vs Alt: TAACAGTTCCTGCATGGG T GGCATGAACCGGAGGCCC

massyah commented 8 years ago

Might be related to issue #18

JRudewicz commented 5 years ago

watch issue #18