Skip to content

Instantly share code, notes, and snippets.

@hyphaltip
Created October 12, 2016 04:59
Show Gist options
  • Save hyphaltip/2d244b9f77a2f63683ea9a128ea4b71c to your computer and use it in GitHub Desktop.
Save hyphaltip/2d244b9f77a2f63683ea9a128ea4b71c to your computer and use it in GitHub Desktop.
# tfastx36 -E 1e-8 sexM.fa ../../ZyGoLife/datasets/genomes/DNA/Phycomyces_blakesleeanus_NRRL_1555.Phybl2.v2.fasta
TFASTX compares a protein to a translated DNA data bank
version 36.3.7a Jan, 2015(preload9)
Please cite:
Pearson et al, Genomics (1997) 46:24-36
Query: sexM.fa
1>>>tr|B0F2G8|B0F2G8_PHYBL Sex minus OS=Phycomyces blakesleeanus GN=sexM PE=4 SV=1 - 188 aa
Library: ../../ZyGoLife/datasets/genomes/DNA/Phycomyces_blakesleeanus_NRRL_1555.Phybl2.v2.fasta
53939167 residues in 80 sequences
Statistics: Altschul/Gish params: n0: 188 Lambda: 0.188 K: 0.034 H: 0.170
statistics sampled from 29 (34) to 29 sequences
Algorithm: TFASTX (3.8 June 2014) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -12/-2, shift: -20
ktup: 2, E-join: 0.5 (0.165), E-opt: 0.1 (0.0411), width: 16
Scan time: 3.580
The best scores are: initn opt bits E(414)
Phybl2|scaffold_1 (4452461) [f] 227 253 74.3 5e-13
>>Phybl2|scaffold_1 (4452461 aa)
Frame: f initn: 227 init1: 227 opt: 253 Z-score: 313.4 bits: 74.3 E(414): 5e-13
Smith-Waterman score: 258; 31.1% identity (57.5% similar) in 212 aa overlap (1-185:3985605-3986237)
10 20 30 40 50
tr|B0F MSLFIFVDESITNNRSSSS--NEEISSHQAD-----TTKEKTPRPKNAFMLYRQEKQKGL
: :.::.: .. . :. .: ::: . .. : ::::::::::: . ..
Phybl2 MESFMFVNEFKNSYQPPSDILTESISSTPVRQKCIISSTVKIRRPKNAFMLYRQAVHPSI
3985620 3985650 3985680 3985710 3985740 3985770
60 70 80 90 100
tr|B0F PLSKTKLLSKDFSKIAAEMWRRESNEVRMYYHRLAEEEKLRHLAKHPGYKYCPK----RQ
:.. . .:..:. :..::. :..::: ::.: :.:::: : : ::: : :. :.
Phybl2 LSSNSTIHNKEISRTAGKMWKNEKEEVRKYYERKADEEKLYHSKKFPGYIYKPQQRKTRR
3985800 3985830 3985860 3985890 3985920 3985950
110 120 130 140 150
tr|B0F SLAHVDKNEYQVLADTTPIVGEMDN-----------QMAFVH---RNVSFDEVSPSSPPS
. : : . .: : . .. : ..: .... ::. ..::.
Phybl2 PQSTVCKPFLRSTSDIQLIYQKTESKSPSEKLSATSQDYIIHCRDESTNEDEIL-NTPPT
3985980 3986010 3986040 3986070 3986100 3986130
160 170 180
tr|B0F SMVQD--LQIFPEGTNLDALSLVDYEFYLLLS
. ..: :. : . ..: .: .: .. : :
Phybl2 QCLDDSILKTFLDDLSFDIFSPIDELIFSLCS
3986160 3986190 3986220
188 residues in 1 query sequences
53939167 residues in 80 library sequences
Tcomplib [36.3.7a Jan, 2015(preload9)] (32 proc in memory [0G])
start: Tue Oct 11 21:58:12 2016 done: Tue Oct 11 21:58:13 2016
Total Scan time: 3.580 Total Display time: 0.030
Function used was TFASTX [36.3.7a Jan, 2015(preload9)]
fasta36 -E 1e-8 sexM.fa ZyGoLife/datasets/genomes/pep/Phycomyces_blakesleeanus_NRRL_1555.Phybl2.v2.aa.fasta | more
# fasta36 -E 1e-8 sexM.fa ZyGoLife/datasets/genomes/pep/Phycomyces_blakesleeanus_NRRL_1555.Phybl2.v2.aa.fasta
FASTA searches a protein or DNA sequence data bank
version 36.3.7a Jan, 2015(preload9)
Please cite:
W.R. Pearson & D.J. Lipman PNAS (1988) 85:2444-2448
Query: sexM.fa
1>>>tr|B0F2G8|B0F2G8_PHYBL Sex minus OS=Phycomyces blakesleeanus GN=sexM PE=4 SV=1 - 188 aa
Library: ../../ZyGoLife/datasets/genomes/pep/Phycomyces_blakesleeanus_NRRL_1555.Phybl2.v2.aa.fasta
5831043 residues in 16528 sequences
Statistics: Expectation_n fit: rho(ln(x))= 6.5022+/-0.00239; mu= 10.1257+/- 0.129
mean_var=59.3897+/-12.015, 0's: 0 Z-trim(93.6): 31 B-trim: 0 in 0/38
Lambda= 0.166425
statistics sampled from 2807 (2816) to 2807 sequences
Algorithm: FASTA (3.8 Nov 2011) [optimized]
Parameters: BL50 matrix (15:-5), open/ext: -10/-2
ktup: 2, E-join: 1 (0.519), E-opt: 0.2 (0.17), width: 16
Scan time: 0.310
The best scores are: opt bits E(16528)
Phybl2|Phybl2_154054 fgenesh1_pm.1_#_199 ( 212) 257 69.6 7.5e-13
>>Phybl2|Phybl2_154054 fgenesh1_pm.1_#_199 (212 aa)
initn: 289 init1: 243 opt: 257 Z-score: 339.0 bits: 69.6 E(16528): 7.5e-13
Smith-Waterman score: 272; 32.4% identity (56.8% similar) in 213 aa overlap (1-185:1-211)
10 20 30 40 50
tr|B0F MSLFIFVDESITNNRSSSSN---EEISSHQAD-----TTKEKTPRPKNAFMLYRQEKQKG
: :.::.: . :. . :. : ::: . .. : ::::::::::: . .
Phybl2 MESFMFVNE-FKNSYQPPSDILTESISSTPVRQKCIISSTVKIRRPKNAFMLYRQAVHPS
10 20 30 40 50
60 70 80 90 100
tr|B0F LPLSKTKLLSKDFSKIAAEMWRRESNEVRMYYHRLAEEEKLRHLAKHPGYKYCPK----R
. :.. . .:..:. :..::. :..::: ::.: :.:::: : : ::: : :. :
Phybl2 ILSSNSTIHNKEISRTAGKMWKNEKEEVRKYYERKADEEKLYHSKKFPGYIYKPQQRKTR
60 70 80 90 100 110
110 120 130 140 150
tr|B0F QSLAHVDKNEYQVLADTTPIVGEMDN-----------QMAFVH-RNVSF--DEVSPSSPP
. . : : . .: : . .. : ..: :. : ::. ..::
Phybl2 RPQSTVCKPFLRSTSDIQLIYQKTESKSPSEKLSATSQDYIIHCRDESTNEDEIL-NTPP
120 130 140 150 160 170
160 170 180
tr|B0F SSMVQD--LQIFPEGTNLDALSLVDYEFYLLLSIMS
.. ..: :. : . ..: .: .: .. : :
Phybl2 TQCLDDSILKTFLDDLSFDIFSPIDELIFSLCS*
180 190 200 210
188 residues in 1 query sequences
5831043 residues in 16528 library sequences
Tcomplib [36.3.7a Jan, 2015(preload9)] (32 proc in memory [0G])
start: Tue Oct 11 21:57:39 2016 done: Tue Oct 11 21:57:39 2016
Total Scan time: 0.310 Total Display time: 0.010
Function used was FASTA [36.3.7a Jan, 2015(preload9)]
>tr|B0F2G8|B0F2G8_PHYBL Sex minus OS=Phycomyces blakesleeanus GN=sexM PE=4 SV=1
MSLFIFVDESITNNRSSSSNEEISSHQADTTKEKTPRPKNAFML
YRQEKQKGLPLSKTKLLSKDFSKIAAEMWRRESNEVRMYYHRLAEEEKLRHLAKHPGY
KYCPKRQSLAHVDKNEYQVLADTTPIVGEMDNQMAFVHRNVSFDEVSPSSPPSSMVQD
LQIFPEGTNLDALSLVDYEFYLLLSIMS
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment