PASA pipeline updating annotation results in no changes
0
0
Entering edit mode
3.0 years ago
rob234king ▴ 610

I'm running these commands leaving the numbers as default in config and cleaned transcripts already.

module load Miniconda3/4.9.2

# create database
$PASAHOME/scripts/create_sqlite_cdnaassembly_db.dbi -c alignAssembly.config -S /home/data/pest_genomics/pasa/opt/pasa-2.4.1/schema/cdna_alignment_sqliteschema

# Upload annotations
$PASAHOME/scripts/Load_Current_Gene_Annotations.dbi -c alignAssembly.config -g Chilo_suppressalis_v2_genome_220620_correct2.fasta -P chilo_transfer2_corrected2.gff3

module load SAMtools/1.12-GCC-10.2.0

$PASAHOME/Launch_PASA_pipeline.pl -c annotCompare.config -A -g Chilo_suppressalis_v2_genome_220620_correct2.fasta -t clean_transcripts.fa --CPU 10

When I do this though no changes are made. I use the same transcripts in maker to update the gene set and lots of UTR are added. I was hoping PASA would do a better job but I'm not getting any differences. I've changed the headers of the transcripts thinking maybe that would have an effect but no luck. Anyone experienced the same? I've posted on PASA google group but no response so far.

PASA • 818 views
ADD COMMENT
0
Entering edit mode

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. If your code has long lines with a single command, break those lines into multiple lines with proper escape sequences so they're easier to read and still run when copy-pasted. I've done it for you this time.
code_formatting

ADD REPLY

Login before adding your answer.

Traffic: 2229 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