Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,195,874 | 0 | A | C | 22.7% | 107.3 / 13.6 | 44 | T318P (ACC→CCC) | aceB | malate synthase |
Reads supporting (aligned to +/- strand): ref base A (29/5); new base C (0/10); total (29/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.21e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAATA > REL606/4195799‑4195949 | cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTaaaa < 2:621855/76‑1 (MQ=255) cgcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTaaaa > 2:38963/1‑76 (MQ=255) gcgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTCTTGATTAAAAc > 2:72997/1‑76 (MQ=255) cgTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAAcc > 1:383441/1‑76 (MQ=255) gTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCt < 2:866089/76‑1 (MQ=255) tCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTg < 2:1079876/76‑1 (MQ=255) tCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTg < 2:1199232/76‑1 (MQ=255) tCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTg > 2:818626/1‑76 (MQ=255) ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc < 2:1610190/76‑1 (MQ=255) ccTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGc > 1:378621/1‑76 (MQ=255) tGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCa > 2:343161/1‑76 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt > 1:1174205/1‑76 (MQ=255) gCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCAt > 2:335074/1‑76 (MQ=255) cagacagaCAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaa > 1:545060/1‑76 (MQ=255) tgacagacaGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATaaa > 2:248530/2‑76 (MQ=255) gacagacaGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAc < 1:1051976/76‑1 (MQ=255) agacagGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAAcgcg > 1:1383081/1‑76 (MQ=255) agGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGc > 2:931903/1‑76 (MQ=255) gCAGGGAGGATGGATAACCCATTCCTGAATTCTCACCCCCGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCtt < 2:756256/76‑1 (MQ=255) aGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtttt > 2:659395/1‑76 (MQ=255) aGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCtttt > 1:1343822/1‑76 (MQ=255) tGACGATGGATAAACCATTCCTGAATGCTTCCTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGc < 1:557549/76‑1 (MQ=255) aTGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATgg > 2:251292/1‑76 (MQ=255) tGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATggg > 2:92410/1‑76 (MQ=255) ggATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggc > 1:1504985/1‑76 (MQ=255) gATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcg > 1:1089551/1‑76 (MQ=255) aTAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGggcgg > 1:1685225/1‑76 (MQ=255) tAAACCATCCCTGAATGCTTACTCCCGCCTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGggcggc < 1:635150/76‑1 (MQ=255) acccATTCCTAAATGCTTCCCCCCCCCTTTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCAt < 1:285638/74‑1 (MQ=255) aCCATTCCTGAATGTTTCCTCCCGCCTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATg < 1:1114760/76‑1 (MQ=255) tCCTGAATGCTTCCCCACGCCTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcggc < 1:854664/76‑1 (MQ=255) tCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATggcggc > 1:1131071/1‑76 (MQ=255) cTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGt > 1:1529996/1‑76 (MQ=255) gAATGCTTACTCACGCCTTTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGttt < 1:166056/76‑1 (MQ=255) aaTGCTTACTCACGCCTTTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTa < 2:128100/76‑1 (MQ=255) aTGCTTCCCCCCGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt < 1:122207/76‑1 (MQ=255) aTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt > 1:583304/1‑76 (MQ=255) aTGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAt > 1:1049228/1‑76 (MQ=255) tGCTTCCTCCCGCCTGTTAATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTAtt < 2:214427/76‑1 (MQ=255) gCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTc > 1:782734/1‑76 (MQ=255) gCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTc > 1:775884/1‑76 (MQ=255) tACTCACGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGa < 1:643284/76‑1 (MQ=255) ccccGCCTGTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCa < 1:1067398/73‑1 (MQ=255) ccGCCTGTTGATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCaaa < 1:161429/75‑1 (MQ=255) cGCCTTTTGATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAg < 1:92410/76‑1 (MQ=255) gCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGa > 1:1425734/1‑76 (MQ=255) ccTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt > 2:148263/1‑76 (MQ=255) ccTGTTAATAAACCCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGAt < 1:1272971/76‑1 (MQ=255) cTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATg > 1:230069/1‑76 (MQ=255) cTGTTAATAAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATg < 2:978521/76‑1 (MQ=255) gTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGaa > 1:139461/1‑76 (MQ=255) gTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGGTTATTCCGAGCAAAGATGaa > 1:502515/1‑76 (MQ=255) gATTAAACCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag < 2:1341079/76‑1 (MQ=255) tAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcac > 2:642772/1‑76 (MQ=255) aaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca < 1:1086847/76‑1 (MQ=255) aaacCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGcaca < 1:596734/76‑1 (MQ=255) aCCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAATa > 2:488686/1‑76 (MQ=255) | CGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAACGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAGCACAATA > REL606/4195799‑4195949 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |