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

quick fix for broken read (mapped to no ref bases) #3662

Merged
merged 1 commit into from
Oct 5, 2017

Conversation

tedsharpe
Copy link
Contributor

there is a read with the cigar 44I107S, which is complete gibberish, in the NA12878 bam.
it's nonsense, but we can't afford to blow up.

@SHuang-Broad
Copy link
Contributor

Is there a way to drop such reads, considering it is in such an incoherent state?

@tedsharpe
Copy link
Contributor Author

Not easy, because then we'll have a read that claims to be part of a pair but which has no mate, for example. It'll mess up other stuff.

@SHuang-Broad
Copy link
Contributor

For documentation purposes, can you post the SAM record of the offending read here?

@codecov-io
Copy link

codecov-io commented Oct 4, 2017

Codecov Report

Merging #3662 into master will decrease coverage by 0.002%.
The diff coverage is 0%.

@@               Coverage Diff               @@
##              master     #3662       +/-   ##
===============================================
- Coverage     79.264%   79.263%   -0.002%     
  Complexity     17394     17394               
===============================================
  Files           1144      1144               
  Lines          62930     62930               
  Branches        9551      9551               
===============================================
- Hits           49881     49880        -1     
  Misses          9241      9241               
- Partials        3808      3809        +1
Impacted Files Coverage Δ Complexity Δ
.../hellbender/tools/spark/sv/utils/SVFastqUtils.java 64.646% <0%> (-0.505%) 5 <0> (ø)

@tedsharpe
Copy link
Contributor Author

tedsharpe commented Oct 4, 2017

HJYFJCCXX160204:4:2122:2646:19786	163	2	233750853	70	44I107S	=	233750920	214	GTGTATGCATGTGCGGTGTGTACGGTGTGATGTATGTGCCATGTTTGTGCATGTGTGCTGTGTATACACGTGTGCAGTGTTGTGTGCCGTGCATATGTGCAGAGTGGTGTGTGTGGGTGTGGTGTGTGTGGGCAGGCGGGGGGTGTGCTGG	=?>???>>=?>>==8=>====>=8=3;==<>>>>>>>:>><<;>;>>>>>>/60><>=?;82><=6;:%%8><2./4;;>???;;>=/&<<############################################################	MC:Z:4S147M	OC:Z:29M4D15M107S	PG:Z:MarkDuplicates.1N	RG:Z:HJYFJ.4	NM:i:44	MQ:i:60	OQ:Z:AAFFFKKKAFKFKKFKKKKKKKKKK<AKKKFFKFKFKAKKA<AF<FKKKKK,<,FAFFKA77FAA7F<,,7FA77,7<AFKKKA<FF7,AA############################################################	UQ:i:0	AS:i:29

EDITED_FOR_EYEBALLING

@tedsharpe tedsharpe force-pushed the tws_quick_fix_for_broken_read branch from 4d88b0e to f39ab50 Compare October 4, 2017 20:54
@tedsharpe tedsharpe force-pushed the tws_quick_fix_for_broken_read branch from f39ab50 to c0f77e5 Compare October 5, 2017 15:33
Copy link
Member

@cwhelan cwhelan left a comment

Choose a reason for hiding this comment

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

:+

@tedsharpe tedsharpe merged commit 9609cb3 into master Oct 5, 2017
@tedsharpe tedsharpe deleted the tws_quick_fix_for_broken_read branch October 5, 2017 17:00
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.

4 participants