Skip to content

Conversation

@lucidtronix
Copy link
Contributor

@yfarjoun @fleharty @takutosato BQSRing with CNNs.

@takutosato takutosato self-requested a review May 25, 2018 13:40
@takutosato takutosato self-assigned this May 25, 2018
bed_dict = bqsr_bed_file_to_dict(args.bed_file)
record_dict = SeqIO.to_dict(SeqIO.parse(args.reference_fasta, "fasta"))

margin = args.window_size/2

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It looks like pysam's fetch() captures all reads that intersect the [start,end] interval. So the current code is capturing many reads that have NIST variants. We should get somewhat better empirical quality if we get rid of the /2.

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

Successfully merging this pull request may close these issues.

3 participants