|
Commit message:
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rnasnp commit 6e51a331dbeab2786d8df5fd379ae3a63eb61d83 |
|
added:
rnasnp.xml test-data/seq1.txt test-data/seq1_snp1_mode1.out test-data/seq2.txt test-data/seq2_snp2_mode2.out test-data/seq3.txt test-data/seq3_mode3.out test-data/snp1.txt test-data/snp2.txt test-data/snp3.txt |
| b |
| diff -r 000000000000 -r 327f5721612f rnasnp.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/rnasnp.xml Tue Dec 13 12:31:56 2016 -0500 |
| [ |
| b'@@ -0,0 +1,176 @@\n+<tool id="rnasnp" name="RNAsnp" version="1.2.0">\n+ <description> Efficient detection of local RNA secondary structure changes induced by SNPs </description>\n+\n+ <requirements>\n+ <requirement type="package" version="1.2">rnasnp</requirement>\n+ </requirements>\n+\n+ <stdio>\n+ <exit_code range="1" level="fatal" description="Error occurred. Please check Tool Standard Error" />\n+ <exit_code range=":-1" level="fatal" description="Error occurred. Please check Tool Standard Error" />\n+ <exit_code range="2:" level="warning" description="Warning: Tool exited with error_code != 0" />\n+ </stdio>\n+\n+ <version_command>RNAsnp --version</version_command>\n+ <command>\n+<![CDATA[\n+\n+ RNAsnp\n+ --seq=\'$seq\'\n+ --mode=$mode.mode_selector\n+ #if $mode.mode_selector == "1"\n+ --snp=$mode.snp\n+ --cutoff=$mode.cutoff\n+ --minLen=$mode.minLen\n+ #end if\n+ #if $mode.mode_selector == "2"\n+ --snp=$mode.snp\n+ --winsize=$mode.winsize\n+ --span=$mode.span\n+ --regionX=$mode.regionX\n+ --regionY=$mode.regionY\n+ #end if\n+ #if $mode.mode_selector == "3"\n+ --pvalue1=$mode.pvalue1\n+ --pvalue2=$mode.pvalue2\n+ --winsizeExt=$mode.winsizeExt\n+ #end if\n+ --edist=$edist\n+ --boltzmannPreFactor=$boltzmannPreFactor\n+ --winsizeFold=$mode.winsizeFold\n+ > \'$outfile\'\n+\n+]]>\n+ </command>\n+ <inputs>\n+ <param argument="--seq" format="txt" type="data" label="Input Alignment File" />\n+ <conditional name="mode">\n+ <param name="mode_selector" type="select" label="Mode of operation">\n+ <option value="1" selected="true">1: Perform global folding by using RNAfold</option>\n+ <option value="2">2: Perform local folding by using RNAplfold</option>\n+ <option value="3">3: Screen putative structure-disruptive SNPs in an RNA sequence</option>\n+ </param>\n+ <when value="1">\n+ <param argument="--snp" format="txt" type="data" label="List of SNP" help="The list of SNPs to be tested have to be provided in separate lines"/>\n+ <param argument="--winsizeFold" type="select" label="Folding windows size" help="Length of flanking sequence on either side of SNP considered for folding. please make sure that the input sequence is at least twice the size of chosen flanking. This restriction is necessary to keep the size of parameter tables for the p-value calculations manageable. In case the input sequence is less than twice the size of chosen flanking, the RNAsnp takes the nts up to the start and end position of the given sequence from the SNP position and perform the analysis.">\n+ <option value="100"></option> <option value="150" ></option>\n+ <option value="200" selected ="true"></option> <option value="250"></option> <option value="300"></option> <option value="350"></option> <option value="400"></option> <option value="450"></option>\n+ <option value="500"></option> <option value="600"></option> <option value="650"></option> <option value="700"></option> <option value="750"></option> <option value="800"></option>\n+ </param>\n+\n+ <param argument="--cutoff" type="float" value="0.01" label="Cut-off" help="Cut-off for the base pair probabilities."/>\n+ <param argument="--minLen" type="integer" value="50" min="1" label="Minimum length of the sequence interval" help="The structural difference between wild-type and mutant is computed for all sequence intervals with the selected minimum length" />\n+ </when>\n+ <when value="2">\n+ <param argument="--snp" format="txt" type="data" label="List of SNP" help="The list of SNPs to be tested have to be provided in separate lines"/>\n+ <par'..b'bel="Windows size" help="size of the flanking region on either side of SNP that includes the local window returned by Mode 2. This subsequence is then passed to Mode 1 for re-computation"/>\n+ </when>\n+ </conditional>\n+ <param argument="--edist" type="integer" value="0" min="0" label="Ensemble Euclidean distance" help="compute ensemble Euclidean distance between the distribution of structures between two sequences"/>\n+ <param argument="--boltzmannPreFactor" type="float" value="1" label="Boltzmann PreFactor" help="Multiply the bolztmann factor with a prefactor alpha"/>\n+\n+ </inputs>\n+ <outputs>\n+ <data name="outfile" format="tabular" />\n+ </outputs>\n+ <tests>\n+ <test>\n+ <param name="seq" value="seq1.txt"/>\n+ <conditional name="mode">\n+ <param name="mode_selector" value="1"/>\n+ <param name="snp" value="snp1.txt"/>\n+ </conditional>\n+ <output name="outfile" file="seq1_snp1_mode1.out"/>\n+ </test>\n+ </tests>\n+\n+ <help>\n+<![CDATA[\n+\n+**What it does**\n+ Efficient detection of local RNA secondary structure changes induced by SNPs\n+\n+\n+ RNAsnp requires an RNA sequence and optionally a list of SNPs to be analyzed. The effect of SNPs on local RNA secondary structure can be detected in three possible modes,\n+\n+ Mode 1: The first mode is designed to compute the effect of SNPs by using global folding. This option should be used only for short input sequences, since the base pair probabilities are calculated using RNAfold.\n+\n+ Mode 2: The second mode is designed to compute the effect of SNPs on large sequences. Here the local base pair probabilities are calculated using RNAplfold (with the parameters -W 200 and -L 120).\n+\n+ Mode 3: The third mode is the combination of the above two. It is intended to determine the positions of putative structure-disruptive SNPs using either transcript or genome sequence.\n+\n+\n+**Input**\n+\n+Sequence file must contain one sequence (preferably in FASTA format).\n+\n+SNP file must contain the list of SNPs that are given in separate lines. The\n+SNPs are described as, wild-type nucleotide followed by nucleotide position\n+followed by mutant nucleotide. In case of multiple SNPs, the SNPs are delimited\n+by the special character "-".\n+\n+Example formats:\n+Single SNP: A201G\n+where, A is the wild-type nucleotide in the given sequence, 201 is the sequence\n+position of wild-type nucleotide and G is the mutant (or SNP).\n+\n+Multiple SNPs: A201G-U257A-C260G\n+The multiple SNPs (which occurs together) are defined next to each other with\n+the delimiter "-" between them.\n+\n+\n+**Output**\n+\n+SNP w\tSlen GC interval d pvalue1 ewin\tinterval d_max\tpvalue2\n+\n+G1C 200\t3344 0.5522 1-46 0.0421 0.0755 200 1-50\t0.1581\t0.0183\n+\n+G7A 200\t3344 0.5556 1-43 0.2236 0.0207 200\t 1-50\t0.1570\t0.0996\n+\n+\n+Output details:\n+\n+Column 1: details of the screened SNPs\n+\n+Column 2: length of the flanking region considered on either side of the SNP to fold\n+\n+Column 3: length of the given input sequence\n+\n+Column 4: GC percent of the sequence interval considered for folding\n+\n+Column 5-7: results of initial screen returned by mode 2\n+\n+Column 8: ewin is the length of the flanking region considered on either of the SNP that also includes the local interval returned by mode 2. This subsequence is then passed to mode 1 to compute the d_max\n+\n+Column 9-11: results of final screen returned by mode 1\n+\n+\n+Note: By default, the ewin(-e) uses the window length of 200. A higher computation speed is achieved if the ewin value is reduced to 100 or 150, however, there will be little difference in the result compare to ewin with 200. In addition, RNAsnp does automatically increase the length of ewin if the defined length is not enough to cover the local interval returned by mode 2.\n+\n+]]>\n+\n+ </help>\n+\n+ <citations>\n+ <citation type="doi">10.1002/humu.22273</citation>\n+ </citations>\n+\n+</tool>\n' |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/seq1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/seq1.txt Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,50 @@ +>gi|1015320|dbj|D32050.1|HUMATS Homo sapiens mRNA for alanyl-tRNA synthetase, complete cds +GGTACAGCTGCGCGTCTGCGGGAATAGGTGCAGCGGGCCCTTGGCGGGGGACTCTGAGGGAGGAGCTGGG +GACGGCGACCCTAGGAGAGTTCTTTGGGGTGACTTTCAAGATGGACTCTACTCTAACAGCAAGTGAAATC +CGGCAGCGATTTATAGATTTCTTCAAGAGGAACGAGCATACGTATGTTCACTCGTCTGCCACCATCCCAT +TGGATGACCCCACTTTGCTCTTTGCCAATGCAGGCATGAACCAGTTTAAACCCATTTTCCTGAACACAAT +TGACCCATCTCACCCCATGGCAAAGCTGAGCAGAGCTGCCAATACCCAGAAGTGCATCCGGGCTGGGGGC +AAACAAAATGACCTGGACGATGTGGGCAAGGATGTCTATCATCACACCTTCTTCGAGATGCTGGGCTCTT +GGTCTTTTGGAGATTACTTTAAGGAATTGGCATGTAAGATGGCTCTGGAACTCCTCACCCAAGAGTTTGG +CATTCCCATTGAAAGACTTTATGTTACTTACTTTGGCGGGGATGAAGCAGCTGGCTTAGAAGCAGATCTG +GAATGCAAACAGATCTGGCAAAATTTGGGGCTGGATGACACCAAAATCCTCCCAGGCAACATGAAGGATA +ACTTCTGGGAGATGGGTGACACGGGCCCCTGTGGTCCTTGCAGTGAGATCCACTACGACCGGATTGGTGG +TCGGGACGCCGCACATCTTGTCAACCAGGACGACCCTAATGTGCTGGAGATCTGGAACCTTGTGTTCATC +CAGTATAACAGGGAAGCTGATGGCATTCTGAAACCTCTTCCCAAGAAAAGCATTGACACAGGGATGGGCC +TGGAACGACTGGTATCTGTGCTGCAGAATAAGATGTCCAACTATGACACTGACCTTTTTGTCCCTTACTT +TGAAGCCATTCAGAAGGGCACAGGTGCCCGACCATACACTGGGAAAGTTGGTGCTGAGGATGCCGATGGG +ATTGACATGGCCTACCGGGTGCTGGCTGACCATGCTCGGACCATCACTGTGGCACTGGCTGATGGTGGCC +GGCCTGACAACACAGGGCGTGGATATGTGTTGAGACGGATTCTCCGCCGAGCTGTCCGATACGCCCATGA +AAAGCTCAATGCCAGCAGGGGCTTCTTTGCTACGTTAGTGGATGTTGTCGTCCAGTCCCTGGGAGATGCA +TTTCCTGAGCTGAAGAAGGACCCAGACATGGTGAAGGACATCATTAATGAAGAAGAGGTGCAGTTTCTCA +AGACTCTCAGCAGAGGGCGTCGCATCCTGGACAGGAAAATTCAGAGCCTGGGAGACAGCAAGACCATTCC +CGGAGACACTGCTTGGCTCCTCTATGACACCTATGGGTTTCCAGTGGATCTGACTGGACTGATTGCTGAA +GAGAAGGGCCTGGTGGTAGACATGGATGGCTTTGAAGAGGAGAGGAAACTGGCCCAGCTGAAATCACAGG +GCAAGGGAGCTGGTGGGGAAGACCTCATTATGCTGGACATTTACGCTATCGAAGAGCTCCGGGCACGGGG +TCTGGAGGTCACAGATGATTCCCCAAAGTACAATTACCATTTGGACTCCAGTGGTAGCTATGTATTTGAG +AACACAGTGGCTACGGTGATGGCTCTGCGCAGGGAGAAGATGTTCGTGGAAGAGGTGTCCACAGGCCAGG +AGTGTGGAGTGGTGCTGGACAAGACCTGTTTCTATGCTGAGCAAGGAGGCCAGATCTATGACGAAGGCTA +CCTGGTGAAGGTGGATGACAGCAGTGAAGATAAAACAGAGTTTACAGTGAAGAATGCTCAGGTCCGAGGA +GGGTATGTGCTACACATTGGAACCATCTACGGTGACCTGAAAGTGGGGGATCAGGTCTGGCTGTTTATTG +ATGAGCCCCGACGAAGACCCATCATGAGCAACCACACAGCTACGCACATTCTGAACTTCGCCCTGCGCTC +AGTGCTTGGGGAAGCTGACCAGAAAGGCTCATTGGTTGCTCCTGACCGCCTCAGATTTGACTTTACTGCC +AAGGGAGCCATGTCCACCCAACAGATCAAGAAGGCTGAAGAGATTGCTAATGAGATGATTGAGGCAGCCA +AGGCCGTCTATACCCAGGATTGCCCCCTGGCAGCAGCGAAAGCCATCCAGGGCCTACGGGCTGTGTTTGA +TGAGACCTATCCTGACCCTGTGCGAGTCGTCTCCATTGGGGTCCCGGTGTCCGAGTTGCTGGATGACCCC +TCTGGGCCTGCTGGCTCCCTGACTTCTGTTGAGTTCTGTGGGGGAACGCACCTGCGGAACTCGAGTCATG +CAGGAGCTTTTGTGATCGTGACGGAAGAAGCCATTGCCAAGGGTATCCGGAGGATTGTGGCTGTCACAGG +TGCCGAGGCCCAGAAGGCCCTCAGGAAAGCAGAGAGCTTGAAGAAATGTCTCTCTGTCATGGAAGCCAAA +GTGAAGGCTCAGACTGCTCCAAACAAGGATGTGCAGAGGGAGATCGCTGACCTTGGAGAGGCCCTGGCCA +CTGCAGTCATCCCCCAGTGGCAGAAGGATGAATTGCGGGAGACTCTCAAATCCCTAAAGAAGGTCATGGA +TGACTTGGACCGAGCCAGCAAAGCCGATGTCCAGAAACGAGTGTTAGAGAAGACGAAGCAGTTCATCGAC +AGCAACCCCAACCAGCCTCTTGTCATCCTGGAGATGGAGAGCGGCGCCTCAGCCAAGGCCCTGAATGAAG +CCTTGAAGCTCTTCAAGATGCACTCCCCTCAGACTTCTGCCATGCTCTTCACGGTGGACAATGAGGCTGG +CAAGATCACGTGCCTGTGTCAAGTCCCCCAGAATGCAGCCAATCGGGGCTTAAAAGCCAGCGAGTGGGTG +CAGCAGGTGTCAGGCTTGATGGACGGTAAAGGTGGTGGCAAGGATGTGTCTGCACAGGCCACAGGCAAGA +ACGTTGGCTGCCTGCAGGAGGCGCTGCAGCTGGCCACTTCCTTCGCCCAGCTGCGCCTCGGGGATGTAAA +GAACTGAGTGGGGAAGGAGGAGGCTCCCACTGGATCCATCCGTCCAGCCAAGAGCTCTTCATCTGCTACA +AGAACATTTGAATCTTGGGACCTTTAAAGAGCCCCTCCTAACCCAGCAGTAACTGGAACACACTTGGGAG +CAGTCCTATGTCTCAGTGCCCCTTAAATTTCTGCCCTGAGCCCTCCACGTCAGTGCCATCGGTCTAGAAC +CACTAACCCCGCATTGCTGTTGATCGTCACGCTCGCATCTATAGATAACGGCTCTCCAGACCTGAGCTTT +CCGCGTCAGCAAGTAGGAATCGTTTTTGCTGCAGAGAATAAAAGGACCACGTGC + |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/seq1_snp1_mode1.out --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/seq1_snp1_mode1.out Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +SNP w Slen GC interval d_max p-value interval r_min p-value +U1013C 200 3344 0.5411 975-1025 0.2432 0.0724 998-1052 0.0615 0.0932 |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/seq2.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/seq2.txt Tue Dec 13 12:31:56 2016 -0500 |
| b |
| b'@@ -0,0 +1,139 @@\n+>gi|5420376|emb|AJ238799.1| Hepatitis C virus type 1b complete genome, isolate Con1\n+GCCAGCCCCCGATTGGGGGCGACACTCCACCATAGATCACTCCCCTGTGAGGAACTACTGTCTTCACGCA\n+GAAAGCGTCTAGCCATGGCGTTAGTATGAGTGTCGTGCAGCCTCCAGGACCCCCCCTCCCGGGAGAGCCA\n+TAGTGGTCTGCGGAACCGGTGAGTACACCGGAATTGCCAGGACGACCGGGTCCTTTCTTGGATCAACCCG\n+CTCAATGCCTGGAGATTTGGGCGTGCCCCCGCGAGACTGCTAGCCGAGTAGTGTTGGGTCGCGAAAGGCC\n+TTGTGGTACTGCCTGATAGGGTGCTTGCGAGTGCCCCGGGAGGTCTCGTAGACCGTGCACCATGAGCACG\n+AATCCTAAACCTCAAAGAAAAACCAAACGTAACACCAACCGCCGCCCACAGGACGTCAAGTTCCCGGGCG\n+GTGGTCAGATCGTCGGTGGAGTTTACCTGTTGCCGCGCAGGGGCCCCAGGTTGGGTGTGCGCGCGACTAG\n+GAAGACTTCCGAGCGGTCGCAACCTCGTGGAAGGCGACAACCTATCCCCAAGGCTCGCCAGCCCGAGGGT\n+AGGGCCTGGGCTCAGCCCGGGTACCCCTGGCCCCTCTATGGCAATGAGGGCTTGGGGTGGGCAGGATGGC\n+TCCTGTCACCCCGTGGCTCTCGGCCTAGTTGGGGCCCCACGGACCCCCGGCGTAGGTCGCGCAATTTGGG\n+TAAGGTCATCGATACCCTCACGTGCGGCTTCGCCGATCTCATGGGGTACATTCCGCTCGTCGGCGCCCCC\n+CTAGGGGGCGCTGCCAGGGCCCTGGCGCATGGCGTCCGGGTTCTGGAGGACGGCGTGAACTATGCAACAG\n+GGAATCTGCCCGGTTGCTCCTTTTCTATCTTCCTTTTGGCTTTGCTGTCCTGTTTGACCATCCCAGCTTC\n+CGCTTATGAAGTGCGCAACGTATCCGGAGTGTACCATGTCACGAACGACTGCTCCAACGCAAGCATTGTG\n+TATGAGGCAGCGGACATGATCATGCATACCCCCGGGTGCGTGCCCTGCGTTCGGGAGAACAACTCCTCCC\n+GCTGCTGGGTAGCGCTCACTCCCACGCTCGCGGCCAGGAACGCTAGCGTCCCCACTACGACGATACGACG\n+CCATGTCGATTTGCTCGTTGGGGCGGCTGCTCTCTGCTCCGCTATGTACGTGGGAGATCTCTGCGGATCT\n+GTTTTCCTCGTCGCCCAGCTGTTCACCTTCTCGCCTCGCCGGCACGAGACAGTACAGGACTGCAATTGCT\n+CAATATATCCCGGCCACGTGACAGGTCACCGTATGGCTTGGGATATGATGATGAACTGGTCACCTACAGC\n+AGCCCTAGTGGTATCGCAGTTACTCCGGATCCCACAAGCTGTCGTGGATATGGTGGCGGGGGCCCATTGG\n+GGAGTCCTAGCGGGCCTTGCCTACTATTCCATGGTGGGGAACTGGGCTAAGGTTCTGATTGTGATGCTAC\n+TCTTTGCCGGCGTTGACGGGGGAACCTATGTGACAGGGGGGACGATGGCCAAAAACACCCTCGGGATTAC\n+GTCCCTCTTTTCACCCGGGTCATCCCAGAAAATCCAGCTTGTAAACACCAACGGCAGCTGGCACATCAAC\n+AGGACTGCCCTGAACTGCAATGACTCCCTCAACACTGGGTTCCTTGCTGCGCTGTTCTACGTGCACAAGT\n+TCAACTCATCTGGATGCCCAGAGCGCATGGCCAGCTGCAGCCCCATCGACGCGTTCGCTCAGGGGTGGGG\n+GCCCATCACTTACAATGAGTCACACAGCTCGGACCAGAGGCCTTATTGTTGGCACTACGCACCCCGGCCG\n+TGCGGTATCGTACCCGCGGCGCAGGTGTGTGGTCCAGTGTACTGCTTCACCCCAAGCCCTGTCGTGGTGG\n+GGACGACCGACCGGTTCGGCGTCCCTACGTACAGTTGGGGGGAGAATGAGACGGACGTGCTGCTTCTTAA\n+CAACACGCGGCCGCCGCAAGGCAACTGGTTTGGCTGTACATGGATGAATAGCACTGGGTTCACCAAGACG\n+TGCGGGGGCCCCCCGTGTAACATCGGGGGGATCGGCAATAAAACCTTGACCTGCCCCACGGACTGCTTCC\n+GGAAGCACCCCGAGGCCACTTACACCAAGTGTGGTTCGGGGCCTTGGTTGACACCCAGATGCTTGGTCCA\n+CTACCCATACAGGCTTTGGCACTACCCCTGCACTGTCAACTTTACCATCTTCAAGGTTAGGATGTACGTG\n+GGGGGAGTGGAGCACAGGCTCGAAGCCGCATGCAATTGGACTCGAGGAGAGCGTTGTAACCTGGAGGACA\n+GGGACAGATCAGAGCTTAGCCCGCTGCTGCTGTCTACAACGGAGTGGCAGGTATTGCCCTGTTCCTTCAC\n+CACCCTACCGGCTCTGTCCACTGGTTTGATCCATCTCCATCAGAACGTCGTGGACGTACAATACCTGTAC\n+GGTATAGGGTCGGCGGTTGTCTCCTTTGCAATCAAATGGGAGTATGTCCTGTTGCTCTTCCTTCTTCTGG\n+CGGACGCGCGCGTCTGTGCCTGCTTGTGGATGATGCTGCTGATAGCTCAAGCTGAGGCCGCCCTAGAGAA\n+CCTGGTGGTCCTCAACGCGGCATCCGTGGCCGGGGCGCATGGCATTCTCTCCTTCCTCGTGTTCTTCTGT\n+GCTGCCTGGTACATCAAGGGCAGGCTGGTCCCTGGGGCGGCATATGCCCTCTACGGCGTATGGCCGCTAC\n+TCCTGCTCCTGCTGGCGTTACCACCACGAGCATACGCCATGGACCGGGAGATGGCAGCATCGTGCGGAGG\n+CGCGGTTTTCGTAGGTCTGATACTCTTGACCTTGTCACCGCACTATAAGCTGTTCCTCGCTAGGCTCATA\n+TGGTGGTTACAATATTTTATCACCAGGGCCGAGGCACACTTGCAAGTGTGGATCCCCCCCCTCAACGTTC\n+GGGGGGGCCGCGATGCCGTCATCCTCCTCACGTGCGCGATCCACCCAGAGCTAATCTTTACCATCACCAA\n+AATCTTGCTCGCCATACTCGGTCCACTCATGGTGCTCCAGGCTGGTATAACCAAAGTGCCGTACTTCGTG\n+CGCGCACACGGGCTCATTCGTGCATGCATGCTGGTGCGGAAGGTTGCTGGGGGTCATTATGTCCAAATGG\n+CTCTCATGAAGTTGGCCGCACTGACAGGTACGTACGTTTATGACCATCTCACCCCACTGCGGGACTGGGC\n+CCACGCGGGCCTACGAGACCTTGCGGTGGCAGTTGAGCCCGTCGTCTTCTCTGATATGGAGACCAAGGTT\n+ATCACCTGGGGGGCAGACACCGCGGCGTGTGGGGACATCATCTTGGGCCTGCCCGTCTCCGCCCGCAGGG\n+GGAGGGAGATACATCTGGGACCGGCAGACAGCCTTGAAGGGCAGGGGTGGCGACTCCTCGCGCCTATTAC\n+GGCCTACTCCCAACAGACGCGAGGCCTACTTGGCTGCATCATCACTAGCCTCACAGGCCGGGACAGGAAC\n+CAGGTCGAGGGGGAGGTCCAAGTGGTCTCCACCGCAACACAATCTTTCCTGGCGACCTGCGTCAATGGCG\n+TGTGTTGGACTGTCTATCATGGTGCCGGCTCAAAGACCCTTGCCGGCCCAAAGGGCCCAATCACCCAAAT\n+GTACACCAATGTGGACCAGGACCTCGTCGGCTGGCAAGCGCCCCCCGGGGCGCGTTCCTTGACACCATGC\n+ACCTGCGGCAGCTCGGACCTTTACTTGGTCACGAGGCATGCCGATGTCATTCCGGTGCGCCGGCGGGGCG\n+ACAGCAG'..b'CACAGCCTCTATCACCAGCCC\n+GCTCACCACCCAACATACCCTCCTGTTTAACATCCTGGGGGGATGGGTGGCCGCCCAACTTGCTCCTCCC\n+AGCGCTGCTTCTGCTTTCGTAGGCGCCGGCATCGCTGGAGCGGCTGTTGGCAGCATAGGCCTTGGGAAGG\n+TGCTTGTGGATATTTTGGCAGGTTATGGAGCAGGGGTGGCAGGCGCGCTCGTGGCCTTTAAGGTCATGAG\n+CGGCGAGATGCCCTCCACCGAGGACCTGGTTAACCTACTCCCTGCTATCCTCTCCCCTGGCGCCCTAGTC\n+GTCGGGGTCGTGTGCGCAGCGATACTGCGTCGGCACGTGGGCCCAGGGGAGGGGGCTGTGCAGTGGATGA\n+ACCGGCTGATAGCGTTCGCTTCGCGGGGTAACCACGTCTCCCCCACGCACTATGTGCCTGAGAGCGACGC\n+TGCAGCACGTGTCACTCAGATCCTCTCTAGTCTTACCATCACTCAGCTGCTGAAGAGGCTTCACCAGTGG\n+ATCAACGAGGACTGCTCCACGCCATGCTCCGGCTCGTGGCTAAGAGATGTTTGGGATTGGATATGCACGG\n+TGTTGACTGATTTCAAGACCTGGCTCCAGTCCAAGCTCCTGCCGCGATTGCCGGGAGTCCCCTTCTTCTC\n+ATGTCAACGTGGGTACAAGGGAGTCTGGCGGGGCGACGGCATCATGCAAACCACCTGCCCATGTGGAGCA\n+CAGATCACCGGACATGTGAAAAACGGTTCCATGAGGATCGTGGGGCCTAGGACCTGTAGTAACACGTGGC\n+ATGGAACATTCCCCATTAACGCGTACACCACGGGCCCCTGCACGCCCTCCCCGGCGCCAAATTATTCTAG\n+GGCGCTGTGGCGGGTGGCTGCTGAGGAGTACGTGGAGGTTACGCGGGTGGGGGATTTCCACTACGTGACG\n+GGCATGACCACTGACAACGTAAAGTGCCCGTGTCAGGTTCCGGCCCCCGAATTCTTCACAGAAGTGGATG\n+GGGTGCGGTTGCACAGGTACGCTCCAGCGTGCAAACCCCTCCTACGGGAGGAGGTCACATTCCTGGTCGG\n+GCTCAATCAATACCTGGTTGGGTCACAGCTCCCATGCGAGCCCGAACCGGACGTAGCAGTGCTCACTTCC\n+ATGCTCACCGACCCCTCCCACATTACGGCGGAGACGGCTAAGCGTAGGCTGGCCAGGGGATCTCCCCCCT\n+CCTTGGCCAGCTCATCAGCTAGCCAGCTGTCTGCGCCTTCCTTGAAGGCAACATGCACTACCCGTCATGA\n+CTCCCCGGACGCTGACCTCATCGAGGCCAACCTCCTGTGGCGGCAGGAGATGGGCGGGAACATCACCCGC\n+GTGGAGTCAGAAAATAAGGTAGTAATTTTGGACTCTTTCGAGCCGCTCCAAGCGGAGGAGGATGAGAGGG\n+AAGTATCCGTTCCGGCGGAGATCCTGCGGAGGTCCAGGAAATTCCCTCGAGCGATGCCCATATGGGCACG\n+CCCGGATTACAACCCTCCACTGTTAGAGTCCTGGAAGGACCCGGACTACGTCCCTCCAGTGGTACACGGG\n+TGTCCATTGCCGCCTGCCAAGGCCCCTCCGATACCACCTCCACGGAGGAAGAGGACGGTTGTCCTGTCAG\n+AATCTACCGTGTCTTCTGCCTTGGCGGAGCTCGCCACAAAGACCTTCGGCAGCTCCGAATCGTCGGCCGT\n+CGACAGCGGCACGGCAACGGCCTCTCCTGACCAGCCCTCCGACGACGGCGACGCGGGATCCGACGTTGAG\n+TCGTACTCCTCCATGCCCCCCCTTGAGGGGGAGCCGGGGGATCCCGATCTCAGCGACGGGTCTTGGTCTA\n+CCGTAAGCGAGGAGGCTAGTGAGGACGTCGTCTGCTGCTCGATGTCCTACACATGGACAGGCGCCCTGAT\n+CACGCCATGCGCTGCGGAGGAAACCAAGCTGCCCATCAATGCACTGAGCAACTCTTTGCTCCGTCACCAC\n+AACTTGGTCTATGCTACAACATCTCGCAGCGCAAGCCTGCGGCAGAAGAAGGTCACCTTTGACAGACTGC\n+AGGTCCTGGACGACCACTACCGGGACGTGCTCAAGGAGATGAAGGCGAAGGCGTCCACAGTTAAGGCTAA\n+ACTTCTATCCGTGGAGGAAGCCTGTAAGCTGACGCCCCCACATTCGGCCAGATCTAAATTTGGCTATGGG\n+GCAAAGGACGTCCGGAACCTATCCAGCAAGGCCGTTAACCACATCCGCTCCGTGTGGAAGGACTTGCTGG\n+AAGACACTGAGACACCAATTGACACCACCATCATGGCAAAAAATGAGGTTTTCTGCGTCCAACCAGAGAA\n+GGGGGGCCGCAAGCCAGCTCGCCTTATCGTATTCCCAGATTTGGGGGTTCGTGTGTGCGAGAAAATGGCC\n+CTTTACGATGTGGTCTCCACCCTCCCTCAGGCCGTGATGGGCTCTTCATACGGATTCCAATACTCTCCTG\n+GACAGCGGGTCGAGTTCCTGGTGAATGCCTGGAAAGCGAAGAAATGCCCTATGGGCTTCGCATATGACAC\n+CCGCTGTTTTGACTCAACGGTCACTGAGAATGACATCCGTGTTGAGGAGTCAATCTACCAATGTTGTGAC\n+TTGGCCCCCGAAGCCAGACAGGCCATAAGGTCGCTCACAGAGCGGCTTTACATCGGGGGCCCCCTGACTA\n+ATTCTAAAGGGCAGAACTGCGGCTATCGCCGGTGCCGCGCGAGCGGTGTACTGACGACCAGCTGCGGTAA\n+TACCCTCACATGTTACTTGAAGGCCGCTGCGGCCTGTCGAGCTGCGAAGCTCCAGGACTGCACGATGCTC\n+GTATGCGGAGACGACCTTGTCGTTATCTGTGAAAGCGCGGGGACCCAAGAGGACGAGGCGAGCCTACGGG\n+CCTTCACGGAGGCTATGACTAGATACTCTGCCCCCCCTGGGGACCCGCCCAAACCAGAATACGACTTGGA\n+GTTGATAACATCATGCTCCTCCAATGTGTCAGTCGCGCACGATGCATCTGGCAAAAGGGTGTACTATCTC\n+ACCCGTGACCCCACCACCCCCCTTGCGCGGGCTGCGTGGGAGACAGCTAGACACACTCCAGTCAATTCCT\n+GGCTAGGCAACATCATCATGTATGCGCCCACCTTGTGGGCAAGGATGATCCTGATGACTCATTTCTTCTC\n+CATCCTTCTAGCTCAGGAACAACTTGAAAAAGCCCTAGATTGTCAGATCTACGGGGCCTGTTACTCCATT\n+GAGCCACTTGACCTACCTCAGATCATTCAACGACTCCATGGCCTTAGCGCATTTTCACTCCATAGTTACT\n+CTCCAGGTGAGATCAATAGGGTGGCTTCATGCCTCAGGAAACTTGGGGTACCGCCCTTGCGAGTCTGGAG\n+ACATCGGGCCAGAAGTGTCCGCGCTAGGCTACTGTCCCAGGGGGGGAGGGCTGCCACTTGTGGCAAGTAC\n+CTCTTCAACTGGGCAGTAAGGACCAAGCTCAAACTCACTCCAATCCCGGCTGCGTCCCAGTTGGATTTAT\n+CCAGCTGGTTCGTTGCTGGTTACAGCGGGGGAGACATATATCACAGCCTGTCTCGTGCCCGACCCCGCTG\n+GTTCATGTGGTGCCTACTCCTACTTTCTGTAGGGGTAGGCATCTATCTACTCCCCAACCGATGAACGGGG\n+AGCTAAACACTCCAGGCCAATAGGCCATCCTGTTTTTTTCCCTTTTTTTTTTTCTTTTTTTTTTTTTTTT\n+TTTTTTTTTTTTTTTTTTCTCCTTTTTTTTTCCTCTTTTTTTCCTTTTCTTTCCTTTGGTGGCTCCATCT\n+TAGCCCTAGTCACGGCTAGCTGTGAAAGGTCCGTGAGCCGCTTGACTGCAGAGAGTGCTGATACTGGCCT\n+CTCTGCAGATCAAGT\n' |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/seq2_snp2_mode2.out --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/seq2_snp2_mode2.out Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +SNP w Slen GC max_k d_max p-value interval d p-value +C9294A-U9296G 200 9605 0.4814 9270 7.0487 0.0624 9270-9298 0.2463 0.0099 |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/seq3.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/seq3.txt Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +>gi|56682960|ref|NM_000146.3| Homo sapiens ferritin, light polypeptide (FTL), mRNA +GCAGTTCGGCGGTCCCGCGGGTCTGTCTCTTGCTTCAACAGTGTTTGGACGGAACAGATCCGGGGACTCTCTTCCAGCCTCCGACCGCCCTCCGATTTCCTCTCCGCTTGCAACCTCCGGGACCATCTTCTCGGCCATCTCCTGCTTCTGGGACCTGCCAGCACCGTTTTTGTGGTTAGCTCCTTCTTGCCAACCAACCATGAGCTCCCAGATTCGTCAGAATTATTCCACCGACGTGGAGGCAGCCGTCAACAGCCTGGTCAATTTGTACCTGCAGGCCTCCTACACCTACCTCTCTCTGGGCTTCTATTTCGACCGCGATGATGTGGCTCTGGAAGGCGTGAGCCACTTCTTCCGCGAATTGGCCGAGGAGAAGCGCGAGGGCTACGAGCGTCTCCTGAAGATGCAAAACCAGCGTGGCGGCCGCGCTCTCTTCCAGGACATCAAGAAGCCAGCTGAAGATGAGTGGGGTAAAACCCCAGACGCCATGAAAGCTGCCATGGCCCTGGAGAAAAAGCTGAACCAGGCCCTTTTGGATCTTCATGCCCTGGGTTCTGCCCGCACGGACCCCCATCTCTGTGACTTCCTGGAGACTCACTTCCTAGATGAGGAAGTGAAGCTTATCAAGAAGATGGGTGACCACCTGACCAACCTCCACAGGCTGGGTGGCCCGGAGGCTGGGCTGGGCGAGTATCTCTTCGAAAGGCTCACTCTCAAGCACGACTAAGAGCCTTCTGAGCCCAGCGACTTCTGAAGGGCCCCTTGCAAAGTAATAGGGCTTCTGCCTAAGCCTCTCCCTCCAGCCAATAGGCAGCTTTCTTAACTATCCTAACAAGCCTTGGACCAAATGGAAATAAAGCTTTTTGATGCA |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/seq3_mode3.out --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/seq3_mode3.out Tue Dec 13 12:31:56 2016 -0500 |
| b |
| b'@@ -0,0 +1,158 @@\n+SNP\tw\tSlen\tGC\tinterval\td\tpvalue1\tewin\tinterval\td_max\tpvalue2\n+C18G\t200\t871\t0.5917\t7-134\t0.1938\t0.0521\t200\t8-136\t0.2894\t0.0300\n+G19A\t200\t871\t0.5890\t10-67\t0.0779\t0.2841\t200\t147-216\t0.2377\t0.0488\n+G19U\t200\t871\t0.5890\t9-67\t0.1096\t0.1738\t200\t14-78\t0.2050\t0.0686\n+U22G\t200\t871\t0.5856\t11-66\t0.2165\t0.0348\t200\t15-64\t0.2482\t0.0518\n+U22A\t200\t871\t0.5856\t17-62\t0.1034\t0.1982\t200\t15-64\t0.2379\t0.0572\n+U24G\t200\t871\t0.5804\t13-66\t0.1306\t0.1342\t200\t15-64\t0.2300\t0.0612\n+U24C\t200\t871\t0.5804\t17-62\t0.2214\t0.0354\t200\t15-64\t0.3698\t0.0165\n+C29G\t200\t871\t0.5764\t17-54\t0.1108\t0.1880\t200\t15-64\t0.1828\t0.0998\n+G47C\t200\t871\t0.5870\t27-53\t0.1553\t0.0924\t200\t13-62\t0.2068\t0.0970\n+A58G\t200\t871\t0.5853\t7-124\t0.2689\t0.0171\t200\t11-134\t0.3584\t0.0286\n+C61G\t200\t871\t0.5862\t7-124\t0.2035\t0.0469\t200\t12-67\t0.2907\t0.0436\n+G119C\t200\t871\t0.5768\t102-134\t0.1579\t0.0634\t200\t106-162\t0.2463\t0.0656\n+C123G\t200\t871\t0.5759\t104-152\t0.0976\t0.2067\t200\t106-162\t0.2816\t0.0527\n+C132G\t200\t871\t0.5753\t101-134\t0.1489\t0.0857\t200\t86-135\t0.2334\t0.0898\n+U173C\t200\t871\t0.5791\t161-205\t0.2468\t0.0151\t200\t157-206\t0.2965\t0.0409\n+G175A\t200\t871\t0.5760\t168-206\t0.1995\t0.0337\t200\t157-206\t0.2359\t0.0722\n+G175C\t200\t871\t0.5760\t159-205\t0.2822\t0.0084\t200\t157-206\t0.2988\t0.0400\n+G175U\t200\t871\t0.5760\t160-205\t0.2448\t0.0158\t200\t157-206\t0.2590\t0.0580\n+U176G\t200\t871\t0.5771\t168-210\t0.2892\t0.0083\t200\t148-211\t0.3160\t0.0330\n+U176A\t200\t871\t0.5771\t161-206\t0.1782\t0.0523\t200\t148-245\t0.2834\t0.0449\n+U177G\t200\t871\t0.5782\t168-210\t0.2112\t0.0275\t200\t148-211\t0.2638\t0.0535\n+U177A\t200\t871\t0.5782\t159-198\t0.2630\t0.0112\t200\t157-206\t0.2751\t0.0480\n+C180A\t200\t871\t0.5816\t161-198\t0.3182\t0.0049\t200\t157-206\t0.2923\t0.0420\n+C180G\t200\t871\t0.5816\t142-196\t0.1680\t0.0591\t200\t150-209\t0.2780\t0.0480\n+C183A\t200\t871\t0.5822\t168-206\t0.2411\t0.0148\t200\t157-206\t0.2633\t0.0557\n+C183G\t200\t871\t0.5822\t168-206\t0.2123\t0.0245\t200\t157-206\t0.2418\t0.0682\n+U185G\t200\t871\t0.5844\t164-196\t0.2181\t0.0224\t200\t157-206\t0.2211\t0.0739\n+C186G\t200\t871\t0.5829\t168-206\t0.2129\t0.0243\t200\t157-206\t0.2409\t0.0598\n+C190A\t200\t871\t0.5821\t168-206\t0.2089\t0.0342\t200\t157-206\t0.2323\t0.0781\n+C190G\t200\t871\t0.5821\t169-206\t0.2357\t0.0223\t200\t157-206\t0.2500\t0.0664\n+A193C\t200\t871\t0.5852\t161-198\t0.1819\t0.0387\t200\t157-206\t0.2047\t0.0970\n+C194A\t200\t871\t0.5838\t161-206\t0.1696\t0.0525\t200\t157-206\t0.2068\t0.0923\n+C194G\t200\t871\t0.5838\t169-205\t0.2459\t0.0140\t200\t157-206\t0.2578\t0.0570\n+C195U\t200\t871\t0.5848\t162-203\t0.2574\t0.0103\t200\t157-206\t0.2341\t0.0684\n+C207A\t200\t871\t0.5835\t168-210\t0.2023\t0.0355\t200\t148-211\t0.2693\t0.0573\n+C208G\t200\t871\t0.5810\t146-222\t0.1694\t0.0608\t200\t146-222\t0.2252\t0.0852\n+G211C\t200\t871\t0.5736\t168-212\t0.2152\t0.0287\t200\t148-211\t0.2651\t0.0594\n+A212G\t200\t871\t0.5736\t145-222\t0.2256\t0.0242\t200\t146-222\t0.2624\t0.0609\n+A230C\t200\t871\t0.5786\t225-245\t0.1066\t0.1712\t200\t227-304\t0.2132\t0.0953\n+G239A\t200\t871\t0.5835\t161-240\t0.1107\t0.1598\t200\t142-245\t0.2323\t0.0799\n+G239C\t200\t871\t0.5835\t154-245\t0.1454\t0.0903\t200\t142-245\t0.3070\t0.0411\n+A240U\t200\t871\t0.5835\t154-245\t0.1253\t0.1258\t200\t142-245\t0.2684\t0.0577\n+A244G\t200\t871\t0.5810\t131-248\t0.0613\t0.3614\t200\t242-441\t0.2532\t0.0662\n+U273G\t200\t871\t0.5786\t237-283\t0.1633\t0.0673\t200\t227-287\t0.2089\t0.0991\n+G277C\t200\t871\t0.5761\t238-284\t0.0631\t0.3506\t200\t317-366\t0.3635\t0.0253\n+A315G\t200\t871\t0.5686\t298-436\t0.0863\t0.2390\t200\t294-382\t0.2233\t0.0868\n+G320C\t200\t871\t0.5661\t315-358\t0.0625\t0.3540\t200\t317-366\t0.2257\t0.0849\n+U327G\t200\t871\t0.5661\t311-359\t0.1266\t0.1231\t200\t301-366\t0.3298\t0.0337\n+U327C\t200\t871\t0.5661\t315-357\t0.2849\t0.0092\t200\t316-366\t0.3349\t0.0323\n+U327A\t200\t871\t0.5661\t311-358\t0.1683\t0.0620\t200\t301-366\t0.3630\t0.0254\n+G328A\t200\t871\t0.5661\t312-359\t0.1532\t0.0795\t200\t317-366\t0.3304\t0.0335\n+G328C\t200\t871\t0.5661\t317-357\t0.4111\t0.0012\t200\t312-446\t0.3102\t0.0400\n+G328U\t200\t871\t0.5661\t317-357\t0.1720\t0.0583\t200\t312-446\t0.2868\t0.0490\n+G329A\t200\t871\t0.5686\t317-356\t0.0696\t0.3149\t200\t312-446\t0.3310\t0.0334\n+G329C\t200\t871\t0.5686\t317-357\t0.1854\t0.0468\t200\t312-361\t0.2975\t0.0447\n+G329U\t200\t871\t0.5686\t322-352\t0.1447\t0.0915\t200\t312-444\t0.3268\t0.0346\n+C330A\t200\t871\t0.5711\t315-358\t0.3260\t0.0047\t200\t317-3'..b'9\t0.0408\n+G503U\t200\t871\t0.5761\t439-504\t0.1370\t0.1037\t200\t375-556\t0.3375\t0.0316\n+C504G\t200\t871\t0.5736\t484-504\t0.1543\t0.0781\t200\t485-536\t0.2287\t0.0826\n+C505U\t200\t871\t0.5736\t462-531\t0.1406\t0.0977\t200\t467-531\t0.2824\t0.0510\n+C506A\t200\t871\t0.5761\t498-545\t0.1543\t0.0781\t200\t500-549\t0.2323\t0.0799\n+U507G\t200\t871\t0.5786\t501-536\t0.1399\t0.0988\t200\t487-536\t0.2751\t0.0544\n+U507A\t200\t871\t0.5786\t498-536\t0.1579\t0.0736\t200\t487-536\t0.2953\t0.0455\n+A510G\t200\t871\t0.5786\t461-531\t0.1746\t0.0559\t200\t467-531\t0.2998\t0.0437\n+G526C\t200\t871\t0.5786\t504-527\t0.1315\t0.1135\t200\t371-543\t0.2739\t0.0550\n+C528G\t200\t871\t0.5786\t484-504\t0.1174\t0.1433\t200\t448-560\t0.2623\t0.0610\n+U538G\t200\t871\t0.5761\t479-550\t0.1485\t0.0858\t200\t448-560\t0.3685\t0.0242\n+C542G\t200\t871\t0.5761\t437-556\t0.1471\t0.0878\t200\t436-556\t0.2684\t0.0577\n+C546G\t200\t871\t0.5786\t442-560\t0.0765\t0.2811\t200\t447-559\t0.2643\t0.0599\n+C546U\t200\t871\t0.5786\t436-555\t0.0950\t0.2071\t200\t437-556\t0.2194\t0.0900\n+A567U\t200\t871\t0.5786\t564-641\t0.1412\t0.0968\t200\t562-641\t0.2266\t0.0842\n+G611C\t200\t871\t0.5711\t587-637\t0.4460\t0.0007\t200\t588-637\t0.4910\t0.0088\n+G611U\t200\t871\t0.5711\t582-631\t0.2850\t0.0091\t200\t588-637\t0.3935\t0.0196\n+U615G\t200\t871\t0.5736\t580-623\t0.3044\t0.0066\t200\t575-631\t0.4179\t0.0160\n+G616C\t200\t871\t0.5711\t580-617\t0.3617\t0.0026\t200\t580-631\t0.3970\t0.0191\n+A623G\t200\t871\t0.5586\t564-636\t0.1137\t0.1523\t200\t562-636\t0.2217\t0.0881\n+A627C\t200\t871\t0.5511\t564-636\t0.1108\t0.1595\t200\t562-636\t0.2152\t0.0935\n+A629G\t200\t871\t0.5511\t565-636\t0.1348\t0.1076\t200\t562-636\t0.2352\t0.0778\n+U637G\t200\t871\t0.5486\t561-641\t0.1374\t0.1030\t200\t562-641\t0.2370\t0.0766\n+A639U\t200\t871\t0.5486\t561-641\t0.1267\t0.1229\t200\t562-641\t0.2487\t0.0689\n+G673C\t200\t871\t0.5338\t632-674\t0.0744\t0.2848\t200\t632-681\t0.2025\t0.0896\n+C678U\t200\t871\t0.5381\t670-745\t0.2001\t0.0314\t200\t632-681\t0.2926\t0.0412\n+G682U\t200\t871\t0.5359\t670-781\t0.1327\t0.1047\t200\t634-683\t0.3026\t0.0371\n+C700A\t200\t871\t0.5323\t687-734\t0.3024\t0.0069\t200\t688-737\t0.2130\t0.0874\n+C700G\t200\t871\t0.5323\t687-734\t0.2348\t0.0200\t200\t688-737\t0.2303\t0.0739\n+A702G\t200\t871\t0.5351\t670-745\t0.2305\t0.0204\t200\t677-745\t0.2007\t0.0987\n+G705A\t200\t871\t0.5313\t688-734\t0.2973\t0.0101\t200\t688-737\t0.3114\t0.0402\n+G705C\t200\t871\t0.5313\t688-734\t0.3418\t0.0052\t200\t688-737\t0.3277\t0.0348\n+G705U\t200\t871\t0.5313\t688-734\t0.3782\t0.0030\t200\t688-737\t0.3589\t0.0265\n+G706C\t200\t871\t0.5301\t688-734\t0.2760\t0.0092\t200\t688-737\t0.2796\t0.0480\n+G706U\t200\t871\t0.5301\t688-734\t0.3926\t0.0013\t200\t688-737\t0.3664\t0.0220\n+C709G\t200\t871\t0.5289\t688-734\t0.2638\t0.0130\t200\t688-737\t0.2477\t0.0659\n+U725G\t200\t871\t0.5331\t688-734\t0.3898\t0.0018\t200\t688-737\t0.3540\t0.0271\n+A729U\t200\t871\t0.5306\t700-734\t0.2789\t0.0089\t200\t681-741\t0.2549\t0.0618\n+G730U\t200\t871\t0.5292\t688-734\t0.3585\t0.0039\t200\t688-737\t0.3104\t0.0448\n+C731A\t200\t871\t0.5279\t688-734\t0.4156\t0.0021\t200\t687-736\t0.3854\t0.0223\n+C731G\t200\t871\t0.5279\t688-734\t0.3453\t0.0057\t200\t688-737\t0.3380\t0.0330\n+C732A\t200\t871\t0.5294\t688-734\t0.3023\t0.0076\t200\t686-736\t0.2990\t0.0407\n+U734C\t200\t871\t0.5325\t688-734\t0.2446\t0.0187\t200\t686-735\t0.2282\t0.0829\n+A767G\t200\t871\t0.5213\t746-771\t0.0963\t0.2207\t200\t737-786\t0.2207\t0.0937\n+A775U\t200\t871\t0.5185\t770-843\t0.0638\t0.3782\t200\t737-786\t0.2391\t0.0796\n+U787G\t200\t871\t0.5193\t770-814\t0.2678\t0.0135\t200\t757-814\t0.2145\t0.0946\n+G811C\t200\t871\t0.5211\t773-814\t0.1850\t0.0578\t200\t673-819\t0.2061\t0.0954\n+C812G\t200\t871\t0.5192\t777-814\t0.3720\t0.0024\t200\t765-814\t0.2643\t0.0577\n+C838U\t200\t871\t0.5342\t824-862\t0.0937\t0.2344\t200\t821-870\t0.2084\t0.0786\n+U855G\t200\t871\t0.5300\t812-861\t0.1019\t0.1781\t200\t813-862\t0.1722\t0.0971\n+A858G\t200\t871\t0.5280\t824-862\t0.1187\t0.1149\t200\t813-862\t0.2035\t0.0669\n+C860G\t200\t871\t0.5283\t826-863\t0.1000\t0.1610\t200\t814-863\t0.1686\t0.0927\n+U862G\t200\t871\t0.5238\t800-870\t0.0646\t0.2581\t200\t764-871\t0.1600\t0.0978\n+U862A\t200\t871\t0.5238\t830-866\t0.1325\t0.0769\t200\t770-866\t0.1756\t0.0819\n+U863G\t200\t871\t0.5215\t824-866\t0.1291\t0.0891\t200\t817-866\t0.2549\t0.0360\n+A867G\t200\t871\t0.5171\t830-871\t0.1078\t0.0925\t200\t817-871\t0.1844\t0.0678\n+G869C\t200\t871\t0.5172\t830-871\t0.0884\t0.0736\t200\t814-871\t0.1601\t0.0522\n+C870G\t200\t871\t0.5149\t823-870\t0.0480\t0.0976\t200\t817-870\t0.1080\t0.0490\n' |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/snp1.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/snp1.txt Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,1 @@ +U1013C |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/snp2.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/snp2.txt Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,1 @@ +C9294A-U9296G |
| b |
| diff -r 000000000000 -r 327f5721612f test-data/snp3.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/snp3.txt Tue Dec 13 12:31:56 2016 -0500 |
| b |
| @@ -0,0 +1,2 @@ +U22G +A56U |