Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,056:1 | +C | 73.7% | 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 | 73.7% | 89.1 / 32.3 | 38 | intergenic (+139/+76) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (5/5); new base C (14/14); total (19/19) | |||||||||||
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 |
CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585926‑4586198 | cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAAt > 5:23782/1‑149 (MQ=39) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg < 6:23782/149‑1 (MQ=39) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGGTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg < 1:135161/149‑1 (MQ=38) cGGCGTGAACAGCTGATATGCGCTGCTACATTGAAAGGTTGAGGATGAAAGGGTGCTTGGCAAGGTGAAGATATTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACAAGGTACC‑GTGCCGGCCGGGAATGCGGCGCt > 8:215360/1‑149 (MQ=25) aGCTGATCTGCGCTGCTCCATTGAAAGGTTGCGGATGACCGGGTGCTTGGCAAGGTGAAGATCATGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa < 7:70556/149‑1 (MQ=38) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg < 3:67362/149‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTgg < 6:191517/145‑1 (MQ=38) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTgg > 5:191517/1‑145 (MQ=38) ctctgctCACTTGAAAGGTTTCGGATGACCGGGTGCTTGTCAAGGTGAAGATCTTGCATTGCTCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCCGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCggg < 6:329813/145‑1 (MQ=37) ctCCATTGCAAGGTTGAGGATGACCGGGGGCTTGGCAAGGGGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc < 7:215360/147‑1 (MQ=37) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc < 6:305149/147‑1 (MQ=38) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg < 7:34025/131‑1 (MQ=38) aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAg > 8:34025/1‑131 (MQ=38) ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGa > 1:203041/1‑142 (MQ=38) ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGa < 2:203041/142‑1 (MQ=38) ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGATAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggcc > 8:180261‑M2/1‑142 (MQ=255) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga > 3:236461‑M2/1‑127 (MQ=255) ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAaga < 4:236461‑M2/130‑4 (MQ=255) ttGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGcccc > 4:217805/1‑98 (MQ=25) ttGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGcccc < 3:217805/98‑1 (MQ=25) gTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg > 5:253859/1‑81 (MQ=21) gTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATgcgg < 6:253859/81‑1 (MQ=21) tcgcggg‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg > 4:28339/6‑149 (MQ=255) gcggg‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCa > 8:134624/4‑149 (MQ=255) cggg‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCgcgc > 5:242260/3‑94 (MQ=255) cggg‑TAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCgcgc < 6:242260/92‑1 (MQ=255) aaCCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCGATCTCCGCCACACTGATCTCACGCATGCGGaa > 7:119830/1‑149 (MQ=255) cGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTc > 1:68536/1‑149 (MQ=255) cTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaa > 7:57117‑M2/1‑64 (MQ=255) cTACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaa < 8:57117‑M2/130‑67 (MQ=255) tACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct < 6:253855‑M2/149‑87 (MQ=255) tACCAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcaacgcagatcagct < 4:311607‑M2/149‑87 (MQ=255) ccAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcaacgcagatcagctg > 3:311607‑M2/1‑61 (MQ=255) cAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCcgc < 2:63358/73‑1 (MQ=255) cAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCcgc > 1:63358/1‑73 (MQ=255) cAGGGACC‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTg > 7:2335/1‑149 (MQ=255) cAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccga > 6:176509‑M2/1‑60 (MQ=255) cAGGGA‑CCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccga < 5:176509‑M2/74‑15 (MQ=255) | CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GG‑TAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTG > NC_002947/4585926‑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: ‑ |