Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,195,874 | 0 | A | C | 28.4% | 153.8 / 12.5 | 67 | T318P (ACC→CCC) | aceB | malate synthase |
Reads supporting (aligned to +/- strand): ref base A (35/13); new base C (0/19); total (35/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.44e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT > REL606/4195799‑4195948 | cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTaaaa < 1:1167002/76‑1 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc > 1:411449/1‑76 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc > 1:1244035/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt < 1:1345132/76‑1 (MQ=255) ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc > 1:505662/1‑76 (MQ=255) cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc > 2:1568895/1‑76 (MQ=255) cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc > 1:1113311/1‑76 (MQ=255) cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc < 1:1687738/76‑1 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCCCGCCTGTTGATTAAAACCTGCCAt > 2:834205/1‑76 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt > 1:4581/1‑76 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt > 2:286872/1‑76 (MQ=255) ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTTATTAAAACCTGCCATa < 2:1378556/76‑1 (MQ=255) ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATa < 1:459872/76‑1 (MQ=255) ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATa > 1:817903/1‑76 (MQ=255) cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa > 1:1162012/1‑76 (MQ=255) cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa > 1:847564/1‑76 (MQ=255) acagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcg < 2:1071380/76‑1 (MQ=255) acagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcg < 1:1683752/76‑1 (MQ=255) cagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcgc < 1:1345071/76‑1 (MQ=255) agacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcgcg > 1:146102/1‑76 (MQ=255) gacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCCTCAACGCgg > 1:1533882/1‑76 (MQ=255) gacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCgg > 1:705405/1‑76 (MQ=255) acagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGt > 2:1681217/1‑76 (MQ=255) aGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCtttt < 1:121602/76‑1 (MQ=255) gTGACGATGGAAAACCCATCCCTGAATGCTTCCTCCCGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTg < 1:1047666/76‑1 (MQ=255) aCGATGGATAAACCATTCCTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGa < 1:892576/76‑1 (MQ=255) aCGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGa > 1:351002/1‑76 (MQ=255) cGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGAt > 1:972058/1‑76 (MQ=255) aTGGATACCCCATCCCTTAATTCTTACTCACGCCTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATgg < 2:668617/76‑1 (MQ=255) ggATAATCCATTCTGGATTGTTTACTCCCGTTTGTGGTTTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc < 2:1190439/76‑1 (MQ=255) ggATAAACCATTCCTAAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggc < 1:1244221/76‑1 (MQ=255) gATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggcn < 1:1474850/76‑2 (MQ=255) aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 2:187534/1‑76 (MQ=255) aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg < 2:650662/76‑1 (MQ=255) aatccATTCCTGATTGCTTACCCCCCCCTTTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa < 1:846019/73‑1 (MQ=255) aaaCCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa < 2:1642823/76‑1 (MQ=255) cccATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg < 1:1569453/75‑1 (MQ=255) aCCATTCCTGAATGCTTCCCCCCGCTTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg < 2:790577/76‑1 (MQ=255) aCCATTCCTGAATGCTTACCCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg < 2:1480125/76‑1 (MQ=255) aTTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcg > 2:1130179/1‑76 (MQ=255) ccTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcg > 1:1118502/1‑76 (MQ=255) cTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGt > 1:642796/1‑76 (MQ=255) tGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGtt < 1:1012805/76‑1 (MQ=255) aaTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa > 1:551756/1‑76 (MQ=255) aaTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa > 2:245153/1‑76 (MQ=255) aTGCTTCCTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt < 1:979984/76‑1 (MQ=255) aTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt > 1:80396/1‑76 (MQ=255) gCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTc > 1:1612259/1‑76 (MQ=255) cTTCCTCCCGCCGGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc < 1:1148108/76‑1 (MQ=255) cTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc < 1:1561521/76‑1 (MQ=255) cTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc > 1:1111890/1‑76 (MQ=255) ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg > 2:681988/1‑76 (MQ=255) ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg > 1:1434645/1‑76 (MQ=255) ttACCCACGCCTGTTAATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg < 1:593019/76‑1 (MQ=255) tACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGa < 1:2047/76‑1 (MQ=255) tACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGa > 1:1133340/1‑76 (MQ=255) ccTCACGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAg < 1:677842/75‑1 (MQ=255) cTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 2:1627181/76‑1 (MQ=255) cTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 1:153530/76‑1 (MQ=255) ccACGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa < 1:289269/75‑1 (MQ=255) cccGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa < 1:1221044/74‑1 (MQ=255) cACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa < 2:1073452/76‑1 (MQ=255) ccGCCTTTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaaa < 1:334546/75‑1 (MQ=255) cGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAg < 2:1423959/76‑1 (MQ=255) ccTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt < 1:1038300/76‑1 (MQ=255) tGTTGATTAAGCCCTGCCATAAACTCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa < 1:120475/76‑1 (MQ=255) tGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa > 1:1247836/1‑76 (MQ=255) gTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGaa < 1:1659688/76‑1 (MQ=255) tGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAaga > 2:206174/1‑76 (MQ=255) tGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAaga > 1:698791/1‑76 (MQ=255) gATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag < 2:218579/76‑1 (MQ=255) gATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag > 2:27627/1‑76 (MQ=255) gATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag > 1:1219953/1‑76 (MQ=255) ataaaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGc < 1:380754/73‑1 (MQ=255) tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac > 1:1151306/1‑76 (MQ=255) aaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAt > 1:1624192/1‑76 (MQ=255) | CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT > REL606/4195799‑4195948 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |