Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,033:1 | +G | 100% | intergenic (+116/+99) | 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 | 97.1% | 103.2 / ‑2.8 | 35 | intergenic (+116/+99) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (0/1); new base G (17/17); total (17/18) | |||||||||||
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.73e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGA > NC_002947/4585887‑4586169 | gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGCCCGGGTGCTTGGCAAGGTTAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGC‑gg > 5:236786/1‑149 (MQ=255) aGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAAccc > 6:184361/1‑149 (MQ=255) gCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTc < 8:162515/149‑1 (MQ=255) aCCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTa > 6:59285/1‑149 (MQ=255) ggTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCa < 6:352665/149‑1 (MQ=255) tGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTTCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑gg < 7:244209/149‑1 (MQ=39) ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGCTCTGCGCTGCTCCATTGAAAGGTGGAGGATGACCGGGTGCTTGGCAAGGTGAAGCGCTTGCATGGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑gggccc > 8:216536/1‑146 (MQ=38) ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑ggggca > 8:272952/1‑146 (MQ=39) ggCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAAACCGCTACTACCA‑ggggcc > 6:82887/1‑146 (MQ=255) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑ggggca < 4:92806/146‑4 (MQ=39) ttAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑ggggca > 3:92806/1‑143 (MQ=39) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCg < 5:57061/149‑1 (MQ=39) ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAAc < 6:249427/123‑1 (MQ=255) ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAAc > 5:249427/1‑123 (MQ=255) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTAcc > 1:108825/1‑127 (MQ=255) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTAcc < 2:108825/127‑1 (MQ=255) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATg‑cg < 5:184361/149‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTGCATTGTAAGGTTGAGGATGTCCGGGTGCTTGGCAAGGTTAAGATCTTGGATTGATCGGGCCGGCCTCTTCGGG‑GGTAAACCCGCTCCTACGT‑GGGACC‑GTGCCGGCGTTGAATG‑CGGCGCTGCCCCTg < 6:236786/149‑1 (MQ=255) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt > 1:80456/1‑130 (MQ=38) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt < 2:80456/130‑1 (MQ=38) gCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAg < 5:59285/149‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAg > 4:225076/1‑148 (MQ=39) aTCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGa < 2:333616/147‑1 (MQ=39) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCgg < 3:225076/146‑1 (MQ=38) ctgctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTAATCGGGCCGGCCTC‑TTAGCGGGTAAAACCCCTCCTAACAGGGGC‑C‑CTGGCGGGCGTGAATT‑AGGCGCTGCCCCGGG‑TAGGAGCGGGgtt > 1:98486/1‑148 (MQ=255) ctCCCTTGAAAGGTTGAGGCTGACCGGGCGCTTGGCAAGGTGAAGATCTTGCATTGCTCGGGCAGGACTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTAc < 7:216536/147‑1 (MQ=37) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTTAAACCGCGCCTAACACGGGC‑C‑GTTGCGGGCGTGTATT‑AGGGGCGGGCCCAGGTTTGGCGCGGGTTTAcc > 3:319391/1‑149 (MQ=255) ccATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TGCGCGGGTAAACCCGCTCCTACCA‑GGGG‑CAGCGCC‑GCATTCAAGGCCGGCACGGTCCCTGG‑TAGGAGCGGGTTTACCcgc < 1:214946/149‑1 (MQ=14) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCgg < 2:229419/100‑1 (MQ=25) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCgg > 1:229419/1‑100 (MQ=25) ggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAggcc < 5:82887/149‑1 (MQ=38) tGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTAGCGGGTAAACCCGCACCTAACA‑GGGACC‑GTGCCGGCCGTGACTG‑CGGGGCTGGCCCGGG‑GATGGTCGGGGTTTCCGGCGACGAggcgggc > 4:221871/1‑149 (MQ=255) cTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGA‑ccgt > 6:189711/1‑70 (MQ=18) ggTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt < 7:320659/78‑1 (MQ=18) ggTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TTCGCGGGTAAACCCGCTCCTACCA‑GGGACC‑GTGCCGGCCTTGAAt > 8:320659/1‑78 (MQ=18) ggccTC‑TTCGCGGGTAAACCCGCTACTACCAGGGGT‑A‑GCGCC‑GCATTCAAGGCCGGCACGGTCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGGGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGa > 6:351747/1‑149 (MQ=9) | GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTC‑TCTGC‑GGTAAACCCGCTCCTACCA‑GGGA‑C‑GTGCCGGCCTTGAATG‑CGGCGCTGCCCCTGG‑TAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGA > NC_002947/4585887‑4586169 |
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: ‑ |