Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,584,941 | G→T | 28.9% | intergenic (+45/‑90) | benC → / → benD | benzoate 1,2‑dioxygenase electron transfer component/1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,584,941 | 0 | G | T | 28.9% | 26.3 / 4.2 | 17 | intergenic (+45/‑90) | benC/benD | benzoate 1,2‑dioxygenase electron transfer component/1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base G (4/8); new base T (2/3); total (6/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGTGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTT > NC_002947/3584793‑3585057 | tgtgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGa < 2:479685/150‑1 (MQ=255) cGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGccc < 2:173200/150‑1 (MQ=255) cGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc < 2:70233/150‑1 (MQ=255) ccACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCca < 1:463791/150‑1 (MQ=255) cAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTc < 1:122633/150‑1 (MQ=255) cAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTc < 2:119059/150‑1 (MQ=255) aGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAagaga < 2:418844/150‑1 (MQ=255) aGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGtt < 1:55810/150‑1 (MQ=255) ttACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGcc > 1:353112/1‑150 (MQ=255) cGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAggg > 2:233590/1‑150 (MQ=255) ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc > 1:503242/1‑150 (MQ=255) cccTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaagg < 2:80174/150‑1 (MQ=255) tCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGAGTTACATATGAACAACAGATTCCAAGGCAAGGtt > 1:158474/1‑150 (MQ=255) ggAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGaa < 1:482504/85‑1 (MQ=255) ggAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGaa > 2:482504/1‑85 (MQ=255) gAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAacaaca > 2:429577/1‑108 (MQ=255) gAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAacaaca < 1:429577/108‑1 (MQ=255) | TGTGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTT > NC_002947/3584793‑3585057 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |