Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,142,011 | T→C | 26.2% | F152L (TTC→CTC) | katE → | hydroperoxidase HPII(III) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,142,011 | 0 | T | C | 26.2% | 65.0 / 27.7 | 42 | F152L (TTC→CTC) | katE | hydroperoxidase HPII(III) |
Reads supporting (aligned to +/- strand): ref base T (13/18); new base C (5/6); total (18/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01 |
GCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGG > minE/1141954‑1142078 | gCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATcc > 1:932440/1‑71 (MQ=255) gCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATcc < 1:50362/71‑1 (MQ=255) ccGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg < 1:2204566/66‑1 (MQ=255) cTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCtctc < 1:659687/61‑1 (MQ=255) aCGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACaaa < 1:1104871/71‑1 (MQ=255) aCGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACaa < 1:1188207/70‑1 (MQ=255) cGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACaaaa > 1:1439121/1‑71 (MQ=255) ggTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt < 1:2027904/71‑1 (MQ=255) aTTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTc > 1:662444/1‑47 (MQ=255) tCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAccc > 1:2143475/1‑69 (MQ=255) cAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAg < 1:2127698/70‑1 (MQ=255) gCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt < 1:1633224/60‑1 (MQ=255) tataAAAGCTTAAGCGATATTACCAAAGCTGATTTCCTCTCAg < 1:398269/43‑1 (MQ=255) tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct < 1:514017/40‑1 (MQ=255) tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct < 1:942588/40‑1 (MQ=255) tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg < 1:776631/43‑1 (MQ=255) tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg < 1:148479/43‑1 (MQ=255) tataAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAg < 1:1778255/43‑1 (MQ=255) taAAAGCTTAAGCGATATTACCAAAGCGGATTTCctct > 1:99136/1‑38 (MQ=255) taAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCAc > 1:1924583/1‑57 (MQ=255) aaaaGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAAc > 1:1349565/1‑48 (MQ=255) aaGCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGtt > 1:1468295/1‑71 (MQ=255) aGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAg < 1:234172/58‑1 (MQ=255) gCTTAAGCGATATTACCAAAGCGGATCTCCTCTCAGATCCGAACa > 1:555137/1‑45 (MQ=255) ttAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt < 1:2086282/47‑1 (MQ=255) gCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTAcc > 1:2087631/1‑71 (MQ=255) cGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCg < 1:358086/71‑1 (MQ=255) atTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGtt > 1:21753/1‑69 (MQ=255) ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt < 1:1111339/37‑1 (MQ=255) ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt < 1:411177/37‑1 (MQ=255) ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAAt < 1:327390/37‑1 (MQ=255) ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg > 1:144684/1‑71 (MQ=255) ttACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAg > 1:2030461/1‑71 (MQ=255) ttACCAAAGCGGATCTCCTCTCAGATCCGAACAAAATCA‑‑CCAGTATTTGTACGTTTCTCTACCGTTCAggg > 1:422037/1‑71 (MQ=255) ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa > 1:64402/1‑36 (MQ=255) ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCa > 1:36124/1‑36 (MQ=255) ccAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTg < 1:115072/48‑1 (MQ=255) cAAAGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtg < 1:1917510/71‑1 (MQ=255) aaGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGggtggt > 1:1815519/1‑71 (MQ=255) aGCGGATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGgtggtg > 1:378540/1‑71 (MQ=255) gATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTgg < 1:1590141/71‑1 (MQ=255) gATCTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTgg < 1:1572166/71‑1 (MQ=255) tttCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTtctc > 1:621141/1‑48 (MQ=255) | GCCGCTCACGGTTATTTCCAGCCATATAAAAGCTTAAGCGATATTACCAAAGCGGATTTCCTCTCAGATCCGAACAAAATCACCCCAGTATTTGTACGTTTCTCTACCGTTCAGGGTGGTGCTGG > minE/1141954‑1142078 |
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: ‑ |