| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,407,343 | T→A | 33.4% | G151G (GGA→GGT) | ydaM ← | diguanylate cyclase, csgD regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,407,343 | 0 | T | A | 33.4% | 90.1 / 75.3 | 78 | G151G (GGA→GGT) | ydaM | diguanylate cyclase, csgD regulator |
| Reads supporting (aligned to +/- strand): ref base T (28/24); new base A (14/12); total (42/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01 | |||||||||||
CCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGA > NC_000913/1407247‑1407442 | cccGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCa > 1:1627175/1‑100 (MQ=255)cccGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCa > 2:1852886/1‑100 (MQ=255)cccGAGCATATTTATTTCCCAGGTATGTTTCTGACACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCa < 2:1326091/100‑1 (MQ=255) gAGCATATTTCTTTCCCGGGTATGTTTCTGGCACATCGTTTCATGGTTATAAACATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCAtct < 2:679335/100‑1 (MQ=255) atatTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTg < 1:2063927/100‑1 (MQ=255) atatTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTg < 1:143370/100‑1 (MQ=255) atatTTATTGCCCAGGTATGTTTCTGGCACATCGTTTCATTGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTg > 2:744314/1‑100 (MQ=255) tatttattTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGa < 1:508660/100‑1 (MQ=255) tttatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACg > 1:1137772/1‑100 (MQ=255) tttatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACg < 1:659256/100‑1 (MQ=255) tatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACggg < 2:1994136/100‑1 (MQ=255) tatttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACggg < 1:1994233/100‑1 (MQ=255) tttCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtc > 1:948633/1‑100 (MQ=255) tCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaa > 2:422030/1‑100 (MQ=255) tCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtcaa > 2:607940/1‑100 (MQ=255) tCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGtcaa > 1:2248593/1‑100 (MQ=255) cccAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaat > 2:629357/1‑100 (MQ=255) cccAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaat < 2:1137769/100‑1 (MQ=255) ccAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGtcaatc < 1:238829/100‑1 (MQ=255) ggTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaac > 2:1567004/1‑100 (MQ=255) gTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaaca < 1:1544126/100‑1 (MQ=255) tATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaacaa > 2:983467/1‑100 (MQ=255) gaTGTTTCTGGCACAGCGTTTCATCGTTATAACCATAGCAAATGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATcaaca < 2:1574426/98‑1 (MQ=255) ttCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTg < 1:880165/100‑1 (MQ=255) tGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTg < 1:1829825/100‑1 (MQ=255) ggCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGc < 2:1804155/100‑1 (MQ=255) cacaTCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGAGTCAATCAACAACATTGGTGCag < 1:727875/100‑1 (MQ=255) acaTCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTGCATCTCGTGACGGGTCAATCAACAACATTGGTGCag > 2:1668373/1‑99 (MQ=255) caTCGTTTCATGGTTATAAGCAGAGAAAGGGAGCGCTGCGGGGTTATCATCGACTTGTGGTCCATCTCCTGACGGGTCAAGCAACAACAGAGGTGCagag > 1:661369/1‑100 (MQ=255) tCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGtt < 1:455901/100‑1 (MQ=255) atttctaTGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTa > 1:325011/7‑100 (MQ=255) ttCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTc < 1:983469/100‑1 (MQ=255) ttCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAg > 2:2084564/1‑98 (MQ=255) accATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTc < 1:824628/98‑1 (MQ=255) cATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAg > 2:719355/1‑100 (MQ=255) cATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAg > 2:1865788/1‑100 (MQ=255) aTGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGa < 2:1416086/100‑1 (MQ=255) aTGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGa > 2:127203/1‑100 (MQ=255) tagTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGaa > 1:1524088/3‑100 (MQ=255) ggTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGata < 1:1522207/100‑3 (MQ=255) ttATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGaaaaa < 2:1627164/100‑1 (MQ=255) tataACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCCATCAACAACATTGGTGCAGAGTTAGTCAGaaaaaa > 1:1688136/1‑100 (MQ=255) ataACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGaaaaatg > 2:2078204/1‑98 (MQ=255) ataACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAg < 1:2096049/100‑1 (MQ=255) taACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgc > 1:396918/1‑100 (MQ=255) aaCCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcg > 2:1708292/1‑100 (MQ=255) cATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAgcgcgc > 1:1767729/1‑100 (MQ=255) cATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATAAACAACCTCGATGCAGAGTTAGTCAGAAAACAgcgcgc > 1:1067012/1‑100 (MQ=255) aTAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCa < 1:1612364/100‑1 (MQ=255) aTAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCa > 2:2065995/1‑100 (MQ=255) tAGAAATTGCGCGCGGCGAGGTTAGCATCGACGATTTGTCAATCTCGTCACGGGGCAATCAACAACATTGGTGCAGAGTTAGCCAGAAAAAAGCGCGCa > 2:1351718/1‑99 (MQ=255) aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATa < 1:849953/100‑1 (MQ=255) aGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATa < 2:1507269/100‑1 (MQ=255) gAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATaa > 1:1054724/1‑100 (MQ=255) aaTTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATaaaa < 2:2184312/100‑1 (MQ=255) tGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAAcc < 1:5100/100‑1 (MQ=255) tGAGCGCGGCGAGGTTAACATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAAcc > 2:1827377/1‑100 (MQ=255) aGCGCGGCGAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCtt > 2:4234/1‑100 (MQ=255) gcgAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGttttt > 2:1801346/1‑100 (MQ=255) gcgAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGtttt > 2:1443477/1‑99 (MQ=255) gcgAGGTTAGCATCGAAGATTTAAACAACTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGttttt > 2:955140/1‑100 (MQ=255) gAGGTTAGCATCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTg < 1:4426/100‑1 (MQ=255) aGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTg < 1:753712/99‑1 (MQ=255) ggTTAGCAATGAAGATTTGTCCATCTCGGCACGAGTCAATCAACAACATTGGTGCAGAGTTAGTCATAATTAAGCGCGCATAAAAACCTTg < 2:1754527/91‑1 (MQ=255) ttAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCt > 1:634600/1‑100 (MQ=255) tAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTg < 1:1446737/100‑1 (MQ=255) tAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTg > 1:1297260/1‑100 (MQ=255) gCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAACAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGAt > 1:1672325/1‑100 (MQ=255) gCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGAt > 1:1551604/1‑100 (MQ=255) gCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGAt > 2:1757680/1‑100 (MQ=255) gCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGAt < 1:959130/100‑1 (MQ=255) aTCGACGATTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATaa > 2:1744850/1‑100 (MQ=255) gacgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAAttt > 1:2186635/1‑100 (MQ=255) gacgaTTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGATGATAAtt > 2:1342531/1‑99 (MQ=255) acgaTTTGTCCATCTCGTGACGGGTTAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAATCCTTGTTTTTTGCGCTGATAATTTg < 2:2109250/100‑1 (MQ=255) cgaTTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGc < 1:2243394/100‑1 (MQ=255) aTTTGACCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGcc < 1:607940/99‑1 (MQ=255) gTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCg > 2:160622/1‑100 (MQ=255) tCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGa < 1:1801358/100‑1 (MQ=255) | CCCGAGCATATTTATTTCCCAGGTATGTTTCTGGCACATCGTTTCATGGTTATAACCATAGAAATTGAGCGCGGCGAGGTTAGCATCGACGATTTGTCCATCTCGTGACGGGTCAATCAACAACATTGGTGCAGAGTTAGTCAGAAAAAAGCGCGCATAAAAACCTTGTTTTTTGCGCTGATAATTTGCCGAGCGA > NC_000913/1407247‑1407442 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |