Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 447,412 | G→A | 100% | R102C (CGC→TGC) | kdpD ← | fused sensory histidine kinase in two‑component regulatory system with KdpE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 447,412 | 0 | G | A | 81.0% | 41.7 / 6.2 | 21 | R102C (CGC→TGC) | kdpD | fused sensory histidine kinase in two‑component regulatory system with KdpE |
Reads supporting (aligned to +/- strand): ref base G (1/3); new base A (5/12); total (6/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TGGGATGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTAACGG > minE/447349‑447481 | tGGGATGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCggcgggc < 1:3773906/70‑1 (MQ=255) tGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGgcagcg < 1:3209299/62‑1 (MQ=255) tGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGgcagcg < 1:3010971/62‑1 (MQ=255) tGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGgcagcg < 1:2455633/62‑1 (MQ=255) aaCCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATc > 1:1282074/1‑70 (MQ=255) aCCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAATATCAGCGCCGGGCGGCGGGCGAGggcggc < 1:2030142/66‑1 (MQ=255) aCCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATc < 1:2001751/69‑1 (MQ=255) cTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCga > 1:3545854/1‑69 (MQ=255) gcgcATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCAGGGCGGCGGGCGAGGGCGGCATCGAGATCa < 1:983885/71‑1 (MQ=255) ttACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGggcagc < 1:426626/43‑1 (MQ=255) aCTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGc > 1:1856608/1‑71 (MQ=255) tgCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGa < 1:3288517/70‑1 (MQ=255) ttCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATgccgcc > 1:364018/1‑71 (MQ=255) gTCCATTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCa < 1:3062706/70‑1 (MQ=255) aTTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCa < 1:1808996/45‑1 (MQ=255) aTTAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCa < 1:1804088/45‑1 (MQ=255) tAAGATCAGCGCCGGGCAGCGGGCGAGGGCGTCATc < 1:709487/36‑1 (MQ=37) tAAGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCa < 1:1723331/43‑1 (MQ=255) aaGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATc > 1:4258932/1‑35 (MQ=255) aaGATCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGc > 1:3468738/1‑71 (MQ=255) aTCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTg < 1:2657086/71‑1 (MQ=255) aTCAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTg < 1:1870240/71‑1 (MQ=255) gcAGCGCCGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTg < 1:3066240/69‑1 (MQ=255) cGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTa < 1:1974320/71‑1 (MQ=255) cGGGCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTa < 1:1751355/71‑1 (MQ=255) gggCAGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTa < 1:2469977/70‑1 (MQ=255) ggcgggcgAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTAACgg > 1:1527779/1‑70 (MQ=255) | TGGGATGACGGGAACCTGGCGCATTACTGTGCGCCAGTTCGTCCATTAAGATCAGCGCCGGGCGGCGGGCGAGGGCGGCATCGAGATCAAACTCGCTGATATGCCGCCCACGGTACGCCTGGCGTTTTAACGG > minE/447349‑447481 |
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: ‑ |