-
-
Notifications
You must be signed in to change notification settings - Fork 13
Open
Description
Hello dear developers,
I'm using smash to compare two public genomes but I'm not getting anything. What parameters should I adjust?
(base) enriquepola@DESKTOP-62E8R2K:~/smash$ smashpp -r fungal_reesei_QM6a_chromosome.fa -t fungal_virens_FT333_masked.fa
====[ PREPARE DATA ]==================================
[+] fungal_reesei_QM6a_chromosome.fa (FASTA) -> fungal_reesei_QM6a_chromosome.seq (seq) finished.
[+] fungal_virens_FT333_masked.fa (FASTA) -> fungal_virens_FT333_masked.seq (seq) finished.
====[ REGULAR MODE ]==================================
[+] Creating model of fungal_reesei_QM6a_chromosome.fa done.
[+] Filtering fungal_virens_FT333_masked.fa done => 0 segments
====[ INVERTED MODE ]=================================
[+] Creating model of fungal_reesei_QM6a_chromosome.fa done.
[+] Filtering fungal_virens_FT333_masked.fa done => 0 segments
Total time: 13 sec.
(base) enriquepola@DESKTOP-62E8R2K:~/smash$ head fungal_reesei_QM6a_chromosome.fa
>CP016232.1 fungal reesei QM6a chromosome I, complete sequence
cctaaccctaaccctaaccctaacctaACCTAACCCTaaaccctaaccctaaccctaaaCCTAACCCTAAACCTAACCCT
TTAAAAAATgctattaaaaattaaaaagaagttttaaaatatagaTTTTTTAATagattttattatattttttatatctA
TTTTAACTTaagatttaaaaaaatttttttctttttttaccTTAACCTTTACCTTATAAGTAAGATATAAaagctagata
tttttatttaagtaaaaaatttaaatttaataagaaatctttatttaagtaaagaatttaagttaagtaagaaatcttta
tttaaataaagaatttaagttaagtaagaaatctttatttaagtaaagaatttaaagtaggtaagttaggctaaatatct
ttatttaagtaaagaatttaagtttagctaaatatttttatttaagtaaagaatttaagtttaataagaaatctttattt
aagtaaaaaatttaagttaagtaagaaatctttatttaagtaaaaaatttaagttaagtaagaaatttttatttaagtaa
agaaaTATTTAAAGTAGGTAAGTTaagctagatatttttatttaagtaaagaatttaagtttaataagaaatctttattt
aagtaaaaaatctaagttaagtaagaaatttttatttaagtaaagaatttaggttaagtaagaaatctttatttaagtaa
(base) enriquepola@DESKTOP-62E8R2K:~/smash$ head fungal_virens_FT333_masked.fa
>Ch_1
ATTGTACTTCGTTAGTTACGTACAAAAGACACCGACAACTTTCTTCAGCA
CCTCCCTAAccctaaccctaaccctaaccctaaccctaaccctaacccta
accctaaccctaaccctaaccctaaccctaaccgtaCCGCTGGATAGCAT
TCCACAGACCTTCTAGAAGAAAGCACGTGATCGGAAAAATCCGGAGATCT
CTGAAGTTAGTAGGTTTTCTCTTTGATTTAGTACGGTTATGTGCCCTCCG
ATAGAAGGAGGGTTCTGGGACCCAAAAAGACCTGAGGATAAACCTCGGGT
TGGCGGATATAATCTTTAAATAAAACTAGCCCCCTTTAAAATCCTATATT
TCAGTTAGATATAGGAGGGCATTAATAGTAATAGGAGAGTAAGCGCGAGG
ATTCCAGCCTTCGTCTGAAATTTGAATTTGCATTGCTGCCTAACTAAGGG
Could you help me? The genomes measure 34Mbp (ref) and the other 41 (target).
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels