Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 4,586,035 | +G | intergenic (+118/+97) | PP_4061 → / ← PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,586,033 | 1 | . | G | 85.7% | 74.6 / 9.0 | 28 | intergenic (+116/+99) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (1/3); new base G (7/17); total (8/20) | |||||||||||
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.27e-01 |
GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑CGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGC > NC_002947/4585887‑4586115 | gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gg < 5:122680/149‑1 (MQ=255) cgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gggt < 2:365213/149‑3 (MQ=255) aGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCCAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gggtaa < 7:193186/149‑5 (MQ=255) ccAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCTAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAAc < 4:218252/149‑1 (MQ=255) ccAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAAc < 5:14148/149‑1 (MQ=255) cAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAAc < 3:108168/148‑1 (MQ=255) gAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGc < 2:27670/149‑1 (MQ=255) aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCACTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCCCTcca > 3:332242/1‑147 (MQ=255) gTGAACGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTTGCAATGTTAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑GGTTAAACCCGCTCCTACCAg < 8:82584/149‑1 (MQ=255) cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGA‑ccgt > 2:117322/1‑146 (MQ=39) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCg < 3:329877/149‑1 (MQ=39) aCCACGGCACCGGCGTGAACGGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTg > 3:297874/1‑148 (MQ=38) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAAt < 3:267345/149‑1 (MQ=39) gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGc < 1:117322/149‑1 (MQ=39) gCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGGGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGc < 8:64901/149‑1 (MQ=38) gTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGccc < 4:124133/149‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTg < 2:112474/149‑1 (MQ=39) gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGAG‑GG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGa < 4:105436/125‑1 (MQ=38) gATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGAAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGAG‑GG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGa > 3:105436/1‑125 (MQ=38) tGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTAcc < 2:285241/101‑1 (MQ=38) tGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTAcc > 1:285241/1‑101 (MQ=38) ctCCATTGATAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAg < 2:229073/135‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAg > 1:229073/1‑136 (MQ=38) cATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgcg < 4:297874/149‑1 (MQ=38) aacggTTGCGCATCACC‑GGGGTCTTGGCACGGGTAATATTTTGCATTGATGGGGCCGGCCTCTTCGGG‑GG‑TAACCCCGCTCCTCCCTGGCACCGTGCCGG‑CCTTTAATGCGGCGCTGCCCCTGGTTGTAGCGGGTTTACCCGCGAaga < 4:332242/145‑1 (MQ=255) aaaGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga < 6:206151/148‑1 (MQ=38) ggtgagCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggc < 6:183771/130‑1 (MQ=37) gCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGGCCCTTGAATGCGGCGCTGCCCCTGGTAGGa < 1:18278/104‑1 (MQ=25) gCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGa > 2:18278/1‑103 (MQ=25) ttGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCt > 7:137434/1‑82 (MQ=21) ttGCATTGATCGGGCCGGCCTC‑TTCGCGGG‑TAAACCCGCTCCTACCAGGGACCGTGCCGG‑CCTTGAATGCGGCGCTGCCCCt < 8:137434/82‑1 (MQ=21) | GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACC‑GGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑CGTGCCGG‑CCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGC > NC_002947/4585887‑4586115 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |