Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 34,492 | T→C | 48.3% | R342R (CGA→CGG) | caiC ← | predicted crotonobetaine CoA ligase:carnitine CoA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 34,492 | 0 | T | C | 48.3% | ‑5.5 / 49.6 | 56 | R342R (CGA→CGG) | caiC | predicted crotonobetaine CoA ligase:carnitine CoA ligase |
Reads supporting (aligned to +/- strand): ref base T (9/20); new base C (18/9); total (27/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.53e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-02 |
TGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAA > minE/34431‑34560 | tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt > 1:3103770/1‑69 (MQ=255) tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt > 1:2005865/1‑69 (MQ=255) tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt > 1:688466/1‑69 (MQ=255) tgtgATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTAt > 1:1673844/1‑69 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:913686/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:2156192/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:2641760/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:3303440/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:534624/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:163544/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:1238546/69‑1 (MQ=255) aTCGTCGCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:698033/69‑1 (MQ=255) tcgtcgCGGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGcc < 1:738024/68‑1 (MQ=255) cgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAg < 1:1791720/69‑1 (MQ=255) cgtcgCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAg < 1:2942670/69‑1 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCa > 1:1441268/1‑64 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAgggct > 1:2083482/1‑66 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:617291/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:39801/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:3307254/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:3080975/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:2724977/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:981021/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:1742977/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:2035523/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAAGCCCGCCCGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCg > 1:2054424/1‑69 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg < 1:2150487/69‑1 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg < 1:858926/69‑1 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg < 1:754091/69‑1 (MQ=255) gcgGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACg < 1:1685169/69‑1 (MQ=255) cTTCGTAGCAAAACCCCACCCGACCAATCGACGTCCAGCGTCGTTTATCGCCAGGACGATCGCCGAt < 1:2875741/67‑1 (MQ=255) cTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGAt < 1:2991562/67‑1 (MQ=255) gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac < 1:1216865/69‑1 (MQ=255) gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac < 1:46654/69‑1 (MQ=255) gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac < 1:2087627/69‑1 (MQ=255) gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac < 1:2113037/69‑1 (MQ=255) gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac < 1:2704419/69‑1 (MQ=255) gCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCcac < 1:2233556/69‑1 (MQ=255) aCCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATgg > 1:1046631/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:1763214/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:525993/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:74124/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:3120966/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:2909198/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:2127314/1‑69 (MQ=255) ccGACCAATCGACGGCCAGCGCCGTTTATCGCCAGGGCGATCGCCGATAATGCCCACAATGGTTTCCGt > 1:2492063/1‑69 (MQ=255) cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt > 1:3317422/1‑55 (MQ=255) cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt > 1:3228199/1‑55 (MQ=255) cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt > 1:2322124/1‑55 (MQ=255) cAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAAt > 1:2145839/1‑55 (MQ=255) aaTCGACGGCCAGCGTCGTTTATCGCCAGGCCGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc < 1:1864592/69‑1 (MQ=255) aaTCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc < 1:892129/69‑1 (MQ=255) aaTCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc < 1:940830/69‑1 (MQ=255) aTCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATccc < 1:1141228/68‑1 (MQ=255) gACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATaa < 1:1710907/69‑1 (MQ=255) gACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATaa < 1:1154894/69‑1 (MQ=255) tcgtTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCaa > 1:1481386/1‑69 (MQ=255) | TGTGATCGTCGCGGATCTCCGCTTCGTAGCAAAACCCCACCCGACCAATCGACGGCCAGCGTCGTTTATCGCCAGGACGATCGCCGATAATGCCCACAATGGTTTCCGTCATCCCATAAGACGTCAGCAA > minE/34431‑34560 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |