Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,244,132 | A→C | 100% | A180A (GCT→GCG) | hemE ← | uroporphyrinogen decarboxylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,244,132 | 0 | A | C | 80.6% | 94.8 / 11.6 | 36 | A180A (GCT→GCG) | hemE | uroporphyrinogen decarboxylase |
Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/29); total (6/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.59e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TGAGCACCGGCTTTAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCGGCATACATCATTTTTTTGATCACGGTGAAG > minE/2244067‑2244183 | tGAGCACCGGCTTTAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCagag < 1:1088013/67‑1 (MQ=255) tAATCTGCGCATTCAGATACAACGTGACGCGTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:805448/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2065273/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:6425/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:470870/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:265784/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:250742/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:645838/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2164805/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2160741/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2130334/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2107875/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2101727/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2092161/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2073282/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2070598/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:2043029/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1987702/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1866633/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1852038/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:181459/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:653433/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1571927/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:707252/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1414869/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1001581/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGGGCa < 1:1663221/60‑1 (MQ=255) tAATCTGCGCATTCAGATACAAAGCGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1645659/60‑1 (MQ=255) aaTCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:360350/59‑1 (MQ=255) cTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa < 1:1072331/56‑1 (MQ=255) ggtACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg > 1:52537/3‑60 (MQ=255) gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg > 1:2215497/1‑60 (MQ=255) gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg > 1:145673/1‑60 (MQ=255) gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg > 1:1422637/1‑60 (MQ=255) gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg > 1:903754/1‑60 (MQ=255) gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGAAGCGCCTGCGGATCgg > 1:1681633/1‑60 (MQ=255) cAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCGGCATACATCATTTTTTTGATCACGGTGAAg < 1:624526/67‑1 (MQ=255) | TGAGCACCGGCTTTAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCGGCATACATCATTTTTTTGATCACGGTGAAG > minE/2244067‑2244183 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |