Repository 'kallisto_quant'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/kallisto_quant

Changeset 9:2568a3b975cb (2023-05-31)
Previous changeset 8:c971db6f0fe5 (2021-07-18) Next changeset 10:4f9c4e6566e5 (2023-06-05)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/kallisto/ commit 168993a4e148506b1d3998c536caa2e501c36ccf
modified:
kallisto_quant.xml
macros.xml
test-data/kallisto_quant_out4.bam
test-data/kallisto_quant_out5.bam
added:
test-data/annotation.gtf.gz
test-data/chromosome_size.tabular
test-data/kallisto_quant_out7.tab
test-data/reads_forward.fastq.gz
test-data/reads_reverse.fastq.gz
test-data/transcripts.fasta
b
diff -r c971db6f0fe5 -r 2568a3b975cb kallisto_quant.xml
--- a/kallisto_quant.xml Sun Jul 18 17:54:19 2021 +0000
+++ b/kallisto_quant.xml Wed May 31 20:09:49 2023 +0000
[
@@ -1,9 +1,11 @@
-<?xml version="1.0"?>
 <tool id="kallisto_quant" name="Kallisto quant" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="20.05">
-    <description>- quantify abundances of RNA-Seq transcripts</description>
+    <description>quantify abundances of RNA-Seq transcripts</description>
     <macros>
         <import>macros.xml</import>
     </macros>
+    <xrefs>
+        <xref type="bio.tools">kallisto</xref>
+    </xrefs>
     <expand macro="requirements" />
     <command detect_errors="exit_code">
         <![CDATA[
@@ -15,7 +17,7 @@
             #set index_path = $reference_transcriptome.index.fields.path
         #end if
         kallisto quant -i '$index_path'
-            $bias --bootstrap-samples $bootstrap_samples --seed $seed $fusion $lib_type $pseudobam
+            $bias --bootstrap-samples $bootstrap_samples --seed $seed $fusion  $pseudobam
             #if $pseudobam:
                 -o .
             #else:
@@ -29,6 +31,7 @@
                 --sd $single_paired.sd
                 '$single_reads'
             #else:
+                $single_paired.lib_type
                 #if str($single_paired.collection.collection_selector) == 'datasets':
                     #set $forward_reads = $single_paired.collection.forward
                     #set $reverse_reads = $single_paired.collection.reverse
@@ -39,6 +42,12 @@
                 #set $reads = "'%s' '%s'" % ($forward_reads, $reverse_reads)
                 $reads
             #end if
+            $single_overhang
+            #if $genomebam_option.selector
+                $genomebam_option.selector
+                --gtf $genomebam_option.gtf
+                --chromosomes $genomebam_option.chromosomes
+            #end if
             #if $pseudobam:
                 && samtools sort --no-PG -O bam -@ \${GALAXY_SLOTS:-1} -T "\${TMPDIR:-.}" -o '$pseudobam_output' pseudoalignments.bam
             #end if
@@ -54,7 +63,7 @@
             </param>
             <when value="single">
                 <param name="reads" type="data" format="fastq,fastq.gz" label="Reads in FASTQ format" />
-                <param name="fragment_length" argument="--fragment-length" type="integer" value="200" label="Average fragment length" help="Illumina typically produces reads of 180-200bp" />
+                <param argument="--fragment-length" type="integer" value="200" label="Average fragment length" help="Illumina typically produces reads of 180-200bp" />
                 <param argument="--sd" type="integer" value="20" label="Estimated standard deviation of fragment length" />
             </when>
             <when value="paired">
@@ -71,18 +80,34 @@
                         <param name="reads" type="data_collection" format="fastq,fastq.gz" collection_type="paired" label="Collection of reads" />
                     </when>
                 </conditional>
+                <param name="lib_type" type="select" label="Library strandness information">
+                    <option value="">Unstranded</option>
+                    <option value="--fr-stranded">Strand specific reads, first read forward</option>
+                    <option value="--rf-stranded">Strand specific reads, first read reverse</option>
+                </param>
             </when>
         </conditional>
-        <param argument="--bias" type="boolean" truevalue="--bias" falsevalue="" label="Perform sequence based bias correction" />
-        <param name="bootstrap_samples" argument="--bootstrap-samples" type="integer" value="0" label="Number of bootstrap samples" help="default: 0" />
-        <param argument="--seed" type="integer" value="42" label="Seed for the bootstrap sampling" help="default: 42" />
-        <param argument="--fusion" type="boolean" truevalue="--fusion" falsevalue="" label="Search for fusions" help="for Pizzly" />
-        <param name="lib_type" type="select" label="Library strandness information">
-            <option value="">Unstranded</option>
-            <option value="--fr-stranded">Strand specific reads, first read forward</option>
-            <option value="--rf-stranded">Strand specific reads, first read reverse</option>
-        </param>
+        <param argument="--bias" type="boolean" truevalue="--bias" falsevalue="" label="Perform sequence based bias correction" help="It allows to learn 
+            parameters for a model of sequences specific bias and corrects the abundances accordlingly"/>
+        <param argument="--bootstrap-samples" type="integer" value="0" label="Number of bootstrap samples" help="Running with bootstraps 
+            is mandatory if you want to perform differential expression analysis of isoforms with Sleuth.Default: 0" />
+        <param argument="--fusion" type="boolean" truevalue="--fusion" falsevalue="" label="Search for fusions" help="It generates the required files for Pizzly. This option does normal quantification, but 
+            additionally looks for reads that do not pseudoalign because they are potentially from fusion genes." />
+        <param argument="--single-overhang" type="boolean" truevalue="--single-overhang" falsevalue="" checked="false" label="Single overhang" help="Include reads where 
+            unobserved rest of fragment is predicted to lie outside a transcript" />
         <param argument="--pseudobam" type="boolean" truevalue="--pseudobam" falsevalue="" label="Output pseudoalignments in BAM format" />
+        <conditional name="genomebam_option">
+            <param name="selector" type="select" label="Project pseudoalignments to genome">
+                <option value="--genomebam">Enabled</option>
+                <option value="" selected="true">Disabled</option>
+            </param>
+            <when value="--genomebam">
+                <param argument="--gtf" type="data" format="gtf" label="GTF file" help="GTF file for transcriptome information" />
+                <param argument="--chromosomes" type="data" format="tabular" label="Chromosome names and lengths"/>
+            </when>
+            <when value=""/>
+        </conditional>
+        <param argument="--seed" type="integer" value="42" label="Seed for the bootstrap sampling" help="Default: 42" />
     </inputs>
     <outputs>
         <data format="h5" name="abundance_h5" from_work_dir="abundance.h5" label="${tool.name} on ${on_string}: Abundances (HDF5)" />
@@ -170,7 +195,25 @@
             <param name="reverse" ftype="fastq" dbkey="hg38" value="hg38_R.fq.gz" />
             <output name="fusion_output" file="fusion.txt" ftype="tabular" />
         </test>
-
+        <test>
+            <param name="reference_transcriptome_source" value="history" />
+            <param name="reference" ftype="fasta" value="transcripts.fasta" />
+            <param name="single_paired_selector" value="paired" />
+            <param name="collection_selector" value="datasets" />
+            <param name="forward" ftype="fastq" value="reads_forward.fastq.gz" />
+            <param name="reverse" ftype="fastq" value="reads_reverse.fastq.gz" />
+            <conditional name="genomebam_option">
+                <param name="selector" value="--genomebam"/>
+                <param name="gtf" value="annotation.gtf.gz"/>
+                <param name="chromosomes" value="chromosome_size.tabular"/>
+            </conditional>
+            <output name="abundance_tab" file="kallisto_quant_out7.tab" ftype="tabular" />
+            <assert_command>
+                <has_text text="--genomebam" />
+                <has_text text="--chromosomes" />
+                <has_text text="--gtf" />
+            </assert_command>
+        </test>
     </tests>
     <help>
  <![CDATA[
b
diff -r c971db6f0fe5 -r 2568a3b975cb macros.xml
--- a/macros.xml Sun Jul 18 17:54:19 2021 +0000
+++ b/macros.xml Wed May 31 20:09:49 2023 +0000
b
@@ -1,11 +1,11 @@
 <?xml version="1.0"?>
 <macros>
-    <token name="@TOOL_VERSION@">0.46.2</token>
+    <token name="@TOOL_VERSION@">0.48.0</token>
     <token name="@VERSION_SUFFIX@">0</token>
     <xml name="requirements">
         <requirements>
             <requirement type="package" version="@TOOL_VERSION@">kallisto</requirement>
-            <requirement type="package" version="1.12">samtools</requirement>
+            <requirement type="package" version="1.14">samtools</requirement>
         </requirements>
     </xml>
     <xml name="citations">
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/annotation.gtf.gz
b
Binary file test-data/annotation.gtf.gz has changed
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/chromosome_size.tabular
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/chromosome_size.tabular Wed May 31 20:09:49 2023 +0000
b
@@ -0,0 +1,1 @@
+chr1    24999771
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/kallisto_quant_out4.bam
b
Binary file test-data/kallisto_quant_out4.bam has changed
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/kallisto_quant_out5.bam
b
Binary file test-data/kallisto_quant_out5.bam has changed
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/kallisto_quant_out7.tab
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/kallisto_quant_out7.tab Wed May 31 20:09:49 2023 +0000
b
@@ -0,0 +1,16 @@
+target_id length eff_length est_counts tpm
+NM_001168316 2283 2105.9 160.693 12587.6
+NM_174914 2385 2207.9 1500.61 112117
+NR_031764 1853 1675.9 102.695 10108.4
+NM_004503 1681 1503.9 330.216 36221.2
+NM_006897 1541 1363.9 664 80309.8
+NM_014212 2037 1859.9 55 4878.16
+NM_014620 2300 2122.9 589.807 45831.5
+NM_017409 1959 1781.9 47 4351.09
+NM_017410 2396 2218.9 42 3122.45
+NM_018953 1612 1434.9 227.994 26211.1
+NM_022658 2288 2110.9 4881 381438
+NM_153633 1666 1488.9 362.124 40121.2
+NM_153693 2072 1894.9 74.8502 6516.13
+NM_173860 849 671.903 962 236185
+NR_003084 1640 1462.9 0.00847856 0.956071
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/reads_forward.fastq.gz
b
Binary file test-data/reads_forward.fastq.gz has changed
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/reads_reverse.fastq.gz
b
Binary file test-data/reads_reverse.fastq.gz has changed
b
diff -r c971db6f0fe5 -r 2568a3b975cb test-data/transcripts.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/transcripts.fasta Wed May 31 20:09:49 2023 +0000
b
b'@@ -0,0 +1,498 @@\n+>NM_001168316\n+TTAGTGAGGTTGGGGAGAGATAACGCTGTAAACTTTTATTTTTCAGGAAATCTGGAAACC\n+TACAGTCTCCAAGCCTGCTCAGCCAAGAAGGAGCTCACTGTGGGCACCAGAGACAGGGAC\n+CCAATGTGGAGACCTGTGAGCCTGTGTCCGGCCCTGAACTCTCAAGCACAGGGCAGGCTT\n+CCTGAGCATTGAAGAGAATATGTGGGAGAACAAAACAGAAACTGAAAGAATATGCAAGGT\n+GTCTTTCTTGGATGTTATTCCATGATAGATAGTAGGGGCAGGAGTGAGAGAGGCTGACTA\n+GGTCTGGACATGGAGGCTGGAAGAGTCAGGGTGTGATTCGGAGAGGCGCATGAGAAGGAA\n+GGTGGATTTTAAGGCTGGAAATCTGAGGGTCAGTGGTCCAAGTCACTCAGAGACAGAATC\n+ACAGCATAGCCCTTGCTGATGGCAAACAAAGGAGGACAAGAGGACTGGAAAGAATTCTGC\n+TAGCAGGCAGGAGCTAGTAAGGATGAATTTGTAGCAAAATTAGCAAGTGGAAAGGATGAT\n+TTTTGGCCATTTTTCCTGTTCTTCAAGAAAACAGGAGATAACTAGAAGGACTAGAGAATG\n+GGGCTGCCAGAACTAGTGGGAAGCTCCCTAGAAATGGTGACATCGCCCACCAAACAGACC\n+CCGCCAAGGCTGCACCTGGCCTTATGTCTCCTTCACCTTTCTGGCCCCACGCAGCCACCA\n+GACAGCCATGCCCAGCAGCTTCCCACAAAGCCATAGAGTCCCCAGAGTGAGCCCCAGCAG\n+AAACACAAAAACGTCGAGCAGGTACTGCTCATGCCAGGGCTGCTGAAAGACATAGGGCTT\n+GAGGTGCGTCGCGCCCCCTGTCTGGAGGACGTGGTCAATCCAGCCCACCAGCCGCTGTGT\n+GGGGCTGAGCGGGTGGGAGCGCAGGATGACACTGGCAGCCACTGCCGCGGACTTGTATCT\n+CTTGTCTTCCATGATTTGTTTCATCTTAAGAGCCAATGTCTCTGCCTTGAGCTTCTTTAA\n+CTGAATAGAAACACCAAACTTTTTGGCTTCTACTCGGACCATGTTTTCAGGCTGGTCTCC\n+AAAGAGAGGGATCCCCACCATGGGCACACCATGCTGGATGGCCTCCATTATGCTATTCTG\n+CCCGCCGTGGGTGACAAACAGACGGATGCTTGGGTGAGCCAGGAGGTCACTCTGAGGAAG\n+CCAGTCCACAATTTTCACATTTGCAGCCAGGTGGACATCTTTGGGCCAATGAGAACACTG\n+ACACTTCCATATCACCCCTTGGGGTAGGTGAGCAAAGGCATTGTTCATCTCCTTGAAGAT\n+TTCCGGATTCTGACAGGTGTTCACCATGGAGCCCAAGGTCACAAGGACAAAACCAGAGTC\n+CCCAAACTTGGCAATGAAGTTCTCCAAGTCTTGTGGTACTGGTTTAATAGGTTTTTCCAT\n+CAAGCCTCCAACATAAACAGTGTTGGGAAGCAGAGGTCGAGCAAAATCAAAGGCAAAGTC\n+AGAGTTAATGAACCACAACTCTGCTTTCAGTAGAAGATGAGACAAAACTGGCCTAGAGCC\n+TTCTGTGAAATGTTCCTTGATGGTGTTGTCAAATGTAGACTGCATGTGCTGTTGCCTCCT\n+GCAGAAACTAAAGAACATCAGAAAATTCTTCACTCGGCCCCAGAAGTCCATGTGATCAGT\n+CAGCAAGGAACGGAATACTGGAACATAAGACAAGGGGATTGGTAGCCCAAATTCCAAAGA\n+GCCGAATGAAGTGGAAAGAATGGCCACAAATGGCTTCCCAAGCTTCTCAGCAATCAGGAA\n+AGGACAGTAGTCAAAAGTTTCAACTATCACCATGTCGAAGTTCTCATTCTTTAAGGAATC\n+CATGATATCCTTTCTATTTAAAAAATGACTGCACTGCAACGCCAAGTATTCTAGAACATT\n+TAATAAGTTTTCAAATTTTCCTCTGCCACCTAAAGTTTCTTCCAGAAAGAAATCAAAACT\n+CTTTTTAAATTCTCTTTGATGATCTTCAGGTGCAAGCCAACTGATAACTTGATATGATTT\n+TTCTTCCTTTTTAAAATCTACTGTAGATATTGTCAGGATTTTGGCAGCCTCTGAGAGCAG\n+GACCCCAGGGAGAAGGAAGCCCACTAGAAGAAGCACTCGCTGCCCAGCCATGCTCACTTC\n+TACGGAAGCCGCGGATCTCAGCCTGGGCTGCGCGCCCTGCGCCCGGCTAAGGGACCCTGT\n+GCACCTCAGTGCGCCAAAGGCACTGGCTGTGGGTAGAGGTAGGAGTAAGGATCCTGCCAG\n+TCT\n+>NM_174914\n+TTAGTGAGGTTGGGGAGAGATAACGCTGTAAACTTTTATTTTTCAGGAAATCTGGAAACC\n+TACAGTCTCCAAGCCTGCTCAGCCAAGAAGGAGCTCACTGTGGGCACCAGAGACAGGGAC\n+CCAATGTGGAGACCTGTGAGCCTGTGTCCGGCCCTGAACTCTCAAGCACAGGGCAGGCTT\n+CCTGAGCATTGAAGAGAATATGTGGGAGAACAAAACAGAAACTGAAAGAATATGCAAGGT\n+GTCTTTCTTGGATGTTATTCCATGATAGATAGTAGGGGCAGGAGTGAGAGAGGCTGACTA\n+GGTCTGGACATGGAGGCTGGAAGAGTCAGGGTGTGATTCGGAGAGGCGCATGAGAAGGAA\n+GGTGGATTTTAAGGCTGGAAATCTGAGGGTCAGTGGTCCAAGTCACTCAGAGACAGAATC\n+ACAGCATAGCCCTTGCTGATGGCAAACAAAGGAGGACAAGAGGACTGGAAAGAATTCTGC\n+TAGCAGGCAGGAGCTAGTAAGGATGAATTTGTAGCAAAATTAGCAAGTGGAAAGGATGAT\n+TTTTGGCCATTTTTCCTGTTCTTCAAGAAAACAGGAGATAACTAGAAGGACTAGAGAATG\n+GGGCTGCCAGAACTAGTGGGAAGCTCCCTAGAAATGGTGACATCGCCCACCAAACAGACC\n+CCGCCAAGGCTGCACCTGGCCTTATGTCTCCTTCACCTTTCTGGCCCCACGCAGCCACCA\n+GACAGCCATGCCCAGCAGCTTCCCACAAAGCCATAGAGTCCCCAGAGTGAGCCCCAGCAG\n+AAACACAAAAACGTCGAGCAGGTACTGCTCATGCCAGGGCTGCTGAAAGACATAGGGCTT\n+GAGGTGCGTCGCGCCCCCTGTCTGGAGGACGTGGTCAATCCAGCCCACCAGCCGCTGTGT\n+GGGGCTGAGCGGGTGGGAGCGCAGGATGACACTGGCAGCCACTGCCGCGGACTTGTATCT\n+CTTGTCTTCCATGATTTGTTTCATCTTAAGAGCCAATGTCTCTGCCTTGAGCTTCTTTAA\n+CTGAATAGAAACACCAAACTTTTTGGCTTCTACTCGGACCATGTTTTCAGGCTGGTCTCC\n+AAAGAGAGGGATCCCCACCATGGGCACACCATGCTGGATGGCCTCCATTATGCTATTCTG\n+CCCGCCGTGGGTGACAAACAGACGGATGCTTGGGTGAGCCAGGAGGTCACTCTGAGGAAG\n+CCAGTCCACAATTTTCACATTTGCAGCCAGGTGGACATCTTTGGGCCAATGAGAACACTG\n+ACACTTCCATATCACCCCTTGGGGTAGGTGAGCAAAGGCATTGTTCATCTCCTTGAAGAT\n+TTCCGGATTCTGACAGGTGTTCACCATGGAGCCCAAGGTCACAAGGACAAAACCAGAGTC\n+CCCAAACTTGGCAATGAAGTTCTCCAAGTCTTGTGGTACTGGTTTAATAGGTTTTTCCAT\n+CAAGCCTCCAACATAAACAGTGTTGGGAAGCAGAGGTCGAGCAAAATCAAAGGCAAAGTC\n+AGAGTTAATGAACCACAACTCTGCTTTCAGTAGAAGATGAGA'..b'G\n+GATCTAATTCCTTTTACCAGGAGAAAGACATGCTCTCAAACTGCAGACAAAACACCTTAG\n+GACATAACACACAGACCTCAATCGCTCAGGATTTTAGTTCTGAGCAGGGCAGGACTGCGC\n+CCCAGGACCAGAAAGCCAGTATCCAGATTTACCCCTGGATGCAGCGAATGAATTCGCACA\n+GTGGGGTCGGCTACGGAGCGGACCGGAGGCGCGGCCGCCAGATCTACTCGCGGTACCAGA\n+CCCTGGAACTGGAGAAGGAATTTCACTTCAATCGCTACCTAACGCGGCGCCGGCGCATCG\n+AGATCGCCAACGCGCTTTGCCTGACCGAGCGACAGATCAAAATCTGGTTCCAGAACCGCC\n+GGATGAAGTGGAAAAAAGAATCTAATCTCACATCCACTCTCTCGGGGGGCGGCGGAGGGG\n+CCACCGCCGACAGCCTGGGCGGAAAAGAGGAAAAGCGGGAAGAGACAGAAGAGGAGAAGC\n+AGAAAGAGTGACCAGGACTGTCCCTGCCACCCCTCTCTCCCTTTCTCCCTCGCTCCCCAC\n+CAACTCTCCCCTAATCACACACTCTGTATTTATCACTGGCACAATTGATGTGTTTTGATT\n+CCCTAAAACAAAATTAGGGAGTCAAACGTGGACCTGAAAGTCAGCTCTGGACCCCCTCCC\n+TCACCGCACAACTCTCTTTCACCACGCGCCTCCTCCTCCTCGCTCCCTTGCTAGCTCGTT\n+CTCGGCTTGTCTACAGGCCCTTTTCCCCGTCCAGGCCTTGGGGGCTCGGACCCTGAACTC\n+AGACTCTACAGATTGCCCTCCAAGTGAGGACTTGGCTCCCCCACTCCTTCGACGCCCCCA\n+CCCCCGCCCCCCGTGCAGAGAGCCGGCTCCTGGGCCTGCTGGGGCCTCTGCTCCAGGGCC\n+TCAGGGCCCGGCCTGGCAGCCGGGGAGGGCCGGAGGCCCAAGGAGGGCGCGCCTTGGCCC\n+CACACCAACCCCCAGGGCCTCCCCGCAGTCCCTGCCTAGCCCCTCTGCCCCAGCAAATGC\n+CCAGCCCAGGCAAATTGTATTTAAAGAATCCTGGGGGTCATTATGGCATTTTACAAACTG\n+TGACCGTTTCTGTGTGAAGATTTTTAGCTGTATTTGTGGTCTCTGTATTTATATTTATGT\n+TTAGCACCGTCAGTGTTCCTATCCAATTTCAAAAAAGGAAAAAAAAGAGGGAAAATTACA\n+AAAAGAGAGAAAAAAAGTGAATGACGTTTGTTTAGCCAGTAGGAGAAAATAAATAAATAA\n+ATAAATCCCTTCGTGTTACCCTCCTGTATAAATCCAACCTCTGGGTCCGTTCTCGAATAT\n+TTAATAAAACTGATATTATTTTTAAAACTTTA\n+>NM_173860\n+ATGGGCGAGCATAATCTCCTGAATCCCGGGTTTGTGGGGCCGCTGGTAAACATCCACACG\n+GGAGACACCTTCTACTTCCCCAACTTCCGCGCGTCCGGGGCGCAGCTTCCCGGGCTGCCT\n+TCGCTGTCCTACCCACGCCGCGACAACGTGTGCTCCCTGTCCTGGCCGTCGGCGGAGCCG\n+TGCAATGGCTACCCGCAGCCCTACCTCGGCAGCCCAGTGTCTCTCAACCCTCCCTTCGGC\n+CGCACGTGCGAGCTGGCGCGCGTGGAGGACGGCAAGGGTTACTACCGCGAGCCGTGCGCC\n+GAGGGTGGCGGCGGGGGCCTGAAGCGTGAGGAGCGCGGGCGCGACCCGGGAGCCGGGCCC\n+GGGGCAGCGCTGCTCCCGCTGGAGCCGTCGGGGCCGCCTGCGCTCGGCTTCAAGTACGAC\n+TACGCGGCGGGCGGCGGCGGTGGCGACGGCGGCGGCGGCGCAGGACCTCCGCACGACCCG\n+CCCTCCTGCCAGTCGCTGGAATCCGACTCCAGTTCGTCCCTGCTCAACGAGGGCAACAAG\n+GGCGCCGGCGCAGGCGACCCCGGCAGCTTGGTATCGCCGTTGAACCCCGGCGGCGGGCTC\n+TCGGCCAGCGGCGCGCCCTGGTACCCGATCAACAGCCGCTCTCGGAAGAAGCGCAAGCCC\n+TATTCGAAGTTGCAACTGGCAGAGCTGGAGGGCGAGTTTCTGGTCAACGAGTTCATCACA\n+CGCCAGCGCCGGAGGGAACTCTCAGACCGCTTGAATCTTAGTGACCAGCAGGTCAAGATC\n+TGGTTTCAGAACCGGAGAATGAAAAAGAAAAGACTTCTGTTGAGGGAGCAAGCTCTCTCC\n+TTCTTTTAA\n+>NR_003084\n+TTATTGTGGTTTGTCCGTTCCGAGCGCTCCGCAGAACAGTCCTCCCTGTAAGAGCCTAAC\n+CATTGCCAGGGAAACCTGCCCTGGGCGCTCCCTTCATTAGCAGTATTTTTTTTAAATTAA\n+TCTGATTAATAATTATTTTTCCCCCATTTAATTTTTTTTCCTCCCAGGTGGAGTTGCCGA\n+AGCTGGGGGCAGCTGGGGAGGGTGGGGATGGGAGGGGAGAGACAGAAGTTGAGGGCATCT\n+CTCTCTTCCTTCCCGACCCTCTGGCCCCCAAGGGGCAGGAGGAATGCAGGAGCAGGAGTT\n+GAGCTTGGGAGCTGCAGATGCCTCCGCCCCTCCTCTCTCCCAGGCTCTTCCTCCTGCCCC\n+CTTCTTGCAACTCTCCTTAATTTTGTTTGGCTTTTGGATGATTATAATTATTTTTATTTT\n+TGAATTTATATAAAGTATATGTGTGTGTGTGTGGAGCTGAGACAGGCTCGGCAGCGGCAC\n+AGAATGAGGGAAGACGAGAAAGAGAGTGGGAGAGAGAGAGGCAGAGAGGGAGAGAGGGAG\n+AGTGACAGCAGCGCTCGAGACGGACGGCAAGCGGTCCCGAACCAGTTACACGCGCTACCA\n+GACTCTGGAACTCGAGAAAGAATTCCACTTTAACCGCTACCTCACTCGCCGCAGGCGCAT\n+AGAGATCGCCAACAACTTGTGTCTCAATGAGAGACAGATCAAGATCTGGTTCCAGAACCG\n+CAGGATGAAGTGGAAGAAAGATTCCAAAATGAAAAGCAAAGAGGCTCTTTAGAGGCAGCG\n+GGGGAGGCCCGCAGAGCGCGCCCCTAGCCGGTTCCTGTCCCTGCGCCTTTCCTTTTCGCC\n+TTTCCTCTCTATATTTCGGGTCGGGGGCAGGTGCTGGAGCACTGGGCTCCCGGGCCCCAC\n+AGACAAAAGCGCTTTTCCTTGGCATTCCGCATCCCTACCGACCCAGGGTTCCCGCGGGGC\n+TGTCGGCGCTGCCCCATCTCCCCTCAGCTCGGCTCAGCTCGGTACCCGGGGCCCAGGGCA\n+AGCTCCGCAGGACTTCCCCGGAGGGCTGCGGCGTACAGGCTGGCGCAGAACGAACCTTGG\n+CCTGGGCCGTATCTCCGGCTCCCAGCCTCAGCGCGGCCCTCCCGAGTTAAGGTGGGCCCG\n+GCCCGCGCCACAGGACCCTCGCCGGACCCTCTAACCTCGCCCTCTCCTTTGTTCCCGGCT\n+GGACGGGTTAGACAGCCAAAGGCTGGCGAGAGTCTGGCCCTAGACTCGGGGTGCTTCCTT\n+GTAGCGACTAAACTAGATTTTCACTTATGAATGATTTGCATATGAAAGGAGAGCATCGGC\n+CTAGGGCCCCCACAGTTGCTCTATGCTTTCCAAACCTTATCTCCACAACCTCTTCCCCCC\n+AAAACCCGGGAACCTCCCCAGCCTGCGCCTGCTGCATGCCCTCTCAGGCCGGCAGCCCCA\n+GCCTGCTAGCTAGCTCAACTAGTGGGGTTTCCTGGCACTGGACCCCAGCAAGTGGTCCTA\n+GAGGCCCTTTGCTGTCCCATAGTCCCTGCCACGAATTTCTGTGCCCTCCTGACCCATTGC\n+TGTTGTCCAACTATTTATTGACTCTGGGTCCTTCCTGAAACTATATTTTGTCATATCAAA\n+TAAAGAGAGAACAGGACTAA\n'