| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,845,328 | A→G | 50.5% | L223L (TTG→CTG) | cysD ← | sulfate adenylyltransferase, subunit 2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,845,328 | 0 | A | G | 50.5% | ‑1.0 / 88.1 | 97 | L223L (TTG→CTG) | cysD | sulfate adenylyltransferase, subunit 2 |
| Reads supporting (aligned to +/- strand): ref base A (19/29); new base G (42/7); total (61/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.66e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01 | |||||||||||
ACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAGCGGAACAATGTCGATATTTTC > minE/1845260‑1845394 | aCGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCaa < 1:2691315/69‑1 (MQ=255) cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg > 1:1602721/1‑69 (MQ=255) cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg > 1:2858576/1‑69 (MQ=255) cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg > 1:3313215/1‑69 (MQ=255) cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg > 1:2704760/1‑69 (MQ=255) cGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTcg > 1:897425/1‑69 (MQ=255) ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:2886547/69‑1 (MQ=255) ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:2895425/69‑1 (MQ=255) ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:749442/69‑1 (MQ=255) ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:2088270/69‑1 (MQ=255) ttttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:19629/69‑1 (MQ=255) tttttAATCACTTCGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:3100896/68‑1 (MQ=255) tCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa < 1:858823/68‑1 (MQ=255) tCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa < 1:388012/68‑1 (MQ=255) tCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa < 1:3202641/68‑1 (MQ=255) cACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGa < 1:1289153/67‑1 (MQ=255) cGCCCGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTcgcg < 1:2146372/55‑1 (MQ=255) cccGCTTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:1335153/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCg > 1:854625/1‑64 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGAc > 1:1788984/1‑67 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:2411541/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:2310490/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:2615273/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:2858479/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:553793/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:601418/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:635553/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:833126/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:846047/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:987563/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:990884/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:101849/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:1339143/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:1965156/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:1521937/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:1264233/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:120868/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:220109/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACg > 1:220660/1‑68 (MQ=255) cccGATTGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTACAGAACCGGACg > 1:1519777/1‑68 (MQ=255) ggTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTc > 1:793767/1‑68 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGtt > 1:2774254/3‑52 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTc > 1:2548974/3‑53 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:578148/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:1866992/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:1659684/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:1709409/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:3011691/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:3391723/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:3428392/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:196943/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:1972634/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:2297671/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:990389/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:2072251/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:2630485/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGtt > 1:2414375/3‑67 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGt > 1:3228848/3‑66 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTCGCGTTCCAGAACCGGACGTTc > 1:432895/3‑68 (MQ=255) gattGCAGGTTGATACGGTTGTCATCAATCATCATCAGCATACCGTAGCGTTCCAGAACCGGAc > 1:2370415/3‑64 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:1816165/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:373128/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:447119/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:886641/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:878464/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:864866/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:1408448/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:841110/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:1336236/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:1911265/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:2295882/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:274074/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:2772302/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:2814777/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:1010045/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:294564/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:3404124/69‑1 (MQ=255) ggTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCga < 1:314461/69‑1 (MQ=255) gTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgag < 1:3161041/69‑1 (MQ=255) tCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga < 1:3202175/69‑1 (MQ=255) tCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga < 1:3925/69‑1 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGCACCGGACGTTCCGCAGCgaga > 1:733206/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACg < 1:3393690/55‑1 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACg < 1:2142513/55‑1 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACg < 1:79632/55‑1 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:412582/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:2302711/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:2439147/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:2002357/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:2295556/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:409437/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:2963415/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:3042495/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga > 1:414398/1‑68 (MQ=255) cGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgag > 1:1297572/1‑67 (MQ=255) atcaatcaTCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCgaga < 1:3197766/55‑1 (MQ=255) aatcaTCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAGCGGAACAATGTc > 1:3187834/1‑69 (MQ=255) atcaACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAGCGGAACAATGTCGATAtt > 1:1706139/1‑68 (MQ=255) aaCATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAGCGGAACAATGTCGATATTTTc > 1:3384954/1‑68 (MQ=255) | ACGCACCATCCGTTTTTTAATCACTTCGCCCGGTTGCAGGTCGATACGGTTGTCATCAATCATCATCAACATACCGTCGCGTTCCAGAACCGGACGTTCCGCAGCGAGATATAGCGGAACAATGTCGATATTTTC > minE/1845260‑1845394 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |