Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,173,295 | G→C | 20.7% | intergenic (‑288/+44) | fadD‑I ← / ← fadD‑II | long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,173,295 | 0 | G | C | 20.7% | 61.2 / 12.8 | 29 | intergenic (‑288/+44) | fadD‑I/fadD‑II | long‑chain‑fatty acid‑‑CoA ligase/long‑chain‑fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base G (11/12); new base C (2/4); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01 |
AAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAT‑GGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTG > NC_002947/5173151‑5173444 | aaaaCAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAgt < 1:97559/150‑1 (MQ=255) cAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAg < 2:24444/150‑1 (MQ=255) ggCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGa < 1:132707/150‑1 (MQ=39) ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTCATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg < 1:536211/150‑1 (MQ=255) ttATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCgg < 1:140174/150‑1 (MQ=255) aaaaaTGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCcgc < 2:336127/150‑1 (MQ=255) cTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGc < 2:244076/150‑1 (MQ=255) cTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg > 2:541397/1‑132 (MQ=38) cTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg < 1:541397/132‑1 (MQ=38) aaTACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTa > 1:275028/1‑150 (MQ=255) aTACAGGTCATTTTTAGAGTAATTATCCAAAATGCCGCAGAGAAGGCTGTATAGAGCCTGTTAGAGCGGTTTCAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTAt < 1:347543/150‑1 (MQ=255) gTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg < 2:94746/115‑1 (MQ=37) gTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACACTGCCCAATGTAGGAGCgg > 1:94746/1‑115 (MQ=37) aaTTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACg < 1:520975/134‑1 (MQ=255) aaTTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACg > 2:520975/1‑134 (MQ=255) tCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGaagaag > 1:223377/1‑113 (MQ=255) tCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGaagaag < 2:223377/113‑1 (MQ=255) aaaaTGTCGCAGATAAGGCTGTCTAGAGCCTGTTAGCGCGGTTTAAGGCATTCTCGGGTAAACCCACTCCTAAATTGGGCCGTGTAGGAGCGGGTTTACCCGAGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCa > 1:4987/1‑150 (MQ=255) gCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCCCCGCAAAATCTTg > 1:400625/1‑150 (MQ=255) aGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCg < 2:275028/150‑1 (MQ=255) gCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGc > 2:60954/1‑150 (MQ=255) tGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATcgc > 1:430319/1‑116 (MQ=255) tGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATcgc < 2:430319/116‑1 (MQ=255) aataGAGCCTTTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTTGTCATTGTAGGCTCGGGTTTACCCGCGAAGAAACCATCCCGGTATCACGCTATATCGCGCAAGTCCCGCCGCAAAATCTTGCCCACCGGCGTCa < 2:189705/149‑1 (MQ=255) cGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAT‑GGGCAACGATTCgcgc > 2:583183/1‑150 (MQ=255) gCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAT‑GGGCAACGATTCGCGCAGCACGAtg < 1:583183/150‑1 (MQ=255) cccGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAT‑GGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg > 1:388015/1‑150 (MQ=255) ccagCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACGGGCGTCAT‑GGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTg > 2:410776/4‑150 (MQ=255) gCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCACGCCGCAAAATCTTGACAACCGGCG‑CATCCAGCAGCGACGCGAGCAGCACTAAGTTCTTGGTTACTTTGTAg > 2:313785/1‑150 (MQ=255) ggCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCGCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAT‑GGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTg > 2:353650/1‑150 (MQ=255) | AAAACAGGCCTTATGGTCACAAAAATGACTCAATGAACAATACAGGTCATTTTTAGAGTAATTATCCAAAATGTCGCAGATAAGGCTGTATAGAGCCTGTTAGAGCGGTTTAAGGCATTCGCGGGTAAACCCGCTCCTACATTGGGCAGTGTAGGAGCGGGTTTACCCGCGAAGAAGCCATCCCGGTATCACGCTATATCGCGCAACTCCCGCCGCAAAATCTTGCCCACCGGCGTCAT‑GGGCAACGATTCGCGCAGCACGATGTGCTTGGGTACTTTGTAGCCGGTGAAGTTG > NC_002947/5173151‑5173444 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |