Drosophila species x D. melanogaster DNA-align stats (dmel-dxxx-dna alignments blastn) Don Gilbert, gilbertd@indiana.edu, Nov 2005 Dmel r4.1 source genome size 2L 22407834 2R 20766785 3L 23771897 3R 27905053 4 1281640 X 22224390 Euch. Total 118357599 # Heterochromatin (not in dmel-chr matches) 2h 1694122 3h 2955737 4h 88110 Xh 359526 Yh 396896 U 8724946 Hetch. Total 14219337 Euch. + Hetr 132576936 Nov. 2005 assemblies Total bases Matched to Pct of Pct of in assembly Dmel DNA Dmel asmbly dmel 118357599 118357599 100.00% 100.00% dsim 118374359 117123159 98.96% 98.94% dsec 166577145 113537874 95.93% 68.16% dyak 168647858 116445353 98.38% 69.05% updated to dyak2 dere 152862534 117010609 98.86% 76.55% dana 231676663 112468944 95.02% 48.55% dper 188374079 109919775 92.87% 58.35% dpse 152738921 110531782 93.39% 72.37% updated to dpse2 dmoj 194270144 103661429 87.58% 53.36% dvir 206998770 92070099 77.79% 44.48% dgri 201006931 93235198 78.77% 46.38% 2004 assemblies (earlier than aug05): Total bases Matched to Pct of Pct of in assembly Dmel DNA Dmel asmbly dere 133575709 102698324 76.88% dana 232189998 48290795 20.80% dmoj 205886621 60834298 29.55% dvir 196204567 53966234 27.51% dgri 226113271 103271421 45.67% dyak 155068951 115882003 97.91% 74.73% << update dyak2 dpse 140943297 109631216 92.63% 77.78% << update dpse2 dyak* wustl asm has some yeast genome contam. in chrU (~25MB); used dyak asm from NCBI w/o yeast dsim assembly here is dgpg (syntenic) lacking heterochromatin and U contigs. Dsim parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dsim-dna1.gff.gz to gff: stdout filters: skipsubseg=3771, joinoverlap=7210, minsizecut=9 Chr cover gaps n_inverts n_segs n_targets prob 2L 22063191 320007 2 90 1 1.00e-10 2R 20520198 239414 0 101 1 1.98e-20 3L 23638389 113900 2 83 1 2.41e-11 3R 27723964 171021 0 94 1 0.00e+00 4 1267656 13487 0 11 1 0.00e+00 X 21909761 313944 0 116 1 6.90e-147 Total 117123159 1171773 4 495 6 2.22e-11 =========================================================================== Dsec parsing sorted dmel-dsec-dna.gff0.gz to gff: stdout filters: skipsubseg=25870, joinoverlap=1977, minsizecut=46 Chr cover gaps n_inverts n_segs n_targets prob 2L 21411657 986873 2 361 182 8.31e-10 2R 19996854 764010 4 398 233 1.11e-11 3L 22896244 870087 5 355 188 5.64e-10 3R 26686482 1218719 8 392 168 2.55e-10 4 1252966 28517 1 83 74 4.82e-08 X 21293671 929186 13 452 260 2.88e-08 Total 113537874 4797392 33 2041 1105 8.63e-09 =========================================================================== chr x targets over 1% of 113537874 cover, from n=704 targets dsec_br0510:scaffold_0 20924877 18.4% 3R 53.0% 3L 47.0% dsec_br0510:scaffold_1 13632915 12.0% 2R 100.0% dsec_br0510:scaffold_2 7365845 6.5% 3L 100.0% dsec_br0510:scaffold_3 7091294 6.2% 2L 100.0% dsec_br0510:scaffold_4 6053796 5.3% X 54.1% 3R 45.8% dsec_br0510:scaffold_5 5387123 4.7% 2L 85.9% 3R 9.5% dsec_br0510:scaffold_6 3967509 3.5% 3R 99.6% dsec_br0510:scaffold_8 3565902 3.1% X 99.0% dsec_br0510:scaffold_7 3507381 3.1% 2L 100.0% dsec_br0510:scaffold_9 3095495 2.7% 2R 99.6% dsec_br0510:scaffold_10 3042853 2.7% X 100.0% dsec_br0510:scaffold_11 2681144 2.4% 3L 99.3% dsec_br0510:scaffold_14 2055043 1.8% 2L 99.7% dsec_br0510:scaffold_12 1963707 1.7% 3R 100.0% dsec_br0510:scaffold_13 1936360 1.7% 3R 100.0% dsec_br0510:scaffold_16 1881261 1.7% 2L 99.2% dsec_br0510:scaffold_15 1876006 1.7% X 100.0% dsec_br0510:scaffold_17 1410172 1.2% X 100.0% dsec_br0510:scaffold_19 1278470 1.1% 3L 100.0% --------------------------------------------------------------------------- Dyak dyak-wu20051101 parsing sorted dmel-dyak2-dna.gff to gff: dmel-dyak2-algn.gff filters: skipsubseg=5916, joinoverlap=2240, minsizecut=11 Chr cover gaps n_inverts n_segs n_targets prob 2L 21783058 613383 10 178 41 5.06e-09 2R 20339876 382679 10 215 95 2.09e-09 3L 23501917 228126 7 381 250 2.89e-11 3R 27686077 212610 16 353 159 8.50e-12 4 1239824 33633 6 36 22 2.78e-10 X 21894601 314440 18 254 86 1.57e-06 Total 116445353 1784871 67 1417 653 2.83e-07 =========================================================================== dyak2 nt=168647858 chr x targets over 1% of 116445353 cover, from n= 17 targets dyak_wu051101:chr3R 26006308 22.3% 3R 99.7% dyak_wu051101:chrX 21349236 18.3% X 99.8% dyak_wu051101:chr2L 20923434 18.0% 2L 77.4% 2R 22.2% dyak_wu051101:chr3L 20566292 17.7% 3L 99.9% dyak_wu051101:chr2R 19898201 17.1% 2R 73.0% 2L 26.8% dyak_wu051101:chr3L_random 2200735 1.9% 3L 97.3% 2L 1.2% dyak_wu051101:chr3R_random 1600312 1.4% 3R 100.0% dyak_wu051101:chrU 1560291 1.3% 2R 39.6% 3L 24.2% --------------------------------------------------------------------------- > dyak asm1 parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dyak-dna1.gff.gz to gff: stdout filters: skipsubseg=6884, joinoverlap=5139, minsizecut=5 Chr cover gaps n_inverts n_segs n_targets prob 2L 21796998 590837 15 172 15 4.07e-15 2R 20108232 610845 12 175 47 1.16e-24 3L 23399322 349815 18 301 132 6.64e-15 3R 27612045 284503 10 243 55 1.65e-28 4 1230442 24579 3 26 20 3.85e-25 X 21734964 472724 18 259 52 3.09e-42 Total 115882003 2333303 76 1176 321 2.30e-15 =========================================================================== Dere parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dere-dna1.gff.gz to gff: stdout filters: skipsubseg=17840, joinoverlap=6369, minsizecut=36 Chr cover gaps n_inverts n_segs n_targets prob 2L 22101188 294787 6 214 73 2.84e-06 2R 20430892 323776 5 209 109 2.40e-06 3L 23563410 186648 13 252 107 3.21e-08 3R 27703242 195279 6 179 42 3.35e-11 4 1251538 27669 3 43 37 9.30e-10 X 21960339 259516 13 310 121 3.25e-07 Total 117010609 1287675 46 1207 489 1.01e-06 =========================================================================== chr x targets over 1% of 117010609 cover, from n=164 targets dere_ag01aug05:scaffold_4784 23165649 19.8% 3L 99.0% dere_ag01aug05:scaffold_4929 21889578 18.7% 2L 74.5% 2R 23.5% dere_ag01aug05:scaffold_4845 20523675 17.5% 2R 72.2% 2L 27.1% dere_ag01aug05:scaffold_4690 18943751 16.2% X 99.9% dere_ag01aug05:scaffold_4770 17508004 15.0% 3R 99.6% dere_ag01aug05:scaffold_4820 10187806 8.7% 3R 100.0% dere_ag01aug05:scaffold_4644 2458832 2.1% X 100.0% --------------------------------------------------------------------------- >> old dere.asm1 parsing sorted /bio/archive-w/eugenes/genomes/dspp.old/dmel-dere-dna.gff.gz to gff: stdout filters: skipsubseg=8046, joinoverlap=107, minsizecut=38 Chr cover gaps n_inverts n_segs n_targets prob 2L 19147819 3245664 0 1345 1338 2.24e-12 2R 17087870 3673897 2 1107 1091 2.71e-10 3L 20610139 3157937 1 1486 1473 5.94e-09 3R 24188468 3711645 2 1742 1723 2.33e-10 4 1206084 73186 0 106 101 2.83e-10 X 20457944 1764258 2 2657 2645 1.51e-08 Total 102698324 15626587 7 8443 8371 5.87e-09 =========================================================================== Dana parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dana-dna1.gff.gz to gff: stdout filters: skipsubseg=30042, joinoverlap=4068, minsizecut=511 Chr cover gaps n_inverts n_segs n_targets prob 2L 21267252 1107150 33 499 118 6.84e-07 2R 19748965 1012039 69 464 106 5.44e-07 3L 22734696 1009917 58 540 95 1.82e-06 3R 27128698 769436 63 586 119 1.84e-07 4 910565 363711 5 78 60 1.05e-06 X 20678768 1532201 60 625 245 1.12e-07 Total 112468944 5794454 288 2792 743 6.59e-07 =========================================================================== chr x targets over 1% of 112468944 cover, from n=178 targets dana_ag01aug05:scaffold_13340 23191576 20.6% 3R 100.0% dana_ag01aug05:scaffold_13337 22154571 19.7% 3L 99.9% dana_ag01aug05:scaffold_13266 18927703 16.8% 2R 99.9% dana_ag01aug05:scaffold_12916 15928326 14.2% 2L 100.0% dana_ag01aug05:scaffold_13117 5211900 4.6% X 100.0% dana_ag01aug05:scaffold_12943 4338716 3.9% 2L 99.1% dana_ag01aug05:scaffold_13248 3257709 2.9% X 100.0% dana_ag01aug05:scaffold_12911 3110910 2.8% 3R 97.6% 3L 1.4% dana_ag01aug05:scaffold_12929 2777774 2.5% X 100.0% dana_ag01aug05:scaffold_13335 2214728 2.0% X 96.4% 4 2.7% dana_ag01aug05:scaffold_13417 2156934 1.9% X 95.4% 2L 2.1% dana_ag01aug05:scaffold_13047 1864879 1.7% X 99.7% dana_ag01aug05:scaffold_13334 1297381 1.2% X 100.0% --------------------------------------------------------------------------- >> old dana1.asm; much poorer parsing sorted /bio/archive-w/eugenes/genomes/dspp.old/dmel-dana-dna.gff.gz to gff: stdout filters: skipsubseg=10566, joinoverlap=24, minsizecut=458 Chr cover gaps n_inverts n_segs n_targets prob 2L 7324773 14959976 3 520 490 6.21e-07 2R 7866492 12890620 2 566 539 1.07e-06 3L 8338995 15395630 0 546 526 1.90e-06 3R 8199364 19638924 2 616 582 5.05e-08 4 725827 548371 0 105 105 3.90e-08 X 15835344 6375293 5 1502 1426 7.95e-08 Total 48290795 69808814 12 3855 3668 5.49e-07 =========================================================================== Dper parsing sorted dmel-dper-dna.gff0.gz to gff: stdout filters: skipsubseg=27428, joinoverlap=3811, minsizecut=757 Chr cover gaps n_inverts n_segs n_targets prob 2L 20997853 1382141 61 537 140 7.29e-08 2R 19030590 1710777 41 606 248 3.40e-07 3L 21783429 1932488 27 777 425 9.73e-07 3R 26855995 1038819 73 676 141 1.46e-08 4 1102648 171702 14 42 19 1.00e-12 X 20149260 2061602 60 748 374 5.76e-08 Total 109919775 8297529 276 3386 1347 3.11e-07 =========================================================================== chr x targets over 1% of 109919775 cover, from n=542 targets dper_br051028:scaffold_0 10546038 9.6% 3R 100.0% dper_br051028:scaffold_2 8088520 7.4% 2R 100.0% dper_br051028:scaffold_1 7909431 7.2% 2L 99.5% dper_br051028:scaffold_4 6875218 6.3% 2R 100.0% dper_br051028:scaffold_3 5682112 5.2% 3R 98.1% 3L 1.9% dper_br051028:scaffold_6 5007713 4.6% 3R 100.0% dper_br051028:scaffold_5 4982522 4.5% 2L 100.0% dper_br051028:scaffold_7 3977208 3.6% 3R 99.6% dper_br051028:scaffold_9 3205734 2.9% 3L 100.0% dper_br051028:scaffold_8 3157799 2.9% 2L 100.0% dper_br051028:scaffold_10 2440257 2.2% 2L 99.9% dper_br051028:scaffold_12 2248344 2.0% 3L 76.4% X 23.6% dper_br051028:scaffold_13 1981579 1.8% X 100.0% dper_br051028:scaffold_11 1969019 1.8% X 99.5% dper_br051028:scaffold_15 1941422 1.8% X 100.0% dper_br051028:scaffold_14 1778054 1.6% X 100.0% dper_br051028:scaffold_19 1610086 1.5% 3R 100.0% dper_br051028:scaffold_20 1568239 1.4% 3L 100.0% dper_br051028:scaffold_17 1523148 1.4% X 100.0% dper_br051028:scaffold_18 1501884 1.4% X 100.0% dper_br051028:scaffold_23 1362058 1.2% 3L 100.0% dper_br051028:scaffold_21 1349260 1.2% X 100.0% dper_br051028:scaffold_24 1347883 1.2% 3L 100.0% dper_br051028:scaffold_22 1293524 1.2% X 100.0% dper_br051028:scaffold_16 1226341 1.1% 2L 97.3% 2R 2.7% --------------------------------------------------------------------------- Dpse > dpse2, 20051111 filters: skipsubseg=2641, joinoverlap=2098, minsizecut=123 Chr cover gaps n_inverts n_segs n_targets prob 2L 20891248 1488594 48 385 138 2.79e-08 2R 19142216 1600552 104 384 43 5.31e-10 3L 22230819 1481879 55 403 124 1.89e-09 3R 26811728 1081031 125 500 63 2.01e-08 4 995900 269991 8 31 11 9.69e-14 X 20459871 1745595 78 497 205 1.83e-09 Total 110531782 7667642 418 2200 584 1.03e-08 ===================================================================== nt=152738921 chr x targets over 1% of 110531782 cover, from n= 22 targets dpse_r201:2 26440716 23.9% 3R 100.0% dpse_r201:3 18573949 16.8% 2R 100.0% dpse_r201:XR_group6 12023065 10.9% 3L 100.0% dpse_r201:XL_group1e 10360267 9.4% X 100.0% dpse_r201:4_group3 9558119 8.6% 2L 100.0% dpse_r201:XR_group8 8019288 7.3% 3L 100.0% dpse_r201:XL_group1a 6912084 6.3% X 99.9% dpse_r201:4_group4 4875786 4.4% 2L 99.4% dpse_r201:4_group1 3748617 3.4% 2L 98.8% 3L 1.2% dpse_r201:XL_group3a 2163093 2.0% X 99.3% dpse_r201:U 1675439 1.5% X 30.8% 2R 24.8% dpse_r201:4_group5 1647222 1.5% 2L 100.0% dpse_r201:XR_group3a 1232444 1.1% 3L 100.0% --------------------------------------------------------------------------- > dpse1 parsing sorted dmel-dpse-dna.gff to gff: stdout filters: skipsubseg=10169, joinoverlap=3879, minsizecut=348 Chr cover gaps n_inverts n_segs n_targets prob 2L 20867315 1512654 51 512 132 1.08e-09 2R 18778093 1962648 116 579 51 4.51e-09 3L 21886963 1793940 54 571 147 1.14e-09 3R 26712923 1181906 134 691 39 5.73e-11 4 1008597 250602 17 44 13 2.27e-09 X 20377325 1833433 91 644 212 2.56e-07 Total 109631216 8535183 463 3041 594 5.54e-08 ========================================================================= Dmoj parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dmoj-dna1.gff.gz to gff: stdout filters: skipsubseg=21049, joinoverlap=4352, minsizecut=911 Chr cover gaps n_inverts n_segs n_targets prob 2L 19882638 2494887 131 555 48 9.33e-07 2R 18085072 2654012 130 553 81 1.11e-06 3L 20941345 2818043 123 595 110 2.75e-07 3R 25347482 2550751 195 728 48 1.62e-07 4 923751 343292 10 43 18 2.90e-09 X 18481141 3720473 104 764 339 2.15e-07 Total 103661429 14581458 693 3238 644 4.87e-07 =========================================================================== chr x targets over 1% of 103661429 cover, from n=136 targets dmoj_ag01aug05:scaffold_6540 25305337 24.4% 3R 99.6% dmoj_ag01aug05:scaffold_6500 19885765 19.2% 2L 99.2% dmoj_ag01aug05:scaffold_6680 18865493 18.2% 3L 99.7% dmoj_ag01aug05:scaffold_6496 17635854 17.0% 2R 99.8% dmoj_ag01aug05:scaffold_6473 10075495 9.7% X 99.6% dmoj_ag01aug05:scaffold_6359 2716293 2.6% X 99.9% dmoj_ag01aug05:scaffold_6328 2648994 2.6% X 100.0% dmoj_ag01aug05:scaffold_6654 1889782 1.8% 3L 99.7% dmoj_ag01aug05:scaffold_6308 1875446 1.8% X 100.0% --------------------------------------------------------------------------- > old dmoj.asm1 parsing sorted /bio/archive-w/eugenes/genomes/dspp.old/dmel-dmoj-dna.gff.gz to gff: stdout filters: skipsubseg=4822, joinoverlap=13, minsizecut=332 Chr cover gaps n_inverts n_segs n_targets prob 2L 11954684 10400188 8 1020 950 1.59e-07 2R 10368447 10239281 5 867 812 2.75e-07 3L 12008324 11681075 5 1187 1122 3.22e-07 3R 13454524 14380421 4 1365 1301 9.12e-07 4 497007 770029 2 36 33 8.33e-06 X 12551312 9642395 3 1506 1478 3.12e-07 Total 60834298 57113389 27 5981 5696 4.68e-07 =========================================================================== Dvir parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dvir-dna.gff0.gz to gff: stdout filters: skipsubseg=13957, joinoverlap=1989, minsizecut=673 Chr cover gaps n_inverts n_segs n_targets prob 2L 17363636 5000714 124 701 150 7.18e-07 2R 15914333 4826536 110 710 151 3.38e-07 3L 19392241 4285270 136 719 87 3.73e-07 3R 21975107 5923295 111 867 211 1.62e-07 4 798862 443037 8 44 28 1.93e-07 X 16625920 5576718 139 881 295 5.36e-08 Total 92070099 26055570 628 3922 922 3.08e-07 =========================================================================== chr x targets over 1% of 92070099 cover, from n=109 targets dvir_ag0508:scaffold_13049 18050132 19.6% 3L 99.5% dvir_ag0508:scaffold_12875 12852583 14.0% 2R 99.7% dvir_ag0508:scaffold_13047 12761962 13.9% 3R 99.7% dvir_ag0508:scaffold_12963 12345017 13.4% 2L 99.5% dvir_ag0508:scaffold_12970 6666235 7.2% X 99.1% dvir_ag0508:scaffold_12855 6522326 7.1% 3R 99.5% dvir_ag0508:scaffold_12928 4002228 4.3% X 100.0% dvir_ag0508:scaffold_12723 3265917 3.5% 2L 99.2% dvir_ag0508:scaffold_12822 2590380 2.8% 3R 100.0% dvir_ag0508:scaffold_13042 2422238 2.6% X 100.0% dvir_ag0508:scaffold_12726 1724746 1.9% X 98.5% 3L 1.5% dvir_ag0508:scaffold_13246 1521823 1.7% 2L 100.0% dvir_ag0508:scaffold_12823 1457836 1.6% 2R 100.0% dvir_ag0508:scaffold_12932 1013201 1.1% X 100.0% --------------------------------------------------------------------------- >> old dvir.asm1 parsing sorted /bio/archive-w/eugenes/genomes/dspp.old/dmel-dvir-dna.gff.gz to gff: stdout filters: skipsubseg=6673, joinoverlap=7, minsizecut=470 Chr cover gaps n_inverts n_segs n_targets prob 2L 10229655 12134669 4 674 651 1.12e-06 2R 8897441 11829383 4 592 577 4.62e-07 3L 10448303 13222211 2 659 630 2.21e-07 3R 11225844 16613219 4 656 622 2.01e-07 4 592787 673061 0 39 38 3.63e-08 X 12572204 9511821 5 1420 1385 7.94e-07 Total 53966234 63984364 19 4040 3903 6.03e-07 =========================================================================== Dgri parsing sorted /bio/bio-grid/shared/chrs//dgri2/gff/dmel-dgri-dna.gff0.gz to gff: stdout filters: skipsubseg=8561, joinoverlap=1944, minsizecut=958 Chr cover gaps n_inverts n_segs n_targets prob 2L 18011830 4356081 113 799 250 8.43e-07 2R 16321684 4415162 117 708 157 2.03e-06 3L 18625534 5053333 151 813 160 1.42e-07 3R 22861891 5033216 81 945 295 2.48e-07 4 750886 496991 14 39 8 1.12e-09 X 16663373 5538327 136 850 287 1.66e-07 Total 93235198 24893110 612 4154 1157 6.27e-07 =========================================================================== chr x targets over 1% of 93235198 cover, from n=238 targets dgri_ag0509:scaffold_15110 17988208 19.3% 3L 99.9% dgri_ag0509:scaffold_15245 12893379 13.8% 2R 99.9% dgri_ag0509:scaffold_15252 11288339 12.1% 2L 99.7% dgri_ag0509:scaffold_14906 9785373 10.5% 3R 99.4% dgri_ag0509:scaffold_15203 7595965 8.1% X 99.5% dgri_ag0509:scaffold_14853 6315918 6.8% X 100.0% dgri_ag0509:scaffold_15126 5671369 6.1% 2L 99.7% dgri_ag0509:scaffold_15074 5171336 5.5% 3R 99.9% dgri_ag0509:scaffold_14830 4580434 4.9% 3R 99.6% dgri_ag0509:scaffold_15112 3178817 3.4% 2R 100.0% dgri_ag0509:scaffold_15081 2762231 3.0% X 99.7% dgri_ag0509:scaffold_14624 2679058 2.9% 3R 99.4% --------------------------------------------------------------------------- >> old dgri.asm1 ; better dna cover ?? parsing sorted /bio/argos/species/web/data/dmel-dspp/gff1.old/dmel-dgri-dna1.gff.gz to gff: stdout filters: skipsubseg=17793, joinoverlap=3570, minsizecut=1201 Chr cover gaps n_inverts n_segs n_targets prob 2L 19712730 2655062 57 678 350 3.08e-07 2R 17881426 2855305 61 593 252 3.79e-07 3L 20890604 2788173 72 723 352 2.81e-08 3R 25367855 2530471 87 782 297 1.70e-07 4 751637 496248 11 47 20 1.09e-05 X 18667169 3534449 111 768 351 1.72e-07 Total 103271421 14859708 399 3591 1622 3.42e-07 =========================================================================== $bg/blast/chrgff2align.pl -in=dmel-dsec-dna.gff0.gz -nogff -stats > & dmel-dsec-algn.stats $bg/blast/chrgff2align.pl -in=dmel-dper-dna.gff0.gz -nogff -stats > & dmel-dper-algn.stats $bg/blast/chrgff2align.pl -in=$sc/dgri2/gff/dmel-dgri-dna.gff0.gz -nogff -stats > & dmel-dgri-algn.stats $bg/blast/chrgff2align.pl -in=$sp/web/data/dmel-dspp/gff1.old/dmel-dana-dna1.gff.gz -nogff -stats > & dmel-dana-algn.stats $bg/blast/chrgff2align.pl -in=$sp/web/data/dmel-dspp/gff1.old/dmel-dere-dna1.gff.gz -nogff -stats > & dmel-dere-algn.stats $bg/blast/chrgff2align.pl -in=$sp/web/data/dmel-dspp/gff1.old/dmel-dmoj-dna1.gff.gz -nogff -stats > & dmel-dmoj-algn.stats $bg/blast/chrgff2align.pl -in=$sp/web/data/dmel-dspp/gff1.old/dmel-dvir-dna.gff0.gz -nogff -stats > & dmel-dvir-algn.stats $bg/blast/chrgff2align.pl -in=$sp/web/data/dmel-dspp/gff1.old/dmel-dsim-dna1.gff.gz -nogff -stats > & dmel-dsim1-algn.stats $bg/blast/chrgff2align.pl -in=$sp/web/data/dmel-dspp/gff1.old/dmel-dyak-dna1.gff.gz -nogff -stats > & dmel-dyak1-algn.stats /bio/archive-w/eugenes/genomes/dspp.old/ -- asm1 dmel-dxxx-dna data $bg/blast/chrgff2align.pl -in=/bio/archive-w/eugenes/genomes/dspp.old/dmel-dana-dna.gff.gz \ -nogff -stats > & dmel-dana1-algn.stats