Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,569,200 | G→A | 12.1% | intergenic (+295/‑32) | pcaK → / → pcaF‑I | 4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,569,200 | 0 | G | A | 12.1% | 80.3 / 3.7 | 33 | intergenic (+295/‑32) | pcaK/pcaF‑I | 4‑hydroxybenzoate transporter/beta‑ketoadipyl‑CoA thiolase |
Reads supporting (aligned to +/- strand): ref base G (15/14); new base A (2/2); total (17/16) | |||||||||||
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.75e-01 |
AGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGAAGGCGCTGATCGAG > NC_002947/1569064‑1569348 | aGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCACTTGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCa > 1:167529/1‑150 (MQ=255) ggTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAg > 2:598530/1‑150 (MQ=255) taccAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCa < 1:598530/150‑1 (MQ=255) accAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAg > 2:569623/1‑150 (MQ=255) cAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGga < 1:553112/150‑1 (MQ=255) tAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACTCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGa < 1:184925/150‑1 (MQ=255) tAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGa < 1:184620/150‑1 (MQ=255) aTCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGga < 2:254090/139‑1 (MQ=255) aTCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGga > 1:254090/1‑139 (MQ=255) tCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGc > 1:570043/1‑150 (MQ=255) ggTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCacga < 1:163565/150‑1 (MQ=255) aaGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTAACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCAt < 2:338871/150‑1 (MQ=255) cAGGAGCCTCCAAATGCCACGCGCATACCATTCACAAACCCACCCAAATCGAGCTGTCCGCCATGCGCGACCTCGACGCCCGTACCGCCACTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATATGTGACGCCAt > 1:386767/1‑150 (MQ=255) aCGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATGTGTGACGCCATCCGTACCCCGATCGGCCg > 1:701235/1‑150 (MQ=255) cAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCt < 1:63733/117‑1 (MQ=255) cAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCt > 2:63733/1‑117 (MQ=255) aGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCATACCCCGATCGGCCGCTTcgg > 1:744808/1‑150 (MQ=255) cATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCggcggc > 1:81659/1‑150 (MQ=255) ttCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGcc < 1:185962/150‑1 (MQ=255) cGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTg > 2:215237/1‑150 (MQ=255) acccaAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCg > 1:109457/1‑150 (MQ=255) ccaAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCACTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTg > 1:210967/1‑150 (MQ=255) aaaaTCGAGCTGTCCGCCATGCGCGACCTCGACGCAGGGACCGCCACTGACCCTGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGc < 2:386767/149‑1 (MQ=255) tCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGAcc > 2:176483/1‑150 (MQ=255) cGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTTGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTg < 2:364099/150‑1 (MQ=255) gCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTgg < 2:109457/150‑1 (MQ=255) gCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTgg > 2:705366/1‑150 (MQ=255) ccATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGc > 2:98951/1‑87 (MQ=255) ccATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGc < 1:98951/87‑1 (MQ=255) tCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCa > 2:480973/1‑77 (MQ=255) tCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCa < 1:480973/77‑1 (MQ=255) gACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGa < 2:740262/150‑1 (MQ=255) cGTACCGCCACTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGAAGgcgc < 2:210967/150‑1 (MQ=255) ccgccgCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGAAGGCGCTGAt > 1:720299/1‑150 (MQ=255) cgcTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGAAGGCGCTGATCGAg < 2:46674/150‑1 (MQ=255) cgcTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGAAGGCGCTGATCGAg < 1:176483/150‑1 (MQ=255) | AGGTACCTACCAGTCTCCTTAATCTGGGTTTACCGAAGCACCCTGACCAGGAGCCTCCAAATGCCACGCGCAGACCATTTCCGAACCCACCCAAATCGAGCTGTCCGCCATGCGTGACCTCGACGCCCGTACCGCCGCTGACCATGGCCAGCTGGCAGGAGAAGCCTGATGCACGACGTATTCATCTGTGACGCCATCCGTACCCCGATCGGCCGCTTCGGCGGCGCCCTGGCCAGCGTGCGGGCCGACGACCTGGCCGCCGTGCCGCTGAAGGCGCTGATCGAG > NC_002947/1569064‑1569348 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |