Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 720,464 | 0 | C | G | 51.2% | 5.3 / 67.9 | 41 | A22G (GCC→GGC) | phrB | deoxyribodipyrimidine photolyase, FAD‑binding |
Reads supporting (aligned to +/- strand): ref base C (3/17); new base G (21/0); total (24/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.52e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 127 total aligned reads displayed.
GAGTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT > REL606/720393‑720539 | gAGTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg‑cc > 1:5456477/1‑76 (MQ=255) gAGTTTTATGATTACCCATCTGGTCTGGTTG‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg‑cc < 1:62057/76‑1 (MQ=255) ggTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg‑cct < 2:6740361/75‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg‑cct < 1:4087897/76‑1 (MQ=255) aGTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg‑cct < 1:3762803/76‑1 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg‑cccg > 2:4814164/1‑74 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgcgg‑ccgc > 2:3881296/1‑74 (MQ=255) gTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctgccg‑cctg < 2:3906749/76‑1 (MQ=255) ttATAATAACCCATCTGGTCTGGTTC‑CGCCAGGTTTTACTTCTCCACAATAATCGCGCCCTGGctcccc‑cctg‑ccg < 2:332883/76‑1 (MQ=255) aTGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg‑cccg‑ccccc > 1:1628506/1‑73 (MQ=255) tGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGctggcg‑cccg‑ccg‑cca > 1:1474421/1‑74 (MQ=255) aTTCCCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCGTAATAATCCCGCACGGGCTCCCG‑CCTG‑CCG‑CAAtt < 2:5798714/76‑1 (MQ=255) ttACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑GCCG‑CCC‑CAAttt > 1:5104896/1‑75 (MQ=255) ttACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAAATc > 1:6802849/1‑76 (MQ=255) tcccccATCTGGTCTGGTTC‑CGCCAGAATTTCCGTCTCCACGATAACCTCGCCCGGCCGCCCC‑CCTG‑CCG‑CAATTc < 2:6952095/73‑1 (MQ=255) tACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGG‑CCCG‑CCG‑CCACTCg > 2:3969474/1‑76 (MQ=255) aCCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTCGt > 1:5465381/1‑76 (MQ=255) aCCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCCG‑CCG‑CAACTCGt > 1:1155687/1‑76 (MQ=255) cccATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCACTCGTc > 2:6437505/1‑76 (MQ=255) ccATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCGCCCGG‑CCG‑GAATTctcc > 1:6368457/1‑73 (MQ=255) ccATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCC‑CCG‑CAAATCGTCt > 2:2696312/1‑76 (MQ=255) ccATCT‑GTTTGGTTTCCGCGGGGTTTTCCGTCTGCACGATACTCTCGCCCTGGCCCCCC‑CCTG‑CCG‑CAATTCGTCt < 2:6917553/76‑1 (MQ=255) cATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCC‑CCG‑CCATTCGTccg > 1:58999/1‑74 (MQ=255) ttctggttgggTTT‑CGCCGGGTTTTCGGTCTCCAGAATATCCTCGCCCGGCCTCCCC‑CCTG‑CCG‑CAATTCGTCTGc < 2:185913/75‑1 (MQ=255) aTCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCTG‑CCC‑CAATTCGTCTGc > 1:1967393/1‑76 (MQ=255) tctggtctggTTC‑CGCCGGTATTGACGTCTGCACAATAGTCTCGCGCGGGCCCCCG‑CCTG‑CCG‑CAATTCGTCTGCa < 2:6973233/76‑1 (MQ=255) ctggtctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTG‑CCC‑CCAATCGTCTGCAc > 1:6522306/1‑76 (MQ=255) ctggtctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAATTCGTCTGCAc > 1:3683079/1‑76 (MQ=255) tgttctggtTT‑CCCCAGGTTTTACGTCTGCACGATAATCTCGCACGGGCGGCCC‑CCTG‑CCG‑CAATTCGTCTGCAcg < 1:1648199/73‑1 (MQ=255) ggtctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTCGTCTGCAcgc > 1:7353727/1‑76 (MQ=255) gtctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGG‑GCCG‑CCC‑CACATCGTCTTCAcccc > 2:6481680/1‑73 (MQ=255) tctgttTT‑CGCCAGGTTTTCCTTCGGCACGATAATCTCCCCCTGGCGCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGt < 1:2234999/76‑1 (MQ=255) tctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑GCCG‑CCG‑CCATTCGCCTTCAcacgg > 2:7778770/1‑75 (MQ=255) tctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAACTCGTCTGCCcgcct > 2:1488831/1‑74 (MQ=255) ctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGG‑CCCG‑GCG‑CAATTAGACTCCAcacgcg > 1:4828091/1‑74 (MQ=255) ctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTCGTCTGCAcgcgcg > 1:4426366/1‑74 (MQ=255) ctggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAACTCGTCTCCCCGCGTg > 1:1269512/1‑76 (MQ=255) tggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCGG‑CCG‑CCATTCGTCTGCCCGCGTGc > 2:5888788/1‑76 (MQ=255) ggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGGG‑CCCG‑CCGCCATTTCCTC‑GCACGCGTGcg > 2:1017563/1‑75 (MQ=255) ggtTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCGC‑CCG‑CGATTCGTCTCCCCGCGGGcc > 1:5603694/1‑75 (MQ=255) gtTT‑CGCCAGGATTTACGTCGGCACGATATCCTCGCCCTGGCGCCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTg < 1:2522464/76‑1 (MQ=255) ttt‑CGCCGTGATTTCCGCCGGCACGATATCCCCCCCCTGGCCCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTgg < 1:4821311/76‑1 (MQ=255) ttt‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑GCG‑CCCTTTGTCTGCCAGCGGGCTgg > 1:1034435/1‑76 (MQ=255) tt‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTG‑CCC‑CCAATTGTCTGCACCCGGGGTggg > 1:6140825/1‑75 (MQ=255) cGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTG‑CCG‑CCATTTGTCTGCACGCGggcgggcgg > 1:6543213/1‑75 (MQ=255) gcgaGGTTTTGCGTCGCCAGGATAACCCCGCCCTGGCCCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGtt < 2:6317960/73‑1 (MQ=255) gCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTG‑GCC‑CCATTTGTCTGGCAGCGGGCTGGCggt > 1:1414290/1‑74 (MQ=255) ccaaGGTTTTCCGTCTGCAGGATACTCTCGCCCTGGCTGCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGtt < 2:3157627/73‑1 (MQ=255) ccAGGTTTTCCGTCTCCACAATAATCCCGCCCTGGCCCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTta < 1:1570561/76‑1 (MQ=255) cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTG‑CCG‑CCATTCGTCTGCCCGCGTGCTGGCGGtat > 1:1066023/1‑76 (MQ=255) cAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCC‑CCG‑CCATTCGTCTTCACGCGGGCTGGCGGtat > 2:4317243/1‑76 (MQ=255) cAGGATTTACGTCTCCACGATAATCGCGCCCTGGCTCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTtat < 1:7175735/76‑1 (MQ=255) aGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCCG‑CCC‑CCATTTCTCTGCACCCGCGCGGGGGGtata > 1:646221/1‑76 (MQ=255) gATTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑GCCC‑CCG‑CAAATTCTCTGTCCGCGTGGTGGCGTgatata > 1:5351427/1‑76 (MQ=255) aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCTG‑CCG‑CAAATCATCTTCCCCCGCGCTGGCGttatttat > 1:51721/1‑76 (MQ=255) aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCCATTGTCTGCACCCGCGCCGGCGTTatatat > 2:3505309/1‑76 (MQ=255) aTTTACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCAATCGTCTTCACCCGTGCTGGCGGtgtattt > 2:6121015/1‑74 (MQ=255) tttACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTCGTCTGCCCGCGTGCCGGGGTTATATATc > 1:7040742/1‑76 (MQ=255) tttACGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAATTCGTCTGCACGCGGGCTGGGGTTATATATc > 1:7397900/1‑76 (MQ=255) tACGTCTGCACGATAATCTCGCACTGGCTGGGG‑CCCG‑CCC‑CCAATTGTCTGCCCGCGCGCCGGGGGTTTATATCGc > 2:2217843/1‑76 (MQ=255) aCGTCTGCAGAATACTCCCGCCGGGGCTGCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCt < 1:2165400/76‑1 (MQ=255) cGTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCC‑CCATTCGCCTTCCCGCGCGGTGGCGGTATATATCGCTa > 1:6476993/1‑76 (MQ=255) gTCTGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAAATTGTCTGCACGCGTGCTGGCGGTTTATATAGCtca > 1:1618321/1‑74 (MQ=255) tCTGCACTATACTCTCGCACGGGCTCCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTaca < 2:3943489/76‑1 (MQ=255) cTGCAGAATAACCTCGCACTGGCTGCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTacac < 1:1933531/76‑1 (MQ=255) tGCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTCGTCTGCACGCGTGCTGGCGTTTTATATCCCTAcacc > 1:6940487/1‑76 (MQ=255) gCACGATAATCTCGCACTGGCTGGCG‑CCTG‑CCG‑CCATTCGTCTGCACGCGGGCTGGGGTTTTATATCGCTAcaccc > 1:4675020/1‑75 (MQ=255) gCACGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CAATTCGTCTGCACCCGCGCTGGCGGTATATATCGCTCcaccc > 2:5902886/1‑75 (MQ=255) gCACGATAATCTCGCACTGGCTGGCG‑CCCC‑CCG‑CCAATCGTCTGCCACCGCGCTGGCGGTTTATATTACTAcacca > 2:7308785/1‑76 (MQ=255) aggATAATCGCGCCCGGGCTGCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg < 2:2291501/74‑1 (MQ=255) aCGATATCCTCGCCCTGCCCGCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg < 1:5384292/76‑1 (MQ=255) aCGATAATCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTTCTCTGCACGCCTGGCGGCGGTTTTTATTGCTAcaccacc > 2:6529214/1‑75 (MQ=255) aCGATAATCCCGCACGGGCCCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACg < 2:2030947/76‑1 (MQ=255) caaTAATCCCCCACTGCCGGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc < 2:6641347/74‑1 (MQ=255) agATAATCTCCCCCGGGCCCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGc < 1:573052/75‑1 (MQ=255) tAATCTCGCACTGGCTCCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg < 2:5170626/76‑1 (MQ=255) ataTCTCGCTCTGCCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg < 2:1348046/74‑1 (MQ=255) aaaTCTCGCACTGGCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAg < 2:2346061/75‑1 (MQ=255) aaTCTCGCACTGGCTGGCG‑CCCG‑CCG‑CCATTCGTCTGCACGCGGGCTGGCGGTATATATCGCTACACCAccccccg > 2:2581693/1‑73 (MQ=255) aaTCTCGCACTGGCTGGCG‑CCCG‑CCC‑CAAATCGTCTGGCCGCGTGCTGGGGTTATATATTGCTACACCACCCCacg > 2:7560845/1‑74 (MQ=255) ctcGCAGTGGCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTggg < 2:6649312/76‑1 (MQ=255) ctcGCACTGGCTGGCG‑CCTGCCCC‑CAATTC‑TTTTCCCCCGTGGTGGCGCTTTATATAGCTACACCACGCCACTggg > 2:5465914/1‑76 (MQ=255) tcGCCCGGGCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGc < 2:4328116/76‑1 (MQ=255) cGCCCTGGCCCCCC‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCg < 1:7825954/76‑1 (MQ=255) cGCACTGGCTGGGG‑CCCG‑CCG‑CAATTCGTCTGCAAGCGGGCTGGCGGTATATATCGCTACACCCCCCCCCTGGGCg > 2:3658054/1‑76 (MQ=255) cGCACTGGCTGGCG‑CCCG‑CCC‑CCATTCGTCTGCACGCGTGCTGGGGGTATATATTGCTACACCACCCCCGTGGGCg > 2:4454638/1‑76 (MQ=255) tGGCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCa > 1:2620211/1‑76 (MQ=255) ggCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt > 2:2572206/1‑76 (MQ=255) ggCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCAt > 2:3589421/1‑76 (MQ=255) ctgccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATaa > 1:4092071/1‑76 (MQ=255) cgcccg‑cctg‑ccg‑cAATTCGTCGGCACGCGTGCTGGCGTTATATATCGCTACGCCACGCCAGTGGGCGACGCATaa < 2:2067915/73‑1 (MQ=255) tgccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 1:6057206/1‑76 (MQ=255) tgccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 1:4626971/1‑76 (MQ=255) tgccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 2:1934252/1‑76 (MQ=255) tgccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 2:4885697/1‑76 (MQ=255) tgccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAAc > 1:2655632/1‑76 (MQ=255) gccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 1:1070014/1‑76 (MQ=255) gccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACa > 1:7129302/1‑76 (MQ=255) ccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt > 1:1466490/1‑76 (MQ=255) ccg‑cctg‑ccg‑cAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAt > 2:5631628/1‑76 (MQ=255) | GAGTTTTATGATTACCCATCTGGTCTGGTTT‑CGCCAGGATTTACGTCTGCACGATAATCTCGCACTGGCTGCCG‑CCTG‑CCG‑CAATTCGTCTGCACGCGTGCTGGCGTTATATATCGCTACACCACGCCAGTGGGCGACGCATAACAT > REL606/720393‑720539 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |