| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,846,899 | 0 | A | . | 48.1% | 90.1 / 102.2 | 54 | coding (654/1038 nt) | iap | aminopeptidase in alkaline phosphatase isozyme conversion |
| Reads supporting (aligned to +/- strand): ref base A (28/0); new base . (0/26); total (28/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.33e-16 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTCAACAGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTGCG > minE/1846884‑1846953 | tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggt < 1:2385091/55‑2 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:560592/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:2996776/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:2708405/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:3060435/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:3075651/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:309818/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:2605252/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:3108310/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:2369416/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:315369/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:3300741/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:3430318/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:2053939/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:567097/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:168753/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:677009/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:1572260/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:150599/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:725225/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:865212/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:1298057/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:1296791/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:1274313/55‑1 (MQ=255)tttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:1192656/55‑1 (MQ=255) ttCAACAGCGGTGT‑AAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAggg < 1:3001041/54‑1 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3473947/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:889055/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:1463819/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:82857/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:799371/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:741856/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:1474591/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:698454/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:698329/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:1681823/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2161915/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:520085/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3593665/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2846626/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3467180/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2221739/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3379177/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2226973/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3209431/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3166065/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:3165024/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:223787/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2447140/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2616852/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2641439/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:2643506/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgcg > 1:1059227/1‑60 (MQ=255) ggTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTgc > 1:1856191/1‑59 (MQ=255) | TTTCAACAGCGGTGTAAAAACCCCTGAGGCAGTAAGGAAATTAACGCGCGACAGGGCGCTGGCAATTGCG > minE/1846884‑1846953 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |