| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,784,184 | T→C | 100% | I223T (ATA→ACA) | yfjB → | NAD kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,784,184 | 0 | T | C | 84.6% | 28.3 / ‑0.0 | 13 | I223T (ATA→ACA) | yfjB | NAD kinase |
| Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (11/0); total (12/1) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.54e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTG > minE/1784141‑1784236 | gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCAc > 1:183133/1‑59 (MQ=255)gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGAt > 1:194819/1‑62 (MQ=255)gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGAt > 1:318790/1‑62 (MQ=255)gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGAt > 1:543332/1‑62 (MQ=255)gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGAt > 1:57420/1‑62 (MQ=255)gTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCACAAACAGCAGCAGCACGAt > 1:6295/1‑62 (MQ=255) gCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGttttt > 1:51476/1‑62 (MQ=255) cACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc > 1:138096/1‑62 (MQ=255) cACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc > 1:205084/1‑62 (MQ=255) cACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc > 1:35606/1‑62 (MQ=255) cACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc > 1:364635/1‑62 (MQ=255) cACGACCACTGGTCACAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc > 1:96031/1‑62 (MQ=255) cACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTg < 1:476444/62‑1 (MQ=255) | GTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTG > minE/1784141‑1784236 |
| 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: ‑ |