Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,195,874 | 0 | A | C | 42.2% | 68.2 / 19.6 | 45 | T318P (ACC→CCC) | aceB | malate synthase |
Reads supporting (aligned to +/- strand): ref base A (19/7); new base C (0/19); total (19/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.93e-07 | |||||||||||
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 > 2:531442/1‑76 (MQ=255) gcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAc < 2:188057/76‑1 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc < 1:1495552/76‑1 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc > 2:1358156/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCACTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt > 2:39879/1‑76 (MQ=255) ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc < 2:1180440/76‑1 (MQ=255) cTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGcc < 1:529744/76‑1 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt < 1:1386557/76‑1 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt < 1:1923809/76‑1 (MQ=255) ccAGACAGACAGGCGGGAAGGATGGATAAACCATTCCTGAAGGCTTACTCACGCCTGTTGATTAAACCCTGCCATa < 2:751547/76‑1 (MQ=255) cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa > 2:127047/1‑76 (MQ=255) acagacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAAcg < 2:1172325/76‑1 (MQ=255) gacagGCAGTGACGATGGATAACCCATTCCTGAATTCTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCgg < 1:1729506/76‑1 (MQ=255) ccagGCAGTGACGATGGATAAACCATTCCTGAATCCTTCCTCGCGCCTGTTGATAAAAACCTGCCATAAACGCGGt < 2:1361220/75‑1 (MQ=255) acagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGt > 1:1150692/1‑76 (MQ=255) ggCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCt > 1:1004257/1‑76 (MQ=255) cAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCttt < 1:1388910/76‑1 (MQ=255) gTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTg > 2:871602/1‑76 (MQ=255) gACGATGGATAACCCATCCCTGAATGCTTCCTCACGCCTGTTAATTAAAACCTGCCATAAACGCGGTGCTTTTGCg < 1:1494567/76‑1 (MQ=255) aCGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGa > 2:1143887/1‑76 (MQ=255) cGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGAt > 2:77206/1‑76 (MQ=255) cGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGAt > 2:951847/1‑76 (MQ=255) ggATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc > 1:1135086/1‑76 (MQ=255) acccATTCCTAAATGCTTCCTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCAt < 1:1071270/74‑1 (MQ=255) cATTCCTGAATGCTTCCCCACGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggc < 1:139523/76‑1 (MQ=255) aTTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcg > 1:2018900/1‑76 (MQ=255) ttCCTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcgg < 2:1979491/76‑1 (MQ=255) tCCTGAATGCTTACTCCCTCCTGTTGATGAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcggc < 1:1862893/76‑1 (MQ=255) ccTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcg > 1:1967440/1‑76 (MQ=255) tGAATGCTTCCCCCCGCCTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGtt < 2:179793/76‑1 (MQ=255) gAATGCTTCCTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt < 1:542718/76‑1 (MQ=255) aaTGCTTACTCCCCCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa < 2:735612/76‑1 (MQ=255) aaTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa > 1:1497945/1‑76 (MQ=255) tGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAtt > 1:1794584/1‑76 (MQ=255) tGCTTACCCACGCTTGTGGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAtt < 1:1692837/76‑1 (MQ=255) gCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTc > 1:1919726/1‑76 (MQ=255) cTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTcc < 2:340987/76‑1 (MQ=255) ttccTCGCGCTTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg < 1:1348384/73‑1 (MQ=255) ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg > 1:1531144/1‑76 (MQ=255) cccACGCCTTTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 2:1952019/74‑1 (MQ=255) cTCGCGCTTTTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 1:1241022/76‑1 (MQ=255) cTCCCGCGTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 1:1076257/76‑1 (MQ=255) cTCCCGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 1:1114318/76‑1 (MQ=255) cTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc > 1:2083994/1‑76 (MQ=255) cccGCCTGTCTATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa < 1:738556/74‑1 (MQ=255) cACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa > 1:1640778/1‑76 (MQ=255) cGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAg < 1:376111/76‑1 (MQ=255) cccTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGa < 1:220461/75‑1 (MQ=255) ccTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt < 2:1967440/76‑1 (MQ=255) ccTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt < 2:1053565/76‑1 (MQ=255) tGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAaga > 1:1288032/1‑76 (MQ=255) taaaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGca < 2:2013287/74‑1 (MQ=255) aaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca < 2:1359234/76‑1 (MQ=255) aaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAt > 2:974334/1‑76 (MQ=255) | CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT > REL606/4195799‑4195948 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |