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

Gaps_to_fasta.py issue #6

Open
yangyxt opened this issue Dec 2, 2020 · 2 comments
Open

Gaps_to_fasta.py issue #6

yangyxt opened this issue Dec 2, 2020 · 2 comments

Comments

@yangyxt
Copy link

yangyxt commented Dec 2, 2020

From my understanding that this python script is used to parse the gaps.fa file into many contigs with skip_length=10(default) between each contig. Which creates a new fasta file.

image

In this for loop command, I realize it only writes the sequence to the new file when encountering the next in gaps.fa file. In my case, I only have one contig, which leads to an empty output file. And it seems in this way, we'll always skip the last contig in gaps.fa file.

Pls let me know your thoughts. Thanks!

@jeremiek
Copy link
Contributor

jeremiek commented Dec 2, 2020

Dear Yangyxt,

You may be working with an old version. Have you tried pulling more recently? There should be another sequence following the for loop to handle that case.

Thank you

@yangyxt
Copy link
Author

yangyxt commented Dec 2, 2020

Dear jeremiek,

Thanks for the info!

I just checked the gaps_to_fasta.py in the latest release, and it seems the same for this part
image

I actually change it to while loop and readline() for my own usage and it works. Thanks!

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

2 participants