| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,976,696 | C→T | 36.3% | intergenic (+24/‑53) | galP → / → sprT | D‑galactose transporter/conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,976,696 | 0 | C | T | 36.3% | 33.9 / 30.4 | 55 | intergenic (+24/‑53) | galP/sprT | D‑galactose transporter/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (15/20); new base T (5/15); total (20/35) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGAAAAC > minE/1976634‑1976756 | aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:3043593/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2332509/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2289386/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2154429/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2472992/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:206264/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2690819/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2981291/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:3293021/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:1641400/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:1454502/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:645463/69‑1 (MQ=255)aaaGGTCGTAAACTGCGCGAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:1167390/69‑1 (MQ=255) gTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGcc < 1:3321251/69‑1 (MQ=255) cGTAAACTGCGCGAAATATGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:2649210/69‑1 (MQ=255) cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:1548149/69‑1 (MQ=255) cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:2625165/69‑1 (MQ=255) cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:3436286/69‑1 (MQ=255) cGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:354378/69‑1 (MQ=255) gTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:2572398/68‑1 (MQ=255) gTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:2456886/68‑1 (MQ=255) aaaCTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:1949984/66‑1 (MQ=255) cgcgAAATCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:3033371/54‑1 (MQ=255) cgcgAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:190218/60‑1 (MQ=255) cgAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTc > 1:575430/1‑69 (MQ=255) cgAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTc > 1:182946/1‑69 (MQ=255) aaaTCGGCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATggagga < 1:2214146/50‑1 (MQ=255) aaaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCAt < 1:3381842/69‑1 (MQ=255) aaaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCAt < 1:1748310/69‑1 (MQ=255) aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc < 1:3144606/69‑1 (MQ=255) aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc < 1:1412990/69‑1 (MQ=255) aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc < 1:3259295/69‑1 (MQ=255) aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc < 1:2460781/69‑1 (MQ=255) aaTAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATc < 1:247650/69‑1 (MQ=255) atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt > 1:934969/4‑69 (MQ=255) atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt > 1:1234768/4‑69 (MQ=255) atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt > 1:596624/4‑69 (MQ=255) atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt > 1:2533900/4‑69 (MQ=255) atcggCGCTCACGATTAATCTCCTCAGGCTTCCTCCCATCATGGAGGAAGCCACCTCTTGCAGTCATCt > 1:3290564/4‑69 (MQ=255) gcgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCAc < 1:1660064/50‑1 (MQ=255) cgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTttct > 1:2857005/1‑69 (MQ=255) cgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTttct > 1:628232/1‑69 (MQ=255) cgcTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTttc > 1:1540356/1‑68 (MQ=255) tAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATc > 1:2868335/1‑69 (MQ=255) aTCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCct > 1:3018651/1‑69 (MQ=255) tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCctct > 1:3073049/1‑68 (MQ=255) tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCctct > 1:76762/1‑68 (MQ=255) tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTg > 1:3232412/1‑69 (MQ=255) tcCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTg > 1:2345227/1‑69 (MQ=255) cccAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGcc > 1:1583595/1‑69 (MQ=255) cccAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGcc > 1:1836475/1‑69 (MQ=255) tcctccCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGaa > 1:2114436/1‑69 (MQ=255) tcctccCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGaa > 1:476875/1‑69 (MQ=255) tccCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGaaaa < 1:3342548/68‑1 (MQ=255) cccATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGAAAAc < 1:1018433/68‑1 (MQ=255) | AAAGGTCGTAAACTGCGCGAAATAGGCGCTCACGATTAATCTCCCCAAGCTTCCTCCCATCGCGGAGGAAGCCACCTCTTGCAGTCATCTTTTCTTCGCTCTATCCTCTGCCGCTATGAAAAC > minE/1976634‑1976756 |
| 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: ‑ |