| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,688,990 | 0 | A | C | 37.5% | 45.0 / 12.3 | 32 | T92P (ACC→CCC) | nth | DNA glycosylase and apyrimidinic (AP) lyase (endonuclease III) |
| Reads supporting (aligned to +/- strand): ref base A (14/6); new base C (0/12); total (14/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.12e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCG > REL606/1688915‑1689065 | cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa < 2:307809/76‑1 (MQ=255)cTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa < 2:643688/76‑1 (MQ=255) ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCg > 2:69130/1‑76 (MQ=255) gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt < 1:852975/76‑1 (MQ=255) gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt < 2:995932/76‑1 (MQ=255) gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt < 2:249088/76‑1 (MQ=255) gggTGAAAACCTATATCAAAACGATGGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCtt < 2:226359/76‑1 (MQ=255) aacccTATATAAAAACGATGGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGa < 2:997053/73‑1 (MQ=255) aCCTATATCAAAACGAGTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc < 2:682037/76‑1 (MQ=255) aCCTATAACAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc < 2:682048/76‑1 (MQ=255) ccTATATCAAAACGATTGGTTTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 2:661433/76‑1 (MQ=255) ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 1:690910/76‑1 (MQ=255) ccTATATAAAAACGATTGGGTTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGagca < 2:856869/76‑1 (MQ=255) tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagcagc > 2:778727/1‑76 (MQ=255) atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa > 1:323104/1‑76 (MQ=255) atCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa < 2:749757/76‑1 (MQ=255) atCAAAACGATTGGGCTTTATAACAGAAAAGAAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATa < 1:126893/76‑1 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt > 1:418909/1‑76 (MQ=255) aCGTTTGGGTTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCg < 2:409953/76‑1 (MQ=255) cGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa > 2:1212234/1‑76 (MQ=255) gATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAg > 1:16244/1‑76 (MQ=255) aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg > 1:441689/1‑76 (MQ=255) ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgat > 1:153438/1‑74 (MQ=255) ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt > 1:412691/1‑76 (MQ=255) tGGGTTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt < 1:12583/76‑1 (MQ=255) tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt > 1:787444/1‑76 (MQ=255) tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt > 2:908447/1‑76 (MQ=255) tGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt < 2:148305/76‑1 (MQ=255) ggCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc > 1:1153414/1‑76 (MQ=255) cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCgg > 1:1168416/1‑76 (MQ=255) tttATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa < 2:1106631/76‑1 (MQ=255) tttATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa < 2:924018/76‑1 (MQ=255) ttAAAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGaa < 1:755149/76‑1 (MQ=255) tttaACAGCAAAGCAGAAATTTTTTTCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg < 1:123872/74‑1 (MQ=255) tataACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg < 2:964726/76‑1 (MQ=255) taACAGCAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt < 2:1024973/76‑1 (MQ=255) cAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct > 1:226329/1‑76 (MQ=255) gCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTgcgc > 1:906619/1‑76 (MQ=255) aaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAg < 1:249539/76‑1 (MQ=255) tcttcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt < 1:342678/73‑1 (MQ=255) aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg < 2:104896/76‑1 (MQ=255) aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCg < 1:907867/76‑1 (MQ=255) aCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCg > 1:195217/1‑76 (MQ=255) | CTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCG > REL606/1688915‑1689065 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |