Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 720,466 | 0 | G | C | 55.0% | 2.5 / 23.9 | 20 | A23P (GCC→CCC) | phrB | deoxyribodipyrimidine photolyase, FAD‑binding |
Reads supporting (aligned to +/- strand): ref base G (8/1); new base C (0/11); total (8/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 110 total aligned reads displayed.
AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATG > REL606/720391‑720540 | aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg > 2:8361971/1‑76 (MQ=255) aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg > 2:8320164/1‑76 (MQ=255) aGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg < 2:6768984/76‑1 (MQ=255) gAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcc < 2:2673380/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct < 2:4372579/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcct < 1:2689819/76‑1 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg < 1:4442986/76‑1 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg < 1:2297103/76‑1 (MQ=255) ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑c > 1:1054547/1‑76 (MQ=255) ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑c > 2:1067235/1‑76 (MQ=255) ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑c > 1:181170/1‑76 (MQ=255) ttttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg‑c < 1:5089141/76‑1 (MQ=255) tttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccgcctg‑cc > 1:2825572/1‑76 (MQ=255) ttATGATTCCCCATCTGGTCTGGTTCCGCCAGAATTTACGTCGGCAGGATAATCTCGCACTGGctgccccctg‑ccg < 2:6710850/76‑1 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggggcccg‑ccc > 2:3581869/1‑75 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggccg‑ccc > 1:3606556/1‑75 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggccg‑ccc > 1:4978749/1‑75 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑ccc > 1:1965067/1‑75 (MQ=255) ttATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑ccc > 1:3709454/1‑75 (MQ=255) tATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑cccc > 2:1950501/1‑74 (MQ=255) aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcggccg‑ccgcc > 2:1795141/1‑75 (MQ=255) aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgccgg‑ccgcc > 1:3870114/1‑75 (MQ=255) aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑ccccc > 1:5760397/1‑73 (MQ=255) aTGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGGCCG‑CCGCa > 1:8332936/1‑76 (MQ=255) tGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑ccgcca > 1:959287/1‑74 (MQ=255) tGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑ACGCac > 1:3460797/1‑75 (MQ=255) gATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcgcccg‑ccgccac > 1:625903/1‑73 (MQ=255) aTTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCGG‑CCGCCAtt > 1:5664496/1‑76 (MQ=255) aTTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCCAtt > 1:5495811/1‑76 (MQ=255) ttACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCCATTc > 1:4255015/1‑76 (MQ=255) tACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCGCCCGCAAtcg > 2:6204402/1‑74 (MQ=255) ccATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCTG‑CCGCCATTCGTCt > 1:3953316/1‑76 (MQ=255) aTCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑GCGCCAATCGTCttc > 2:7039889/1‑74 (MQ=255) ctggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCCCCATTCGTCTGCAc > 1:321780/1‑76 (MQ=255) tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCG‑CCCCCCATCGTCTGCAcc > 2:3765643/1‑75 (MQ=255) tggtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCCCCCATCCTCTGCAcg > 1:6206513/1‑76 (MQ=255) gtctggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCG‑CCGCCCCTCGTCTTCAcgcc > 2:8350115/1‑75 (MQ=255) ggtTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCC‑CCGCCCATCGTCTTCACCCGTggt > 1:3988015/1‑74 (MQ=255) ggtTCCGCCAGGTTTTCCTTCAGCAGAATAATCCCGCACTGGCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCt < 2:7610845/76‑1 (MQ=255) tttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGGCCG‑CCGCAATTCGTCTGCCCGCGGGCTgg > 1:5043198/1‑76 (MQ=255) tttCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCCCCCATCGTCTTCTCGCGCGGTgg > 2:7644549/1‑76 (MQ=255) ttgcGCCAGTATTTAGGTCTGCAGAATACCCCCGCCTTGCCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTgg < 2:2244740/73‑1 (MQ=255) ttgcGCCAGGTTTTACGTCGGCGCGATAATCCCGCCCGGCCTCCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTgg < 1:3934289/73‑1 (MQ=255) cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCCATTCGTCTTCCCCCCGGCTGGCgg > 1:6555278/1‑75 (MQ=255) cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGt > 1:1798654/1‑76 (MQ=255) ccAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCCCCATTCGTCTTCACGCGGGGTGGCggtt > 1:813186/1‑73 (MQ=255) cAGGATTTCCTTCTGCACGATAATCCCGCCCGGGCGGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTtat < 1:910335/76‑1 (MQ=255) cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCAATTCGTCTGCACGCGCGCTGGCGCtat > 1:1213439/1‑76 (MQ=255) gtttttACGTCTGCAAGATACTCTCGCCGTGGCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTtatat < 2:271451/73‑1 (MQ=255) ggttttCCGTCTGCAGGATAATCCCGCACTGGCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTtatat < 2:4659464/73‑1 (MQ=255) ggttttCCGTCGGCACAATAATCGCGCCGTGGCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTtatat < 1:7668027/73‑1 (MQ=255) tttttCCGTCTGCACAATAATCTCGCCCTGGCCCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTtatata < 1:6713449/74‑1 (MQ=255) taTTTCCGTCTGCACAATATCCCCCCCCTGCCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTtatata < 1:4621101/75‑1 (MQ=255) ttttGCGCCTCCACAATAATCTCGCCCGGGCTCCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTatatat < 1:4066959/75‑1 (MQ=255) tttCCGTCTGCCAGATACTCCCGCCCTGCCCCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATc < 2:6478063/76‑1 (MQ=255) tttACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGGCAATCCTCTGCACGCGCGCTGGCGTTatatata > 1:1038078/1‑75 (MQ=255) ttACGTCTGCACGATAATCTCGCACTGGCTGGGGCCCG‑CCGCCAATCGTCTTCCCGCGCGCTGGGGGTatttatat > 2:2650637/1‑74 (MQ=255) tctgTCTGAACAATAATCGCGCCCTGCCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGc < 1:6604242/73‑1 (MQ=255) tcGTCCGCACGATAATCTCGCCGGGGCTGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt < 1:5633250/75‑1 (MQ=255) cGTCTGCCCGATACTCCCGCCCTGGCTGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTa < 2:3636592/76‑1 (MQ=255) tCTGCACGATAATCTCGCACTGGCTGCCGCCCG‑CCGCCATTCGTCTGCACCCCCGCTGGGGGTTTATATCGCTaca > 1:5586474/1‑76 (MQ=255) cTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCCATTCGTCTGCACGCGTGCTGGCGGTATATATCGCTacac > 1:1497442/1‑76 (MQ=255) cTGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCCCCAATCGTCTGCACCCGGGCTGGCGGTTTATATCGCTacac > 2:7923937/1‑76 (MQ=255) cTGCACGATAATCTCGCACTGGCTGCCGCCTC‑CCGCCATTCCCCTGCCCCCGTGCTGGGGTTTTATATTGCCacac > 1:4522078/1‑76 (MQ=255) tccACGATACCCGCGCTCTCGCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc < 1:4419467/74‑1 (MQ=255) tGCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCCCTTCGTCTGCACGCGGGCTGGCGGTATATATTTCTAcacc > 2:4318216/1‑76 (MQ=255) ggCACGATACCCCCGCCCTGCCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc < 2:6340055/75‑1 (MQ=255) gccacaATAATCCCGCACTGCCGCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc < 1:6899750/74‑1 (MQ=255) gccACGATATCCTCGCCGGGCCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacc < 2:791374/74‑1 (MQ=255) gCACGATAATCTCGCACTGGCTGGCGCCCG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcacca > 1:1434100/1‑76 (MQ=255) gCACGATAATCTCGCACTGGCTGCCGCCTG‑CCGCAATTCGTCTTCCCCCGTGCTGGCGGTTTATTTCGCTAcacca > 1:1199857/1‑76 (MQ=255) cACTATAATCCCGCCCTGGCTGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTAcaccac < 1:3829329/76‑1 (MQ=255) aCGATAATCTCGCCCTGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg < 2:2567412/76‑1 (MQ=255) aCGATAATCTCGCACTGGCTGGCGCCCG‑CCCCCAATCGTCTTCTCCCCCGCTGGCGGTTTATATTTCTTcaccacc > 2:3348617/1‑75 (MQ=255) caaTAGTCTCGCCGTGGCGGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc < 2:4929666/74‑1 (MQ=255) cGATAATCTCGCACTGGCTGGCGCCTG‑CCGCAATTCGTCTGCACGCGTGCCGGCGTTATATATTGCTAcaccaacc > 1:707895/1‑73 (MQ=255) gATAATCTCGCACTGGCTGGCGCCCG‑CCGCCAATCGTCTGCCCCCGCGCTTGGGGTTTATATCGCTACACCAcccc > 2:4372231/1‑76 (MQ=255) aTAATCTCGCACTGGCTGGCCCCCG‑CCCCCATTTCTCTGCACGCCGGCCGGGGGTATATATTGCTTCACCACAccc > 2:6661725/1‑75 (MQ=255) tAATCTCGCACTGGCTGGCGCCCG‑CCGCCATTCGTCTTCACCCGGGCTGGCGGTTTATATCGCTTCACCCCCCCAg > 2:2871819/1‑76 (MQ=255) aaTCTCGCACTGGCTGGCGCCCG‑CCCCCATTTGTCTGCACCCGTGCTTGGGGTTTATATCGCTACACCCCCCCaca > 2:8152880/1‑74 (MQ=255) aTCTCGCCCGGGCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTg < 2:8320882/76‑1 (MQ=255) ccccgcgCTGGCTGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTgg < 2:5227659/73‑1 (MQ=255) ccgcTGGCCCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGa < 2:2923935/73‑1 (MQ=255) cTGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCCGTGGGCGACGc > 1:3169068/1‑76 (MQ=255) cTGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGc > 2:110568/1‑76 (MQ=255) cTGGCTGACCCCTG‑CCGCAATTCGTCAGCACGCCTGCTGGCGTTATATATCGCTGAAACACGCAAGCGGGCGACGc > 2:2030324/1‑76 (MQ=255) tGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACCCCACGCCAGGGGGCGACGCa > 2:2857306/1‑76 (MQ=255) ggCTCCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt < 2:6202407/76‑1 (MQ=255) ggCCGCCCCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt < 2:4123238/76‑1 (MQ=255) gCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATa > 1:1632217/1‑76 (MQ=255) tgccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 2:6025175/1‑76 (MQ=255) tgccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 2:825339/1‑76 (MQ=255) ccccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc < 1:6285822/74‑1 (MQ=255) gccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 2:6367300/1‑76 (MQ=255) gccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 1:5817405/1‑76 (MQ=255) gccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 1:3725496/1‑76 (MQ=255) ccgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt > 1:3662379/1‑76 (MQ=255) cgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATg > 2:5697159/1‑76 (MQ=255) cgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATg > 2:4480653/1‑76 (MQ=255) cgcctg‑ccgcAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATg > 1:6665703/1‑76 (MQ=255) | AGGAGTTTTATGATTACCCATCTGGTCTGGTTTCGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCGCCTG‑CCGCAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACATG > REL606/720391‑720540 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |