Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,558,439 | A→G | Y518H (TAT→CAT) | dosP ← | oxygen‑sensing cyclic‑di‑GMP phosphodiesterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,558,439 | 0 | A | G | 93.3% | 88.2 / ‑1.0 | 30 | Y518H (TAT→CAT) | dosP | oxygen‑sensing cyclic‑di‑GMP phosphodiesterase |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (11/17); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAA > NZ_CP009273/1558353‑1558519 | tGCGCCAAAAACCAAACGCTCTTTTCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc < 2:715501/90‑1 (MQ=255) ggcgcATAAAACCAACCTTTTTTTTCCCTTTTCTTTCATCGCGGGGCTGAAGAACGGCCAGCCGTTACCGCCATTCTTGCGAATATGATc < 1:662666/89‑1 (MQ=255) gcgcCTAAAACCAAACGCTCTTTTTCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc > 1:313272/1‑90 (MQ=255) gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATcc > 1:250818/1‑90 (MQ=255) gcCTAAACCAAACCGCTTTTTTCCCATTTCGTCCACCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcat < 2:655104/90‑1 (MQ=255) gcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcat > 2:940822/1‑90 (MQ=255) gcCAAAAACCAACCGCTTTTTTCCCTTTTCGTTCATCCCGGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcat < 1:810575/90‑1 (MQ=255) acccAACCGCTTTTTTCCCTTTTCTTTCATCGCGGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcatt < 1:427197/88‑1 (MQ=255) cccAACCGTTTTTTTCCCTTTTCTTTCCCCGCGGGGCTGAAAACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg < 1:911660/89‑1 (MQ=255) aCCAACCGCTTTTTTCCCTTTTCTTTCACCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCcattgcattg < 1:704115/90‑1 (MQ=255) ccGCTTTTTTCCCTTTTCTTTCACCGCGGGGTTGAAGACCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc < 1:940822/89‑1 (MQ=255) ccGCTTTTTTCCCTTTCCTTTCACGGCGGGGTTAAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc < 2:554618/89‑1 (MQ=255) tttCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGgag < 2:941367/90‑1 (MQ=255) aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcatt > 1:578319/1‑66 (MQ=255) aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcatt < 2:578319/66‑1 (MQ=255) aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg > 1:607771/1‑90 (MQ=255) aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg > 1:939906/1‑90 (MQ=255) tCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt > 1:392525/1‑90 (MQ=255) tggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTa > 1:78635/2‑90 (MQ=255) ttaaagacCCGCCGGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc < 1:813087/87‑1 (MQ=255) cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc < 2:351761/90‑1 (MQ=255) gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc > 1:276198/1‑42 (MQ=255) gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc < 2:276198/42‑1 (MQ=255) gCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGAt < 1:447923/90‑1 (MQ=255) gCCATTCTTGCGAACATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACt < 1:911325/90‑1 (MQ=255) ttCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCa < 1:79014/64‑1 (MQ=255) ttCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCa > 2:79014/1‑64 (MQ=255) ttCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAg > 1:759097/1‑90 (MQ=255) tGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTaa > 2:579742/1‑90 (MQ=255) tGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTaa < 2:617345/90‑1 (MQ=255) | TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAA > NZ_CP009273/1558353‑1558519 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |