| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,016,202:1 | +T | 100% | intergenic (+19/‑121) | metC → / → yghB | cystathionine beta‑lyase, PLP‑dependent/conserved inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,016,202 | 1 | . | T | 75.8% | 81.6 / 25.0 | 33 | intergenic (+19/‑121) | metC/yghB | cystathionine beta‑lyase, PLP‑dependent/conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base . (7/1); new base T (25/0); total (32/1) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
GGACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAG > minE/2016156‑2016255 | ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:357709/1‑66 (MQ=255)ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:479575/1‑66 (MQ=255)ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:1210698/1‑66 (MQ=255)ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:1276695/1‑66 (MQ=255)ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:540335/1‑66 (MQ=255)ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:717894/1‑66 (MQ=255)ggACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATtg > 1:190252/1‑66 (MQ=255) taACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCg < 1:1136969/67‑1 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAAttgtt > 1:620571/1‑52 (MQ=39) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCaa > 1:754364/1‑66 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:864890/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:695764/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:505514/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:498470/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:79304/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:415135/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1057651/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:229783/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:225697/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:2195593/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:2194244/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:2186466/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:2137405/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:211851/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1964842/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1937419/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1716300/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1690914/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1652287/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1602369/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1191708/1‑67 (MQ=255) cACTTTTGGACAATTTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:1185854/1‑67 (MQ=255) cACTTTTGGACAATTTTTGAAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAg > 1:836861/1‑67 (MQ=255) | GGACGCCGGTTTTGCGCGAATTGTATAACATTGCCACTTTTGGACAA‑TTTTGCAGACATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAG > minE/2016156‑2016255 |
| 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: ‑ |