Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,554,830 | 0 | A | G | 39.5% | 118.2 / 24.2 | 86 | intergenic (+30/+27) | yjiV/mcrC | b4486; conserved protein/5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
Reads supporting (aligned to +/- strand): ref base A (25/27); new base G (34/0); total (59/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.48e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 161 total aligned reads displayed.
AGAGAGCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATACAC > REL606/4554755‑4554905 | agagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGa < 2:7366723/76‑1 (MQ=255) agagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGa < 2:6704195/76‑1 (MQ=255) gagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAg < 2:4921297/76‑1 (MQ=255) gagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAg < 2:3096105/76‑1 (MQ=255) agagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGc < 2:6281645/76‑1 (MQ=255) agagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGc < 1:3489634/76‑1 (MQ=255) agagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGc < 1:983539/76‑1 (MQ=255) ggCTGGGCAGCGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCgg < 2:7295111/75‑1 (MQ=255) agCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCgg > 2:315141/1‑76 (MQ=255) agCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCgg > 2:1590035/1‑76 (MQ=255) gCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGt < 2:7587664/76‑1 (MQ=255) gCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGt < 2:2543530/76‑1 (MQ=255) cTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCgggg > 2:1475293/1‑74 (MQ=255) cTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTa < 2:4883975/76‑1 (MQ=255) gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 1:1191198/76‑1 (MQ=255) gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 2:5395781/76‑1 (MQ=255) gggCAACGGATGTGCGGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 2:4864828/76‑1 (MQ=255) ggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGTTGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTg < 2:5831407/76‑1 (MQ=255) ggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTg < 2:3683587/76‑1 (MQ=255) ggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTg < 2:3553179/76‑1 (MQ=255) gCAGCGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGt < 2:414754/76‑1 (MQ=255) gCAACGGATGTGCGGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGt < 1:1882864/76‑1 (MQ=255) aaCGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCt > 1:4838605/1‑76 (MQ=255) cGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGCTTTGTCtat > 2:7628453/1‑76 (MQ=255) cGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGGTTGTCtat > 2:4988502/1‑76 (MQ=255) ggtgtgTTGGGGGGGGTGTCCCCCTCCTGAGTATGCGCCCCCCCCGCCGTCGTGGGGGGAGGGTTATTGTCTATaa < 2:4296494/74‑1 (MQ=255) tgtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGTTTGTCTGTtaaa > 1:5262857/1‑76 (MQ=255) tgtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATaaaa > 2:1286483/1‑76 (MQ=255) gtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGGTTTTaaaaa > 1:6206598/1‑75 (MQ=255) gcgggtggAGTGGACCGCCTCCTGATGATGCGCCCCCCCCGACGTCGGGGCGGGAGCGGTATTGTCTATAAAACtt < 2:4633894/73‑1 (MQ=255) gCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGGTTGTTTTTaaaaaat > 1:4411345/1‑73 (MQ=255) gCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTTTaaaaact > 2:7544765/1‑73 (MQ=255) cgggtggAGGTGACCGCCTCCTGATGATGCGCCGCCCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTa < 2:7072847/74‑1 (MQ=255) cTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTa < 1:7533293/76‑1 (MQ=255) tggtggAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGGTTTTCTTTAAAAATTAc > 1:5820004/1‑76 (MQ=255) ggggTGACCGCCTCCTGATGATGGGCCCCTCCCGACGTCGTGTGGGGGGCGTTATTGTCTATAAAACTTACCGCtt < 2:3672229/74‑1 (MQ=255) gAGGTGGTTCCCTCCTGAGGAGGTGCCGCCCCCGCGGTGGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCtt < 1:4237670/76‑1 (MQ=255) tGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCTATAAAAATTTACGCtttttt > 1:5930959/1‑76 (MQ=255) tGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTAttt > 2:6513658/1‑76 (MQ=255) tGATCCCCTCCTGATGAGGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTAttt < 2:1717714/76‑1 (MQ=255) gATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTAtttt > 1:5494202/1‑75 (MQ=255) cGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGGCTTATAAACATACCGCCTGTTTtaga > 2:6504720/1‑76 (MQ=255) cGCCCCCTGATGAGGGCCCGCCCCCGACGTCGTGCGGGGAGGGGTATTGTCTATAAAACTTACCGCTTATTTgaga < 2:1217994/76‑1 (MQ=255) gCCTCGTGAGGGTGTGCCGCCCCCGACGTCGTGTGGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGat < 1:6463694/76‑1 (MQ=255) ccccctGAGGATGGCCCGCCCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGata < 2:6250520/73‑1 (MQ=255) ctcctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGTATTGTGTATAAAAATTACCGCTTTTTTTTGatat > 2:6927689/1‑76 (MQ=255) ctcctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGGTTGTCTATAAAAACTACCGCTTTTTTGAGatat > 2:615231/1‑76 (MQ=255) tcttgaggaTGAGCCGCTCCCGACGTCGGGTAGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATAtt < 1:1541553/73‑1 (MQ=255) tcctGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTATTTTAGATAtt > 2:6547086/1‑76 (MQ=255) tcctGATGATGCGCCGCTCCCGACGTCGGGGGGGGAGGGGTATTGTCTATAAAACTTACCGCTTATTTGAGATAtt < 1:457113/76‑1 (MQ=255) cctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCTTTAAAACTTACCGCTTTTTTTAGATAttt > 1:5974375/1‑75 (MQ=255) cctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTATTTGAGATATTc > 2:7565731/1‑76 (MQ=255) cctGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTc < 1:7520056/76‑1 (MQ=255) cctGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTc > 1:2420854/1‑76 (MQ=255) ctGATGATGTCCCGCCCCCGCCGTGGGGACGGGAGTGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCa < 1:4234106/76‑1 (MQ=255) tgatgaTGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCTTTAAAAATTACCGCTTTTTTGTGATATTTAt > 1:2031494/1‑76 (MQ=255) tgatgaTGTGCCGCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCAt > 1:7380459/1‑76 (MQ=255) tgatgaTGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTATTTTAGATATTCAt > 2:5360362/1‑76 (MQ=255) gatgatGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTTTATAAAACTTACCGCTTATTTGAGATATTCATc > 2:1908022/1‑76 (MQ=255) gatgatGTGCCGCCCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATc < 1:124134/76‑1 (MQ=255) gatgatGTCCCCCCCCCGACGTGGTGTCGGGAGGGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATc < 1:7329509/76‑1 (MQ=255) atgatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAAATTACCGCTTATTTTTGATATTAATCg > 2:1779448/1‑76 (MQ=255) atgatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTATTTGAGATATTCATCg > 2:2101704/1‑76 (MQ=255) atgatgTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAAATTACCCCTTATTTTAGATATTTATCg > 2:6810528/1‑76 (MQ=255) tgtgCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCCCTTATTTGAGATATTCATCGaaaa > 2:7208361/1‑76 (MQ=255) gtgCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAAt > 1:1156612/1‑76 (MQ=255) tgCCGCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAAATTAACGCTTATTTGAGATATTCATCGaaaaag > 2:3158319/1‑74 (MQ=255) gCCGCTCCCGCCGTCGGGGGGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATgg < 2:3752139/76‑1 (MQ=255) ccGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGc > 2:4726825/1‑76 (MQ=255) ccGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTATTTTAGATATTCATCGAAAATGGc > 2:7057471/1‑76 (MQ=255) cGCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATAGAAAATGGCg > 2:3372912/1‑76 (MQ=255) gCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGa > 2:1819167/1‑76 (MQ=255) gCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGa > 1:7207543/1‑76 (MQ=255) cTCCCGACGTGGTGTGGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAg < 1:5192026/76‑1 (MQ=255) cTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAg < 2:6916006/76‑1 (MQ=255) ccccGCCGTCGTGGCGGGAGGGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGt < 2:3464168/75‑1 (MQ=255) cccGAGGTCGGGTCGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTa < 2:4042273/76‑1 (MQ=255) cccGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAgga > 2:6312057/1‑74 (MQ=255) cccGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAgga > 2:5165731/1‑74 (MQ=255) cGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGGAAt > 2:2299078/1‑76 (MQ=255) gACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAAtt > 2:1431094/1‑76 (MQ=255) ccgtcggGTGGGGGGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTc < 2:4314340/75‑1 (MQ=255) aCGTGGTGTCGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTc < 2:3590421/76‑1 (MQ=255) cgtggtGTGGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCt < 2:3399380/76‑1 (MQ=255) cgtcgtGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCt > 2:5240929/1‑76 (MQ=255) gtcgtGTGGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCtt < 1:1796039/76‑1 (MQ=255) gtgtggCGGGGGGGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTg < 2:1048590/74‑1 (MQ=255) cgtGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGGAATTCTTGa > 2:4971068/1‑76 (MQ=255) tgtCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATg > 2:2872864/1‑76 (MQ=255) tgtCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATg > 2:5156185/1‑76 (MQ=255) gtggggAGGGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt < 1:5079655/73‑1 (MQ=255) gtCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt > 2:3972615/1‑76 (MQ=255) gtCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt > 1:7152121/1‑76 (MQ=255) tCGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGta < 1:7057236/76‑1 (MQ=255) cGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGtat > 2:5822400/1‑76 (MQ=255) gggAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGtata > 2:1864401/1‑76 (MQ=255) ggAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATac < 1:315141/76‑1 (MQ=255) ggAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATac < 1:785832/76‑1 (MQ=255) gAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATaca < 1:6617454/76‑1 (MQ=255) aGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATacac > 2:122858/1‑76 (MQ=255) | AGAGAGCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATACAC > REL606/4554755‑4554905 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |