Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,711,796 | 0 | A | C | 32.0% | 47.2 / 14.2 | 25 | T92P (ACC→CCC) | nth | endonuclease III |
Reads supporting (aligned to +/- strand): ref base A (17/0); new base C (0/8); total (17/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.25e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGC > NC_000913/1711732‑1711853 | aGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa > 1:2666380/1‑65 (MQ=255) aGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCaaaa > 1:1837852/1‑65 (MQ=255) ggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt > 1:2467434/1‑65 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTAt > 1:252489/1‑65 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 1:2698897/1‑65 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 1:40262/1‑65 (MQ=255) aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt > 1:2316162/1‑65 (MQ=255) ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATc > 1:340493/1‑65 (MQ=255) ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATc > 1:2862710/1‑65 (MQ=255) atCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCt > 1:1527954/1‑65 (MQ=255) tCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg > 1:1489825/1‑65 (MQ=255) tCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTg > 1:1633924/1‑65 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg > 1:2596923/1‑65 (MQ=255) cAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTgg > 1:951057/1‑65 (MQ=255) aaaaCGATTGGGTTTTAAAACAACAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGa < 1:1587391/65‑1 (MQ=255) aaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagc > 1:481267/1‑65 (MQ=255) ggttGGGGTTTATAACAGCAAAGCAGAAAATATCATCAACCCCTGCCGTATCTTGCTGGagcagc < 1:1442785/63‑1 (MQ=255) gagtGGGTTATATAACACAAAAGCAGAAAATAGCACCAAAACCTGCCGTATCTTGCTGGagcagc < 1:1467529/62‑1 (MQ=255) gATTGGGTTTTAAAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcagc < 1:2034860/65‑1 (MQ=255) gATTGGGTTTTAAAACAGAAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGagcagc < 1:942797/65‑1 (MQ=255) gATTGGGCTTTAAAACAGAAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGagcagc < 1:1183367/65‑1 (MQ=255) gATTGGGCTTTAAAACACAAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGagcagc < 1:2119108/65‑1 (MQ=255) gATTGGCCTTAACAACAACAAGCCGCAAAATTTATTAAAACCCTGCCGTATCTTGCTGGagcagc < 1:273042/65‑1 (MQ=255) ggCTTTATAACAAAAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAAt < 1:913397/65‑1 (MQ=255) ataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGctgt > 1:1293659/1‑62 (MQ=255) taACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg > 1:2220333/1‑65 (MQ=255) taACAACAAAGCAAAAAATATAATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg < 1:653461/65‑1 (MQ=255) aaacaacAAAGGAGAAAATAACATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg < 1:551216/64‑1 (MQ=255) aaaCAGAAAAGAAGAAAATCCCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg < 1:2974767/64‑1 (MQ=255) actagcaaagcaaAAAATAGCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt < 1:2045421/62‑1 (MQ=255) aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt > 1:1039052/1‑65 (MQ=255) tgCAAAGAAAAAAATACCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc < 1:1564874/64‑1 (MQ=255) gcaaaaGCAAAAAATAAAAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTcc < 1:1352009/62‑1 (MQ=255) cAGAAAATAGCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt < 1:1554439/65‑1 (MQ=255) aaaaaaaTACCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt < 1:842505/62‑1 (MQ=255) aGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc > 1:1074215/1‑65 (MQ=255) aaaaTATCAAAAACCCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGt < 1:68730/65‑1 (MQ=255) atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGANGTTCCGGAAGATCGggnt > 1:2321327/1‑65 (MQ=255) atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgcc > 1:2246402/1‑64 (MQ=255) atatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGCTCCGGAAGATCGtgct > 1:529472/1‑65 (MQ=255) catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTgcgc < 1:298796/65‑1 (MQ=255) | AGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGC > NC_000913/1711732‑1711853 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |