| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,228,052 | C→T | 100% | A152A (GCG→GCA) | yebZ ← | predicted inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,228,052 | 0 | C | T | 80.8% | 90.1 / 13.3 | 62 | A152A (GCG→GCA) | yebZ | predicted inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (6/6); new base T (50/0); total (56/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.50e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
GCATACACCAGAGCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCTGTTACCCCTTCATTCAGCGTCGCATGCCCCACT > minE/1227991‑1228112 | gCATACACCAGAGCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGAt < 1:2809539/69‑1 (MQ=255)gCATACACCAGAGCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGAt < 1:2867578/69‑1 (MQ=255) cacCAGAGCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGt > 1:1672728/1‑68 (MQ=255) agagCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGAt < 1:951497/60‑1 (MQ=255) agCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATgg < 1:726966/69‑1 (MQ=255) aGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGc > 1:2783190/1‑68 (MQ=255) aGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCt > 1:322011/1‑69 (MQ=255) aGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCt > 1:238528/1‑69 (MQ=255) aTCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCTGTTa > 1:2142872/1‑69 (MQ=255) ccccAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCTGTTAc > 1:1350125/1‑67 (MQ=255) ccccAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCTGTTAc < 1:988418/67‑1 (MQ=255) ccccAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCTGTTAc < 1:1746118/67‑1 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:3079915/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:3088244/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:3117263/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2499283/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2375165/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2372834/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2246683/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2801104/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2197423/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:3274164/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:2188215/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1850900/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:681573/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1578828/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:149958/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:140448/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:833955/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1389540/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1367242/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:912195/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1216701/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:920761/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1153215/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGAttt > 1:1020541/1‑40 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCa > 1:1022571/1‑69 (MQ=255) gcACAGATCAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCa > 1:1057461/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTAcc > 1:2722538/1‑43 (MQ=38) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGcccga > 1:536261/1‑66 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:1664605/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2712531/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:983980/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:1112286/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:1171816/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:122900/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:1392221/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:833340/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:79986/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:750605/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:282873/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:604701/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2827135/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:402936/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:169532/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2220655/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2418210/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2437585/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:3088046/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2490423/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATTCAACGTCGCATGCCCGAc > 1:2535103/1‑69 (MQ=255) cAGGTGAATTGCGTGATTAGTCTGATGGATTTTCGCCGTTACCCCTTCATT‑AACGTCGCATGCCCGACt > 1:1485329/1‑69 (MQ=255) | GCATACACCAGAGCACCGGCAGTAATCCCCCAAACCACGCCGCCGCACAAATCAGGTGAATCGCGTGATTAGTCTGATGGATTTTCGCTGTTACCCCTTCATTCAGCGTCGCATGCCCCACT > minE/1227991‑1228112 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |