Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | REL606 | 239,002 | A→C | 100% | A4A (GCA→GCC) | dnaQ → | DNA polymerase III subunit epsilon |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 239,002 | 0 | A | C | 92.8% | 138.0 / 3.6 | 42 | A4A (GCA→GCC) | dnaQ | DNA polymerase III subunit epsilon |
Reads supporting (aligned to +/- strand): ref base A (2/1); new base C (18/21); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.98e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCGTCGCCTGACGGGCAATAACTTCCATGTTTATCTCA > REL606/238858‑239147 | gcCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTac < 2:175407/150‑1 (MQ=255) ccGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTaca < 1:565504/150‑1 (MQ=255) cccAGGTCTTGGATTGGCCAGACACGAACCGTCGGTGAAAATTTTTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCTAACTGTAACACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACCCCAGAt > 1:262682/1‑150 (MQ=255) ccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATc < 2:246220/150‑1 (MQ=255) ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCCTTACACGCCAGATCGtt > 1:440341/1‑150 (MQ=255) ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGtt > 1:27455/1‑150 (MQ=255) ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGtt > 2:577312/1‑150 (MQ=255) ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGtt > 1:441300/1‑150 (MQ=255) gCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGGTAAGTCTGGCATAAATGCCCGCTATGAGCACTGCCATTACACGCCAGATCGTTCCCGATACCgg > 1:653266/1‑149 (MQ=255) gCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGa > 2:594078/1‑150 (MQ=255) gCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGa > 2:401884/1‑150 (MQ=255) gACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAacca > 1:174982/1‑150 (MQ=255) gACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATCACACGCCAAATCGTTTTCGATACCGAAacca > 1:398190/1‑150 (MQ=255) acacGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGCACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTGCTCGATACCGAAaccac > 1:198383/1‑150 (MQ=255) aaCCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa < 2:468790/150‑1 (MQ=255) tCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc < 1:592335/150‑1 (MQ=255) cGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACAGGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc < 1:258031/150‑1 (MQ=255) gAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt < 2:263431/150‑1 (MQ=255) aaaaTTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt < 2:511765/150‑1 (MQ=255) aaTTTCTACCTGTTTAAGCATCTCTGGTAGACGTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg < 1:574973/150‑1 (MQ=255) gTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg > 1:571412/1‑150 (MQ=255) aaGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg < 2:263800/150‑1 (MQ=255) aaGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg < 2:336350/150‑1 (MQ=255) aaGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg < 2:175481/150‑1 (MQ=255) aaGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg < 2:573081/150‑1 (MQ=255) tctGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg < 2:27538/150‑1 (MQ=255) ctGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa > 1:456129/1‑150 (MQ=255) ctGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa > 2:571714/1‑150 (MQ=255) ctGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa < 2:258100/150‑1 (MQ=255) ctGGTAGACTTCCAGTAATTGAATCGAACTGTAAAACGCTAAGTCTGACATAAATGACCCCTATGACCACTCCCATTACCCGCCAGATCGTTCTCGATACCGAAACCCCCCGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAgc < 1:620436/150‑2 (MQ=255) gTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCCCTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGGTTGGt > 1:314755/1‑150 (MQ=255) cTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAgtggt > 2:280823/1‑150 (MQ=255) gTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGa > 1:410990/1‑150 (MQ=255) aaaaCGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAAc < 2:441588/150‑1 (MQ=255) aaaaCGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAAc < 2:269109/150‑1 (MQ=255) aaCGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCg < 2:442550/150‑1 (MQ=255) aaCGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCg < 2:550356/150‑1 (MQ=255) gTCTGACATAAATGACCGCTATGAGCACTGCCATTACAAGCCAGAACGTTCTCGATACCGAAACCAACGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGATGCCGTTGTAGTGATTAACCGTCGCGTGa > 1:324161/1‑150 (MQ=255) aCATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCGTCGCCTGACGGGc > 2:607437/1‑150 (MQ=255) aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCGTCGCCTGACGGGCAATAAc < 1:280023/150‑1 (MQ=255) aCTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCGTCGCCTGACGGGCAATAACTTCCATGTTTAtctc < 2:396888/150‑1 (MQ=255) cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCGTCGCCTGACGGACAATAACTTCCATGTTTATCTCa > 2:503590/1‑150 (MQ=255) | GCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGAAGTGGTGAACCGTCGCCTGACGGGCAATAACTTCCATGTTTATCTCA > REL606/238858‑239147 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |