Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 38,036 | G→T | F360L (TTC→TTA) | caiB ← | gamma‑butyrobetainyl‑CoA:carnitine CoA transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 38,036 | 0 | G | T | 89.5% | 52.6 / 0.6 | 19 | F360L (TTC→TTA) | caiB | gamma‑butyrobetainyl‑CoA:carnitine CoA transferase |
Reads supporting (aligned to +/- strand): ref base G (0/2); new base T (9/8); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.74e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01 |
GCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTGAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTGATTGATTCGCGAGCCA > NC_000913/37958‑38112 | gCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTaa > 2:7602/1‑81 (MQ=255) cAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCatga < 2:7577/80‑1 (MQ=255) aaaaTGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCatgat < 1:7578/80‑1 (MQ=255) aaTGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGtt < 1:7579/81‑1 (MQ=255) aaTGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGtt > 2:7608/1‑81 (MQ=255) aaTGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGtt > 2:7605/1‑81 (MQ=255) gCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCt > 2:7615/1‑81 (MQ=255) catgccatgTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGt < 1:7580/81‑1 (MQ=255) aGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTGAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGAccatcc < 1:7581/80‑1 (MQ=255) ttCCGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGttt < 2:7582/81‑1 (MQ=255) ccGCGCCAGATTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGc > 2:7612/1‑81 (MQ=255) aTTTGTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGtg < 2:7583/81‑1 (MQ=255) tttGTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGtga > 1:7614/1‑81 (MQ=255) gTCCGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGtgatt > 1:7617/1‑80 (MQ=255) tCCGGGGTTATTTTTGAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGtgattg < 2:7584/80‑1 (MQ=255) cGGGGTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTgattgatt > 2:538423/1‑81 (MQ=255) gggTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTGATTGATTcg < 2:7585/81‑1 (MQ=255) gggTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTGATTGATTcg > 2:7610/1‑81 (MQ=255) gTTATTTTTTAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTGATTGATTcgcg < 1:7586/81‑1 (MQ=255) ttttAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTGATTGATTCGCGAGCca < 2:7587/80‑1 (MQ=255) | GCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTTTTGAATTTCGGCATGATGTTCGGCCCTTTGCAGGTGCGACCATCCATCGTTTGCCACTGAGTGATTGATTCGCGAGCCA > NC_000913/37958‑38112 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |