| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 4,195,874 | 0 | A | C | 29.5% | 87.6 / 13.0 | 44 | T318P (ACC→CCC) | aceB | malate synthase |
| Reads supporting (aligned to +/- strand): ref base A (28/3); new base C (0/13); total (28/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT > REL606/4195799‑4195948 | cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTaaaa > 2:220615/1‑76 (MQ=255)cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTaaaa > 2:1099887/1‑76 (MQ=255)cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTaaaa > 1:1528690/1‑76 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAAcc > 1:263179/1‑76 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAAcc > 2:1352113/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCt > 1:1642427/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCt > 1:1742427/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCt > 1:1677351/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCt > 2:1050043/1‑76 (MQ=255) tCCTGCCAGACAGACAGGCAGTGACGATGGATAATCCATTCCTG‑AATGCTTACTCACGCCTGTTGATTCAAACCTg > 1:656760/1‑76 (MQ=255) ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGc < 1:1396614/76‑1 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt > 2:1148003/1‑76 (MQ=255) ccAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATa > 1:1247538/1‑76 (MQ=255) agacagacGGGCAGGGACGATGGATAAACCATCCCTG‑AATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATaaa < 1:739499/76‑1 (MQ=255) gacagacaGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAc > 2:1218657/1‑76 (MQ=255) gacagGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCgg > 2:386760/1‑76 (MQ=255) gacagGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCgg > 2:1277734/1‑76 (MQ=255) gCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtt > 1:313228/1‑76 (MQ=255) cAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCttt > 2:1205070/1‑76 (MQ=255) aCGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGa > 1:590768/1‑76 (MQ=255) aaTGGAT‑AATCATTTCTGAAATCCTTGCTCACCCTATTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATg < 2:504445/75‑1 (MQ=255) aTGGATAAACCATTCCTG‑AATGCTTACTCGCGCTTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATgg < 2:52938/76‑1 (MQ=255) ggATAAACCATTCCTG‑AATGCTTCCTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggc < 1:82791/76‑1 (MQ=255) gATAACCCATCCCTG‑AATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggcg < 2:365166/76‑1 (MQ=255) aTAACCCATTCCTG‑AATGCTTACCCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggcgg < 1:176547/76‑1 (MQ=255) aTAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 1:1468164/1‑76 (MQ=255) aTAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 1:1374318/1‑76 (MQ=255) aTAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 2:1166566/1‑76 (MQ=255) aaaaCCATGCCGG‑AATGCTCCCCCCCGCTTTTTAATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcggc < 2:1625949/75‑1 (MQ=255) aaaCCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa > 1:91620/1‑76 (MQ=255) aCCATTCCTA‑AATGCTTCCCCACGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg < 1:1298088/76‑1 (MQ=255) cATTCCTG‑AATGCTTCCTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggc < 2:728018/76‑1 (MQ=255) ttCCTG‑AATGCTTCCCCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcgg < 2:171307/76‑1 (MQ=255) tCCTG‑AATGCTTCCCCCCGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcggc < 2:1544426/76‑1 (MQ=255) ccTG‑AATGCTTACTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcg < 2:95025/76‑1 (MQ=255) tG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcggt > 2:728438/1‑74 (MQ=255) tG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGtt > 1:1842397/1‑76 (MQ=255) g‑AATGCTTACTCGCGCGTGTTGATAAAACCCTGCCATAAACGCGGTGCTTGTGCGATGGGCGGCATGGCGGCGttt < 1:217858/76‑1 (MQ=255) aaTGCTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa < 2:201880/76‑1 (MQ=255) aaTGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa < 1:1050043/76‑1 (MQ=255) tcTTACCCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTc < 1:386760/75‑1 (MQ=255) ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg > 2:182890/1‑76 (MQ=255) ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg > 2:1408143/1‑76 (MQ=255) cTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGc < 1:39044/76‑1 (MQ=255) tCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa > 1:1392440/1‑76 (MQ=255) ccACGCCTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa < 1:798423/75‑1 (MQ=255) gCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGa > 1:502060/1‑76 (MQ=255) tGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa > 2:1275942/1‑76 (MQ=255) gattGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa < 1:641730/74‑1 (MQ=255) gTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGaa < 1:1819669/76‑1 (MQ=255) tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac > 2:1626898/1‑76 (MQ=255) tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac > 1:1075740/1‑76 (MQ=255) aaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACaa > 2:263227/1‑76 (MQ=255) aaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAt > 2:700728/1‑76 (MQ=255) | CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTG‑AATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAT > REL606/4195799‑4195948 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |