https://stackoverflow.com/questions/27551921/how-to-extend-ambiguous-dna-sequence from Bio import Seq from itertools import product def extend_ambiguous_dna(seq): """return list of all possible sequences given an ambiguous DNA input""" d = Seq.IUPAC.IUPACData.ambiguous_dna_values r = [] for i in product(*[d[j] for j in seq]): r.append("".join(i)) return r In [1]: extend_ambiguous_dna("AV") Out[1]: ['AA', 'AC', 'AG'] # 2^18= 262144 combinations extend_ambiguous_dna("TCTRGGGTCGTATGATCTTACCTATTCTCAAACTATAAATGGGTAAGGTAATGGGCAGCATGCTCGGATGATGCTGCCCTCAAAGCGGTTGAAAAGCAAACAGTGCCCGTGAGGGTGCTAGCTAGATATCGGTGTGCTTAGTGGGCCAAGTTTTGGTAAGCAGAACTGGTGCTGTGGGATGAACCAAACGTGATGTTACGGCGCCTAAATAAACGACGCATCATAGATACCATGAAAGGTGTTGATTGCTAAAGACAGCAGGACGGTGGACATGGAAGTTGTCATCCGCTAAGGAGTGTGTAACAACTCACCTGCCGAAGCAATTAGCCCTTAAAATGGATGGCGCTCAAGTCGTTTGCCTATACATCGCCGCTAGCGGCAGAATCTTGGTGGGAAGCCGGCGCGCTGTCCCACTTTGAGGCCCTAGTGAGTAGGAGGGTACGGTGGTGGCGTTGAAGTGTTTGGCGCAAGCCGGCATGGAGCCGCCACTGGCACAGATCTTGGTGGTAGTAGCAAATATTCGAATGAGATCTTGGATGACTGAAGTGGAGGAGGGTTTCGTGTCAACAGCAGTTGAACCGAAARGTTAGCCASTTCCGGTTTTCCAARATTSCWTTCAAATWTTTTCTGTTMCSCCMGGGTGTGGGCCCGTTGGGGGGGCGYTGGGCGSGGGGGCGGGCSTGTGGCGGGGGCGGYGYGTCTCTGTCTGWGGGGYGGCSGGGGGGGGGGGGGGGGGGGGGGGGGSG") ~180Mb Nucleic acid code Meaning Mnemonic A A Adenine C C Cytosine G G Guanine T T Thymine U U Uracil R A or G Purine Y C, T or U Pyrimidines K G, T or U Bases which are ketones M A or C Bases with amino groups S C or G Strong interaction W A, T or U Weak interaction B Not A (i.e. C, G, T or U) B comes after A D Not C (i.e. A, G, T or U) D comes after C H Not G (i.e., A, C, T or U) H comes after G V Neither T nor U (i.e. A, C or G) V comes after U N A C G T U Nucleic acid X Masked - Gap of indeterminate length #gsnap gmap_build --help -D, --dir=STRING Destination directory for installation (defaults to gmapdb directory specified at configure time) -d, --db=STRING Genome name -k, --kmer=INT k-mer value for genomic index (allowed: 15 or less, default is 15) -q INT sampling interval for genomoe (allowed: 1-3, default 3) reczko@max:/data/images/proton2/run375/YK$ gmap_build -D . -d IK28species -q 1 IK28species.fa ../bin/MosaikBuild -fr /data/results/tools/align/mosaik/MOSAIK/test/pipiwhite.fa -oa ik28/pipiwhite.dat &> log1 ../bin/MosaikBuild -q /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1IonXpress_017.fastq -out t.mkb -st 454 reczko@max:/data/results/tools/align/mosaik/MOSAIK/test$ ../bin/MosaikAligner -mm 30 -hs 11 -in t.mkb -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out t.mka -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log ../bin/MosaikBuild -q /data/images/proton2/run374/R_2017_12_07_15_50_07_user_IONAS-374-DrKobold_PHlab_YK_KELPNO_171207.KLPNC4-trf5504p.IonXpress_017.fastq -out t.mkb -st 454 ../bin/MosaikAligner -mm 30 -hs 11 -in t.mkb -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out t.mka -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log #0% ../bin/MosaikAligner -mm 30 -hs 11 -in t.mkb -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out t.mka -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log #0% ../bin/MosaikBuild -q /data/images/proton2/run374/R_2017_12_07_15_50_07_user_IONAS-374-DrKobold_PHlab_YK_KELPNO_171207.YKR1-RNAMosq.IonXpress_016.fastq -out t.mkb -st 454 #0% ../bin/MosaikBuild -q /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1IonXpress_017.fastq -out t.mkb -st 454 #0% ../bin/MosaikBuild -q /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1.YKR1-RNAMosq.IonXpress_016.fastq -out t.mkb -st 454 #0% for i in /data/images/proton2/run374/*fastq do echo $i /data/results/tools/align/bowtie2-2.2.9/bowtie2 -x IK28species-unambig $i --no-unal --local --very-sensitive-local -S 375_017.sam done for i in /data/images/proton2/run375/*fastq do echo $i /data/results/tools/align/bowtie2-2.2.9/bowtie2 -x IK28species-unambig $i --no-unal --local --very-sensitive-local -S 375_017.sam done /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1.PH3R42-MCas12-Con2.IonXpress_032.fastq 7487246 reads; of these: 7487246 (100.00%) were unpaired; of these: 7486535 (99.99%) aligned 0 times 550 (0.01%) aligned exactly 1 time 161 (0.00%) aligned >1 times 0.01% overall alignment rate /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1.MT3R22-iP13.IonXpress_023.fastq 4778294 reads; of these: 4778294 (100.00%) were unpaired; of these: 4776932 (99.97%) aligned 0 times 194 (0.00%) aligned exactly 1 time 1168 (0.02%) aligned >1 times 0.03% overall alignment rate Dear Pantelis+Vagelis searching for the sample of Yannis Karakasiliotis, I've aligned all samples in runs 374,5,6,7,8 to the 13 seqs YK sent with very sensitive local alignment. There is no sample seems to contain his material: Most unique alignments: /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1.PH3R42-MCas12-Con2.IonXpress_032.fastq 7487246 reads; of these: 7487246 (100.00%) were unpaired; of these: 7486535 (99.99%) aligned 0 times 550 (0.01%) aligned exactly 1 time 161 (0.00%) aligned >1 times 0.01% overall alignment rate Most multi alignments: /data/images/proton2/run374/R_2017_12_07_15_50_07_user_IONAS-374-DrKobold_PHlab_YK_KELPNO_171207.KLPNC4-trf5504p.IonXpress_017.fastq 28782052 reads; of these: 28782052 (100.00%) were unpaired; of these: 28780352 (99.99%) aligned 0 times 53 (0.00%) aligned exactly 1 time 1647 (0.01%) aligned >1 times 0.01% overall alignment rate Pls check. BW Martin All alignments: #@ nobarcode: for i in nom* > do > /data/results/tools/formats/bam2fastq-1.1.0/bam2fastq -o $i.fastq -q $i > done run374/nomatch_rawlib.basecaller.bam 6227610 reads; of these: 6227610 (100.00%) were unpaired; of these: 6096984 (97.90%) aligned 0 times 1658 (0.03%) aligned exactly 1 time 128968 (2.07%) aligned >1 times 2.10% overall alignment rate run375/nomatch_rawlib.basecaller.bam.fastq 10263218 reads; of these: 10263218 (100.00%) were unpaired; of these: 10098014 (98.39%) aligned 0 times 3791 (0.04%) aligned exactly 1 time 161413 (1.57%) aligned >1 times 1.61% overall alignment rate Dear all, analysis of the nobarcodes in run374,375 and 338 as neg. control suggests that ~250k reads of 374,375 map to YK's reference (vs 180 reads in the neg control). The mapping quality distribution indicates an enrichment of high-qual alingments in 374,375 (vs no enrichm. in control). The detected species distribution is very similar in 374 and 375 (and both are dissimilar to the control. Thus about 250k sequences of YK have no barcode in 374,375. BW Martin Details: #mapping ../bin/MosaikBuild -q /data/images/proton2/run374/nomatch_rawlib.basecaller.bam.fastq -out t.mkb -st 454 ../bin/MosaikAligner -p 5 -mmp 0.5 -act 15 -hs 10 -in t.mkb -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out r374c -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log2c # unaligned mates(X): 6084617 ( 97.7 %) # filtered out(F): 27159 ( 0.4 %) # uniquely aligned mates(U): 8518 ( 0.1 %) # multiply aligned mates(M): 107316 ( 1.7 %) ================================================ total aligned:  115834 ( 1.9 %) total: 6227610 ../bin/MosaikBuild -q /data/images/proton2/run375/nomatch_rawlib.basecaller.bam.fastq -out r375 -st 454 ../bin/MosaikAligner -p 5 -mmp 0.5 -act 15 -hs 10 -in r375 -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out r375c -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log1c # unaligned mates(X): 10083565 ( 98.2 %) # filtered out(F): 32940 ( 0.3 %) # uniquely aligned mates(U): 11228 ( 0.1 %) # multiply aligned mates(M): 135485 ( 1.3 %) ================================================ total aligned:  146713 ( 1.4 %) total: 10263218 ##neg control /data/images/proton2/run338/nomatch_rawlib.basecaller.bam ../bin/MosaikBuild -q /tmp/r338nobarcode.fq -out r338 -st 454 ../bin/MosaikAligner -p 5 -mmp 0.5 -act 15 -hs 10 -in r338 -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out r338 -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log338 # unaligned mates(X): 5153227 ( 99.9 %) # filtered out(F): 3859 ( 0.1 %) # uniquely aligned mates(U): 73 ( 0.0 %) # multiply aligned mates(M): 109 ( 0.0 %) ================================================ total aligned:  182 ( 0.0 %) total: 5157268 #map qual samtools view r374c.bam -F 4 | awk -F "\t" '{o[$5]++;n++;}END{for (i=0;i<256;i++){x=o[i];if(x>0){print i,x,100.0*x/n}}}' 0 1456 1.25697 1 1295 1.11798 2 9244 7.98039 3 1980 1.70934 4 1853 1.5997 5 1441 1.24402 6 797 0.688054 7 2368 2.0443 8 1945 1.67913 9 328 0.283164 10 419 0.361725 11 500 0.431652 12 418 0.360861 13 700 0.604313 14 792 0.683737 15 1301 1.12316 16 1474 1.27251 17 1848 1.59539 18 1702 1.46934 19 2812 2.42761 20 2263 1.95366 21 2738 2.36373 22 1692 1.46071 23 1100 0.949635 24 1011 0.872801 25 1403 1.21122 26 1735 1.49783 27 5425 4.68343 28 10006 8.63822 29 4833 4.17235 30 5727 4.94414 31 9085 7.84312 32 10199 8.80484 33 17811 15.3763 34 5959 5.14443 35 173 0.149352 36 1 0.000863304 samtools view r375c.bam -F 4 | awk -F "\t" '{o[$5]++;n++;}END{for (i=0;i<256;i++){x=o[i];if(x>0){print i,x,100.0*x/n}}}' 0 1886 1.2855 1 2069 1.41024 2 12347 8.41575 3 2459 1.67606 4 2377 1.62017 5 1744 1.18872 6 1103 0.751808 7 2851 1.94325 8 2123 1.44704 9 479 0.326488 10 563 0.383742 11 623 0.424639 12 676 0.460764 13 1020 0.695235 14 1234 0.841098 15 1942 1.32367 16 2335 1.59154 17 2634 1.79534 18 2194 1.49544 19 2881 1.9637 20 2294 1.5636 21 2930 1.9971 22 1929 1.31481 23 1361 0.927661 24 1421 0.968558 25 1832 1.2487 26 2283 1.5561 27 7282 4.96343 28 13084 8.91809 29 6754 4.60355 30 7706 5.25243 31 10705 7.29656 32 12937 8.8179 33 20925 14.2625 34 7487 5.10316 35 243 0.165629 ##neg control samtools view r338.bam -F 4 | awk -F "\t" '{o[$5]++;n++;}END{for (i=0;i<256;i++){x=o[i];if(x>0){print i,x,100.0*x/n}}}' 0 15 8.24176 1 16 8.79121 2 23 12.6374 3 9 4.94505 4 8 4.3956 5 2 1.0989 6 4 2.1978 7 5 2.74725 8 1 0.549451 9 9 4.94505 10 4 2.1978 13 2 1.0989 14 3 1.64835 15 1 0.549451 16 2 1.0989 17 4 2.1978 18 6 3.2967 19 9 4.94505 20 8 4.3956 21 3 1.64835 22 6 3.2967 25 1 0.549451 26 1 0.549451 27 3 1.64835 28 3 1.64835 29 1 0.549451 30 3 1.64835 31 8 4.3956 32 1 0.549451 33 13 7.14286 34 8 4.3956 #species dist reczko@max:/data/results/tools/align/mosaik/MOSAIK/test$ samtools view r374c.bam -F 4 | awk -f /data/results/tools/align/get-sam-refid-stats1.awk detri 3655 3.15538 richi 8404 7.25521 caspius 3241 2.79797 marti 4979 4.29839 pipiwhite 3268 2.82128 longi 13789 11.9041 F10_PSEUDO.ab1 12574 10.8552 pipiblack 572 0.49381 vexans 1860 1.60575 albo 4084 3.52574 B10_MAC.ab1 39587 34.1756 pulcri 5023 4.33638 annu 14798 12.7752 samtools view r375c.bam -F 4 | awk -f /data/results/tools/align/get-sam-refid-stats1.awk richi 10827 7.37971 detri 5165 3.52048 caspius 3946 2.6896 marti 5919 4.03441 pipiwhite 4590 3.12856 longi 18170 12.3847 F10_PSEUDO.ab1 14933 10.1784 pipiblack 780 0.53165 vexans 2664 1.81579 albo 5227 3.56274 pulcri 6171 4.20617 B10_MAC.ab1 50382 34.3405 annu 17939 12.2273 ##neg control samtools view r338.bam -F 4 | awk -f /data/results/tools/align/get-sam-refid-stats1.awk richi 3 1.64835 detri 4 2.1978 caspius 1 0.549451 marti 1 0.549451 pipiwhite 10 5.49451 longi 5 2.74725 F10_PSEUDO.ab1 49 26.9231 vexans 1 0.549451 pipiblack 2 1.0989 pulcri 2 1.0989 albo 5 2.74725 B10_MAC.ab1 92 50.5495 annu 7 3.84615 cd /data/results/tools/align/mosaik/MOSAIK/test reczko@max:/data/results/tools/align/mosaik/MOSAIK/test$ ../bin/MosaikBuild -q /data/images/proton2/run375/R_2017_12_18_17_11_17_user_IONAS-375-DrTokameh_PHlab_YK_171218_MT3R19-23_PH3R41-46_YKR1.YKR1-RNAMosq.IonXpressRNA_016.fastq -out r375b -st 454 ../bin/MosaikAligner -p 5 -mmp 0.5 -act 15 -hs 10 -in r375b -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out r375b -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log375b # too short: 10449 ( 0.3 %) # failed hash: 3740425 ( 95.3 %) ------------------------------------------------ # unaligned mates(X): 3750874 ( 95.6 %) # filtered out(F): 20980 ( 0.5 %) # uniquely aligned mates(U): 5403 ( 0.1 %) # multiply aligned mates(M): 147046 ( 3.7 %) ================================================ total aligned:  152449 ( 3.9 %) total: 3924303 ../bin/MosaikBuild -q /data/images/proton2/run375/nomatch_rawlib.basecaller.bam.fastq -out r375 -st 454 ../bin/MosaikAligner -p 5 -mmp 0.5 -act 15 -hs 10 -in r375 -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out r375c -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log1c # unaligned mates(X): 10083565 ( 98.2 %) # filtered out(F): 32940 ( 0.3 %) # uniquely aligned mates(U): 11228 ( 0.1 %) # multiply aligned mates(M): 135485 ( 1.3 %) ================================================ total aligned:  146713 ( 1.4 %) total: 10263218 ../bin/MosaikBuild -q /data/images/proton2/run374/nomatch_rawlib.basecaller.bam.fastq -out r374 -st 454 ../bin/MosaikAligner -p 5 -mmp 0.5 -act 15 -hs 10 -in r374 -ia /data/results/tools/align/mosaik/MOSAIK/test/ik28/pipiwhite.dat -out r374c -annse /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.se.100.005.ann -annpe /data/results/tools/align/mosaik/MOSAIK/src/networkFile/2.1.26.pe.100.0065.ann &> log374c # unaligned mates(X): 6084617 ( 97.7 %) # filtered out(F): 27159 ( 0.4 %) # uniquely aligned mates(U): 8518 ( 0.1 %) # multiply aligned mates(M): 107316 ( 1.7 %) ================================================ total aligned:  115834 ( 1.9 %) total: 6227610 samtools view -F 4 r374c.bam | awk -v th=20 -f /data/results/tools/ngs/sam/sam2aligned-fasta1.awk > r374.fa samtools view -F 4 r375c.bam | awk -v th=20 -f /data/results/tools/ngs/sam/sam2aligned-fasta1.awk > r375.fa #@ vs Ref /data/results/tools/align/bowtie2-2.2.9/bowtie2 -x /data/results/reference/insect/aedes_albopictus/MNAF02 /data/images/proton2/run375/nomatch_rawlib.basecaller.bam.fastq --no-unal --local --very-sensitive-local -S 375.sam /data/results/tools/align/bowtie2-2.2.9/bowtie2 -x /data/results/reference/insect/aedes_albopictus/MNAF02 /data/images/proton2/run374/nomatch_rawlib.basecaller.bam.fastq --no-unal --local --very-sensitive-local -S 374.sam #neg /data/results/tools/align/bowtie2-2.2.9/bowtie2 -x /data/results/reference/insect/aedes_albopictus/MNAF02 /tmp/r338nobarcode.fq --no-unal --local --very-sensitive-local -S 338.sam run374 nomatch 6227610 reads; of these: 6227610 (100.00%) were unpaired; of these: 3572288 (57.36%) aligned 0 times 75752 (1.22%) aligned exactly 1 time 2579570 (41.42%) aligned >1 times 42.64% overall alignment rate run375 nomatch 10263218 reads; of these: 10263218 (100.00%) were unpaired; of these: 6283187 (61.22%) aligned 0 times 252984 (2.46%) aligned exactly 1 time 3727047 (36.31%) aligned >1 times 38.78% overall alignment rate 5157268 reads; of these: 5157268 (100.00%) were unpaired; of these: 4958699 (96.15%) aligned 0 times 78941 (1.53%) aligned exactly 1 time 119628 (2.32%) aligned >1 times 3.85% overall alignment rate #@ samtools view r374c.bam -F 4 > /tmp/foo samtools view r375c.bam -F 4 >> /tmp/foo cat /tmp/foo | awk -f sam2aligned-fasta1.awk > r375-all.fa cat /tmp/foo | awk -v th=10 -f sam2aligned-fasta1.awk > r375-MQ10.fa #Passed 210402 mappings of 262547 80.1388 % cat /tmp/foo | awk -v th=20 -f sam2aligned-fasta1.awk > r375-MQ20.fa #Passed 182334 mappings of 262547 69.4481 % cat /tmp/foo | awk -v th=30 -f sam2aligned-fasta1.awk > r375-MQ30.fa #Passed 108958 mappings of 262547 41.5004 % for i in /data/results/tools/align/mosaik/MOSAIK/test/*r375-*fa do ln -s $i done ]0;/data/results/tools/align/mosaik/MOSAIK/testreczko@max:/data/results/tools/align/mosaik/MOSAIK/test$ cat /tmp/foo | awk -v th=30 -f /data/results/tools/align/get-sam-refid-stats1.awk |sort -rn -k3,3 MQ cutoff30 Passed 108958 mappings of 262547 41.5004 % longi 26041 23.9 B10_MAC.ab1 25187 23.1162 annu 25051 22.9914 F10_PSEUDO.ab1 8821 8.09578 richi 7512 6.8944 pulcri 5182 4.75596 marti 4332 3.97584 caspius 2920 2.67993 albo 2644 2.42662 vexans 791 0.725968 pipiblack 238 0.218433 detri 200 0.183557 pipiwhite 39 0.0357936 ]0;/data/results/tools/align/mosaik/MOSAIK/testreczko@max:/data/results/tools/align/mosaik/MOSAIK/test$ cat /tmp/foo | awk -v th=20 -f /data/results/tools/align/get-sam-refid-stats1.awk |sort -rn -k3,3 MQ cutoff20 Passed 182334 mappings of 262547 69.4481 % B10_MAC.ab1 59639 32.7087 longi 30368 16.6551 annu 28120 15.4222 F10_PSEUDO.ab1 16572 9.08882 richi 15766 8.64677 pulcri 9450 5.1828 marti 8082 4.43252 caspius 5327 2.92156 albo 5094 2.79377 vexans 2913 1.59762 pipiblack 520 0.285191 detri 383 0.210054 pipiwhite 100 0.0548444 ]0;/data/results/tools/align/mosaik/MOSAIK/testreczko@max:/data/results/tools/align/mosaik/MOSAIK/test$ cat /tmp/foo | awk -v th=10 -f /data/results/tools/align/get-sam-refid-stats1.awk |sort -rn -k3,3 MQ cutoff10 Passed 210402 mappings of 262547 80.1388 % B10_MAC.ab1 68556 32.5833 annu 31713 15.0726 longi 31366 14.9077 F10_PSEUDO.ab1 21087 10.0222 richi 17530 8.33167 pulcri 11053 5.25328 marti 10168 4.83265 caspius 6156 2.92583 albo 6081 2.89018 vexans 4273 2.03087 detri 1227 0.583169 pipiblack 1005 0.477657 pipiwhite 187 0.0888775