Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,416,838 | 0 | T | G | 26.2% | 128.5 / 17.3 | 65 | intergenic (+1/‑4) | ECB_01334/ECB_01335 | putative Lysozyme/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (26/22); new base G (0/17); total (26/39) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.05e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCG > REL606/1416765‑1416913 | ccTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTc < 1:727237/76‑1 (MQ=255) cTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCa > 1:179696/1‑76 (MQ=255) tGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCaa < 1:1323756/76‑1 (MQ=255) gTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAAt > 2:1476712/1‑76 (MQ=255) gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCa > 1:508219/1‑76 (MQ=255) gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCa < 1:609870/76‑1 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCaa < 2:552960/76‑1 (MQ=255) aTCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaaa > 1:540069/1‑76 (MQ=255) tCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaaaa < 2:784046/76‑1 (MQ=255) cTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAAc > 2:941371/1‑76 (MQ=255) cTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAAc < 1:388823/76‑1 (MQ=255) tGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACt < 1:563459/76‑1 (MQ=255) aCAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTAt > 1:672888/1‑76 (MQ=255) cAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATg < 2:1807559/76‑1 (MQ=255) aaaTCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGa > 2:640645/1‑76 (MQ=255) aaaTCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGa < 2:653311/76‑1 (MQ=255) aaaTCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGa < 2:975759/76‑1 (MQ=255) tCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACt > 2:1381874/1‑76 (MQ=255) tCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACt < 1:806437/76‑1 (MQ=255) cTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTg > 1:1306325/1‑76 (MQ=255) cTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTg < 1:1213147/76‑1 (MQ=255) tttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGa > 2:838071/1‑76 (MQ=255) ttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAt > 1:143763/1‑76 (MQ=255) ttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAt > 1:599897/1‑76 (MQ=255) ttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAt > 1:1707021/1‑76 (MQ=255) ttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGAt < 2:92011/76‑1 (MQ=255) tATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAtt > 1:1353370/1‑76 (MQ=255) tttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgt < 1:272708/76‑1 (MQ=255) tttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgt < 1:358750/76‑1 (MQ=255) ttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgtg < 2:1545454/76‑1 (MQ=255) aGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAg < 1:1759637/76‑1 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt < 1:552415/76‑1 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt > 2:1275598/1‑76 (MQ=255) aagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa > 2:729097/1‑76 (MQ=255) aagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa < 2:35658/76‑1 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt > 2:357284/1‑76 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt > 2:194882/1‑76 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt > 2:1185909/1‑76 (MQ=255) aaGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtttt > 2:596336/1‑76 (MQ=255) gTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAAGTGATTGTGAAAGGGATCCGCAATTTTgaga < 1:583850/76‑1 (MQ=255) gTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgaga < 1:54263/76‑1 (MQ=255) aaGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAAt < 1:530811/76‑1 (MQ=255) aaGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAAt < 1:460230/76‑1 (MQ=255) ggCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAGTaa < 2:258250/76‑1 (MQ=255) gCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaaa < 2:1597213/76‑1 (MQ=255) gCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaaa < 1:252971/76‑1 (MQ=255) tGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGt > 1:1325803/1‑76 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt > 2:1070415/1‑76 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt < 1:1250990/76‑1 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt < 2:1279527/76‑1 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt > 2:1358398/1‑76 (MQ=255) gCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt < 1:1308566/76‑1 (MQ=255) gCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt < 1:388484/76‑1 (MQ=255) aCGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTAc < 1:1612688/76‑1 (MQ=255) ggggAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACgg < 2:31485/76‑1 (MQ=255) gggAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGt < 1:901668/76‑1 (MQ=255) ggAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTa < 2:1034448/76‑1 (MQ=255) ggAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTa < 1:1476712/76‑1 (MQ=255) gagaGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTaa < 1:237875/76‑1 (MQ=255) agTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACtt < 1:328382/76‑1 (MQ=255) gTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt < 1:1392203/76‑1 (MQ=255) gTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt > 2:1671196/1‑76 (MQ=255) tAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTa > 2:1204467/1‑76 (MQ=255) aaTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAg > 1:161752/1‑76 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc < 1:1791105/76‑1 (MQ=255) ttCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCg > 1:705417/1‑76 (MQ=255) | CCTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCG > REL606/1416765‑1416913 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |