Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,806,160 | A→T | 12.1% | S593T (TCG→ACG) | PP_4234 ← | bifunctional aldehyde oxidase/xanthine dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,806,160 | 0 | A | T | 12.1% | 83.5 / 5.5 | 33 | S593T (TCG→ACG) | PP_4234 | bifunctional aldehyde oxidase/xanthine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (15/14); new base T (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 = 2.43e-01 |
CCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAA > NC_002947/4806013‑4806301 | ccGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAAc > 2:377019/1‑150 (MQ=255) ccGACGATGGCGACGCCAGCTTGCCGACCCAGAACGGGTTGGGCGGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAAc > 1:550733/1‑150 (MQ=255) gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCg < 1:470462/150‑1 (MQ=255) cccAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTAAACTCGTCGGCCGAAAGGCCGAGGAAATCGCTGACAAACAAAg > 1:450520/1‑150 (MQ=255) gggTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAg > 2:368808/1‑150 (MQ=255) gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGGCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt > 1:76015/1‑150 (MQ=255) gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt < 2:470741/150‑1 (MQ=255) tGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTgg < 1:368808/150‑1 (MQ=255) gcgTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTGCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAGGCCTGGGAGGTGGACATGCCGGTg < 2:632743/150‑1 (MQ=255) tCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacga < 2:672601/106‑1 (MQ=255) tCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGacga > 1:672601/1‑103 (MQ=255) tCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGccc < 1:76045/150‑1 (MQ=255) ctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCACT‑GCCATTCGGGTTACTGCATTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGata < 1:377019/150‑1 (MQ=255) ggTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGGCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCa < 2:76015/150‑1 (MQ=255) gTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGGGGACATGCCGGTGCCCAGCTCGGTGCCAATATGGCGCAGGCTGag > 2:44768/1‑149 (MQ=255) gCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa < 2:726617/62‑1 (MQ=255) gCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa > 1:726617/1‑62 (MQ=255) aaTTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTg > 1:611032/1‑150 (MQ=255) gCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg > 2:194950/1‑150 (MQ=255) cATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGcc > 2:250002/1‑150 (MQ=255) aTTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATACGGCCATCGGCGCTGAACTCGACACTGGCCa > 1:83821/1‑150 (MQ=255) ttC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCAt > 1:703789/1‑150 (MQ=255) tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACAt < 2:450190/69‑1 (MQ=255) tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACAt > 1:450190/1‑69 (MQ=255) tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACAAAAGCCTGGGAGGTGGACATGCCGGAGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCAATACGCGCCTc > 2:224049/1‑150 (MQ=255) ggcTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGaa < 2:110226/148‑1 (MQ=255) gTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGaa < 2:450520/150‑1 (MQ=255) tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAc < 1:44768/150‑1 (MQ=255) ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt < 1:194950/150‑1 (MQ=255) ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt > 1:402066/1‑150 (MQ=255) gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATACGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa < 2:83821/150‑1 (MQ=255) ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTgg > 2:585633/2‑50 (MQ=255) ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTgg < 1:585633/49‑1 (MQ=255) | CCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTC‑GGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAA > NC_002947/4806013‑4806301 |
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: ‑ |