Answers.Body: <p>The page you link to is a GenBank Flat file, the default format used by GenBank. So, to get that format, which includes the <code>join</code> lines, use:</p>
  
  <pre><code>./efetch -db nuccore -id CP003820.1 -format gb
  </code></pre>
  
  <p>The <code>join()</code> lines give the different ways in which a sequence can be built from this gene. For example, one mRNA sequence produced by the locus you mention (CP003820) is:</p>
  
  <pre><code>  mRNA            complement(join(100..1497,1556..1767,1823..2274,
                   2322..2782,2846..3126,3206..3890,3958..5263,5322..5422,
                   5494..&gt;5645))
  </code></pre>
  
  <p>The corresponding coding sequence is for the mRNA above is:</p>
  
  <pre><code> CDS             complement(join(168..1497,1556..1767,1823..2274,
                   2322..2782,2846..3126,3206..3890,3958..5263,5322..5422,
                   5494..5645))
  </code></pre>
  
  <p>These aren't necessarily splicing "modifications", and aren't "variants", they are simply the way that the exons are joined to make an mRNA. </p>
  
  <p>Note that since CP003820.1 is "Cryptococcus neoformans var. grubii H99 chromosome 1, complete sequence", there are multiple genes described in there and each will have its own (or multiple, if there is alternative splicing) <code>join</code> statements. Specifically, I see 1014 genes described in the file:</p>
  
  <pre><code>$ grep -Pc '^\s+gene\s+' CP003820.gb 
  1014
  </code></pre>
  
  <p>950 mRNAs:</p>
  
  <pre><code>$ grep -c 'mRNA.*join(' CP003820.gb 
  950
  </code></pre>
  
  <p>And 955 proteins:</p>
  
  <pre><code>$ grep -Pc 'protein_id=' CP003820.gb 
  955
  </code></pre>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2019-04-26T12:14:57.327Z
Answers.LastActivityDate: 2019-04-26T12:14:57.327Z
Answers.UserDisplayName: terdon
Body: <p>I am trying to fetch genome assemblies that include splicing information.<br>
  E.g. <a href="https://www.ncbi.nlm.nih.gov/nuccore/CP003820.1" rel="nofollow noreferrer">Here</a>, is a genome assembly of a yeast, where genes show variants like this: <code>join(5928..7626,7685..7785,7857..7982)</code>.</p>
  
  <p>I am retrieving the information from NCBI with the following commands:<br>
  <code>efetch -db nuccore -id CP003820.1 -format tabular</code>. I tried different formats and it changes but I don't get the <code>join(...)</code></p>
  
  <p>1) Am I correct and the <code>join(...)</code> indicates "splicing-like" modifications in yeast?<br>
  2) How can I get the <code>join(...)</code> (any other format is fine too) with efetch?</p>
  
Tags: <phylogenetics><data-retrieval>
Title: NCBI retrieve genome assembly with splicing or "modifications"
ThreadText: # NCBI retrieve genome assembly with splicing or "modifications"
  
  (Asked by: username_0 on 2019-04-26)
  
  I am trying to fetch genome assemblies that include splicing information.  
  E.g. [Here](https://www.ncbi.nlm.nih.gov/nuccore/CP003820.1), is a genome assembly of a yeast, where genes show variants like this: `join(5928..7626,7685..7785,7857..7982)`.
  
  I am retrieving the information from NCBI with the following commands:  
  `efetch -db nuccore -id CP003820.1 -format tabular`. I tried different formats and it changes but I don't get the `join(...)`
  
  1) Am I correct and the `join(...)` indicates "splicing-like" modifications in yeast?  
  2) How can I get the `join(...)` (any other format is fine too) with efetch?
  
  **Tags:** phylogenetics, data-retrieval
  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2019-04-26. Score: 3 | Accepted answer
  
  The page you link to is a GenBank Flat file, the default format used by GenBank. So, to get that format, which includes the `join` lines, use:
  
  ```
  ./efetch -db nuccore -id CP003820.1 -format gb
  
  ```
  
  The `join()` lines give the different ways in which a sequence can be built from this gene. For example, one mRNA sequence produced by the locus you mention (CP003820) is:
  
  ```
    mRNA            complement(join(100..1497,1556..1767,1823..2274,
                   2322..2782,2846..3126,3206..3890,3958..5263,5322..5422,
                   5494..>5645))
  
  ```
  
  The corresponding coding sequence is for the mRNA above is:
  
  ```
   CDS             complement(join(168..1497,1556..1767,1823..2274,
                   2322..2782,2846..3126,3206..3890,3958..5263,5322..5422,
                   5494..5645))
  
  ```
  
  These aren't necessarily splicing "modifications", and aren't "variants", they are simply the way that the exons are joined to make an mRNA.
  
  Note that since CP003820.1 is "Cryptococcus neoformans var. grubii H99 chromosome 1, complete sequence", there are multiple genes described in there and each will have its own (or multiple, if there is alternative splicing) `join` statements. Specifically, I see 1014 genes described in the file:
  
  ```
  $ grep -Pc '^\s+gene\s+' CP003820.gb 
  1014
  
  ```
  
  950 mRNAs:
  
  ```
  $ grep -c 'mRNA.*join(' CP003820.gb 
  950
  
  ```
  
  And 955 proteins:
  
  ```
  $ grep -Pc 'protein_id=' CP003820.gb 
  955
  
  ```
  
  
  ---
  
  