Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,416,838 | 0 | T | G | 22.5% | 167.2 / 11.6 | 71 | intergenic (+1/‑4) | ECB_01334/ECB_01335 | putative Lysozyme/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (33/22); new base G (0/16); total (33/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.31e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGC > REL606/1416763‑1416912 | ggCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAt > 2:1045896/1‑76 (MQ=255) gCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAtt < 2:1488494/76‑1 (MQ=255) tGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCaa < 2:916739/76‑1 (MQ=255) gTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAAt < 1:1347892/76‑1 (MQ=255) cGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTc > 1:715621/1‑76 (MQ=255) cGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTc > 1:423126/1‑76 (MQ=255) gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCa > 1:1456813/1‑76 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa > 1:1301511/1‑76 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa > 1:1382834/1‑76 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa > 1:1301523/1‑76 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa > 1:1301514/1‑76 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCaa < 1:1129632/76‑1 (MQ=255) aTCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaaa > 2:785541/1‑76 (MQ=255) tCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaaaa < 1:179918/76‑1 (MQ=255) tCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCaaaa < 1:1210937/76‑1 (MQ=255) cAAATCTTTATTTCAAGATTAAAGAAGATACGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATg < 2:500947/76‑1 (MQ=255) aaaTCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGa < 2:294161/76‑1 (MQ=255) aaaTCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGa < 2:24705/76‑1 (MQ=255) aTCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAAc > 2:103621/1‑76 (MQ=255) tCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACt > 2:21584/1‑76 (MQ=255) cTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTg < 1:444704/76‑1 (MQ=255) tATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAtt > 2:945803/1‑76 (MQ=255) tCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGa < 1:1799805/76‑1 (MQ=255) aGATTAAAGAAGATAAGCGTAAGGCTGCACGGGTAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAg < 2:1535417/76‑1 (MQ=255) aGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAg < 1:252294/76‑1 (MQ=255) aTTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAggg > 2:1573624/1‑76 (MQ=255) aTTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAggg < 1:1332585/76‑1 (MQ=255) ttAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGa < 1:1821315/76‑1 (MQ=255) ttAAAGAAGAGAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGa < 2:1175391/76‑1 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt > 1:918788/1‑76 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt > 2:918460/1‑76 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt < 2:157191/76‑1 (MQ=255) aaaGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATc > 2:1627208/1‑76 (MQ=255) agaagaTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCg < 1:657198/76‑1 (MQ=255) aagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa > 1:331945/1‑76 (MQ=255) aagaTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa < 1:532339/76‑1 (MQ=255) agaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCaa > 1:110013/1‑76 (MQ=255) agaTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCaa < 1:1835556/76‑1 (MQ=255) gaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAt > 1:538445/1‑76 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt < 1:156289/76‑1 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt < 1:139987/76‑1 (MQ=255) tAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCCAAACTATGAAATGATTGTGAAAGGGATCCGCAAttt > 2:130924/1‑76 (MQ=255) tAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAttt > 2:1529678/1‑76 (MQ=255) aGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTg > 2:52588/1‑76 (MQ=255) cGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgag < 1:1213100/76‑1 (MQ=255) gTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgaga > 2:1370309/1‑76 (MQ=255) gTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgaga > 1:382907/1‑76 (MQ=255) tAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGaa > 2:963628/1‑76 (MQ=255) tAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGaa < 1:1486935/76‑1 (MQ=255) aGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATa < 1:1105784/76‑1 (MQ=255) ggCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaa < 1:1414875/76‑1 (MQ=255) ggCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaa < 1:1414885/76‑1 (MQ=255) gCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaaa < 1:1004536/76‑1 (MQ=255) cTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAg < 2:1005707/76‑1 (MQ=255) cTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAg < 1:711078/76‑1 (MQ=255) tGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGt > 1:1604420/1‑76 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa > 1:366309/1‑76 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa > 2:1064227/1‑76 (MQ=255) cACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa < 1:1390625/76‑1 (MQ=255) cgtggAGAGTAATTCAATGACTCAAAACTATTAACTTATTGTGAAAGGTATCCGCAATTTTGAGAATAAAGTTACg < 1:359349/73‑1 (MQ=17) gggAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGt > 2:961644/1‑76 (MQ=255) gagagaGTAAGTCAATGACGCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGt < 1:844951/74‑1 (MQ=18) ggAGAGTAATTTAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTa > 2:316040/1‑76 (MQ=255) ggAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTa > 1:1476462/1‑76 (MQ=255) gTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt < 1:1257515/76‑1 (MQ=255) gTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt < 2:316429/76‑1 (MQ=255) tAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTa > 2:511144/1‑76 (MQ=255) tAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTa < 1:1796418/76‑1 (MQ=255) aaTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAg < 1:460267/76‑1 (MQ=255) aaTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAg < 2:516845/76‑1 (MQ=255) aaTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAg < 1:1215664/76‑1 (MQ=255) aaTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAg < 1:103621/76‑1 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc > 1:1485875/1‑76 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc > 1:1273775/1‑76 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc < 1:1066072/76‑1 (MQ=255) | GGCCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGC > REL606/1416763‑1416912 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |