NAL-i5K / general_issues

for issues and discussions not tied to a specific repository
2 stars 0 forks source link

New organism: Schistocerca piceifrons #191

Closed mpoelchau closed 2 years ago

mpoelchau commented 2 years ago

See https://gitlab.com/i5k_Workspace/workspace_roadmap/-/wikis/Adding-an-organism-CWL-update for full description of each task (requires gitlab login)

Tasks on Ceres:

Tasks on the NAL servers:

amcooksey commented 2 years ago

functional annotation CERES file path: /project/nal_genomics/amanda.cooksey/protein_sets/Schistocerca_piceifrons/NCBI Annotation Release 100 functional annotation

childers commented 2 years ago
cp /app/data/other_species/schpic/iqSchPice1.1/scaffold/GCF_021461385.2_iqSchPice1.1_genomic.fna /usr/local/i5k/media/blast/db/

cp /app/data/other_species/schpic/iqSchPice1.1/scaffold/analyses/NCBI\ Schistocerca\ piceifrons\ Annotation\ Release\ 100/GCF_021461385.2_iqSchPice1.1_cds_from_genomic.fna   /usr/local/i5k/media/blast/db/

cp /app/data/other_species/schpic/iqSchPice1.1/scaffold/analyses/NCBI\ Schistocerca\ piceifrons\ Annotation\ Release\ 100/GCF_021461385.2_iqSchPice1.1_rna_from_genomic.fna   /usr/local/i5k/media/blast/db/

cp /app/data/other_species/schpic/iqSchPice1.1/scaffold/analyses/NCBI\ Schistocerca\ piceifrons\ Annotation\ Release\ 100/GCF_021461385.2_iqSchPice1.1_translated_cds.faa /usr/local/i5k/media/blast/db/

1. Add organism: 
    1. python manage.py addorganism Schistocerca piceifrons 
2. Add each of the fasta files:
    * genome fasta: 
python manage.py addblast Schistocerca piceifrons -t nucleotide Genome Assembly -f GCF_021461385.2_iqSchPice1.1_genomic.fna -d 'Schistocerca piceifrons genome assembly, iqSchPice1.1' 
    * RNA fasta file: 
python manage.py addblast Schistocerca piceifrons -t nucleotide Transcript -f GCF_021461385.2_iqSchPice1.1_rna_from_genomic.fna -d 'Schistocerca piceifrons NCBI Annotation release 100, transcripts' 
    * Peptide fasta file: 
python manage.py addblast Schistocerca piceifrons -t peptide Protein -f GCF_021461385.2_iqSchPice1.1_translated_cds.faa -d 'Schistocerca piceifrons NCBI Annotation release 100, translated CDS' 
    * CDS fasta file: 
python manage.py addblast Schistocerca piceifrons -t nucleotide Transcript -f GCF_021461385.2_iqSchPice1.1_cds_from_genomic.fna -d Schistocerca piceifrons NCBI Annotation release 100, CDS 
1. Process each fasta file:
    * Make blast database: 

    python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_genomic.fna -m 
    python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_rna_from_genomic.fna -m 
    python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_cds_from_genomic.fna -m 
    python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_translated_cds.faa -m 

    * populate sequence table: 
        python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_genomic.fna -p 
        python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_rna_from_genomic.fna -p 
        python manage.py blast_utility GCF_021461385.2_iqSchPice1.1_cds_from_genomic.fna -p 
        python manage.py blast_utilityGCF_021461385.2_iqSchPice1.1_translated_cds.faa -p 

    * display in UI: 
        python manage.py blast_shown GCF_021461385.2_iqSchPice1.1_genomic.fna --shown true 
        python manage.py blast_shown GCF_021461385.2_iqSchPice1.1_rna_from_genomic.fna --shown true 
        python manage.py blast_shown GCF_021461385.2_iqSchPice1.1_cds_from_genomic.fna --shown true 
        python manage.py blast_shown GCF_021461385.2_iqSchPice1.1_translated_cds.faa --shown true 

1. For the genome assembly fasta file:
python manage.py addjbrowse GCF_021461385.2_iqSchPice1.1_genomic.fna https://apollo.nal.usda.gov/apollo/Schistocerca_piceifrons/jbrowse/

1. Set up hmmer for peptide fasta file: 
python manage.py addhmmer Schistocerca piceifrons -f /usr/local/i5k/media/blast/db/GCF_021461385.2_iqSchPice1.1_translated_cds.faa -d 'Schistocerca piceifrons NCBI Annotation release [release number], translated CDS'
amcooksey commented 2 years ago

Added organism page. Could not add photo. When I try to upload I get this error: "AJAX HTTP request terminated abnormally" Image is only 1.8 MB.

mpoelchau commented 2 years ago

@amcooksey for S. gregaria I saved the file in png format and reduced the size, and was able to upload it. Could you try this for this species?

amcooksey commented 2 years ago

@mpoelchau That seems to have worked.

mpoelchau commented 2 years ago

Out - closing for now, will re-open once the RNA-Seq pipeline is updated.