Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 899,672 | T→C | 21.9% | L15P (CTG→CCG) | pspE → | thiosulfate:cyanide sulfurtransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 899,672 | 0 | T | C | 21.9% | 195.3 / 55.3 | 96 | L15P (CTG→CCG) | pspE | thiosulfate:cyanide sulfurtransferase |
Reads supporting (aligned to +/- strand): ref base T (41/34); new base C (10/11); total (51/45) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.26e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
AATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTC > minE/899602‑899740 | aaTTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACt > 1:2853038/1‑71 (MQ=255) aaTTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACt > 1:1674951/1‑71 (MQ=255) aTTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTg > 1:4110586/1‑71 (MQ=255) ttttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGc > 1:2253094/1‑71 (MQ=255) ttCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGccc > 1:1187676/1‑71 (MQ=255) tCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCg > 1:1339472/1‑71 (MQ=255) tCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCg > 1:2802077/1‑71 (MQ=255) tCATACATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGccc > 1:1258939/1‑70 (MQ=255) aTCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGtt < 1:1831720/71‑1 (MQ=255) cATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGtttt < 1:282394/70‑1 (MQ=255) gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTTCCCGTTTTcgcc < 1:999460/70‑1 (MQ=255) gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTTCCCGTTTTcgcc < 1:3750951/70‑1 (MQ=255) gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgcc < 1:3867457/70‑1 (MQ=255) gAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgcc < 1:3643922/70‑1 (MQ=255) ggACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTg < 1:1385180/71‑1 (MQ=255) ggACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTg < 1:30263/71‑1 (MQ=255) gACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGa < 1:2393291/71‑1 (MQ=255) cGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAac < 1:656659/71‑1 (MQ=255) tttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCg < 1:2500917/70‑1 (MQ=255) ttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGaa < 1:3667855/59‑1 (MQ=255) ttAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGAt > 1:3879977/1‑71 (MQ=255) tAAAAAGGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:542243/1‑71 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:3945242/1‑71 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:3184012/1‑71 (MQ=255) tAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:1705725/1‑71 (MQ=255) aaaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTcgc > 1:3875342/1‑50 (MQ=255) aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGc < 1:4214504/38‑1 (MQ=255) aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGc < 1:2452887/38‑1 (MQ=255) aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTc > 1:3929391/1‑71 (MQ=255) aaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTc > 1:3121831/1‑71 (MQ=255) aaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTc > 1:31702/1‑44 (MQ=255) aaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt < 1:3589835/71‑1 (MQ=255) cTTACTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGtt > 1:2112625/1‑70 (MQ=255) aCTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAaca > 1:1905512/1‑48 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGaa > 1:1446593/1‑45 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGaa > 1:1242015/1‑45 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGaa > 1:1488700/1‑45 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCag < 1:3740876/70‑1 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATCTTCGTGTTCCag < 1:58908/70‑1 (MQ=255) cTTGCTCTGGCACTGGTGTTTTCACCGCCCGTCTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCag < 1:865147/70‑1 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:1694324/53‑1 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATc < 1:2951732/53‑1 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:2692740/1‑71 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:2547711/1‑71 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:3402059/1‑71 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:2666025/1‑71 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:964083/1‑71 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCaga > 1:869599/1‑70 (MQ=255) ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCag > 1:4222468/1‑69 (MQ=255) tGCTCTGGCACTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:2592626/1‑70 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag < 1:1770664/69‑1 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:1460377/1‑69 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:4155046/1‑69 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag < 1:2695875/69‑1 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCagag > 1:589137/1‑69 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGc < 1:2110549/70‑1 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGc < 1:585265/70‑1 (MQ=255) gCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGc < 1:3085268/70‑1 (MQ=255) tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc > 1:1588973/1‑63 (MQ=255) tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc > 1:2221842/1‑63 (MQ=255) tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc > 1:3868698/1‑63 (MQ=255) tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc > 1:1780545/1‑63 (MQ=255) tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTcc > 1:1354717/1‑63 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:3951591/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt > 1:467954/1‑69 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt > 1:2904664/1‑69 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:577329/69‑1 (MQ=255) gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTAt < 1:2805925/69‑1 (MQ=255) aCTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCg < 1:3267146/51‑1 (MQ=255) cTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCa < 1:1891399/71‑1 (MQ=255) cTGGTGTTTTCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCa < 1:3899413/71‑1 (MQ=255) tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:1829885/1‑45 (MQ=255) tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:3142369/1‑45 (MQ=255) tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:1959078/1‑45 (MQ=255) tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATg > 1:2682987/1‑45 (MQ=255) tGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgtg > 1:21085/1‑51 (MQ=255) gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAag < 1:3288631/70‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt < 1:1706569/45‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt < 1:3290527/45‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCgt < 1:407799/45‑1 (MQ=255) ttttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCa < 1:3032676/71‑1 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:2937631/1‑70 (MQ=255) ttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg > 1:146253/1‑71 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:206264/1‑69 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:593895/1‑69 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:3675365/1‑69 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:3200008/1‑69 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:4022009/1‑69 (MQ=255) tCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCAc > 1:3354023/1‑69 (MQ=255) tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg < 1:3169548/70‑1 (MQ=255) tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg < 1:395561/70‑1 (MQ=255) tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg < 1:3279131/70‑1 (MQ=255) tCACCGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACg < 1:3353920/70‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc < 1:1450144/71‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc < 1:710735/71‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc < 1:3495540/71‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc < 1:3706235/71‑1 (MQ=255) aCTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTc < 1:116511/71‑1 (MQ=255) | AATTTTCATCCATAGAAGGACGCTTACATGTTTAAAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTC > minE/899602‑899740 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |