| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,843,423 | C→A | 9.2% | D79Y (GAT→TAT) | cysC ← | adenosine 5'‑phosphosulfate kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,843,423 | 0 | C | A | 9.2% | 226.6 / 8.8 | 97 | D79Y (GAT→TAT) | cysC | adenosine 5'‑phosphosulfate kinase |
| Reads supporting (aligned to +/- strand): ref base C (54/34); new base A (5/4); total (59/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.34e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
CGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACG > minE/1843360‑1843482 | cGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTaa > 1:1284678/1‑69 (MQ=255) aCCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCgag < 1:702799/68‑1 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATAGCTAAAACCGAGATCGc > 1:55242/1‑68 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATAGCTAAAACCGAGATCGc > 1:2117159/1‑68 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATAGCTAAAACCGAGATCGc > 1:2321828/1‑68 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATAGCTAAAACCGAGATCGc > 1:3366866/1‑68 (MQ=255) cAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATAGCTAAAACCGAGATCGc > 1:1203070/1‑68 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:1756292/1‑59 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:1642213/1‑59 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:2397121/1‑59 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:13524/1‑59 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:1306606/1‑59 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:1216892/1‑59 (MQ=255) aaaTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAAcc > 1:544286/1‑59 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:2759643/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:2850381/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:2458946/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:2451985/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:2872612/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:357911/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:998608/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:2672406/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGc > 1:1496419/1‑68 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCa < 1:3070891/69‑1 (MQ=255) aTTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCa < 1:2054324/69‑1 (MQ=255) tCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTg > 1:2289635/1‑68 (MQ=255) tCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTg > 1:1827107/1‑68 (MQ=255) tCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTg > 1:2703553/1‑68 (MQ=255) tCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTg > 1:1301106/1‑68 (MQ=255) tCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTg > 1:790520/1‑68 (MQ=255) cACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTgg > 1:1389887/1‑68 (MQ=255) cACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTgg > 1:3397100/1‑68 (MQ=255) cACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTg > 1:109883/1‑67 (MQ=255) cgacgcgacgGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGaa < 1:534161/69‑1 (MQ=255) cgacgcgacgGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGaa < 1:581767/69‑1 (MQ=255) cgacgcgacgGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGaa < 1:3122296/69‑1 (MQ=255) gacgcgacgGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAAc < 1:718726/69‑1 (MQ=255) atatTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCa < 1:465036/69‑1 (MQ=255) tatTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCa < 1:2912737/68‑1 (MQ=255) ttCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:3171765/1‑69 (MQ=255) ttCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:3023898/1‑69 (MQ=255) ttCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:1172698/1‑69 (MQ=255) ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:1466778/1‑67 (MQ=255) ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:1819672/1‑67 (MQ=255) ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:940077/1‑67 (MQ=255) ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:939623/1‑67 (MQ=255) ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:1013476/1‑67 (MQ=255) ctctTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATcc > 1:3165308/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAAGATTGTCGCCATCcagcag > 1:1233050/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:2054300/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:1068444/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:891727/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:803352/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:3307411/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:2720017/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:679153/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:2586149/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:363635/1‑67 (MQ=255) tACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCcagcag > 1:575406/1‑67 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGc > 1:1175164/1‑44 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:1665423/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:3257025/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:3417557/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:1847666/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:190461/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:1248488/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:2430586/1‑63 (MQ=255) aCGGTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACGTTGTCGCCATCcag > 1:2721515/1‑63 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:1347066/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:1710476/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:967882/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:812223/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2033866/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2425636/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:438450/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2598113/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:3043549/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:1198057/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2902602/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:3361416/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2886382/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2724468/68‑1 (MQ=255) cGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2842834/68‑1 (MQ=255) gATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:783636/67‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATa < 1:2740558/66‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:2156962/68‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:3105650/68‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:1900620/68‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:2224913/68‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:2233108/68‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:110771/68‑1 (MQ=255) aTCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:2573162/68‑1 (MQ=255) aTCGGCATAGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:2027576/68‑1 (MQ=255) aTCGGCATAGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:28902/68‑1 (MQ=255) aTCGGCATAGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:2904845/68‑1 (MQ=255) aTCGGCATAGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:1322744/68‑1 (MQ=255) tCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACg < 1:1951617/67‑1 (MQ=255) | CGGCTTCAACCATCAAATTCGCCACTTCACCGACGCGACGGATATTCTCTTTACGATCGGCATCGCTAAAACCGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACG > minE/1843360‑1843482 |
| 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: ‑ |