Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,195,874 | 0 | A | C | 29.4% | 109.2 / 11.8 | 52 | T318P (ACC→CCC) | aceB | malate synthase |
Reads supporting (aligned to +/- strand): ref base A (26/10); new base C (1/14); total (27/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-05 | |||||||||||
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:630687/76‑1 (MQ=255) ggTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc < 2:1202856/75‑1 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc < 2:1202866/76‑1 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt > 2:287359/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt < 1:556890/76‑1 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt < 1:872312/76‑1 (MQ=255) tCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTg < 2:652482/76‑1 (MQ=255) ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc > 2:920841/1‑76 (MQ=255) tGCCAGATAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa < 2:420181/76‑1 (MQ=255) tGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa > 1:120258/1‑76 (MQ=255) gccagccaGACAGGCAGTGACGATGGATAACCCATTCCTGAATGCTTTCCCACGCCTGTTAATTAAAGCCTGCCAt < 2:243526/76‑1 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt < 2:579922/76‑1 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt > 1:1154118/1‑76 (MQ=255) cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa > 2:363104/1‑76 (MQ=255) cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa < 1:153291/76‑1 (MQ=255) agacagacAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCGNTTGATTAAAACCTACCATaac > 2:270463/1‑75 (MQ=255) gacagacaGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAAc < 1:989408/76‑1 (MQ=255) cagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTg < 2:597175/76‑1 (MQ=255) cagGCAGTGACGATGGAAAAACCATCCCTGAATGCTTACTCACGCCTGTTAATTAAAACCTGCCATAAACGCGGTg < 2:810559/76‑1 (MQ=255) agGCAGTGACGATGGATAAACCATACCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGc > 2:1059390/1‑76 (MQ=255) gCAGTGACGATGGATAAACCATTCCTGAATGCTTANNNACGCCTGTTGATTAAAANCTGCCATAAACGCGGTGCtt > 1:843379/1‑76 (MQ=255) gCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTTTTGATTAAAACCTGCCATAAACGCGGTGCtt < 1:1559927/76‑1 (MQ=255) gCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtt > 1:1164655/1‑76 (MQ=255) gACGATGGATAAACCATTCCTGAATGCTTACCCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCg < 1:617145/76‑1 (MQ=255) cGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGAt > 1:1200911/1‑76 (MQ=255) aTGGATAAACCATTCCTGAATGCTTCCTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATgg < 2:1102012/76‑1 (MQ=255) aTGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAACACTGCCATACACGCGGTGCTTTTGCGACgg > 2:394907/1‑76 (MQ=255) aTGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATgg > 2:715764/1‑76 (MQ=255) tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg > 1:1039567/1‑76 (MQ=255) tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg > 1:147436/1‑76 (MQ=255) tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg > 1:1651774/1‑76 (MQ=255) ggATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc > 2:1479459/1‑76 (MQ=255) ggATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc > 1:267376/1‑76 (MQ=255) aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 1:1207377/1‑76 (MQ=255) aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 1:318126/1‑76 (MQ=255) aaaCCATCCCTGAATGCTTACCCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCa < 1:586976/76‑1 (MQ=255) aaCCATTCCTAAATGCTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCAt < 1:264905/76‑1 (MQ=255) cATCCCTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggc < 2:957404/76‑1 (MQ=255) ccTGAATGCTTACTCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGgcggcg < 1:154307/76‑1 (MQ=255) cTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGt > 1:1128885/1‑76 (MQ=255) tGAATGCTTCCCCACGCCTGTTAATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGtt < 2:82264/76‑1 (MQ=255) gAATGCTTCCTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt < 1:1507467/76‑1 (MQ=255) gAATGCTTACTCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt < 1:515969/76‑1 (MQ=255) gAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt > 1:940771/1‑76 (MQ=255) aTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt > 1:255470/1‑76 (MQ=255) tGCTTCCCCCCGCCTTTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAtt < 2:1267867/76‑1 (MQ=255) ttACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCg > 2:594484/1‑76 (MQ=255) tCACGCATTTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa < 1:1369850/76‑1 (MQ=255) ccACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa < 1:699750/75‑1 (MQ=255) cgcGCCTGTTGATAAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaa < 1:1500226/74‑1 (MQ=255) cGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAg > 1:1182502/1‑76 (MQ=255) cttGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt < 2:1049775/74‑1 (MQ=255) tGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATgn < 1:463590/76‑2 (MQ=255) tGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGa > 1:151703/1‑76 (MQ=255) ttGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag < 1:1386653/76‑1 (MQ=255) ttGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag < 1:1096885/76‑1 (MQ=255) ttGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag > 1:10103/1‑76 (MQ=255) ttGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGTGGCATGGCGGCGTTTATTCCGAGCAAAGATGAag < 1:928609/76‑1 (MQ=255) gATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag < 1:842989/76‑1 (MQ=255) tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac > 1:1464723/1‑76 (MQ=255) aaaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac < 1:715764/75‑1 (MQ=255) aaaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca > 2:1552211/1‑76 (MQ=255) aaaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca < 1:1520167/76‑1 (MQ=255) aaaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACaa > 2:6277/1‑76 (MQ=255) aaCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAAt > 1:1301386/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 < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |