| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,042,145 | C→G | 64.0% | T26T (ACC→ACG) | gst → | glutathionine S‑transferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,042,145 | 0 | C | G | 64.0% | 62.8 / 58.8 | 75 | T26T (ACC→ACG) | gst | glutathionine S‑transferase |
| Reads supporting (aligned to +/- strand): ref base C (4/23); new base G (31/17); total (35/40) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.34e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTA > minE/1042078‑1042212 | tCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTAcc < 1:618759/68‑1 (MQ=255) cAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTc > 1:1942731/1‑67 (MQ=255) cAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTc > 1:1001528/1‑67 (MQ=255) aaCCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTCGTc < 1:2077872/68‑1 (MQ=255) aaCCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTCGTc < 1:3356055/68‑1 (MQ=255) aaCCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTCGTc < 1:2626448/68‑1 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGCGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:1626645/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAgtgt > 1:1090662/1‑61 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:3232183/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:1944847/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:2079462/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:2396642/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:2608405/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:281697/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:2312146/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:46417/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:541824/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:685652/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:744160/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:745460/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:943441/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:1307022/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:1381126/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:1037866/1‑68 (MQ=255) gCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTa > 1:1169739/1‑68 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:2967488/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:2249843/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:2184626/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1356649/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1978658/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:3230842/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1774536/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:3291675/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1174206/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1216874/69‑1 (MQ=255) atCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1317526/69‑1 (MQ=255) tCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:157106/68‑1 (MQ=255) tCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:2698535/68‑1 (MQ=255) tCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:2129016/68‑1 (MQ=255) tCACCCTGCGTGAGAGCGGAAAGGATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:3411396/68‑1 (MQ=255) tCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaaa > 1:130533/1‑69 (MQ=255) tGAGAGCGTAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2614311/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGTTTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:664904/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2790434/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2894585/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:1609699/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2400184/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2347846/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:1671231/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2231517/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:212678/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:1744453/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:143317/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:1797620/68‑1 (MQ=255) tGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACg < 1:2002716/68‑1 (MQ=255) agagCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGa > 1:1484849/1‑64 (MQ=255) ggATTTTACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:634822/46‑1 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGtctc > 1:432178/1‑40 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:848286/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:133997/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:606798/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:430840/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:683555/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:1710155/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:1758707/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:2041137/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:1277463/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:2978241/1‑69 (MQ=255) ttttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa > 1:2949520/1‑69 (MQ=255) ttttACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:1227219/69‑1 (MQ=255) ttttACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:2897298/69‑1 (MQ=255) ttttACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:3262950/69‑1 (MQ=255) ttttACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:408442/69‑1 (MQ=255) ttttACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:2013913/69‑1 (MQ=255) ttttACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:2311071/69‑1 (MQ=255) tttACGCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGaa < 1:1985973/42‑1 (MQ=255) cccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTaa < 1:679772/64‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATTAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:1005190/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACTTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:914300/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:1197609/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:543127/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:1199097/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:3434361/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:3391100/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:326667/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:2835778/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:2735904/68‑1 (MQ=255) ccTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTa < 1:158732/68‑1 (MQ=255) | TCTACAAACCGGGTGCCTGCTCTCTCGCTTCCCATATCACCCTGCGTGAGAGCGGAAAGGATTTTACCCTCGTCAGTGTGGATTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTA > minE/1042078‑1042212 |
| 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: ‑ |