| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 351,261 | T→C | 23.3% | N137D (AAT→GAT) | purK ← | N5‑carboxyaminoimidazole ribonucleotide synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 351,261 | 0 | T | C | 23.3% | 119.1 / 30.7 | 94 | N137D (AAT→GAT) | purK | N5‑carboxyaminoimidazole ribonucleotide synthase |
| Reads supporting (aligned to +/- strand): ref base T (17/55); new base C (17/5); total (34/60) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.79e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 | |||||||||||
AATGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACA > minE/351202‑351316 | aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTaa > 1:3088289/1‑69 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTaa > 1:980309/1‑69 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTaa > 1:2012606/1‑69 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTaa > 1:1852840/1‑69 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTaa > 1:1532825/1‑69 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTaa > 1:1440737/1‑69 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTa > 1:1008406/1‑68 (MQ=255)aaTGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTa > 1:1476385/1‑68 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:2451658/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:2507179/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:1930832/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:182839/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:279492/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:796941/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGACATTGACCGCGAc > 1:2051044/1‑69 (MQ=255) tACATTCGCCGTAACACTCTGCCGGTAACAGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc > 1:2706188/1‑69 (MQ=255) cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt < 1:3193566/69‑1 (MQ=255) cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt < 1:2780012/69‑1 (MQ=255) cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt < 1:1596286/69‑1 (MQ=255) cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt < 1:2757272/69‑1 (MQ=255) cGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCAt < 1:3239612/69‑1 (MQ=255) gCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt < 1:2293477/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTTCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:2532033/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:737158/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:1085674/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:139280/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:3397177/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:2465022/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:1894580/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:3251314/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:187543/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:429586/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:1432436/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:3224990/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:3134605/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:461115/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:1643377/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:879667/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:806346/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGACATTGACCGCGACCGTCATa < 1:614350/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGATTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:219233/68‑1 (MQ=255) ccGTAACACTCTGCCGGTAACTGTTCGGATTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:2657158/68‑1 (MQ=255) cGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa < 1:1718491/67‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgct < 1:1474050/54‑2 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:3163527/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:2481855/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:2097213/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:2437233/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:1550911/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:525506/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:928379/54‑1 (MQ=255) acacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:1714375/54‑1 (MQ=255) cacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:3220416/53‑1 (MQ=255) cacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:2007146/53‑1 (MQ=255) cacTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACcgcg < 1:1634559/53‑1 (MQ=255) ctctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGt > 1:2359285/1‑69 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCTCGTAAACGCCATTg < 1:2098348/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:998588/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:415782/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:3150099/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:3409942/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:3182828/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:3298530/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:3205514/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:264274/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:1880056/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:1239433/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:2360762/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:2250230/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTg < 1:1235162/43‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc < 1:1329288/69‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc < 1:401565/69‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc < 1:2847139/69‑1 (MQ=255) ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc < 1:2381452/69‑1 (MQ=255) gCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc < 1:106759/67‑1 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACcgcg > 1:1065995/1‑37 (MQ=38) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGt > 1:1684451/1‑42 (MQ=37) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGGTTAAc > 1:2506270/1‑67 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTaa > 1:1562334/1‑66 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAAc > 1:2664765/1‑67 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAAc > 1:3140078/1‑67 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:1234713/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:1259342/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:1272807/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:3417717/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:1829828/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:1833971/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:1917195/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:2122539/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:2195853/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:2503493/1‑68 (MQ=255) gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa > 1:274864/1‑68 (MQ=255) gTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc < 1:2550347/52‑1 (MQ=255) tttCATTTGCGCGTAAACGCCATTGACCGCGACCGTGATa < 1:742460/40‑1 (MQ=255) | AATGCCCTGCTCGACAATACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACA > minE/351202‑351316 |
| 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: ‑ |