| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,820,792 | G→C | 51.2% | A55A (GCG→GCC) | nrfE → | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,820,792 | 0 | G | C | 51.2% | ‑0.8 / 97.7 | 80 | A55A (GCG→GCC) | nrfE | heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE |
| Reads supporting (aligned to +/- strand): ref base G (10/29); new base C (22/19); total (32/48) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.29e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01 | |||||||||||
ATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGG > minE/2820752‑2820857 | aTTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTggg < 1:803165/69‑1 (MQ=255)aTTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTggg < 1:319786/69‑1 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTggg > 1:1556816/1‑68 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTggg < 1:3092764/68‑1 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa > 1:125451/1‑69 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa > 1:2399471/1‑69 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa > 1:124258/1‑69 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa > 1:2747446/1‑69 (MQ=255) ttGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGa > 1:220929/1‑69 (MQ=255) tctcATCAGCGATTTCTCGGTCTTTTATGTCGCCCAACATAGCTATAGCctg < 1:1382022/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:1765338/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:1849740/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:1872256/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:1876127/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:244720/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:3394812/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:3296499/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:3411256/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:2494389/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:2096305/52‑1 (MQ=255) tctcATCAGCGATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:1171941/52‑1 (MQ=255) aGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa > 1:192125/1‑69 (MQ=255) aGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa > 1:2843610/1‑69 (MQ=255) aGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCa > 1:640294/1‑69 (MQ=255) gATTTCTCTGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:1617446/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCg > 1:3321801/1‑69 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTc > 1:2455973/1‑57 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:2418535/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:2442189/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:1850271/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:3089418/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:254084/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:2879573/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:283406/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:2771468/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:3190511/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:3221534/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:2074033/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:1914076/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:1671521/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:3347204/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:427817/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:1310477/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:440769/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:1216747/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:690134/1‑58 (MQ=255) gATTTCTCGGTCATTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCa > 1:850551/1‑58 (MQ=255) aTTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCgg < 1:1609248/69‑1 (MQ=255) tttCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCgg < 1:1465115/68‑1 (MQ=255) tttCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCgg < 1:2830726/68‑1 (MQ=255) tttCTCGGTCATTTATGTCGCCCAACATAGCTATAGCctg < 1:1065550/40‑1 (MQ=255) aTTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaaa < 1:1115644/48‑2 (MQ=255) aTTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaaa < 1:1204134/48‑2 (MQ=255) aTTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaaa < 1:610081/48‑2 (MQ=255) aTTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaaa < 1:2553356/48‑2 (MQ=255) aTTTATGTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaaa < 1:1825445/48‑2 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:969818/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:1230293/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:594857/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:1259104/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:1500209/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:1531418/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:1705143/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:325651/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:2134168/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:231553/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:2111882/68‑1 (MQ=255) tGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:2321918/68‑1 (MQ=255) gTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATg < 1:1670406/67‑1 (MQ=255) gTCGCCCAACATAGCTATAGCCTGCTGTCGTGGGAACTCaaa < 1:285152/42‑2 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:2724541/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:2716593/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:266393/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:675308/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:287457/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:298773/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:302938/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:1620210/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:3097799/68‑1 (MQ=255) tCGCGCAACATAGCTATAGCCCGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGaa < 1:747233/68‑1 (MQ=255) gcgcAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAgg < 1:1047710/68‑1 (MQ=255) gcgaAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAgg < 1:2346153/68‑1 (MQ=255) | ATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGG > minE/2820752‑2820857 |
| 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: ‑ |