Entering edit mode
3.9 years ago
543541656
▴
20
perl annotate_variation.pl --geneanno --dbtype refGene --buildver msuv7 --outfile /home/riceUsers/fyx/qtl-bm4-mw/5.compare/20.awk_custom2/mut_index_2/msu_result /home/riceUsers/fyx/qtl-bm4-mw/5.compare/20.awk_custom2/mut_index_2/candidatesnp_final.avinput msudb/
NOTICE: Output files are written to /home/riceUsers/fyx/qtl-bm4-mw/5.compare/20.awk_custom2/mut_index_2/msu_result.variant_function, /home/riceUsers/fyx/qtl-bm4-mw/5.compare/20.awk_custom2/mut_index_2/msu_result.exonic_variant_function
NOTICE: Reading gene annotation from msudb/msuv7_refGene.txt ... Done with 66338 transcripts (including 0 without coding sequence annotation) for 55986 unique genes
NOTICE: Processing next batch with 267 unique variants in 267 input lines
NOTICE: Variants with invalid input format are written to /home/riceUsers/fyx/qtl-bm4-mw/5.compare/20.awk_custom2/mut_index_2/msu_result.invalid_input
after running the 'perl annotate_variation.pl' command,why the message do not display read fasta like this 'NOTICE: Reading FASTA sequences from ricedb/gcf_refGeneMrna.fa ... Done with 52 sequences'
