Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,416,838 | 0 | T | G | 27.7% | 171.5 / 21.9 | 83 | intergenic (+1/‑4) | ECB_01334/ECB_01335 | putative Lysozyme/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (41/19); new base G (0/23); total (41/42) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.54e-09 | |||||||||||
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 | ccTGTTGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTc > 1:379711/1‑76 (MQ=255) ttGACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATg > 2:1812333/1‑76 (MQ=255) gACGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGAc > 1:92452/1‑76 (MQ=255) cGTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTc > 1:1784309/1‑76 (MQ=255) gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCa > 1:1828525/1‑76 (MQ=255) gTATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCa < 1:442364/76‑1 (MQ=255) tATCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaa < 1:921992/76‑1 (MQ=255) aTCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCaaa < 2:563734/76‑1 (MQ=255) tCTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCaaaa > 1:147262/1‑76 (MQ=255) cTGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAAc < 1:261542/76‑1 (MQ=255) tGACAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACt > 1:944069/1‑76 (MQ=255) aCAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTAt < 1:1798374/76‑1 (MQ=255) aCAAATCTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTAt < 1:559747/76‑1 (MQ=255) cTTTATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTg < 1:402066/76‑1 (MQ=255) tttATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGa < 2:732530/76‑1 (MQ=255) tATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGAtt > 2:997314/1‑76 (MQ=255) tATTTCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGAtt < 1:935838/76‑1 (MQ=255) tttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATtgt < 1:362699/76‑1 (MQ=255) tttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGCAATTCAATGACTCAAAACTATGAACTGATtgt > 1:39377/1‑76 (MQ=255) ttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATtgtg < 2:497790/76‑1 (MQ=255) ttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATtgtg > 2:20269/1‑76 (MQ=255) ttCAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATtgtg < 1:794411/76‑1 (MQ=255) cAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGaa > 1:478853/1‑76 (MQ=255) cAAGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGaa > 1:358562/1‑76 (MQ=255) aaGATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGaaa > 2:1785534/1‑76 (MQ=255) gATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAgg > 1:787672/1‑76 (MQ=255) gATTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAgg > 2:1613037/1‑76 (MQ=255) aTTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAggg > 1:645029/1‑76 (MQ=255) aTTAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAggg < 2:1193784/76‑1 (MQ=255) ttAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGa < 1:821713/76‑1 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt < 1:1050943/76‑1 (MQ=255) tAAAGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGAt < 1:1507722/76‑1 (MQ=255) aaaGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATc > 2:228791/1‑76 (MQ=255) aaaGAAGATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATc > 2:504187/1‑76 (MQ=255) aagaagATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATcc > 2:706291/1‑76 (MQ=255) aagaagATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATcc > 2:1290089/1‑76 (MQ=255) aagaagATAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATcc > 2:1604271/1‑76 (MQ=255) agaagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCg > 1:826991/1‑76 (MQ=255) aagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa > 2:222308/1‑76 (MQ=255) aagaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCa > 1:297787/1‑76 (MQ=255) agaTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCaa > 2:204673/1‑76 (MQ=255) agaTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCaa < 1:574455/76‑1 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt > 1:454486/1‑76 (MQ=255) aTAAGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtt < 1:379382/76‑1 (MQ=255) tAAGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAttt > 1:1823179/1‑76 (MQ=255) aaGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtttt > 2:520457/1‑76 (MQ=255) aaGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtttt > 2:1664303/1‑76 (MQ=255) aaGCGTAAGGCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtttt < 2:1227748/76‑1 (MQ=255) aaGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAAtttt < 1:165269/76‑1 (MQ=255) aGCGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTg < 1:465459/76‑1 (MQ=255) cGTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgag < 1:1727081/76‑1 (MQ=255) gTAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTgaga < 1:669427/76‑1 (MQ=255) tAAGGCTGCATGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGaa < 1:960492/76‑1 (MQ=255) tAAGGCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGaa < 1:196880/76‑1 (MQ=255) ggCTGCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaa < 1:980058/76‑1 (MQ=255) gCTGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATaaa < 1:1528276/76‑1 (MQ=255) tGCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGt < 1:1348216/76‑1 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt > 1:635470/1‑76 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt > 2:102588/1‑76 (MQ=255) gCACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt > 2:246515/1‑76 (MQ=255) gCACGGGGAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGtt < 1:1048835/76‑1 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa < 2:16561/76‑1 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa > 2:896787/1‑76 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa > 1:621314/1‑76 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa > 1:1668897/1‑76 (MQ=255) cACGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTa < 2:299217/76‑1 (MQ=255) aCGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTAc < 2:1391501/76‑1 (MQ=255) cGGGGAGAGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACg < 1:315365/76‑1 (MQ=255) ggAGAGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTa < 1:64945/76‑1 (MQ=255) gagaGTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTaa < 1:1345486/76‑1 (MQ=255) gagaGTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTaa < 1:1448767/76‑1 (MQ=255) agagTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAAc < 2:1444698/76‑1 (MQ=255) agagTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAAc > 1:707817/1‑76 (MQ=255) gTAATTCAATGACTCAAACCTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt > 1:792151/1‑76 (MQ=255) gTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt < 1:103827/76‑1 (MQ=255) gTAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt < 1:1763362/76‑1 (MQ=255) gTAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACttt < 1:1813725/76‑1 (MQ=255) tAATTCAATGACTCAAAACTATGAACTGATTTTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTa > 1:739613/1‑76 (MQ=255) tAATTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTa > 2:94566/1‑76 (MQ=255) tAAGTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTa < 1:764429/76‑1 (MQ=255) aaTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAg < 1:157439/76‑1 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc > 1:24193/1‑76 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc > 1:743341/1‑76 (MQ=255) aTTCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGc > 2:835497/1‑76 (MQ=255) ttCAATGACTCAAAACTATGAACTGATTGTGAAAGGGATCCGCAATTTTGAGAATAAAGTTACGGTAACTTTAGCg > 2:1709944/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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |