Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,925,446 | A→T | 19.4% | intergenic (‑151/+414) | PP_2563 ← / ← PP_2564 | antibiotic biosynthesis protein/quercetin 2,3‑dioxygenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,925,446 | 0 | A | T | 19.4% | 66.0 / 11.3 | 31 | intergenic (‑151/+414) | PP_2563/PP_2564 | antibiotic biosynthesis protein/quercetin 2,3‑dioxygenase |
Reads supporting (aligned to +/- strand): ref base A (11/14); new base T (3/3); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01 |
TCGAGTGCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAGCCGCCGAA > NC_002947/2925301‑2925595 | tCGAGTGCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCggg > 2:475616/1‑150 (MQ=255) gTGCTCCGCAGGAACTTTGTGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCACAAACTTAATGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAg < 1:195346/150‑1 (MQ=255) gCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCGGGTTAGca < 1:475616/150‑2 (MQ=255) tCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGgt > 1:283121/1‑150 (MQ=255) cGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAAACCTGCTAACCCGACGGGTTAGCggggt < 2:329125/150‑4 (MQ=255) cGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAAACCTGCTAACCCGACGGGTTAGCggggt > 1:329125/1‑147 (MQ=255) gAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCa > 2:158703/1‑150 (MQ=255) gAACGTTGGGCTGGCTGGTAGCACAACGCTAGGGCATGCTGGCTGATGGCTCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCa > 1:13908/1‑150 (MQ=255) ttGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGga < 2:13908/150‑1 (MQ=255) aGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCt > 1:495856/1‑150 (MQ=255) caACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCg < 2:247175/150‑1 (MQ=255) tGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGgttg < 1:242693/113‑1 (MQ=255) tGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGgttg > 2:242693/1‑113 (MQ=255) gCAGGCTGGCTAATCCAGCGAATATTGGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTg > 2:55529/1‑150 (MQ=255) cAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGt < 2:561828/150‑1 (MQ=255) ggctggctGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCGGGTTAGcag > 2:495125/1‑104 (MQ=255) ggctggctGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGTCGGGTTAGcag < 1:495125/106‑3 (MQ=255) gctgATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGcaca < 1:368832/149‑1 (MQ=255) gctgATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGcaca > 2:368832/1‑149 (MQ=255) aTGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCACATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCaaa < 2:25793/150‑1 (MQ=255) gacAATATTTGAACGCTGCCTTGCGATAGCCGGGCAGGTCCGGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGa > 2:24498/3‑150 (MQ=255) atatTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATgcg > 1:157238/1‑150 (MQ=255) aGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGt < 2:157238/150‑1 (MQ=255) cGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATc < 1:362153/150‑1 (MQ=255) gggCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCCCCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGaaa < 1:568107/150‑1 (MQ=255) cAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATgg < 1:476653/55‑1 (MQ=255) cAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATgg > 2:476653/1‑55 (MQ=255) aaGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCCAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTc < 2:398328/150‑1 (MQ=255) ccgCTAACCCGTCGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTgc < 2:454498/58‑1 (MQ=37) ccgCTAACCCGTCGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTgc > 1:454498/3‑60 (MQ=37) aCCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAgccg < 2:495856/150‑1 (MQ=255) aCGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAGCCGCCGaa < 1:158703/150‑1 (MQ=255) | TCGAGTGCTCCGCAGGAACGTTGGGCTGGCTGGTAGCACAACGCTAGTGCAGGCTGGCTGATGGCGCCAATATTTGAAGGCTGCCTTGCGATAGCCGGGCAGGTCCAGGGCCGCAAGCGTGCCCCCATAAAACCCTGCTAACCCGACGGGTTAGCGGTTGTTCCATGGGAGATCACCCAGCCCTTCCGGGCTGCGCTGTCGCACAGCAAACAGGATGCGCATGACGTCAGATCTGCACTGACATCGAGATCTAGAAAAGCATTCGCCGACTGATCCGGGTAGGTCAGCCGCCGAA > NC_002947/2925301‑2925595 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |