Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 729,603 | A→C | 25.0% | P207P (CCA→CCC) | fabD → | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 729,603 | 0 | A | C | 25.0% | 62.8 / 36.6 | 52 | P207P (CCA→CCC) | fabD | malonyl‑CoA‑[acyl‑carrier‑protein] transacylase |
Reads supporting (aligned to +/- strand): ref base A (19/20); new base C (7/6); total (26/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGGCGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT > minE/729535‑729672 | gcggcgGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACcag > 1:1406939/1‑70 (MQ=255) gcggcgGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACcag > 1:3338842/1‑70 (MQ=255) cAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACaa < 1:3268154/70‑1 (MQ=255) aaaaCGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAAc > 1:2450902/1‑71 (MQ=255) aaCGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTg < 1:2129514/71‑1 (MQ=255) cgcgcgCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTgg > 1:1044125/1‑70 (MQ=255) cgcgcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAg < 1:3306606/71‑1 (MQ=255) tGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACaa < 1:864027/58‑1 (MQ=255) gacGTTACCAGTGAGCGTACCGTGTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAAt > 1:1392662/3‑71 (MQ=255) ccGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAAtt < 1:2746857/71‑1 (MQ=255) ccGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAAtt < 1:2122452/71‑1 (MQ=255) cGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTa < 1:1812753/71‑1 (MQ=255) aCCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCg > 1:1238798/1‑70 (MQ=255) cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa < 1:1333608/71‑1 (MQ=255) cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa < 1:2352206/71‑1 (MQ=255) aGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa > 1:2593168/1‑71 (MQ=255) gTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa < 1:2515254/70‑1 (MQ=255) tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAAtt < 1:2169469/61‑1 (MQ=255) gAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGaa > 1:3487844/1‑58 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATTAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc < 1:3452148/71‑1 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc < 1:3280600/71‑1 (MQ=255) aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc > 1:3185462/1‑71 (MQ=255) cGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa < 1:2991794/64‑1 (MQ=255) cGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCa < 1:1384179/68‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCCGCCGCCGACaaa > 1:2925660/1‑41 (MQ=38) tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGc > 1:39123/1‑58 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa < 1:1199241/62‑1 (MQ=255) tACCGTCTCACTGTGCGCTGATGAAACCAG‑AGCCGACAAACTGGCAGTAGAATTAGCGAAAAt < 1:2549133/63‑1 (MQ=255) tcACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt < 1:787061/64‑1 (MQ=255) tcACTGTGCGCTGATGAAACCAGCAGCCGACAAACTgg < 1:4068101/38‑1 (MQ=255) tcACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAAt > 1:65853/1‑57 (MQ=255) aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa > 1:2497391/1‑64 (MQ=255) aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa > 1:1289678/1‑64 (MQ=255) aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCAcc > 1:3199168/1‑70 (MQ=255) aCTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa > 1:1772401/1‑54 (MQ=255) gtgCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa > 1:3947650/1‑69 (MQ=255) tgCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTa < 1:1417049/37‑1 (MQ=255) tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCa > 1:914748/1‑64 (MQ=255) tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAAc < 1:3863467/69‑1 (MQ=255) tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACa < 1:3969633/70‑1 (MQ=255) tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAg > 1:73272/1‑71 (MQ=255) gcTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa > 1:1555058/1‑65 (MQ=255) cTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTcc > 1:209615/1‑71 (MQ=255) atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtg > 1:2966419/1‑70 (MQ=255) tgaAACCCGCAGCCGACAAACTGGCAGTAGAATTAg > 1:778030/1‑36 (MQ=255) gaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttg > 1:1321206/1‑71 (MQ=255) aaaCCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt < 1:3945791/71‑1 (MQ=255) aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa > 1:2130742/1‑50 (MQ=255) aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt > 1:872711/1‑71 (MQ=255) aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtg < 1:3028040/67‑1 (MQ=255) aaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt < 1:2063715/70‑1 (MQ=255) aaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtggtgt > 1:4127857/1‑70 (MQ=255) aCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGa < 1:919787/71‑1 (MQ=255) ccAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGaa > 1:2105943/1‑71 (MQ=255) cagcagCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAt < 1:2810698/71‑1 (MQ=255) | GCGGCGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT > minE/729535‑729672 |
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: ‑ |