NanoSim icon indicating copy to clipboard operation
NanoSim copied to clipboard

reads headers empty in fasta file

Open jlipovac opened this issue 3 years ago • 3 comments

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.

jlipovac avatar Jun 28 '22 13:06 jlipovac

Hi @jlipovac ,

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

Ka Ming

kmnip avatar Jun 28 '22 16:06 kmnip

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

kmnip avatar Jun 28 '22 18:06 kmnip

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

jlipovac avatar Jun 28 '22 18:06 jlipovac