Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,586,030 | 2 bp→TC | 100% | 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 | 97.4% | 101.5 / ‑5.6 | 36 | intergenic (+113/+102) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (18/16); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
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. | |||||||||||
* | NC_002947 | 4,586,031 | 0 | T | C | 98.9% | 95.2 / ‑6.9 | 36 | 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); major base C (16/17); minor base A (2/0); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
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: ‑ |