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

reads headers empty in fasta file #168

Open
jlipovac opened this issue Jun 28, 2022 · 3 comments
Open

reads headers empty in fasta file #168

jlipovac opened this issue Jun 28, 2022 · 3 comments
Labels

Comments

@jlipovac
Copy link

Hi,

I'm trying to simulate metagenome reads (including chimeric reads) using pre-trained model and I succeeded in it but I got fasta file of simulated aligned reads with reads that have no headers. Something like this:

>
CAAAATCATAATACGGACGTTCAATTCAACAGCACCAACAAGAAAAACCAGCAGGCCGGCACCTTCGGCTGTTGGCTAAGTACCCCAAGCAAGCCTGATCGTCGTTGGGATGACCCTCGGGGGACGATGGTGTTCTACACTTTTACCACCTACTTGCAACAGTTCATGATCAACACGACCGGCCTGTCTAAGGAAAGTGCGCGCAGATGTTTTCTTCGTCGAGCTGCTGGTTGGAGTCCAACCGGTTCGGTCACGCATCCGATAGATTGCCGCAGCGCTCTTGATCTCTTCGGGGTCGGCGGAACCATCTTGACCGTGCCGCTGCCTT...
>
TTTTACGCGAGTCCTTCTATTTTAATTACATAAAAAAAAGCAGTACTTAATTAGCACTACCCTATGTCAAATGTTCATTAAAATTACAATAGAATTCTATTTAAATTTACCCACCTTTTAGATCTGTTACCACTAGCTTTCATATCAATCAAAAGTTGCTG...

In the file with simulated unaligned reads there are normal headers.

I'm wondering did I maybe make a mistake in genome list, abundance or dna type tsv files?

Genome list looks like:

genome1_name   path1
genome2_name   path2
genome3_name   path3
genome4_name   path4

Abundance:

Size   10000
genome1_name   40
genome2_name   20
genome3_name   10
genome4_name   30

Dna type:

genome1_name   genome1_complete_genome   circular
genome2_name   genome2_complete_genome   circular
genome3_name   genome3_complete_genome   circular
genome4_name   genome4_complete_genome   circular

Thanks.

@kmnip
Copy link
Collaborator

kmnip commented Jun 28, 2022

Hi @jlipovac ,

It looks like this is related to pull-request #167, which is pending review now.

Ka Ming

@kmnip kmnip added the bug label Jun 28, 2022
@kmnip
Copy link
Collaborator

kmnip commented Jun 28, 2022

Hi @jlipovac ,

I have merged the PR. The changes should be in the master branch now. Let us know if you still have any issues.

Ka Ming

@jlipovac
Copy link
Author

Hi @kmnip, now it works! Thank you very much! :)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants