Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,965,106 | A→G | 7.5% | K440E (AAG→GAG) | yegE → | phosphodiesterase domain‑containing diguanylate cyclase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,965,106 | 0 | A | G | 7.5% | 148.9 / 2.6 | 53 | K440E (AAG→GAG) | yegE | phosphodiesterase domain‑containing diguanylate cyclase |
Reads supporting (aligned to +/- strand): ref base A (26/23); new base G (2/2); total (28/25) | |||||||||||
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.95e-01 |
GGCAAGCAACCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTGCCC > NC_002947/1964972‑1965251 | ggCAAGCAACCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGac > 8:58325/1‑149 (MQ=255) gcaagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACGTGGCGCGCCTGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGaca < 4:287694/149‑1 (MQ=255) caagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGTCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCTACGacac < 5:246926/149‑1 (MQ=255) aagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcga < 4:325492/142‑3 (MQ=255) aagcaaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcga > 3:325492/1‑140 (MQ=255) aaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCtgtg < 5:11982/149‑1 (MQ=255) aaCCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGGCAACGACACCtgtg > 6:267278/1‑149 (MQ=255) caggttGCGCTTACCTGGTGGGCCAACCGCCGGTATTTGTAATACCGCTCAAGTAGGAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGc < 1:111800/144‑1 (MQ=255) gATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcgag > 2:238984/1‑129 (MQ=255) gATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGCTTGAACTGGTcgag < 1:238984/132‑4 (MQ=255) cGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGAt > 4:69498/1‑148 (MQ=255) gggCTGGCCAACCGCCGCGAGTTCGAATACGGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGa < 2:304064/149‑1 (MQ=255) gggCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGa < 3:69498/149‑1 (MQ=255) gggCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGa < 7:206543/149‑1 (MQ=255) gCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGt < 7:58325/149‑1 (MQ=255) gCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTtgc < 1:273943/148‑1 (MQ=255) ccAACCGGCGCGAGTTCGTTTACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTGGACCTGTTGAAGCTGGTCAACGACACCTGTTGCCACGCTACCGGCGATGAGTtgctg < 4:99563/149‑1 (MQ=255) ccgccgCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCgg > 4:178409/1‑148 (MQ=255) cgAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCtgtg > 2:273850/1‑113 (MQ=255) cgAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCtgtg < 1:273850/113‑1 (MQ=255) aGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTgc < 3:178409/149‑1 (MQ=255) aTACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGtgc < 3:149270/149‑1 (MQ=255) ccGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTgcgc > 7:323622/1‑141 (MQ=255) ccGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTgcgc < 8:323622/141‑1 (MQ=255) ccGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGtgctg > 6:304460/1‑148 (MQ=255) gAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCgg > 7:181261/1‑149 (MQ=255) ggCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGc < 8:181261/149‑1 (MQ=255) cgcTCACCGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTg < 2:284299/71‑1 (MQ=255) cTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTTCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATAc > 6:227844/1‑149 (MQ=255) cTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACAt < 6:58768/110‑1 (MQ=255) cTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACAt > 5:58768/1‑110 (MQ=255) aGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCgg > 3:214901/1‑149 (MQ=255) cGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTgg > 2:177998/1‑149 (MQ=255) ggcggcATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGggc > 1:120461/1‑149 (MQ=255) aTTCGCTGATGTTCCTCGACCTCGACCAGTTAAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATATGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGAt > 3:133779/1‑149 (MQ=255) cGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGAt < 5:255793/68‑1 (MQ=255) cGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGAt > 6:255793/1‑68 (MQ=255) gCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATg > 4:213605/1‑146 (MQ=255) gCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATg < 3:213605/146‑1 (MQ=255) aTGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTg < 2:120461/148‑1 (MQ=255) gTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCg > 1:129848/1‑149 (MQ=255) ttCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCg < 5:227844/148‑1 (MQ=255) cctcgacctcgaCCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACTTCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGtgc > 4:9102/1‑149 (MQ=255) ctcgacctcgacCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGtgct > 3:335599/1‑149 (MQ=255) tcgacctcgaccAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGtgctg > 3:79073/1‑149 (MQ=255) gacctcgaccAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTgctgct > 6:83523/1‑149 (MQ=255) ctcgcccAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTgg > 3:305924/1‑148 (MQ=255) cAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTg > 6:41281/1‑149 (MQ=255) cAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTg > 6:36115/1‑149 (MQ=255) aGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTgc < 4:79073/149‑1 (MQ=255) ttCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGc > 2:218275/1‑129 (MQ=255) ttCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGc < 1:218275/129‑1 (MQ=255) ttCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTgccc > 3:215123/1‑149 (MQ=255) | GGCAAGCAACCCACGATGCGCTGACCGGGCTGGCCAACCGCCGCGAGTTCGAATACCGC‑C‑TGGAGCAGGCGCTCAACGACCTGGCGCGCCAGGCCGGGCGGCATTCGCTGATGTTCCTCGACCTCGACCAGTTCAAGCTGGTCAACGACACCTGTGGCCACGCTGCCGGCGATGAGTTGCTGCGGCACATCTGCGCGGTGCTGCAATCCGGCCTGCGTGAGGGCGATACGCTGGCGCGGCTGGGCGGCGATGAATTTGGCGTGCTGCTGGAAAACTGCCC > NC_002947/1964972‑1965251 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |