| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,857,716 | C→G | 29.1% | E27Q (GAA→CAA) | cysH ← | 3'‑phosphoadenosine 5'‑phosphosulfate reductase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,857,716 | 0 | C | G | 29.1% | 101.0 / 48.6 | 79 | E27Q (GAA→CAA) | cysH | 3'‑phosphoadenosine 5'‑phosphosulfate reductase |
| Reads supporting (aligned to +/- strand): ref base C (35/21); new base G (12/11); total (47/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.54e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01 | |||||||||||
ACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCG > minE/1857649‑1857774 | aCATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCg > 1:3123947/1‑69 (MQ=255)aCATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCg > 1:2484376/1‑69 (MQ=255) atTCACCGGGCAGTTTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:1025312/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCTGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:2353284/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:1527963/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:1795892/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:1889294/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:2049122/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:113925/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:2525620/69‑1 (MQ=255) atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:843859/69‑1 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTGGGCGt < 1:2224393/68‑1 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:2289534/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:265941/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:37392/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:2856797/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:2567707/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:145066/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:1127491/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:1290357/1‑69 (MQ=255) ttCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGtt > 1:1561430/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:1632841/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:2965364/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:737889/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:3068065/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:1824385/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:337245/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:2212220/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:407680/1‑69 (MQ=255) cACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAg > 1:2314132/1‑69 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:3311122/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:3127893/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:2349952/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:1202903/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:3060980/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:1398953/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:660157/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:578889/69‑1 (MQ=255) cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCcgccagcgccag < 1:588846/69‑1 (MQ=255) ccAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCg < 1:510731/69‑1 (MQ=255) ccAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCg < 1:1906247/69‑1 (MQ=255) ccAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCg < 1:1258612/69‑1 (MQ=255) gCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:1274148/69‑1 (MQ=255) gCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:1118871/69‑1 (MQ=255) cTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:2581264/68‑1 (MQ=255) cTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:2604295/68‑1 (MQ=255) cTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:2705967/68‑1 (MQ=255) cTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:2712962/68‑1 (MQ=255) cTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:573379/68‑1 (MQ=255) cTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:2994847/68‑1 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCTCCAGAATCCGGTCTACCttt > 1:1060952/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:978959/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:3391108/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:628852/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:3286394/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:1004973/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:2544470/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:1493338/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:1687885/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:1780323/1‑65 (MQ=255) ccTTCAGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCtt > 1:765543/1‑65 (MQ=255) ccTTCAGCGTCCAGCTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATCCGGTCTACCttt > 1:206733/1‑65 (MQ=255) tCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:2834175/57‑1 (MQ=255) aGCGTCCAGTTTTTCCAGTTGGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCt < 1:3099248/55‑1 (MQ=255) cGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcg > 1:902312/1‑68 (MQ=255) cGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcg > 1:3041478/1‑68 (MQ=255) cGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcg > 1:2794481/1‑68 (MQ=255) cGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAgttcg > 1:3425607/1‑68 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCa > 1:3099234/1‑68 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCa > 1:1073496/1‑68 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCa > 1:345721/1‑68 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAg > 1:1698183/1‑69 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAg > 1:1795447/1‑69 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAg > 1:453488/1‑69 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAg > 1:770155/1‑69 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAg > 1:3097188/1‑69 (MQ=255) cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAg > 1:1007734/1‑69 (MQ=255) ttttCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAggggg > 1:1080028/1‑67 (MQ=255) ttttCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCg > 1:2338160/1‑69 (MQ=255) | ACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCG > minE/1857649‑1857774 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |