Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,093,980 | G→A | 20.8% | D369N (GAC→AAC) | ycdT → | diguanylate cyclase, membrane‑anchored |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,093,980 | 0 | G | A | 20.8% | 155.8 / 34.8 | 77 | D369N (GAC→AAC) | ycdT | diguanylate cyclase, membrane‑anchored |
Reads supporting (aligned to +/- strand): ref base G (24/37); new base A (8/8); total (32/45) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.70e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01 |
CAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAA > NC_000913/1093881‑1094077 | cAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGTGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACgg < 1:1128481/100‑1 (MQ=255) cAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACgg < 1:1006057/100‑1 (MQ=255) cAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACgg < 2:573511/100‑1 (MQ=255) ggTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGAc < 1:2335455/100‑1 (MQ=255) gTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACa < 2:2169839/100‑1 (MQ=255) gATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGTGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATc < 2:1784753/100‑1 (MQ=255) gATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATc < 1:2023026/100‑1 (MQ=255) aTAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCg > 2:2503301/1‑100 (MQ=255) tAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGa > 1:390209/1‑100 (MQ=255) tAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGa > 2:1188709/1‑100 (MQ=255) aaaaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGACGAGTTTGGCGTCTTGCTGACGGACATCGAt > 1:2344264/1‑100 (MQ=255) aaaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGCAGAGTTTGGCGTCTTGCTGACGGACATCGATa > 2:858463/1‑100 (MQ=255) aaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATAc < 2:1984550/100‑1 (MQ=255) aaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATAc < 1:1983011/100‑1 (MQ=255) aaaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATAc < 1:191754/100‑1 (MQ=255) aaCAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACt < 1:1108829/100‑1 (MQ=255) aCAGTGGTGAACATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTg < 1:1904340/100‑1 (MQ=255) aGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGaa < 2:2207362/100‑1 (MQ=255) ggtgAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGt > 2:1422098/1‑100 (MQ=255) gAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCg < 2:1920349/100‑1 (MQ=255) gAATATCATTCGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGTGTGTTGCTGACGGACATCGATACTGAACGTGCg < 1:1700952/100‑1 (MQ=255) cATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCt < 1:2409030/100‑1 (MQ=255) aTTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTGGGCGGCGAAGAGTATGGCGTCTTGCTGACGGACATCGATACTTACCGTGCGAAAGCtt > 1:1262976/1‑100 (MQ=255) aTTGGGAAAAGCAGACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGc < 1:2147105/98‑1 (MQ=255) ttGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCttt < 2:1610021/100‑1 (MQ=255) tGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTa < 2:2165504/100‑1 (MQ=255) tGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTa < 1:1163190/100‑1 (MQ=255) ggAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGAGAAAGCTTTAGc < 2:1705257/100‑1 (MQ=255) gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCg > 2:1705986/1‑100 (MQ=255) gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCg > 1:1727678/1‑100 (MQ=255) gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCg > 1:640625/1‑100 (MQ=255) gAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCg > 1:1018770/1‑100 (MQ=255) aaaaGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCgg < 2:769254/100‑1 (MQ=255) aaGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGaa < 1:1018447/100‑1 (MQ=255) cATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAgg > 2:944493/1‑100 (MQ=255) aTACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGa > 1:937874/1‑100 (MQ=255) tACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGAt > 2:138294/1‑100 (MQ=255) cGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTc > 2:2576345/1‑100 (MQ=255) cGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTc > 2:978094/1‑100 (MQ=255) gACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGTAAGCTTTAGCGGAAAGGATTCg < 1:1505793/100‑1 (MQ=255) gACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCg > 1:751107/1‑100 (MQ=255) aCCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCgg > 2:2038610/1‑100 (MQ=255) cAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGa < 2:896586/100‑1 (MQ=255) gATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGaaa > 2:1279471/1‑100 (MQ=255) tGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAAt < 1:1498952/100‑1 (MQ=255) gACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATg > 2:768327/1‑100 (MQ=255) ccTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGtt > 2:1039948/1‑100 (MQ=255) cTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTg > 1:2575269/1‑100 (MQ=255) aGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGt > 1:1713006/1‑100 (MQ=255) gcgcgcGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGtt > 2:1697956/1‑100 (MQ=255) cgcgcgACGGCGGCGAAGAGTTTGGCGTATTGCTGAAGGACATCGATACTGAACGTGCGAAAGCTTCAGCGGAAAGGATGCGGGAAAAGGTTGAGCGttt > 1:459657/1‑100 (MQ=255) cggcggCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGc < 2:1240158/100‑1 (MQ=255) cggcggCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGc > 2:184905/1‑100 (MQ=255) ggcggcGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCg > 1:2141814/1‑100 (MQ=255) cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt < 2:2470385/100‑1 (MQ=255) cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt < 1:1188705/100‑1 (MQ=255) cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt < 1:2149078/100‑1 (MQ=255) cggcgAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt < 1:657632/100‑1 (MQ=255) cggcgAAGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGAt < 2:443975/100‑1 (MQ=255) gAAGAGTTTGGGGTCTTGCTGACGGACATCGATATTGAACGTGCGAATGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAAGTGGGGATAATc < 1:1362255/100‑1 (MQ=255) aaGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATcc < 2:1980980/100‑1 (MQ=255) aaGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATcc < 1:1352557/100‑1 (MQ=255) aaGAGTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATcc < 2:617535/100‑1 (MQ=255) agagTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCt < 2:1836376/100‑1 (MQ=255) agagTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCt > 1:2590025/1‑100 (MQ=255) gagTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTg < 2:1227334/100‑1 (MQ=255) gagTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTg < 1:1521745/100‑1 (MQ=255) gagTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTg < 2:2407638/100‑1 (MQ=255) agTTTGGCGTCTTGCTGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGa < 1:1260534/100‑1 (MQ=255) gTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGaa < 2:1944657/100‑1 (MQ=255) tttGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAAt < 2:2557079/100‑1 (MQ=255) aggtgttGCGGACGGACGCGGGTGCTGGACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGc < 1:2270095/95‑1 (MQ=255) ggcTTGCTGGCGGACATCGATAGTGAAGGTGCGAATGCTTTAGTTGAAAGGATTCGTCAAAATTTTGAGCGTTTAACTGTCGATAATCCTGAATATGc < 1:2082859/96‑1 (MQ=255) cTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCtata > 1:795462/1‑100 (MQ=255) cTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCtata > 2:1204589/1‑100 (MQ=255) cagACGGACATCGATACTGAACGTGCGAAAGATTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACcac > 2:1745706/3‑100 (MQ=255) cTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACcac < 1:952305/100‑1 (MQ=255) tGACGAACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACcaca > 1:1743741/1‑100 (MQ=255) gACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACaa < 1:2183165/100‑1 (MQ=255) cGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACaaaa < 2:1378263/100‑1 (MQ=255) | CAGGTGATAAAAACAGTGGTGAATATCATTGGGAAAAGCATACGACCAGATGACCTTTTAGCGCGCGTCGGCGGCGAAGAGTTTGGCGTCTTGCTGACGGACATCGATACTGAACGTGCGAAAGCTTTAGCGGAAAGGATTCGGGAAAATGTTGAGCGTTTAACTGGCGATAATCCTGAATACGCTATACCACAAAA > NC_000913/1093881‑1094077 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |