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

Print reversed-complemented sequences #196

Merged
merged 1 commit into from
Dec 8, 2024

Conversation

AndreaGuarracino
Copy link
Contributor

Real read sequences:

>A00132:53:HFHJKDSXX:2:1335:29622:10520_R2
TTATGGTTCACCAACAGTTATTTCACTTTAAGGCAGTCTAAATAGAGGATTTCATGGAGTTAAATCAGACCGTTTGTGGGAAAATTGTTTCTTAAAAAAAAAAAACATTTTAAAGGAAAGGATGATGATTTGTTGTTTTTGGTTTTTTTC
>A00132:53:HFHJKDSXX:2:1335:29622:10520_R1
ATTTTGCCACGTTGTCCAGCCTGGTCTCGAACTCCTGGCCTCAAGTGATCAGCCCTCCTCAGCCTCCCAAAGTGCTGGCATTACAGGCATGAGTCACTGCGCCCAGCCTGTTGGGTTTTTTTGGAAAAAAACCAAAAACAACAAATCATC

Input alignments:

 cat c.gaf | column -t
A00132:53:HFHJKDSXX:2:1335:29622:10520_R2  150  0  150  +  >23>24>25>27>28>31>32>33>35>36     175  1   151  150  150  0  cg:Z:150M
A00132:53:HFHJKDSXX:2:1335:29622:10520_R1  150  0  150  +  <45<44<42<41<40<39<37<36<35<33<32  209  30  180  150  150  0  cg:Z:150M

Before this PR:

fgfa -i x.fgfa gaf c.gaf -s

A00132:53:HFHJKDSXX:2:1335:29622:10520_R2       TTATGGTTCACCAACAGTTATTTCACTTTAAGGCAGTCTAAATAGAGGATTTCATGGAGTTAAATCAGACCGTTTGTGGGAAAATTGTTTCTTAAAAAAAAAAAACATTTTAAAGGAAAGGATGATGATTTGTTGTTTTTGGTTTTTTTC
A00132:53:HFHJKDSXX:2:1335:29622:10520_R1       GCGAGTTCGAGACCAGGCTGGACAACGTGGCAAAGTGAGGAGGGCTGATCACTTGAGGCCACTCATGCCTGTAATGCCAGCACTTTGGGAGGCACAGTGTTTTTTTCCAAAAAAACCCAACAGGCTGGGCGGGATTTGTTGTTTTTGAAA

After this PR:

A00132:53:HFHJKDSXX:2:1335:29622:10520_R2       TTATGGTTCACCAACAGTTATTTCACTTTAAGGCAGTCTAAATAGAGGATTTCATGGAGTTAAATCAGACCGTTTGTGGGAAAATTGTTTCTTAAAAAAAAAAAACATTTTAAAGGAAAGGATGATGATTTGTTGTTTTTGGTTTTTTTC
A00132:53:HFHJKDSXX:2:1335:29622:10520_R1       ATTTTGCCACGTTGTCCAGCCTGGTCTCGAACTCCTGGCCTCAAGTGATCAGCCCTCCTCAGCCTCCCAAAGTGCTGGCATTACAGGCATGAGTCACTGCGCCCAGCCTGTTGGGTTTTTTTGGAAAAAAACCAAAAACAACAAATCATC

@sampsyo
Copy link
Collaborator

sampsyo commented Dec 8, 2024

This is great; thank you for catching this!!

@sampsyo sampsyo merged commit 53c0908 into cucapra:main Dec 8, 2024
12 checks passed
@sampsyo
Copy link
Collaborator

sampsyo commented Dec 8, 2024

While this approach was perfectly good, it made me finally get around to trying a "lazy" reverse-complement strategy in #197.

sampsyo added a commit that referenced this pull request Dec 29, 2024
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.

2 participants