| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,863,770 | A→G | 27.2% | Q31Q (CAA→CAG) | ygcP → | predicted anti‑terminator regulatory protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,863,770 | 0 | A | G | 27.2% | 55.8 / 20.7 | 55 | Q31Q (CAA→CAG) | ygcP | predicted anti‑terminator regulatory protein |
| Reads supporting (aligned to +/- strand): ref base A (2/38); new base G (8/7); total (10/45) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.79e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCTGCCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAG > minE/1863717‑1863830 | gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:2559151/69‑1 (MQ=255)gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:785049/69‑1 (MQ=255)gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:3055751/69‑1 (MQ=255)gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:2496776/69‑1 (MQ=255)gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg > 1:673765/1‑69 (MQ=255)gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:3256080/69‑1 (MQ=255)gctgcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:1598110/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:1113089/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:2851111/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:317133/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:3173831/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:2112353/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:3266747/69‑1 (MQ=255) gcCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGACTCTGAATGCCAGTTTATTTCCGTGTTGTAc < 1:1687723/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:2848306/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:28469/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:1566723/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:1093495/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:872135/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:1231591/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:2269696/69‑1 (MQ=255) cTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGt < 1:308717/69‑1 (MQ=255) gCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGtgttg < 1:1795959/43‑1 (MQ=255) tGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATc > 1:3117238/1‑68 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:402625/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:701878/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:1487874/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:1514307/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:371351/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:38516/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:2902794/69‑1 (MQ=255) cGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:3019867/69‑1 (MQ=255) gATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:2585612/68‑1 (MQ=255) gATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTaaaaa < 1:129679/68‑1 (MQ=255) gATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAg < 1:222461/69‑1 (MQ=255) gATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAg < 1:1983357/69‑1 (MQ=255) gATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAg < 1:1716844/69‑1 (MQ=255) aTTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAg < 1:2160353/68‑1 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAg > 1:3324657/1‑66 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:2500392/1‑67 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:2428772/1‑67 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:2950155/1‑67 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:3335347/1‑67 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:1175066/1‑67 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:395232/1‑67 (MQ=255) tCTGAATGCCAGTTTATTTCCGTGTTGTACGGCAATATCTGCACCATTAGTAATATCGTTAAAAAGa > 1:675337/1‑67 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:906213/66‑1 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:1116039/66‑1 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:1175146/66‑1 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:425790/66‑1 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:1320626/66‑1 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:1886467/66‑1 (MQ=255) tGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:2272457/66‑1 (MQ=255) tGCCAATTTATTTCCGTGCTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:346104/66‑1 (MQ=255) gCCAATTTATTTCCGTGTTGTACGTCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAg < 1:1845722/65‑1 (MQ=255) | GCTGCCGTTAAAGATAATGCCAGCCTGCAACTGGCAATCGATTCTGAATGCCAATTTATTTCCGTGTTGTACGGCAATATCTGCACCATCAGTAATATCGTTAAAAAGATTAAG > minE/1863717‑1863830 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |