Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_002947 | 4,586,030 | 2 bp→TC | intergenic (+113/+101) | 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,030 | 0 | C | T | 96.3% | 77.8 / ‑4.7 | 27 | intergenic (+113/+102) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (7/19); total (8/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.96e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
* | NC_002947 | 4,586,031 | 0 | T | C | 100.0% | 86.1 / NA | 27 | intergenic (+114/+101) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (8/19); total (8/19) |
GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑CTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGA > NC_002947/4585887‑4586109 || gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGC‑gg < 2:204809/149‑1 (MQ=255) ccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGC‑ggg < 5:260940/149‑2 (MQ=255) cgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGC‑gggt < 5:371991/149‑3 (MQ=255) gAGTCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGC‑gggt > 3:49481/1‑146 (MQ=255) gCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGC‑gggtaa > 4:180357/1‑144 (MQ=255) catgAGCGACCGGTGAAGTCTTAATGGTAGACCACGGCACCGGGGTTAACAGATGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGTGTGCTTGGCAAGGTGAAGATCTTTCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAAc < 4:49481/145‑1 (MQ=255) gAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCg > 3:112082/1‑148 (MQ=255) gCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCt > 1:356758/1‑148 (MQ=255) gCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCt < 2:356758/148‑1 (MQ=255) ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGA‑cc < 7:2455/149‑2 (MQ=255) cTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGA‑ccgt > 2:9516/1‑146 (MQ=39) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAg < 2:172645/141‑1 (MQ=255) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAg > 1:172645/1‑141 (MQ=255) aaGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGcc > 2:437233/1‑149 (MQ=39) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCg < 7:167201/149‑1 (MQ=39) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGaa < 5:77137/149‑1 (MQ=39) cGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTAcc < 2:94325/125‑1 (MQ=255) gCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCCACCAGGGACC‑GTGCCGGCCTTGAATg‑cgg < 1:413894/149‑1 (MQ=39) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGcccc < 3:257625/149‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTg < 1:79964/149‑1 (MQ=39) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAg < 1:9516/149‑1 (MQ=39) accTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGGTAGGAGc < 5:396513/147‑1 (MQ=14) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGtt < 1:446701/149‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTAccc < 8:429139/149‑1 (MQ=38) cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAAt < 2:55098/113‑1 (MQ=37) aTTGAAAGGGTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑TCGCGGGTAAACCCGCTCCTACCCGGGACC‑GGGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGa < 3:180357/149‑1 (MQ=37) ggATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATACATCGGGCCGGACTCTACCCG‑GGTAAAACCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATg‑cc > 7:84599/1‑101 (MQ=255) || GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCT‑CTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGA > NC_002947/4585887‑4586109 |
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: ‑ |