e:0.002479 screening
>Contig1 with 5 hitsequence seqlen:1288
UniRef100_O82531 Score:346 Escore:0.000000 Evalue:8.7e-93 querystart:168 querystop:836 sbjctstart:1
UniRef100_A7PY95 Score:334 Escore:0.000000 Evalue:3.9e-89 querystart:168 querystop:836 sbjctstart:1
UniRef100_P42742 Score:311 Escore:0.000000 Evalue:2.1e-82 querystart:168 querystop:836 sbjctstart:1
UniRef100_Q2V3M5 Score:307 Escore:0.000000 Evalue:3.7e-81 querystart:168 querystop:836 sbjctstart:1
UniRef100_Q93X30 Score:295 Escore:0.000000 Evalue:1.4e-77 querystart:192 querystop:767 sbjctstart:1
Contig1 orf: 3 orfshift:0
Contig1 with right orf:3 querystart:168
start in hit region
orf:start:168 stop:836 sureforstart:1 length of the gene:668 suregene
contig_1 with 8 realsnp,with 4 truesnp
snp_278_CT is synonymous snp protein:RMSTGYN
snp_456_CT is non-synonymous snp protein1:MGQLLSN protein2:MGQFLSN
snp_591_CT is synonymous snp protein:RVGYSSQ
snp_669_GT is synonymous snp protein:PLLLPAV
There are 0 bi-indelsnp
There are 4 bi-snp without indel in hit region, 1 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region not hit region,0 of that non-synonymoussnp
There are 4 bi-snp without indel in coding region,1 of that non-synonymoussnp
>Contig5 with 5 hitsequence seqlen:898
UniRef100_Q3HRY7 Score:322 Escore:0.000000 Evalue:8.5e-86 querystart:108 querystop:725 sbjctstart:1
UniRef100_Q6TKT4 Score:306 Escore:0.000000 Evalue:5.8e-81 querystart:108 querystop:707 sbjctstart:1
UniRef100_A5CB91 Score:291 Escore:0.000000 Evalue:1.6e-76 querystart:105 querystop:725 sbjctstart:1
UniRef100_A7PRE6 Score:290 Escore:0.000000 Evalue:3.6e-76 querystart:105 querystop:725 sbjctstart:1
UniRef100_P41127 Score:286 Escore:0.000000 Evalue:4.7e-75 querystart:108 querystop:725 sbjctstart:1
Contig5 orf: 3 orfshift:0
Contig5 with right orf:3 querystart:108
start in hit region
orf:start:108 stop:725 sureforstart:1 length of the gene:617 suregene
contig_5 with 9 realsnp,with 6 truesnp
snp_182_TC is synonymous snp protein:TWFNQPA
snp_344_GC is synonymous snp protein:LEELKAA
snp_368_GA is synonymous snp protein:IPKKLAP
snp_447_CA is synonymous snp protein:NVQRLKT
snp_575_TC is synonymous snp protein:LPITREQ
snp_605_AG is synonymous snp protein:DFVKVTD
There are 0 bi-indelsnp
There are 6 bi-snp without indel in hit region, 0 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region not hit region,0 of that non-synonymoussnp
There are 6 bi-snp without indel in coding region,0 of that non-synonymoussnp
>Contig6 with 5 hitsequence seqlen:990
UniRef100_Q307X7 Score:267 Escore:0.000000 Evalue:2.4e-69 querystart:112 querystop:630 sbjctstart:1
UniRef100_Q6RK02 Score:258 Escore:0.000000 Evalue:1.4e-66 querystart:112 querystop:630 sbjctstart:1
UniRef100_Q947H2 Score:256 Escore:0.000000 Evalue:4.8e-66 querystart:112 querystop:636 sbjctstart:1
UniRef100_A7PHG3 Score:247 Escore:0.000000 Evalue:3.9e-63 querystart:109 querystop:651 sbjctstart:33
UniRef100_A5B412 Score:245 Escore:0.000000 Evalue:1e-62 querystart:112 querystop:627 sbjctstart:1
Contig6 orf: 1 orfshift:0
Contig6 with right orf:1 querystart:112
start in hit region
orf:start:112 stop:630 sureforstart:0 length of the gene:518 suregene
contig_6 with 46 realsnp,with 22 truesnp
snp_79_-G Indelsnp
snp_80_-A Indelsnp
snp_81_A- Indelsnp
snp_82_G- Indelsnp
snp_83_G- Indelsnp
snp_84_A- Indelsnp
snp_85_A- Indelsnp
snp_86_G- Indelsnp
snp_87_C- Indelsnp
snp_88_A- Indelsnp
snp_89_G- Indelsnp
snp_98_AG not in hit region,in UTR
snp_103_TA not in hit region,in UTR
snp_324_GA is synonymous snp protein:FCRGVGR
snp_413_GA is synonymous snp protein:ILDLLKN
snp_435_CT is synonymous snp protein:AESNAEV
snp_510_GA is synonymous snp protein:QKQRRRT
snp_531_CT is synonymous snp protein:YRAHGRI
snp_586_TC is synonymous snp protein:LILSEKE
snp_616_GT is synonymous snp protein:KKEPETQ
snp_661_A- Indelsnp
snp_679_AT not in hit region,in UTR
There are 12 bi-indelsnp
There are 7 bi-snp without indel in hit region, 0 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region not hit region,0 of that non-synonymoussnp
There are 7 bi-snp without indel in coding region,0 of that non-synonymoussnp
>Contig8 with 5 hitsequence seqlen:1201
UniRef100_Q38JI1 Score:471 Escore:0.000000 Evalue:2.9e-130 querystart:177 querystop:1186 sbjctstart:1
UniRef100_Q8VXC9 Score:471 Escore:0.000000 Evalue:2.9e-130 querystart:177 querystop:1186 sbjctstart:1
UniRef100_Q8VXD0 Score:470 Escore:0.000000 Evalue:6.2e-130 querystart:177 querystop:1186 sbjctstart:1
UniRef100_Q9ZRR5 Score:470 Escore:0.000000 Evalue:8.4e-130 querystart:177 querystop:1186 sbjctstart:1
UniRef100_Q8VXD1 Score:469 Escore:0.000000 Evalue:1.5e-129 querystart:177 querystop:1186 sbjctstart:1
Contig8 orf: 3 orfshift:2
ORF shift loci in sequence(not exact):1122 1150 1122 uncorrected 1150 uncorrected 
Contig8 with wrong seqence not because Indel
>Contig12 with 5 hitsequence seqlen:1143
UniRef100_Q38JI1 Score:489 Escore:0.000000 Evalue:1.1e-135 querystart:65 querystop:1125 sbjctstart:96
UniRef100_Q8VXC9 Score:489 Escore:0.000000 Evalue:1.1e-135 querystart:65 querystop:1125 sbjctstart:96
UniRef100_Q2MIR4 Score:487 Escore:0.000000 Evalue:3.8e-135 querystart:62 querystop:1119 sbjctstart:95
UniRef100_Q8VXD1 Score:487 Escore:0.000000 Evalue:6.1e-135 querystart:65 querystop:1122 sbjctstart:96
UniRef100_Q9FRU1 Score:486 Escore:0.000000 Evalue:8.2e-135 querystart:65 querystop:1134 sbjctstart:96
Contig12 orf: 2 orfshift:3
ORF shift loci in sequence(not exact):923 1083 1105 923 uncorrected 1083 uncorrected 1105 uncorrected 
Contig12 with wrong seqence not because Indel
>Contig13 with 5 hitsequence seqlen:972
UniRef100_Q1SN35 Score:213 Escore:0.000000 Evalue:5.7e-53 querystart:96 querystop:587 sbjctstart:1
UniRef100_Q2XPV7 Score:212 Escore:0.000000 Evalue:6.8e-53 querystart:96 querystop:606 sbjctstart:1
UniRef100_O65743 Score:207 Escore:0.000000 Evalue:3.2e-51 querystart:96 querystop:584 sbjctstart:1
UniRef100_A7PPN3 Score:205 Escore:0.000000 Evalue:1.3e-50 querystart:96 querystop:587 sbjctstart:1
UniRef100_Q9FUL4 Score:205 Escore:0.000000 Evalue:1.5e-50 querystart:96 querystop:616 sbjctstart:1
Contig13 orf: 3 orfshift:0
Contig13 with right orf:3 querystart:96
start in hit region
orf:start:96 stop:587 sureforstart:1 length of the gene:491 suregene
contig_13 with 12 realsnp,with 8 truesnp
snp_80_-G Indelsnp
snp_81_-C Indelsnp
snp_82_-T Indelsnp
snp_83_-C Indelsnp
snp_84_-G Indelsnp
snp_85_-C Indelsnp
snp_86_-T Indelsnp
snp_87_-A Indelsnp
There are 8 bi-indelsnp
There are 0 bi-snp without indel in hit region, 0 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region not hit region,0 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region,0 of that non-synonymoussnp
>Contig15 with 5 hitsequence seqlen:873
UniRef100_Q1SN35 Score:211 Escore:0.000000 Evalue:1.5e-52 querystart:100 querystop:591 sbjctstart:1
UniRef100_Q2XPV7 Score:211 Escore:0.000000 Evalue:1.8e-52 querystart:100 querystop:610 sbjctstart:1
UniRef100_Q9FUL4 Score:206 Escore:0.000000 Evalue:5e-51 querystart:100 querystop:632 sbjctstart:1
UniRef100_O65743 Score:205 Escore:0.000000 Evalue:8e-51 querystart:100 querystop:588 sbjctstart:1
UniRef100_A7PPN3 Score:203 Escore:0.000000 Evalue:3.3e-50 querystart:100 querystop:591 sbjctstart:1
Contig15 orf: 1 orfshift:0
Contig15 with right orf:1 querystart:100
start in hit region
orf:start:100 stop:591 sureforstart:1 length of the gene:491 suregene
contig_15 with 7 realsnp,with 3 truesnp
snp_81_C- Indelsnp
snp_82_G- Indelsnp
snp_83_G- Indelsnp
There are 3 bi-indelsnp
There are 0 bi-snp without indel in hit region, 0 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region not hit region,0 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region,0 of that non-synonymoussnp
>Contig16 with 5 hitsequence seqlen:1150
UniRef100_Q2V991 Score:439 Escore:0.000000 Evalue:1.2e-120 querystart:133 querystop:1108 sbjctstart:1
UniRef100_Q9SXS0 Score:436 Escore:0.000000 Evalue:6.4e-120 querystart:133 querystop:1108 sbjctstart:1
UniRef100_Q40403 Score:422 Escore:0.000000 Evalue:9.6e-116 querystart:133 querystop:1108 sbjctstart:1
UniRef100_Q9SXR9 Score:419 Escore:0.000000 Evalue:9.1e-115 querystart:133 querystop:1108 sbjctstart:1
UniRef100_P93340 Score:414 Escore:0.000000 Evalue:3.1e-113 querystart:133 querystop:1108 sbjctstart:1
Contig16 orf: 1 orfshift:2
ORF shift loci in sequence(not exact):1003 1073 1003 uncorrected 1073 uncorrected 
Contig16 with wrong seqence not because Indel
>Contig23 with 5 hitsequence seqlen:1720
UniRef100_A7P5C8 Score:190 Escore:0.000000 Evalue:2e-45 querystart:676 querystop:1359 sbjctstart:136
UniRef100_Q4VT47 Score:189 Escore:0.000000 Evalue:2.3e-45 querystart:676 querystop:1359 sbjctstart:136
UniRef100_A5AI54 Score:188 Escore:0.000000 Evalue:5.2e-45 querystart:685 querystop:1359 sbjctstart:144
UniRef100_Q94IC5 Score:173 Escore:0.000000 Evalue:1.8e-40 querystart:676 querystop:1359 sbjctstart:119
UniRef100_A7PEQ3 Score:171 Escore:0.000000 Evalue:6.3e-40 querystart:604 querystop:1347 sbjctstart:84
Contig23 orf: 1 orfshift:0
Contig23 with right orf:1 querystart:676
start in hit region
orf:start:214 stop:1383 sureforstart:1 length of the gene:1169 suregene
contig_23 with 8 realsnp,with 3 truesnp
snp_767_-A Indelsnp
snp_1155_CG is non-synonymous snp protein1:PYALHFC protein2:PYAVHFC
snp_1231_AT is synonymous snp protein:GTKVNAV
There are 1 bi-indelsnp
There are 2 bi-snp without indel in hit region, 1 of that non-synonymoussnp
There are 0 bi-snp without indel in coding region not hit region,0 of that non-synonymoussnp
There are 2 bi-snp without indel in coding region,1 of that non-synonymoussnp
**********
There are total:9 contig for blast:9 with hit protein and 0 nohit,3 contig with orfshift
There are total:6 contig with right orf,0 contigs with almost right orf and 3 contig with wrong orf
There are toal:90 snp(from realsnp with sequence errors),total 46 truesnp without sequence errors, total:24 bi-indel-snp,total:19 bi-snp without indel with in hit region and 2 bi-snp is non-synonymous
There are total:0 bi-snp not in hit region but in coding region and 0 of that non-synonymous
There are total:19 bi-snp in coding region and 2 of that non-synonymous
There are total:3 bi-snp in 5'UTR and 0 bi-snp in 3'UTR(not including indel-snp)
There are total:46 bi-snp and other snp:44
***add some statistic****
Noorf(orf==0)contig:0, contigwith uncorrectorfshrift:3, M in hitregion,but gene shorter than and start far to hit start 33:0, M not(* in coding) in hit region:0,M in hit region:0
Results analysis******************
SNP in coding region:0.41 (19/46),of which,0.11 (2/19)are nsSNP
SNP in UTR:0.07 (3/46)
Indel SNP in hitcontig:0.52 (24/46)
SNP in fasty hit region:0.41 (19/46)of which,0.11 (2/19)are nsSNP
SNP in coding region not hit region:0.00 (0/46)of which,nan (0/0)are nsSNP
SNP in hit region:1.00 (19/19) 0f all in coding regionsof whic,:1.00 (2/2)are nsSNP in hit region of all nsSNP in coding regions
