Skip to content

No results in two genomes at the chromosome level #18

@enriquesanMx

Description

@enriquesanMx

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
No labels

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions