| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 190,093 | T→G | 21.1% | P530P (CCT→CCG) | ldcC → | lysine decarboxylase 2, constitutive |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 190,093 | 0 | T | G | 21.1% | 88.4 / 36.5 | 57 | P530P (CCT→CCG) | ldcC | lysine decarboxylase 2, constitutive |
| Reads supporting (aligned to +/- strand): ref base T (27/18); new base G (6/6); total (33/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.43e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01 | |||||||||||
ATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTG > minE/190028‑190157 | aTCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATaa > 1:3132968/1‑71 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc > 1:397263/1‑71 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAcc > 1:1751251/1‑71 (MQ=255) cccGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:3751403/1‑70 (MQ=255) ggcggcGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgc < 1:3767392/71‑1 (MQ=255) gcggagCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACctgct > 1:2416242/1‑71 (MQ=255) gcTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtc < 1:2875556/71‑1 (MQ=255) ggTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:3534410/1‑58 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTtctc < 1:4049635/69‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:1680610/71‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:2393807/71‑1 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt > 1:2480545/1‑71 (MQ=255) gTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt < 1:1466335/71‑1 (MQ=255) aGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCtt > 1:1335423/1‑69 (MQ=255) gCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg < 1:738377/71‑1 (MQ=255) aaTACCTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAtc > 1:989448/1‑70 (MQ=255) ataccTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCAt < 1:1156892/66‑1 (MQ=255) accTCGACGAACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcctg < 1:2034573/58‑2 (MQ=19) ccTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCa > 1:2442605/1‑71 (MQ=255) cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt > 1:3564507/1‑71 (MQ=255) cTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCAt > 1:3512026/1‑71 (MQ=255) tCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc > 1:3450749/1‑71 (MQ=255) gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTta > 1:2721051/1‑38 (MQ=255) gacgaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATc < 1:1918597/69‑1 (MQ=255) gacgaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGa > 1:1369840/1‑71 (MQ=255) gaACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAg < 1:178726/56‑1 (MQ=255) gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct > 1:2411829/1‑50 (MQ=16) gaACATGGCATCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTcct > 1:1916936/1‑50 (MQ=16) cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAAc > 1:1728748/1‑36 (MQ=255) cGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaaa < 1:2821963/71‑2 (MQ=255) gTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt < 1:2450675/53‑1 (MQ=255) gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTa > 1:244781/1‑49 (MQ=255) gggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt > 1:3267638/1‑63 (MQ=255) ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa > 1:3703444/4‑70 (MQ=255) ggcaTCGTTGTTGAGAAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGATAAGACCa > 1:1480316/4‑70 (MQ=255) ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTAt > 1:3199325/1‑52 (MQ=255) ggATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa > 1:1260951/1‑66 (MQ=255) tCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCa < 1:3279586/71‑1 (MQ=255) cGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa > 1:4054129/1‑71 (MQ=255) tagtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGt < 1:3837901/44‑1 (MQ=255) agtagAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGAt > 1:4042529/1‑55 (MQ=255) gagaAAACCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt < 1:1978557/51‑1 (MQ=255) agaAAACCGGCCCTTATAACCTGCTGTTTCTCttt < 1:2328990/35‑1 (MQ=255) aaaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt < 1:1400862/48‑1 (MQ=255) aaaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCGAt < 1:123114/48‑1 (MQ=255) aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattatt > 1:2056157/1‑71 (MQ=255) aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattatt > 1:3405240/1‑71 (MQ=255) aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGattat > 1:1062013/1‑70 (MQ=255) aaaaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGatta > 1:3646824/1‑69 (MQ=255) aaaCCGGTCCGTATAACCTGCTGTTCCTGTTCAGCATCGGTATCg < 1:793907/45‑1 (MQ=255) aaCCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg > 1:3524962/1‑70 (MQ=255) cGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTg < 1:1231794/71‑1 (MQ=255) ccTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCaa < 1:3204363/52‑1 (MQ=255) cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa > 1:1858961/1‑42 (MQ=255) cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaaaa > 1:454097/1‑42 (MQ=255) cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat < 1:528183/57‑1 (MQ=255) cTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGat < 1:1104122/57‑1 (MQ=255) ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATaa < 1:601685/39‑1 (MQ=255) ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa < 1:526092/44‑1 (MQ=255) ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCa < 1:2296245/44‑1 (MQ=255) ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg < 1:413607/61‑1 (MQ=255) ttATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTg < 1:2585473/61‑1 (MQ=255) | ATCCCGGCGGCGCTGGTAGCAAAATTCCTCGACGAACGTGGGATCGTAGTAGAGAAAACCGGCCCTTATAACCTGCTGTTTCTCTTTAGTATTGGCATCGATAAAACCAAAGCAATGGGATTATTGCGTG > minE/190028‑190157 |
| 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: ‑ |