MCScanX result problem
2
0
Entering edit mode
4.2 years ago
HZZ0036 ▴ 20

Hello,

I want to use MCScanX to see the synteny blocks for my weed draft genome and maize genome. I downloaded maize genome from NCBI. Then I prepared .blast and .gff file and run ./MCScanX according to the manual, however, in the collinearity result file, it showed:

############### Parameters ###############
# MATCH_SCORE: 50
# MATCH_SIZE: 5
# GAP_PENALTY: -1
# OVERLAP_WINDOW: 5
# E_VALUE: 1e-05
# MAX GAPS: 25
############### Statistics ###############
# Number of collinear genes: 0, Percentage: 0.00
# Number of all genes: 1
##########################################

Besides, there is no .tandem file and html is empty. Does anybody can tell me the reason about this? Thanks.

Zhang

genome sequence • 4.5k views
ADD COMMENT
2
Entering edit mode
4.2 years ago
kashiff007 ★ 1.6k

Can you check your gff and blast.output file? The chr name should be same in both. And the blast.output should contain weed-weed, weed-maize, maize-weed and maize-maize while your gff should have weed the maize.

ADD COMMENT
0
Entering edit mode

My draft genome has no chr name, only scaffold name. I creat a link for .blast and .gff files. Could you please have a look? Thanks. https://filemover.auburn.edu/files/files/1505337598_goosemaize.blast https://filemover.auburn.edu/files/files/1505337635_goosemaize.gff

ADD REPLY
1
Entering edit mode

As you can see, you have ".t1" is present with gene name in blast output file, but not in gff file. Either remove it from blast output or add to gff file.

ADD REPLY
0
Entering edit mode

I have removed it from blast output. Now the gene name in gff file and blast output file is same, but it still doesn't work. Could you please have a look again? Here is the link:

https://filemover.auburn.edu/files/files/1505917638_goosemaize.blast https://filemover.auburn.edu/files/files/1505917678_goosemaize.gff

Thank you so much!

ADD REPLY
1
Entering edit mode

what command you have used for blastall?

ADD REPLY
1
Entering edit mode

Here is the script for blastall: blastp -db /mcScanX/goose_maize -query /mcScanX/goose_maize.fa -out /mcScanX/goose_maize.blast -evalue 1e-10 -outfmt 6 -num_threads 20 -num_alignments 5

ADD REPLY
2
Entering edit mode

remove this: -num_alignments 5

ADD REPLY
0
Entering edit mode

Thanks. I have tried, but it still doesn't work. The log file showed:

Reading BLAST file and pre-processing
Generating BLAST list
0 matches imported (7239387 discarded)
0 pairwise comparisons
0 alignments generated
Pairwise collinear blocks written to /mcScanX_gg_maize_analysis/goose_maize.collinearity [29.246 seconds elapsed]
Writing multiple syntenic blocks to HTML files
mz1 GRMZM5G884466 301344228 301363304^M.html
Done! [0.706 seconds elapsed]
ADD REPLY
1
Entering edit mode

I solved this by making sure that both the .blast and the .gff files were TAB-delimited. I had the same "0 matches imported" problem when the .gff file was space-delimited.

The tr command can be used to change it easily:

```bash < file.gff tr ' ' '\t' > file_tab.gff

ADD REPLY
0
Entering edit mode

Hi,

I have solved this problem, but I got a very strange circle. Could you please see the new post: Very strange circle from MCScanX

Thanks.

ADD REPLY
0
Entering edit mode

Hi, I run into the same problem, how do you solve it? Thanks in advance

ADD REPLY
0
Entering edit mode

whats the solution? facing the same problem.

ADD REPLY
0
Entering edit mode

Hii Anjali,

Did you find solution for this issue. Can you please help me to fix this.

ADD REPLY
0
Entering edit mode

can you share your files so that we can go through them to know the problem?

ADD REPLY
0
Entering edit mode

Hii

I am also struggling with same problem. Can you help me out how to fix this issue. I tried checking blast file and gff file ids . All are same. So please help me with your experience to fix this.

ADD REPLY
0
Entering edit mode

Hi! I used SyMAP instead of MCScan.

ADD REPLY
0
Entering edit mode
3.0 years ago
nareshmvr • 0

Hii Liu,

I am struggling with same issue. If you fixed it please can you help me to sort this out.

ADD COMMENT
0
Entering edit mode

Did you find a way ?

ADD REPLY
0
Entering edit mode

No

发自我的 iPhone

在 2019年10月23日,下午7:25,venura on Biostar <mailer@biostars.org<a rel="nofollow" href="mailto:mailer@biostars.org">mailer@biostars.org> 写道:

Activity on a post you are following on Biostarhttp://www.biostars.org

User venuravenura wrote Comment: MCScanX result problemC: MCScanX result problem:

Did you find a way ?

ADD REPLY
0
Entering edit mode

"Solution is to put .blast and .gff file inside MCScanX folder without creating a separate folder." it woreked!

ADD REPLY

Login before adding your answer.

Traffic: 2258 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6