| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,606,216 | T→C | 33.3% | A269A (GCT→GCC) | amiA → | N‑acetylmuramoyl‑l‑alanine amidase I |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,606,216 | 0 | T | C | 33.3% | 51.7 / 35.6 | 69 | A269A (GCT→GCC) | amiA | N‑acetylmuramoyl‑l‑alanine amidase I |
| Reads supporting (aligned to +/- strand): ref base T (13/33); new base C (10/13); total (23/46) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.80e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01 | |||||||||||
GAAGAAGAACGGCTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGT > minE/1606157‑1606281 | gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:1097043/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:2126809/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:280359/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:2135385/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:3426666/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:1610934/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:424868/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:683693/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:1462051/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:1375169/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:2220520/69‑1 (MQ=255)gaagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:767103/69‑1 (MQ=255) aagaagagCGGCTGTTAGGCACTGCAGCGTTTCGCCAGAAAATCGCCACAGCGATTGCCGAAGGCGTg < 1:843079/68‑1 (MQ=255) cGGCTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTAt < 1:1651215/69‑1 (MQ=255) cttttAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:2887113/66‑1 (MQ=255) cTGTTAGGCACTGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:1206383/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:2161391/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:2362804/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:1952160/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:1934844/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:2471217/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:1873731/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:2718280/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:2943919/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:1455800/69‑1 (MQ=255) cTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTc < 1:724536/69‑1 (MQ=255) aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTAtt > 1:2602653/1‑62 (MQ=255) aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTAtt > 1:2335241/1‑62 (MQ=255) aGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTAtt > 1:1739961/1‑62 (MQ=255) cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATGAGTTATTTCCACTGGt > 1:1262747/1‑68 (MQ=255) cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGt > 1:1120475/1‑68 (MQ=255) cACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGt > 1:3153275/1‑68 (MQ=255) aCGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTc > 1:364425/1‑69 (MQ=255) aCGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTc > 1:1076923/1‑69 (MQ=255) gcggcgTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGAc > 1:3104958/1‑69 (MQ=255) tCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGaa > 1:3051865/1‑69 (MQ=255) tCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGaa > 1:605899/1‑69 (MQ=255) tCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGaa > 1:441381/1‑69 (MQ=255) tCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGaa > 1:2818231/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAg > 1:2743721/1‑36 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAg > 1:2739644/1‑36 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:78765/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:2315583/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:2438242/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:537681/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:1614838/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:3150000/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcaca > 1:3052060/1‑69 (MQ=255) gAAAATCGCCACAGCGATTGCCGAAGGCGTGATCAGTTATTTCCACTGGTTTGACAACCAGAAAGcac > 1:1525185/1‑68 (MQ=255) aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGcaca < 1:779856/68‑1 (MQ=255) aaaaTCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGcaca < 1:1736028/68‑1 (MQ=255) cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:313405/69‑1 (MQ=255) cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:2802082/69‑1 (MQ=255) cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:3359583/69‑1 (MQ=255) cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:2582861/69‑1 (MQ=255) cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:1910305/69‑1 (MQ=255) cGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:1638799/69‑1 (MQ=255) cGATTGCTGAACGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:3076937/69‑1 (MQ=255) gATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGAt < 1:3187376/68‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:944122/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:740420/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:1366595/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:708536/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:147396/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:574779/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:2974145/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:1796640/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:2498404/69‑1 (MQ=255) tGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGt < 1:235721/69‑1 (MQ=255) | GAAGAAGAACGGCTGTTAGGCACGGCGGCGTTTCGTCAGAAAATCGCCACAGCGATTGCTGAAGGCGTGATCAGTTATTTCCACTGGTTCGACAACCAGAAAGCACATTCGAAAAAGCGATAAGT > minE/1606157‑1606281 |
| 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: ‑ |