Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,565,689 | Δ1 bp | coding (45/2400 nt) | dosP ← | oxygen‑sensing c‑di‑GMP phosphodiesterase DosP |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,565,686 | 0 | G | . | 96.5% | 107.9 / ‑1.2 | 30 | coding (48/2400 nt) | dosP | oxygen‑sensing c‑di‑GMP phosphodiesterase DosP |
Reads supporting (aligned to +/- strand): ref base G (1/0); major base . (15/13); minor base T (0/1); total (16/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATA > NC_000913/1565567‑1565814 | ttATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGTAAAACATCACTTCATCATTTTTGTTAATTAACACCGCACCCATCATATTTTTCTCCAGGGCGGGaaaaaaaagcc > 2:131296/1‑130 (MQ=255) aTTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCAt > 2:53023/1‑133 (MQ=255) gCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATcgg < 1:120950/133‑1 (MQ=255) gCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATcgg > 1:531877/1‑133 (MQ=255) ccAATGACTTCTTCACGTTTGTATACCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATcggc < 2:105821/133‑1 (MQ=255) cttcACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGGTGAAAAACATCCCTTCATCATTTTCGTTAATTAACACCGCACCCCTCATATTTTGCTCAAGAGC‑GAGAACAAAATGCCCTCGGCGGCATTATcc > 2:232436/1‑133 (MQ=255) ttGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGt < 2:531877/133‑1 (MQ=255) ttGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTCCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGt < 1:375372/133‑1 (MQ=255) tATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAAAATCACTTCATCATTTTTGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAg > 2:413632/1‑133 (MQ=255) tCCCCAGATCTTCTCTGCGGCGGGGTTGAATAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCt < 1:245147/133‑1 (MQ=255) agagCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATg > 1:438719/1‑133 (MQ=255) ctctGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACct > 1:595950/1‑133 (MQ=255) ctctGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAAGGC‑GGGAAAAAAATGCCATCGGCGGGCTTATCCCCATCGGTTAGCTTCATGATTACct > 2:60023/1‑133 (MQ=255) ctctGCGGCGGGGGTGAAAAACATCACTTCCTCATTTTCGTTAATTAACACCGCACCCACCATATTTTGCTCAAAGGC‑GGGAAAAAAATGCACTCGGCGGAATTCTCCACCTCGGTTCGCTTTATTATTACct > 2:243165/1‑133 (MQ=255) ctGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACctct < 1:60023/133‑1 (MQ=255) ttGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGgcg < 2:438719/133‑1 (MQ=255) catTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTc > 1:118238/1‑133 (MQ=255) aattaaCACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGAGCATCTAGAGA‑GGGTTTTCCATAGTTCAac > 1:200246/1‑133 (MQ=255) taaCACCGCACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAacac < 1:53023/133‑1 (MQ=255) cACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTa > 1:127000/1‑133 (MQ=255) cACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTa > 2:462792/1‑133 (MQ=255) cACCCATCATATTTTGCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTACCATAGTTCAACAAGGTTTCTa > 1:246953/1‑133 (MQ=255) atTTTGCTCAAGGGCTGGG‑AAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt < 1:232436/133‑1 (MQ=255) atTTTGCTCAAGGGC‑AGGAAAACAATGCCACCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCttt < 1:183951/133‑1 (MQ=255) tttGCTCAAGGGC‑GGAAAAAAAATGCCATCTGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTg < 1:243165/133‑1 (MQ=255) gCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg > 2:147700/1‑133 (MQ=255) gCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg < 2:612442/133‑1 (MQ=255) gCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCg > 2:73437/1‑133 (MQ=255) gCTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGCCTTTCCAGTGTTGAACGCGGTTTCTGCCTCGTCTTTTGGCg > 2:183951/1‑133 (MQ=255) cTCAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGa > 1:22243/1‑133 (MQ=255) cAAGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGata < 1:413632/133‑1 (MQ=255) aGGGC‑GGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGAtata < 2:22243/133‑1 (MQ=255) | TTATTGCCAATGACTTCTTCACGTTTGTATCCCCAGAGCTTCTCTGCGGCGGGGTTGAAAAACATCACTTCATCATTTTCGTTAATTAACACCGCACCCATCATATTTTGCTCAAGGGCGGGGAAAAAAATGCCATCGGCGGCATTATCCGCATCGGTTAGCTTCATGATTACCTCTGCATCCTGGCGCATCTAAAGACTGGCTTTCCAGAGTTCAACACGGTTTCTACCTCGTCTTTTGGCGATATA > NC_000913/1565567‑1565814 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |