Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,056:1 | +C | 51.6% | 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 | 51.6% | ‑5.6 / 56.9 | 31 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (6/9); new base C (4/12); total (10/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01 |
ACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585927‑4586198 | aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg > 7:274807/1‑149 (MQ=39) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgc < 6:2032/149‑1 (MQ=39) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcg < 5:385823/149‑1 (MQ=39) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg < 8:274807/149‑1 (MQ=39) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg < 5:244508/149‑1 (MQ=39) caaggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGG‑TAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTg > 5:243003‑M1/95‑149 (MQ=255) aCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTgg < 1:89105/149‑1 (MQ=39) ggtgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGG‑TAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa > 7:207330‑M1/92‑149 (MQ=255) gCTGATCTGCGCTGCTCCATTGACAGGTTGAGGATGACCGGGGGCTT‑GGCAAGGTGAAGATCTTGCATTGCTCGGGCCGGCCGCTTCGCGGG‑TACACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg < 3:249431/149‑1 (MQ=255) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg < 4:147157/149‑1 (MQ=39) tgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGG‑TAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg > 3:209498‑M1/90‑145 (MQ=255) tgcagaagttccgcatgcgtgagatcagtgtggcggagattgcggctgcatctgcgggttgatgtgctggcgcggccggcctcttcgcgGG‑TAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTgg < 4:209498‑M1/56‑1 (MQ=255) aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGa < 1:236007/147‑1 (MQ=39) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGtt > 2:139983/1‑149 (MQ=38) ctgctCCATTGAAAGGTTGAGGAAGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCTTTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTa < 7:251884/149‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc < 1:139983/147‑1 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc > 7:28983/1‑148 (MQ=38) aaaGGTGGAGTATCACCGGG‑GCTTGGGCAAGGTTAAGATGTTGAATTGATCGGGCCGGCCTCTTCGCGGG‑AAAACCCGCTCCTACCCGGGACC‑GGGCCGGCTTTTAATGCGGCGCTGCCCCTGGTTGGAGCGGTTTTACCCGCGAaga < 8:28983/148‑1 (MQ=255) ggTTGAGGATGCCCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACCCGGGA‑CCGGGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggcc > 2:47261‑M2/1‑142 (MQ=255) tt‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt < 2:184932/99‑1 (MQ=25) tt‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCt > 1:184932/1‑99 (MQ=25) ctcttcgcgGG‑TAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACACGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCCCCCTGAtct < 7:238018‑M1/140‑1 (MQ=255) AACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGAGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTACGCCACACTGATCTCACGCATGCGGaa < 5:266261‑M1/149‑1 (MQ=255) cTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagc < 8:129052‑M2/149‑86 (MQ=255) tACCAGGGA‑CCGTTCCGGCCTTTAATGCGCCGCTGCACCTGTTAGTAGCGGTTTAACCCGCTAagaggccgtcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct > 1:285616‑M2/1‑63 (MQ=255) tACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct < 5:107995‑M2/149‑87 (MQ=255) tACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct < 7:201363‑M2/149‑87 (MQ=255) aCCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagctg < 1:57906‑M2/149‑88 (MQ=255) cAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagctgtt < 3:86669‑M2/149‑90 (MQ=255) CAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTg > 2:133346‑M1/1‑142 (MQ=255) CAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGACGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTg < 1:133346‑M1/142‑1 (MQ=255) | ACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTT‑GGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585927‑4586198 |
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: ‑ |