-
Couldn't load subscription status.
- 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).
Metadata
Metadata
Assignees
Labels
No labels