Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 5,988,809 | G→A | 46.6% | intergenic (+103/‑443) | 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,809 | 0 | G | A | 46.6% | 4.3 / 42.9 | 30 | intergenic (+103/‑443) | kefB‑III/PP_mr64 | glutathione‑regulated potassium/H+ antiporter/C4 ncRNA |
Reads supporting (aligned to +/- strand): ref base G (3/13); new base A (4/10); total (7/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.75e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01 |
TTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAG > AE015451/5988579‑5988924 | ttACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACTGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGtt < 2:507242/233‑1 (MQ=255) tcgcggtgtgtatagggcaatgggcggaattcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTcctgcgggatttgcgcg < 2:286225‑M2/67‑9 (MQ=255) cgcggtgtgtatagggcaatgggcggaattcgcgctttggcgggtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATtcctgcgggatttgcgcgc < 1:428360‑M2/68‑10 (MQ=255) aGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGgcg > 1:2265/1‑234 (MQ=255) atgaCCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGc < 1:247778/210‑1 (MQ=255) gtaggtgtttggcgggaggttttcagcgcctgggaggtgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTc < 2:691375‑M2/117‑1 (MQ=255) aGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑c > 1:552725/1‑177 (MQ=255) gCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc < 1:778127/235‑1 (MQ=255) tgttcgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc < 1:659358‑M2/123‑1 (MQ=255) ttGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGa < 1:643403/227‑1 (MQ=255) cgccgggggattttcagcgcctgggaggtgttcgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCg > 1:44919‑M2/77‑190 (MQ=255) tCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATcg < 1:578156/169‑1 (MQ=255) tCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGcc < 2:2265/199‑1 (MQ=255) ggACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTg < 2:334213/143‑1 (MQ=255) gcgcGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTcg < 1:87766/167‑1 (MQ=255) cgccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcgcgcgc < 2:44919‑M2/123‑1 (MQ=255) gccgggagagtttcagcgcctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCgg < 2:674072‑M2/115‑1 (MQ=255) cctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTAt < 2:55428‑M2/49‑1 (MQ=255) cctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTAt > 1:730358‑M2/26‑74 (MQ=255) cctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTAt < 2:730356‑M2/49‑1 (MQ=255) cctgggagatcgagcgccgcccgcgCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTAt > 1:55428‑M2/26‑74 (MQ=255) gcgccgccctcgCGGCGCTTCGCGGGTAACCCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTa < 1:16011‑M2/159‑1 (MQ=255) gcgcgGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCgg > 2:346999/1‑118 (MQ=14) gcgcgGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGA‑‑TT‑TGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCgg < 1:347001/118‑1 (MQ=9) ggcgCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGTCGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTcg < 1:193627/120‑1 (MQ=12) ttCGCGGGTAAACCCGCTCCTCCGTTTGTTTCGGGCCAGTTATTCCCGGGTGA‑‑CG‑CGCACGCGAGCCCCCTGGCGGACGCCTCGATATCGGGTCGTACAAACACGGCGGTcgcgcgc < 2:94761/117‑1 (MQ=25) aaCCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGGCGTACAACCAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaa < 2:552723/138‑1 (MQ=255) cTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg < 1:803225/138‑1 (MQ=255) tACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAggg < 2:524761/113‑1 (MQ=255) tCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaa > 2:657864/1‑114 (MQ=255) | TTACCCAGCATGACGAAGAGGTGCTGGAGGCCCAGGGGCAAGTGCGTGATGACCGGGCCAAGGTGATGCAGACGGCCAAGGAGGCGCGGGTGGAGCTGGAGCGGTTGTTTGATTCGGATGCGGACTGAGGCCTTTTGGCGCGTGCCTTGAGGGTTGTAGTGCCTGTGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGA‑‑CG‑CGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAG > AE015451/5988579‑5988924 |
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: ‑ |