| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,695,249 | G→A | 100% | intergenic (‑432/+20) | iscR ← / ← yfhQ | DNA‑binding transcriptional activator/predicted methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,695,249 | 0 | G | A | 75.7% | 35.1 / 11.3 | 41 | intergenic (‑432/+20) | iscR/yfhQ | DNA‑binding transcriptional activator/predicted methyltransferase |
| Reads supporting (aligned to +/- strand): ref base G (9/1); new base A (0/31); total (9/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.85e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTC > minE/1695202‑1695301 | ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:853504/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:408647/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:340974/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:3388693/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:3032490/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:3031566/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:1929240/1‑69 (MQ=31)ccGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGtt > 1:2802147/1‑69 (MQ=31) tacaGGCCTGATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTGTTACCTTTa < 1:615990/65‑1 (MQ=21) tGTAGGCTTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTa < 1:1032084/68‑1 (MQ=255) caGGCCTGATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTGTTACCTTTa < 1:2833809/65‑1 (MQ=21) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:317334/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2653540/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2654506/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2690918/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2715147/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2543605/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:3259026/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:3310691/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:3316248/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:930065/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:160275/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1138923/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1230982/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1326143/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1408128/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1429944/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1477899/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1586761/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2592519/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1678907/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1753942/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1826847/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:1873190/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2048805/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2103978/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2285692/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2394728/50‑1 (MQ=17) gATAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:2542761/50‑1 (MQ=17) aTAAGACGCGTTAAGCATCGCATCAGGCAACCACTGTTATTCGGCCTTg < 1:475813/49‑1 (MQ=17) aGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc > 1:1275469/1‑66 (MQ=255) | CCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTC > minE/1695202‑1695301 |
| 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: ‑ |