Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,378,978 | A→T | 15.4% | intergenic (‑34/+29) | ampC ← / ← frdD | penicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,378,978 | 0 | A | T | 15.4% | 67.4 / 5.7 | 26 | intergenic (‑34/+29) | ampC/frdD | penicillin‑binding protein; beta‑lactamase, intrinsically weak/fumarate reductase (anaerobic), membrane anchor subunit |
Reads supporting (aligned to +/- strand): ref base A (15/7); new base T (2/2); total (17/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.91e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01 |
CTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACG > NC_000913/4378836‑4379123 | ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCgg > 1:543724/1‑148 (MQ=255) tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTAT‑GGCGGGc > 3:824092/1‑149 (MQ=255) ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTAT‑GGCGGGcc < 4:824092/149‑1 (MQ=255) tattcacaATATCGTTGATTTGTTGATGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAAGAAGGCGCAGAGCGTCGTTTTGAACATAGTGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GTCGGGCCCGATTTACATTGGCGATg < 1:492381/144‑1 (MQ=255) gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGc > 3:511278/1‑149 (MQ=255) gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACTGCCCGCCTAT‑GGCGGGcc < 3:904751/117‑1 (MQ=255) gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACTGCCCGCC‑TATGGCGGGcc > 4:904751/1‑117 (MQ=255) gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGtt > 4:886924/1‑139 (MQ=255) gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGtt < 3:886924/139‑1 (MQ=255) ttGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAAc > 2:259030/1‑148 (MQ=255) gggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg > 4:425725/1‑149 (MQ=255) gggCAGCAAATGTAGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGc > 1:685916/1‑148 (MQ=255) ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg < 2:543724/148‑1 (MQ=255) tgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCTTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa < 4:511278/149‑1 (MQ=255) ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATagcag < 1:259030/148‑1 (MQ=255) gcgcAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAc < 2:685916/148‑1 (MQ=255) tGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAg > 4:3204/1‑149 (MQ=255) gTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGAt > 3:610387/1‑149 (MQ=255) ggCCCGCCTAT‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt > 3:729644/1‑149 (MQ=255) gCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGTCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt > 1:589530/1‑148 (MQ=255) ccGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa < 3:71003/74‑1 (MQ=255) ccGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCa > 4:71003/1‑74 (MQ=255) ccGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa > 4:49964/1‑149 (MQ=255) ccGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa > 4:811996/1‑149 (MQ=255) ccGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa > 3:350544/1‑149 (MQ=255) cATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTAACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACg > 3:750089/1‑149 (MQ=255) | CTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATA‑GGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACG > NC_000913/4378836‑4379123 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |