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 | 96.0% | 71.1 / ‑2.9 | 25 | 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 (8/16); total (8/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-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. |
GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑G‑GTAAACCCGCTCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGG > NC_002947/4585887‑4586097 | gccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑g‑g < 3:149333/149‑1 (MQ=255) ccAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑GTGTAAAc < 1:250645/149‑1 (MQ=255) cAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTTAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAAcc > 6:66224/1‑149 (MQ=255) gAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGc < 8:2421/149‑1 (MQ=255) cACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTcc > 7:281742/1‑148 (MQ=255) cACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCt < 8:322154/149‑1 (MQ=255) cACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCt < 8:263300/149‑1 (MQ=255) ggTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGAGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCa > 7:242595/1‑149 (MQ=39) ggTGAAGGCTTAAGGGTAGACCACGGCACCCGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTAc > 1:35505/1‑147 (MQ=255) aGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCg < 2:35505/149‑1 (MQ=39) gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCgg > 5:366781/1‑149 (MQ=39) gACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCGTGGCAAGGTGATGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCtt < 5:259280/148‑1 (MQ=38) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGaa < 6:222156/149‑1 (MQ=39) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑g‑ggt > 2:266078/1‑111 (MQ=255) ccACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTCGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGC‑g‑ggt < 1:266078/113‑3 (MQ=255) ctcGGCACCGGCGGGAACAGCTGATCTGCGCTGCGCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAg < 8:307231/127‑1 (MQ=38) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAAt > 6:372586/1‑149 (MQ=39) cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAAt < 8:281742/149‑1 (MQ=39) gCACCGGCGTTATCATCTGATCTGCGCTGATCCATTGAAAGGTTGAGGATGACGGGGTGCTTGTCAAGGTGAATATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGGTC‑CGTGCCTGCCTTGAATgcgg < 8:242595/149‑1 (MQ=25) tGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGcccc < 6:366781/149‑1 (MQ=39) atcTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa < 1:361553/147‑1 (MQ=38) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTa < 2:107204/147‑1 (MQ=39) cTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCATTGAATGCGGCGCTGCCCCTGGTa > 1:107204/1‑147 (MQ=38) cgcTGCTCCATTGAAAGGTTTAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 5:66224/146‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGc < 5:372586/144‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg < 5:350049/146‑1 (MQ=38) gATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt < 7:60157/78‑1 (MQ=18) gATCGGGCCGGCCTCTTCGCGG‑GTAAACCCGCTCCTACCAGGG‑ACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGt > 8:60157/1‑78 (MQ=18) | GCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑G‑GTAAACCCGCTCCTACCAGGG‑A‑CGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGG > NC_002947/4585887‑4586097 |
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: ‑ |