Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | CP000730 | 849,783 | (A)6→7 | coding (909/912 nt) coding (1/996 nt) |
yvcJ → yvcK → |
ATP‑binding protein hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 849,777 | 1 | . | A | 89.5% | 52.0 / 1.5 | 19 | coding (903/912 nt) | yvcJ | ATP‑binding protein |
Reads supporting (aligned to +/- strand): ref base . (1/1); new base A (10/7); total (11/8) | |||||||||||
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.52e-01 |
TTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAG‑AAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTAC > CP000730/849640‑849907 | ttAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAG‑aaa < 1:147638/141‑1 (MQ=255) ttAGTAATTGACATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAG‑aaa < 2:353414/141‑1 (MQ=255) aaTTGCTATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGTAGTATTTGAATATAATGTTTATGTGAATCATAGGTATGCACATGTTGAAATTG‑GCGAG‑aaaaaaat > 2:92088/1‑139 (MQ=255) taggTTGTACGGGTGAACAACATAGATATGTAGCATTAGCAGAAAGACTAGGTAATTATCTAGATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAATA‑AAAAGAAAAAAATGAGACaaa < 2:1691/139‑1 (MQ=255) ttGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATaaa > 1:15495/1‑141 (MQ=255) gTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTAc < 2:15495/141‑1 (MQ=255) gTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTAc > 1:178849/1‑141 (MQ=255) tGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACt < 2:217612/141‑1 (MQ=255) gATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGgtggtg > 1:75985/1‑141 (MQ=255) ttAGCAGAACGACGAGGTAATAAACTAAATGAAGTATTTGAATATAATGTTAATGTGCATCATAGGGACGCACATATTGAAAAAAAAAAAA‑AAAAAACGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTa < 2:307222/141‑1 (MQ=255) tATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTa > 2:166423/1‑141 (MQ=255) tAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAagag > 2:63312/1‑141 (MQ=255) tGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTT‑GCGAG‑AAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGCTGGCACTGGCTTCTCAGCTATGCCTACGGGATTAAGAGACTTc > 2:225279/1‑141 (MQ=255) tgtgCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGt < 2:178849/141‑1 (MQ=255) gtgCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTa > 1:21631/1‑141 (MQ=255) tgactcatAGGCACGACCATATTGAAAGGA‑GCTAGAAAAAAATGAGACTAATATAATTTGTACTTATTGGTGGTGGCACTGGCTTATCAGTTATGGCTAAAAGAATAAGAGAATTCCCAAGTCATATTACGGCGATTGTaa < 2:271205/137‑1 (MQ=255) gCATCATAGGGACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAAc < 1:368420/141‑1 (MQ=255) gggACGCACATATTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTg > 2:207318/1‑141 (MQ=255) ggACACACATATTGAAAGTG‑TCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGa < 2:354354/141‑1 (MQ=255) tatTGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAAtggtgg > 2:211374/1‑141 (MQ=255) ttGAAAGTG‑GCGAGAAAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGa > 1:50594/1‑141 (MQ=255) aaGTG‑GCGAGAAAAAAATGAGACAAACAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTAc > 1:356367/1‑141 (MQ=255) | TTAGTAATTGCCATCGGTTGTACGGGTGGACAACATCGATCTGTAGCATTAGCAGAACGACTAGGTAATTATCTAAATGAAGTATTTGAATATAATGTTTATGTGCATCATAGGGACGCACATATTGAAAGTG‑GCGAG‑AAAAAATGAGACAAATAAAAGTTGTACTTATCGGTGGTGGCACTGGCTTATCAGTTATGGCTAGGGGATTAAGAGAATTCCCAATTGATATTACGGCGATTGTAACAGTTGCTGATAATGGTGGGAGTAC > CP000730/849640‑849907 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |