Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,099,555 | C→A | 12.5% | intergenic (‑116/‑153) | acsA‑I ← / → PP_4488 | acetyl‑CoA synthetase/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,099,555 | 0 | C | A | 12.5% | 81.7 / 3.0 | 32 | intergenic (‑116/‑153) | acsA‑I/PP_4488 | acetyl‑CoA synthetase/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (13/15); new base A (2/2); total (15/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 = 7.25e-01 |
TCCGGCAAAACGCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAGGCAAGCGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGG > NC_002947/5099413‑5099697 | tCCGGCAAAACGCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAg < 1:582231/150‑1 (MQ=255) gCAAAACGCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGAGCAATAAAATTCACCTCGCCCGCacc > 2:79794/1‑148 (MQ=255) cAAAACGCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACg < 1:13314/150‑1 (MQ=255) aaaCGCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGc > 1:22232/1‑150 (MQ=255) ctctttGGCGCGGGCATTACAGCCTTCCCGAGGAAAAAACGAGGGGGTTACCACAGAATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCTACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAg < 2:22232/145‑1 (MQ=255) gCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAgg > 1:630088/1‑150 (MQ=255) ggCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGAtt < 1:588938/150‑1 (MQ=255) gCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACt < 1:556825/150‑1 (MQ=255) gCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACt > 2:484984/1‑150 (MQ=255) cATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTg < 2:629546/150‑1 (MQ=255) ttACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGt < 2:501399/150‑1 (MQ=255) aCAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTaa < 1:484984/150‑1 (MQ=255) aaaaaaCGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAgg < 1:329205/150‑1 (MQ=255) cacaGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGc < 1:570830/121‑1 (MQ=255) cacaGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGc > 2:570830/1‑121 (MQ=255) aGGATTGTTACAGATTCTGTCAAAACACACTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTCACACCCCCCACGAAGGCAAGCGTTAATCTCCATCTaa > 2:355130/1‑150 (MQ=255) ttACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAg > 2:732770/1‑89 (MQ=255) ttACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAg < 1:732770/89‑1 (MQ=255) aTTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACATCGCCAGCAACg < 1:668191/78‑1 (MQ=255) aTTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACATCGCCAGCAACg < 1:667724/78‑1 (MQ=255) aTTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACATCGCCAGCAACg > 2:667724/1‑78 (MQ=255) aTTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCAAGCGCCATAAAATTCACATCGCCAGCAACg > 2:668191/1‑78 (MQ=255) gTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAGGCAAGCGTTAATCTCCATCTAATCCCGATAGTTAAAACTa < 2:630088/150‑1 (MQ=255) tCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCAGCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAGGCACGCGTTAATCTCCATCTAATCCCGATAGTTAAAACTAg > 1:174790/1‑150 (MQ=255) tttATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGAt > 2:473524/1‑73 (MQ=255) tttATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGAt < 1:473524/73‑1 (MQ=255) ccACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGa < 1:616833/59‑1 (MQ=255) ccACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGa > 2:616833/1‑59 (MQ=255) cACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGcccc < 2:274154/79‑1 (MQ=255) cACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGcccc > 1:274154/1‑79 (MQ=255) aTAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAGGCAAGCGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGa > 2:21046/1‑150 (MQ=255) aaTTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAGGCAAGCGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAgg > 1:669959/1‑150 (MQ=255) | TCCGGCAAAACGCTCTGGCGCGGGCATTACAGCCTTGCCGACGAAAAAACGAGGGGTTTGCCACAGGATTGTTACAGATTCTGTCAAAACACTTTATCAAAACACCCACTGTTTCGGCCTTGGCCAGCGCCATAAAATTCACCTCGCCAGCAACGGCAAGGCGATTAACAACTGGTAACAGCCCCCACGAAGGCAAGCGTTAATCTCCATCTAATCCCGATAGTTAAAACTAGCAATACTCGCGGCGCCATAAGCGCCGCGTGCCCCCTCACGACCTTAGACAGG > NC_002947/5099413‑5099697 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |