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

Look into possibility to drop --no-discordant for Bowtie 2 mode to support dovetailing of reads #14

Closed
FelixKrueger opened this issue Jan 19, 2016 · 1 comment

Comments

@FelixKrueger
Copy link
Owner

Can we explore how easy it would be to disable the Bowtie 2 option --no-discordant since this would allow dovetailing reads to align straight away. This would be very useful for PBAT-type runs that require clipping from their 5'-ends.

@FelixKrueger
Copy link
Owner Author

Having looked at this in some detail we decided not to drop --no-discordant because this would then all of a sudden start to report really weird alignments, e.g. reads facing away from each other or reads on different chromosomes. Instead we enabled the new option --dovetail (which will be turned on by default for --pbat libraries) which will now allow dovetailing reads to be reported:

 --dovetail        It is possible, though unusual, for the mates to 'dovetail', with the mates
                   seemingly extending 'past' each other as in this example:

  Mate 1:                 GTCAGCTACGATATTGTTTGGGGTGACACATTACGC
  Mate 2:            TATGAGTCAGCTACGATATTGTTTGGGGTGACACAT
  Reference: GCAGATTATATGAGTCAGCTACGATATTGTTTGGGGTGACACATTACGCGTCTTTGAC

  By default, dovetailing is considered inconsistent with concordant alignment, but setting 
  --dovetail causes Bowtie 2 to consider dovetailing alignments as concordant. This becomes 
  relevant whenever reads are clipped from their 5' end prior to mapping, e.g. because of
  quality or bias issues. --dovetail is set automatically for PBAT libraries.

This new option does indeed seem to work nicely (see this QC Fail article). Added in 395a793 and 8662516.

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