Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,688,990 | 0 | A | C | 42.4% | 40.0 / 12.6 | 33 | T92P (ACC→CCC) | nth | DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III) |
Reads supporting (aligned to +/- strand): ref base A (16/3); new base C (0/14); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.77e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG > REL606/1688918‑1689063 | ggCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt > 1:241457/1‑76 (MQ=255) ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg > 1:109011/1‑76 (MQ=255) ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg > 2:1209532/1‑76 (MQ=255) tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt < 1:1002568/76‑1 (MQ=255) aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt < 1:917685/76‑1 (MQ=255) aGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATc > 1:610605/1‑76 (MQ=255) gggTAAAANCGTATCTCAANACGATTGGGCTTTATAACAGCANAGCAGAAAATATCATCAAACCCTGCCGTATCtt < 2:1165848/76‑1 (MQ=255) gTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGc > 2:1080853/1‑76 (MQ=255) taaaaaCCTATAACAAAACGATTGGGCTTTATAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCt < 2:250665/74‑1 (MQ=255) gAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTg < 2:731753/76‑1 (MQ=255) aaCCTATATCAAAACGATGGGGTTTTATAAAAAAAAAAAAAAAAATTTCACCAAACCCTGCCGTATCTTGCTGGag < 1:688399/76‑1 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAATGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc < 2:39529/76‑1 (MQ=255) ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 2:967646/76‑1 (MQ=255) ccTATATAAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 1:1270331/76‑1 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt > 2:3112/1‑76 (MQ=255) aaaaaCGATGGGGTTTTAAAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAAt < 1:1068848/75‑1 (MQ=255) aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATNTCATCAAAACGTGCCCTATCTTGCTGGAGCAGCANAATg > 1:162713/1‑76 (MQ=255) aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATg > 2:752790/1‑76 (MQ=255) tttGGGTTTTATAACAGAAAAGAAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg < 2:858391/75‑1 (MQ=255) tGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt < 2:1247707/76‑1 (MQ=255) tGGGCTCTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt < 1:1382691/76‑1 (MQ=255) gggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTc > 1:317901/1‑76 (MQ=255) gggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTCATGGAGCAGCATAATGGCGAGGTTc > 1:481529/1‑76 (MQ=255) ggCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc < 2:1140230/76‑1 (MQ=255) tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa > 1:927120/1‑76 (MQ=255) tataACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg < 1:604557/76‑1 (MQ=255) tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg > 2:131347/1‑76 (MQ=255) tacaaCAGCAAAGNAAAAAATATCTCCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg < 1:180939/73‑1 (MQ=255) taACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt > 2:1222411/1‑76 (MQ=255) aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc > 2:915576/1‑76 (MQ=255) aaaGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctg < 2:982670/76‑1 (MQ=255) gAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGNTGCGCTTg > 1:1233750/1‑76 (MQ=255) aaaaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCNGAAGATCGTGCTGCGCTTg < 1:1305182/75‑1 (MQ=255) aaaaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg < 2:969853/75‑1 (MQ=255) aaaaaTATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTg < 1:556762/75‑1 (MQ=255) aaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGaa < 1:826462/76‑1 (MQ=255) atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc > 2:660086/1‑76 (MQ=255) catcaCACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg < 2:416194/76‑1 (MQ=255) aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg < 2:241457/76‑1 (MQ=255) aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg > 1:211359/1‑76 (MQ=255) aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCgg > 2:1295747/1‑76 (MQ=255) | GGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGG > REL606/1688918‑1689063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |