Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,056:1 | +C | 38.2% | intergenic (+139/+76) | 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,056 | 1 | . | C | 38.2% | 37.2 / 54.3 | 34 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (8/13); new base C (2/11); total (10/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 |
GACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑G‑GTAAACCCGCTCCTACCAGGGA‑C‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACC > NC_002947/4585923‑4586195 | cccgatgaccgtgaccggcaaggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCtt < 4:168988‑M1/36‑1 (MQ=255) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATgc‑gg < 1:82555/149‑1 (MQ=39) gaccggcaaggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTg < 5:378340‑M1/49‑1 (MQ=255) caaggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTg > 6:378340‑M1/95‑149 (MQ=255) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GTCGCCGGCCTTGAATGC‑GGCGCTGcccc < 1:3512/148‑1 (MQ=38) aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTgg < 4:127632/149‑1 (MQ=39) ggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTa > 7:137716‑M1/92‑149 (MQ=255) ggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTa > 3:242865‑M1/92‑149 (MQ=255) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTa < 3:235960/149‑1 (MQ=39) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAg > 8:155275/1‑149 (MQ=39) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAg < 7:155275/149‑1 (MQ=39) tgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAgg > 3:168988‑M1/90‑149 (MQ=255) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCggg < 4:121726/147‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCggg < 7:139934/147‑1 (MQ=38) cctgctCCATTGAAAGGCTGACGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GCGTAAACCCGCTCCTACCAGCGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGttt < 4:99872/148‑1 (MQ=38) ctgctCCATTGAAAGGTTGAGGATGTCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTa < 5:66604/149‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATTTCCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTAc < 3:274528/147‑1 (MQ=37) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATgc‑g > 1:271275/1‑93 (MQ=25) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGACC‑GT‑GCCGGCCTTGAATgc‑g < 2:271275/93‑1 (MQ=25) tcgcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg < 5:96432‑M1/144‑1 (MQ=255) cgcgG‑GTAAACCGGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTTAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGc > 3:65671‑M1/5‑149 (MQ=255) gcgG‑GTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAcgg < 4:65671‑M1/145‑2 (MQ=255) g‑gTAAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggag < 6:230591‑M2/148‑72 (MQ=255) AAACCCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGCTGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCgg < 3:92397‑M1/148‑1 (MQ=255) cccGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcg < 4:314996‑M2/147‑77 (MQ=255) CCCGCTCCTACCAGGGG‑CAGC‑GCCG‑CATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAc < 2:349667‑M1/137‑1 (MQ=255) CCCGCTCCTACCAGGGG‑CAGC‑GCCG‑CATTCAAGGCCGGCACGGTCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCAc > 1:349667‑M1/1‑137 (MQ=255) CCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTCGC‑‑CAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCt < 8:140143‑M1/149‑1 (MQ=255) CCGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACtt < 4:218181‑M1/149‑1 (MQ=255) CGCTCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTc > 3:218181‑M1/1‑149 (MQ=255) tCCTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagat < 7:326185‑M2/147‑82 (MQ=255) CTACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAc < 3:167320‑M1/149‑1 (MQ=255) tACCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct < 6:106505‑M2/149‑87 (MQ=255) aCCAGGGA‑CCGT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctc > 1:147340‑M2/1‑62 (MQ=255) | GACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑G‑GTAAACCCGCTCCTACCAGGGA‑C‑GT‑GCCGGCCTTGAATGC‑GGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACC > NC_002947/4585923‑4586195 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |