| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 471,243 | 0 | G | C | 20.0% | 108.5 / 19.5 | 45 | A203A (GCG→GCC) | sucD | succinyl‑CoA synthetase, NAD(P)‑binding, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base G (22/14); new base C (4/5); total (26/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.61e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAAGAGCACG > minE/471182‑471310 | tCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGt < 1:629255/67‑1 (MQ=255) gggCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCg > 1:2846628/1‑71 (MQ=255) tAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTga < 1:718979/69‑1 (MQ=255) aCTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCg > 1:1799506/1‑63 (MQ=255) aCTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGt > 1:731676/1‑65 (MQ=255) aCTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGt > 1:3130304/1‑65 (MQ=255) tCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATcggcgg < 1:890972/70‑1 (MQ=255) tctcGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCt < 1:1446633/70‑1 (MQ=255) gAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTga < 1:3029599/68‑1 (MQ=255) aaTGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTgaagaa < 1:2205313/70‑1 (MQ=255) gTTCGAAAAAGATCCGCAGACCGAAGCGATCGtgat > 1:2697846/1‑36 (MQ=255) cGAAAAAGATCCGCAGACAGAAGCGATCGTGATGATCGGTgag > 1:1531318/1‑43 (MQ=255) aaaaaGATCCGCAGACCGAAGCGATCGTGATGATCgg > 1:1040141/1‑37 (MQ=255) aaaaaGATCCGCAGACCGAAGCGATCGTGATGATCGGt > 1:198761/1‑38 (MQ=255) aaaaaGATCCGCAGACCGAAGCGATCGTGATGATCGGt > 1:1683648/1‑38 (MQ=255) aaaaGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTg < 1:1102142/71‑1 (MQ=255) aaaaGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTg < 1:3092842/71‑1 (MQ=255) aaaGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAagca > 1:1972715/1‑66 (MQ=255) aGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCg > 1:1454291/1‑70 (MQ=255) aGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCg > 1:1539881/1‑70 (MQ=255) gATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTa < 1:1964901/71‑1 (MQ=255) gATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTa < 1:1476066/71‑1 (MQ=255) tCCGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaa < 1:2693177/55‑1 (MQ=19) tCCGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaa < 1:1754086/55‑1 (MQ=19) tCCGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaa < 1:309939/55‑1 (MQ=19) tCCGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaa < 1:1939974/55‑1 (MQ=19) ccGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACag > 1:11616/1‑70 (MQ=255) ccGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaa < 1:2610377/54‑1 (MQ=16) cGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaaga > 1:2733009/1‑55 (MQ=19) cGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaaga > 1:378391/1‑55 (MQ=19) cGCAGACCGAAGCCATCGTCATGATCGGCGAGATCGGTGGTTCCGCCgaagaaga > 1:2435297/1‑55 (MQ=19) cGCAGACCGAAGCCATCGTCATGAT‑GGCGAGATCGGTGGTTCCGCCgaagaaga > 1:1623428/1‑54 (MQ=12) gCAGACCGAAGCGATCGTGATGATCGGTGAGATcgg < 1:1605067/36‑1 (MQ=255) cAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAagca > 1:1344181/1‑57 (MQ=255) aGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCa > 1:1187160/1‑69 (MQ=255) gACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAag < 1:1825280/71‑1 (MQ=255) cGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTg < 1:2032037/41‑1 (MQ=255) aaGCGATCTTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAAGAGCAc > 1:1001196/1‑71 (MQ=255) aaGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAAGAGCAc < 1:804325/71‑1 (MQ=255) aaGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAAGAGCAc > 1:1868225/1‑71 (MQ=255) aGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCt > 1:2661015/1‑37 (MQ=255) aGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCt > 1:241005/1‑37 (MQ=255) aGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCt > 1:2284670/1‑37 (MQ=255) aGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCt > 1:1943752/1‑37 (MQ=255) aGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAAGAGCACg > 1:995906/1‑71 (MQ=255) | TCCCGGGCTCTAACTTTATCGACATTCTCGAAATGTTCGAAAAAGATCCGCAGACCGAAGCGATCGTGATGATCGGTGAGATCGGCGGTAGCGCTGAAGAAGAAGCAGCTGCGTACATCAAAGAGCACG > minE/471182‑471310 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |