Multiple sequence alignment - TraesCS2A01G256100

Loading Multiple Alignment...

Download Multiple Sequence Alignment: (Fasta) (Clustal)



BLAST Results


Shown here are the results of a BLAST search against the genome using the input sequence. The pipeline uses this to extract flanking regions for cloning. If the result of the multiple sequence alignment is not as expected, these results could be used for further investigation.

qseqid sseqid percentage.identical alignment.length no.mismatch no.gap.openings qstart qend sstart send evalue bitscore
0 TraesCS2A01G256100 chr2A 100.000 6457 0 0 1 6457 392723405 392729861 0.000000e+00 11924.0
1 TraesCS2A01G256100 chr2D 93.583 3226 99 51 2677 5858 325363069 325359908 0.000000e+00 4711.0
2 TraesCS2A01G256100 chr2D 89.860 2495 103 56 262 2679 325365596 325363175 0.000000e+00 3068.0
3 TraesCS2A01G256100 chr2D 88.276 145 10 5 71 210 325365759 325365617 4.000000e-37 167.0
4 TraesCS2A01G256100 chr2B 95.320 2372 61 26 3495 5858 394361954 394359625 0.000000e+00 3720.0
5 TraesCS2A01G256100 chr2B 90.604 1937 72 45 381 2286 394365006 394363149 0.000000e+00 2468.0
6 TraesCS2A01G256100 chr2B 85.676 1117 27 40 2296 3357 394363108 394362070 0.000000e+00 1053.0
7 TraesCS2A01G256100 chr2B 81.215 181 26 7 5959 6135 109330857 109331033 8.730000e-29 139.0
8 TraesCS2A01G256100 chr2B 93.478 92 2 2 104 192 394365574 394365484 4.060000e-27 134.0
9 TraesCS2A01G256100 chr2B 89.286 56 5 1 5961 6015 779309040 779308985 1.160000e-07 69.4
10 TraesCS2A01G256100 chr6B 91.618 346 21 4 5973 6318 132949625 132949962 7.580000e-129 472.0
11 TraesCS2A01G256100 chr6B 82.464 211 37 0 1156 1366 574460723 574460933 1.110000e-42 185.0
12 TraesCS2A01G256100 chr6B 76.610 295 67 2 3544 3837 61224420 61224713 1.860000e-35 161.0
13 TraesCS2A01G256100 chr3A 93.706 286 16 2 5973 6257 159996537 159996821 1.660000e-115 427.0
14 TraesCS2A01G256100 chr3A 95.192 208 3 3 6256 6457 159997059 159997265 8.080000e-84 322.0
15 TraesCS2A01G256100 chr5B 83.878 459 54 10 6009 6457 498834209 498833761 2.780000e-113 420.0
16 TraesCS2A01G256100 chr6A 83.412 211 35 0 1156 1366 526115974 526116184 5.110000e-46 196.0
17 TraesCS2A01G256100 chr6A 77.365 296 63 4 3544 3837 33584217 33584510 8.610000e-39 172.0
18 TraesCS2A01G256100 chr6A 78.210 257 56 0 3544 3800 33541690 33541434 1.440000e-36 165.0
19 TraesCS2A01G256100 chr6D 82.464 211 37 0 1156 1366 384218667 384218877 1.110000e-42 185.0
20 TraesCS2A01G256100 chr6D 78.327 263 57 0 3538 3800 28999646 28999908 3.100000e-38 171.0
21 TraesCS2A01G256100 chr6D 78.138 247 52 2 3544 3789 28941829 28941584 8.670000e-34 156.0
22 TraesCS2A01G256100 chr7B 78.889 180 32 5 5958 6135 66969847 66970022 4.090000e-22 117.0
23 TraesCS2A01G256100 chr3B 91.228 57 3 2 5972 6026 42602669 42602613 6.940000e-10 76.8
24 TraesCS2A01G256100 chr3B 90.909 55 5 0 5962 6016 130471079 130471133 2.500000e-09 75.0
25 TraesCS2A01G256100 chr1A 95.556 45 2 0 5971 6015 137874078 137874122 8.980000e-09 73.1
26 TraesCS2A01G256100 chr4D 90.566 53 4 1 5964 6015 51993122 51993174 1.160000e-07 69.4


These results show the homologous regions that are presented in the multiple sequence alignment. HSPs from the BLAST search that are within 2000 bases of each other have been grouped together. Homologous regions that are purely upstream or downstream may be removed if they don't overlap both primers.

query scaffold start end length rev.comp avg.bitscore max.bitscore avg.percent.identical query.start query.end num_hsp groupid homo_length
0 TraesCS2A01G256100 chr2A 392723405 392729861 6456 False 11924.000000 11924 100.0000 1 6457 1 chr2A.!!$F1 6456
1 TraesCS2A01G256100 chr2D 325359908 325365759 5851 True 2648.666667 4711 90.5730 71 5858 3 chr2D.!!$R1 5787
2 TraesCS2A01G256100 chr2B 394359625 394365574 5949 True 1843.750000 3720 91.2695 104 5858 4 chr2B.!!$R2 5754
3 TraesCS2A01G256100 chr3A 159996537 159997265 728 False 374.500000 427 94.4490 5973 6457 2 chr3A.!!$F1 484


AutoCloner calculated primer pairs

These primer pairs were selected by Autocloner to minimize primer penalty (calculated by Primer3), whilst remaining within the specified product range where possible.

Download AutoCloner chosen primers as CSV: (Forward) (Reverse)
Position MSA Position Penalty Sequence TM GC Self any TH Self end TH Hairpin End Stability Orientation
217 222 0.104304 CCCGACTGAATCTTACCCCG 59.896 60.0 0.00 0.00 0.00 5.73 F
993 1305 0.112606 ATCTAGCCGGCTGGTAGCTA 59.887 55.0 38.98 19.88 43.35 3.32 F
1484 1813 0.588233 CGTACGTACGCTGCACCTAG 60.588 60.0 32.36 5.42 43.14 3.02 F
1680 2010 0.724453 GCCGTACGTGCATGTGTTTG 60.724 55.0 21.92 8.35 0.00 2.93 F
2933 3481 0.523072 ACGCAGTGATGAAATGGTGC 59.477 50.0 0.00 0.00 42.51 5.01 F
3524 4115 0.038526 GACTGACGGACGGATGTTGT 60.039 55.0 0.00 0.00 0.00 3.32 F
4565 5156 0.173708 GGCTCTTCTACACGGACCTG 59.826 60.0 0.00 0.00 0.00 4.00 F
Position MSA Position Penalty Sequence TM GC Self any TH Self end TH Hairpin End Stability Orientation
1375 1702 0.460987 AGCGAGGAAGGAAATCTGCG 60.461 55.0 0.00 0.0 0.00 5.18 R
2566 2972 0.043334 AGAGACCTACCTGCCCACAT 59.957 55.0 0.00 0.0 0.00 3.21 R
2579 2985 0.173708 CGGCTGGACTACAAGAGACC 59.826 60.0 0.00 0.0 40.14 3.85 R
3505 4096 0.038526 ACAACATCCGTCCGTCAGTC 60.039 55.0 0.00 0.0 0.00 3.51 R
4602 5193 1.123077 TGGATGGACGGGATCTTCTG 58.877 55.0 0.00 0.0 0.00 3.02 R
4926 5517 1.208642 TGCATCGTCATCGTCATCGC 61.209 55.0 0.00 0.0 38.33 4.58 R
5969 6575 0.037605 GCGCGGGTATGGCTAGTTAT 60.038 55.0 8.83 0.0 0.00 1.89 R


All possible primers

Listed here are all possible primers based on the multiple sequence alignment and SNPs generated previously. A lower penalty score indicates a better primer. If the AutoCloner primers did not work during PCR or sequencing, these could be used as alternatives.

Download all possible primers: (Forward) (Reverse)



Based at the University of Bristol with support from BBSRC. BBSRC icon Bristol icon

AutoCloner maintained by Alex Coulton.