Skip to content

Craig mfa extract#669

Draft
craiglowe wants to merge 5 commits intomainfrom
craigMfaExtract
Draft

Craig mfa extract#669
craiglowe wants to merge 5 commits intomainfrom
craigMfaExtract

Conversation

@craiglowe
Copy link
Copy Markdown
Contributor

@craiglowe craiglowe commented Apr 20, 2022

This adds some minor functionality to the code. It will now reverse complement a sequence if the bed file denotes the negative strand

@codecov
Copy link
Copy Markdown

codecov Bot commented Apr 20, 2022

Codecov Report

Merging #669 (6c85b27) into main (82ddeca) will decrease coverage by 0.91%.
The diff coverage is 28.57%.

@@            Coverage Diff             @@
##             main     #669      +/-   ##
==========================================
- Coverage   55.38%   54.47%   -0.92%     
==========================================
  Files         279      288       +9     
  Lines       26960    28043    +1083     
==========================================
+ Hits        14933    15276     +343     
- Misses      10965    11691     +726     
- Partials     1062     1076      +14     
Impacted Files Coverage Δ
fasta/modify.go 28.57% <0.00%> (-6.22%) ⬇️
cmd/multiFaExtract/multiFaExtract.go 25.97% <37.50%> (-5.54%) ⬇️
cmd/vcfWebAnnotate/vcfWebAnnotate.go 0.00% <0.00%> (-68.92%) ⬇️
browser/multiFa.go 79.66% <0.00%> (-3.98%) ⬇️
cmd/gtfToBed/gtfToBed.go 44.44% <0.00%> (-3.18%) ⬇️
cmd/globalAlignmentAnchor/globalAlignmentAnchor.go 82.87% <0.00%> (-2.71%) ⬇️
sam/bamRead.go 85.54% <0.00%> (-1.51%) ⬇️
vcf/header.go 48.48% <0.00%> (-1.30%) ⬇️
cmd/multFaVisualizer/multFaVisualizer.go 5.88% <0.00%> (-1.27%) ⬇️
bed/bed.go 53.21% <0.00%> (-1.26%) ⬇️
... and 23 more

Help us with your feedback. Take ten seconds to tell us how you rate us.

rimangan
rimangan previously approved these changes Apr 20, 2022
Comment thread cmd/multiFaExtract/multiFaExtract.go
Copy link
Copy Markdown
Contributor

@RavenYL RavenYL left a comment

Choose a reason for hiding this comment

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

If following strand could be written as a option that would be great!

Copy link
Copy Markdown
Contributor

@cshoben cshoben left a comment

Choose a reason for hiding this comment

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

Not sure if we are waiting to review this one, but happy to approve when it's time.

ans = fasta.RemoveGaps(ans)
ans = extractMultiHelper(records, b.ChromStart, b.ChromEnd, s.RemoveGaps, b.Strand == bed.Negative)
if b.FieldsInitialized >= 4 {
fasta.Write(fmt.Sprintf("%s.fa", b.Name), ans)
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Can you help me understand why we would name the output file differently if the fields initialized are less than or equal to 4? I'm assuming it would be to not lose information, but then we do not seem to include the chrom start end anymore?

@craiglowe craiglowe marked this pull request as draft March 9, 2023 14:03
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