Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 4,518,633 | C→T | 79.3% | G296D (GGC→GAC) | clpA ← | ATP‑dependent serine protease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 4,518,633 | 0 | C | T | 79.3% | 91.8 / 17.5 | 43 | G296D (GGC→GAC) | clpA | ATP‑dependent serine protease |
Reads supporting (aligned to +/- strand): ref base C (2/7); new base T (15/19); total (17/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.81e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01 |
CGAAACGCCCTTTCAGGCCGCGCAGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCACGCCTG > AE015451/4518412‑4518861 | cGAAACGCCCTTTCAGGCCGCGCAGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGa > 1:265857/1‑236 (MQ=255) cTTTCAGGCCGCGCAGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAAtt < 2:265856/230‑1 (MQ=255) tCAGGCCGCGCAGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGa < 2:593232/234‑1 (MQ=255) aGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGAt < 2:731467/233‑1 (MQ=255) aGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGc > 1:673800/1‑236 (MQ=255) gCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGag > 1:163144/1‑235 (MQ=255) gCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCgca > 2:642251/1‑213 (MQ=255) gCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGACTCCATCACCCCGCCCGACGCCGCACCGGCGTCgat > 2:531789/1‑175 (MQ=255) aaGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTtctc > 2:690874/1‑233 (MQ=255) gACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAg < 2:332219/221‑1 (MQ=255) gcgcCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacc < 1:390393/235‑1 (MQ=255) cgcCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacc < 2:454503/234‑1 (MQ=255) cgcCAGGGCACGATCCTTCCCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacc < 1:227866/234‑1 (MQ=255) gcCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAgg > 1:104480/1‑171 (MQ=255) tcGAAGATGCCGCGAAACTCCGGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTtcg < 2:208043/141‑1 (MQ=255) gAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGc < 1:53703/235‑1 (MQ=255) aTGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGgc < 1:214810/236‑1 (MQ=255) tGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTtcg > 1:742905/1‑134 (MQ=255) tGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGgcg > 1:165893/1‑236 (MQ=255) gcgAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGa < 2:742904/134‑1 (MQ=255) gcgAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCa < 1:13080/236‑1 (MQ=255) aaaCTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacc < 1:93892/201‑1 (MQ=255) tGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacc < 1:690875/194‑1 (MQ=255) aCGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGa < 1:314557/224‑1 (MQ=255) tAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTcacc < 2:558290/187‑1 (MQ=255) tCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAg < 1:751155/223‑1 (MQ=255) aCCCGACGACAGCGGCGGCTGGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGt < 2:38193/174‑1 (MQ=255) agcGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGctct > 2:695156/1‑205 (MQ=255) tataagagacagCATCCATCACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTg < 1:125662/225‑1 (MQ=255) atccatcACCCCGCCCGGCGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGa < 2:660567/160‑1 (MQ=255) atcACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAgg > 2:329233/1‑189 (MQ=255) atcACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAgg < 1:329234/189‑1 (MQ=255) atcACCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTGGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGAt > 2:325601/1‑212 (MQ=255) tcACCCCGCCCGACGCCGCACCGTCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGc > 1:227768/1‑220 (MQ=255) aCCCCGCCCGACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGa < 2:671045/232‑1 (MQ=255) cccgACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAg > 2:125662/1‑146 (MQ=255) cgACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGt > 1:748169/1‑205 (MQ=255) gACGCCGCACCGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGcc > 1:435223/1‑236 (MQ=255) gACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGa > 2:587244/1‑40 (MQ=255) gACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGa < 1:587245/40‑1 (MQ=255) cgccgcACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGa < 2:559435/101‑1 (MQ=255) agcGGCGTCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCACgcctg < 1:562637/234‑1 (MQ=255) cGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAgg < 2:246470/172‑1 (MQ=255) | CGAAACGCCCTTTCAGGCCGCGCAGAATGCCCACGGTGTCTTCGACCGATGGCTCGCTGACATCCACCTTCTGGAAGCGACGCGCCAGGGCACGATCCTTCTCGAAGATGCCGCGAAACTCCTGGAACGTAGTCGAACCGATGCAGCGGATTTCACCCGACGACAGCAGCGGCTTGAGCAGGTTGGAGGCATCCATCACCCCGCCCGACGCCGCACCGGCGCCGATGATGGTGTGAATTTCGTCGATGAACAGGATGGCCTGCGGACGCTTGCGCAGCTCACCGAGCAGCGCCTTGAAGCGCTTCTCGAAGTCACCACGGTATTTGGTACCGGCCAGCAAGGCGCCCAGATCAAGGGAGTACACCACGCTCTGCGCCAGCAGGTCAGGCACCTGACCATCGACGATGCGTTTGGCCAGGCCTTCGGCGATGGCCGTCTTGCCCACGCCTG > AE015451/4518412‑4518861 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |