Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 5,988,894 | C→T | 32.1% | intergenic (+188/‑358) | kefB‑III → / → PP_mr64 | glutathione‑regulated potassium/H+ antiporter/C4 ncRNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 5,988,894 | 0 | C | T | 32.1% | 33.8 / 23.7 | 28 | intergenic (+188/‑358) | kefB‑III/PP_mr64 | glutathione‑regulated potassium/H+ antiporter/C4 ncRNA |
Reads supporting (aligned to +/- strand): ref base C (13/6); new base T (2/7); total (15/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.18e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01 |
GCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGA > AE015451/5988662‑5989106 | gcgcGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGcac < 1:585045/235‑1 (MQ=255) tGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTa < 1:9486/236‑1 (MQ=255) agtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTggccga < 1:698372‑M2/148‑3 (MQ=255) tttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCa > 1:164484/1‑236 (MQ=255) tgtgtCGGGCCAATTATGCCTGCGGGAGTTGCGCGCGGCCGCCGTGGGGGATGTGTCGATATCGCGTCGTACAAACAAGGCGGGCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTa < 2:36938/126‑1 (MQ=255) tCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCCCAGGAGGGATTGGCCAGAAACACACGTACGAGCgg > 2:27605/1‑130 (MQ=255) cAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAgg > 1:573399/1‑119 (MQ=255) aGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCAc > 2:460207/1‑224 (MQ=255) gACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTTGccc > 2:345457/1‑89 (MQ=255) gcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAACCGTAGGAGCGGGTTTAcccgc < 2:59830‑M2/114‑6 (MQ=255) cgcACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGaa > 1:691859/1‑207 (MQ=255) cgcAAGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCa < 1:422677/187‑1 (MQ=255) ccTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCCGAAACAAACGTAGGAGCGGGGTTAcccgcgaagcgccgcgcgggcgg < 2:173370‑M2/109‑24 (MQ=255) cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGATCGGGTTTAcccgcgaagcgacgcgcgggcggcgct < 2:141508‑M2/112‑28 (MQ=255) cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattc < 1:371951‑M2/157‑74 (MQ=255) cGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTaa > 2:669176/1‑226 (MQ=255) tatCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGt > 1:173306/1‑236 (MQ=255) atCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTg > 2:422676/1‑167 (MQ=255) cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCa < 1:591811/221‑1 (MQ=255) cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCa > 2:591810/1‑221 (MQ=255) cgtACAAACAAGGCGGTCGCGCGCGACTGGTGCaggag < 1:172084/38‑1 (MQ=14) cgtACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcccgcgcggcgc > 2:445849‑M2/1‑71 (MQ=255) cgtACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcccgcgcggcgc < 1:445850‑M2/97‑27 (MQ=255) acaaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaat < 1:679242‑M2/161‑94 (MQ=255) acaaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaat > 2:679241‑M2/1‑68 (MQ=255) aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAACACCATCACAAATc > 2:284714/1‑204 (MQ=255) aaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaac < 1:662464‑M2/214‑149 (MQ=255) gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGTCGGCGCTCGATCTCATGTGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCctcgc > 1:185560/1‑227 (MQ=255) gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGAcg > 2:583767/1‑168 (MQ=255) | GCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGCGAACTCACAAAAACCATCACAAATCGTCTAATTCTCAGCTGGGTCCTACTTCCCCGCGA > AE015451/5988662‑5989106 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |