Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 457,813 | C→T | 33.3% | Q41* (CAG→TAG) | ybgL → | predicted lactam utilization protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 457,813 | 0 | C | T | 33.3% | 68.4 / 49.8 | 57 | Q41* (CAG→TAG) | ybgL | predicted lactam utilization protein |
Reads supporting (aligned to +/- strand): ref base C (22/16); new base T (9/10); total (31/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 |
AACGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGAC > minE/457752‑457878 | aaCGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGc > 1:3835427/1‑62 (MQ=255) aaCGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTtgc > 1:3272348/1‑70 (MQ=255) gTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGcgtgcgtg < 1:2213861/70‑1 (MQ=255) tCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAg > 1:3983823/1‑70 (MQ=255) tctGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTtg < 1:1598323/56‑1 (MQ=255) tGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGGGTGAAGCAATa > 1:4112101/1‑70 (MQ=255) tGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATa > 1:3934280/1‑70 (MQ=255) ccAATATTGCCTGTGGATTTCATGCAGGCTATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATa < 1:2061929/68‑1 (MQ=255) cAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAg < 1:1119477/46‑1 (MQ=255) aaTATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAg > 1:1786738/1‑45 (MQ=255) atatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAAt < 1:3014271/71‑1 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCtt < 1:2508663/47‑1 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa < 1:1820242/69‑1 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaa < 1:3835156/66‑1 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaa < 1:3839800/66‑1 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaa < 1:3540747/66‑1 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATg > 1:3016769/1‑71 (MQ=255) tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATg > 1:214500/1‑71 (MQ=255) atTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTtgcg > 1:3723583/1‑49 (MQ=255) atTGACTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGc > 1:2681523/1‑59 (MQ=255) ttGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGgt > 1:2005888/1‑71 (MQ=255) ttGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCtt > 1:3367990/1‑45 (MQ=255) cTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATaaaa < 1:2064656/61‑1 (MQ=255) cTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGAGTGAAGCAATAAAAAATGGTGTcg > 1:3955938/1‑71 (MQ=255) tgtgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaa < 1:4056036/60‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgcg < 1:1724838/70‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgcg < 1:2541379/70‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgc < 1:1854743/69‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGcgtgcgtg < 1:3044354/47‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTtgc < 1:1780446/69‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgc < 1:3212319/69‑1 (MQ=255) tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGa > 1:1371059/1‑71 (MQ=255) gATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGAtt > 1:2041976/1‑71 (MQ=255) gATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGAtt > 1:1809768/1‑71 (MQ=255) aTTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaa < 1:44830/53‑1 (MQ=255) ttCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa > 1:2675740/1‑55 (MQ=255) ttCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa > 1:2497920/1‑55 (MQ=255) cAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCAc > 1:987112/1‑71 (MQ=255) gCGATGCGCAAATCATGTAGGCTTGCATGCGTGAAGCAATAAAAAATg < 1:1945943/48‑1 (MQ=39) cGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCAccc > 1:333278/1‑69 (MQ=255) gATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa > 1:2525955/1‑44 (MQ=255) gcgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCTCTCACCCGAGtt > 1:1979731/1‑70 (MQ=255) gcgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGttt > 1:555552/1‑71 (MQ=255) gcgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGtt > 1:617474/1‑70 (MQ=255) cgcAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCCCCCGAGtttt > 1:427799/1‑71 (MQ=255) cgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGttt < 1:2792444/70‑1 (MQ=255) cAAATCATGTAGGCTTGCGTGCGTGAAGCAATaaaaa > 1:3270021/1‑37 (MQ=255) cAAATCATGTAGGCTTGCGTGCGTGAAGCAATaaaaa > 1:935573/1‑37 (MQ=255) cAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgc < 1:3222428/47‑1 (MQ=255) cAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa < 1:3273294/38‑1 (MQ=255) cAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa < 1:806176/38‑1 (MQ=255) cAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTc < 1:104190/70‑1 (MQ=255) aaaTCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCa > 1:1088487/1‑58 (MQ=255) aaaTCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCa > 1:2678156/1‑58 (MQ=255) tCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTccc > 1:1578746/1‑68 (MQ=255) tCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGAc < 1:1764423/71‑1 (MQ=255) tCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGa > 1:1249778/1‑45 (MQ=255) tCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGAt > 1:1597783/1‑46 (MQ=255) | AACGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGAC > minE/457752‑457878 |
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: ‑ |