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' |