| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 875,783 | C→T | 6.3% | S333S (AGC→AGT) | acnA → | aconitate hydratase 1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 875,783 | 0 | C | T | 6.3% | 203.8 / 3.1 | 63 | S333S (AGC→AGT) | acnA | aconitate hydratase 1 |
| Reads supporting (aligned to +/- strand): ref base C (19/40); new base T (2/2); total (21/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01 | |||||||||||
TGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAAC > minE/875714‑875853 | tGGTGCCACCTGTGTCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGc > 1:1805054/1‑70 (MQ=255) ggTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCg < 1:1577475/70‑1 (MQ=255) gCCACCTGGGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGTGAAg < 1:1821976/70‑1 (MQ=255) aCCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCa < 1:1530661/71‑1 (MQ=255) ggCTTCTTCCCAAACGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTc > 1:2702001/1‑69 (MQ=255) gCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGa < 1:2715393/70‑1 (MQ=255) cttcttCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCg < 1:2863363/56‑1 (MQ=255) cttcttCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGt < 1:3978807/71‑1 (MQ=255) tcttcCCAATCGATGCTGTAACCCTCGATTACATGCGGTTAAGCGGGCGCAGCGAAGATCAGGTCGa < 1:811032/67‑1 (MQ=255) tcCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTc > 1:4138695/1‑71 (MQ=255) cccAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCa > 1:3505275/1‑57 (MQ=255) cccAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGt < 1:3457752/69‑1 (MQ=255) cccAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTc > 1:3317530/1‑70 (MQ=255) ccAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGa < 1:2887591/71‑1 (MQ=255) ccAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGa < 1:3059584/71‑1 (MQ=255) cAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCg < 1:2444777/48‑1 (MQ=255) aaTCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAg < 1:2989619/50‑1 (MQ=255) aaTCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGGCGaaa < 1:3335562/71‑1 (MQ=255) tCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaaa < 1:3091739/71‑1 (MQ=255) tCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaaa < 1:2844211/71‑1 (MQ=255) cGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaaa > 1:1300632/1‑70 (MQ=255) gATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGtt > 1:2026809/1‑59 (MQ=255) gATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGt > 1:138628/1‑58 (MQ=255) gATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAata > 1:1879308/1‑71 (MQ=255) gATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAat < 1:3683331/70‑1 (MQ=255) gATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAat < 1:836099/70‑1 (MQ=255) tGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAg < 1:305871/44‑1 (MQ=255) tGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAg < 1:2954953/44‑1 (MQ=255) gCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGTGAAGATCa > 1:2046536/1‑47 (MQ=255) gCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGTGAAGATCAGGTCGAGTTGGTCGAAAAATATGc > 1:1599820/1‑71 (MQ=255) cTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGa > 1:2173286/1‑43 (MQ=255) cTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGa < 1:2679865/52‑1 (MQ=255) cTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaaaa < 1:3854451/65‑1 (MQ=255) tGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGcc > 1:3404139/1‑70 (MQ=255) tAACCCTCGATTACATGCGTTTAAGCGGGCGCAGTGAAGATCAGGTCGAGTTGGTCGaaaaa < 1:1059464/62‑1 (MQ=255) aaCCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTgg < 1:3867430/53‑1 (MQ=255) aaCCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTgg < 1:156021/53‑1 (MQ=255) aCCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAg < 1:3554558/70‑1 (MQ=255) cccTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGAt < 1:1390249/38‑1 (MQ=255) ccTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAgcgc < 1:3756059/71‑1 (MQ=255) cGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCa > 1:1845760/1‑69 (MQ=255) gATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAgg < 1:381018/70‑1 (MQ=255) aTTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaa < 1:2808027/50‑1 (MQ=255) aTTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAgcgc < 1:2927553/66‑1 (MQ=255) tACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAg > 1:245436/1‑39 (MQ=255) tACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTc > 1:4091412/1‑45 (MQ=255) aCATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtg > 1:3658595/1‑71 (MQ=255) cATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtgt < 1:1059875/71‑1 (MQ=255) cATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtg < 1:3089210/70‑1 (MQ=255) aTGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaa < 1:1385947/46‑1 (MQ=255) aTGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGaaa < 1:1745331/46‑1 (MQ=255) aTGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGCTCGAAAAATATGCCAAAGCGCAGGGCAtgtg < 1:3226403/71‑1 (MQ=255) gCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAat > 1:3503046/1‑47 (MQ=255) gTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATAAGCCAAAGCGCAGGGCATGTGGCGt < 1:1178012/71‑1 (MQ=255) tAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAAc > 1:2295356/1‑71 (MQ=255) aaGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAAc < 1:1415985/70‑1 (MQ=255) aGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAAccc > 1:2037745/1‑71 (MQ=255) gCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtg < 1:3475583/59‑1 (MQ=255) gCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtg < 1:3864301/59‑1 (MQ=255) gCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtg < 1:3900011/59‑1 (MQ=255) ggCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGc < 1:2391415/71‑1 (MQ=255) gCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCa < 1:2883487/36‑1 (MQ=255) gCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCa > 1:554015/1‑48 (MQ=255) gCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCAtg > 1:230718/1‑50 (MQ=255) cGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAAc > 1:3472505/1‑71 (MQ=255) cGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAAc > 1:129120/1‑71 (MQ=255) cGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAAc > 1:3484958/1‑71 (MQ=255) | TGGTGCCACCTGTGGCTTCTTCCCAATCGATGCTGTAACCCTCGATTACATGCGTTTAAGCGGGCGCAGCGAAGATCAGGTCGAGTTGGTCGAAAAATATGCCAAAGCGCAGGGCATGTGGCGTAACCCGGGCGATGAAC > minE/875714‑875853 |
| 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: ‑ |