Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,554,830 | 0 | A | G | 41.5% | 58.4 / 30.0 | 65 | intergenic (+30/+27) | yjiV/mcrC | b4486; conserved protein/5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC |
Reads supporting (aligned to +/- strand): ref base A (17/21); new base G (27/0); total (44/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.00e-07 | |||||||||||
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 124 total aligned reads displayed.
AGAGAGCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATACAC > REL606/4554755‑4554905 | agagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGa < 2:2685887/76‑1 (MQ=255) gagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAg < 2:4230876/76‑1 (MQ=255) gagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAg < 2:1361475/76‑1 (MQ=255) gagagCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAg < 2:4746518/76‑1 (MQ=255) agagCTGGGCAACGGATTTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCCCTCCCGACGTCGTGTCGGGAGc < 1:3864632/76‑1 (MQ=255) gCTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGt < 2:4484924/76‑1 (MQ=255) cTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCgggt > 2:4564914/1‑74 (MQ=255) cTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTa < 2:2969775/76‑1 (MQ=255) cTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTa < 1:2282961/76‑1 (MQ=255) cTGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTa < 2:1137491/76‑1 (MQ=255) tGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCgggtt > 1:5114722/1‑73 (MQ=255) tGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAt < 2:4337514/76‑1 (MQ=255) tGGGCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAt < 2:3068140/76‑1 (MQ=255) tggCAACGGATGTGCTGGGGAAGTGGATCGCCGCCTTCTGACGTGCCACTCCCGGCGTCGTGTCGGGCGCGGTAtt < 2:3102569/75‑1 (MQ=255) gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 2:4286495/76‑1 (MQ=255) gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 2:517715/76‑1 (MQ=255) gggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 2:5064538/76‑1 (MQ=255) gggCAACGGATGTGCGGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTAtt < 2:3349387/76‑1 (MQ=255) ggCAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTg < 2:1216435/76‑1 (MQ=255) cAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGGTTGtt > 2:2686259/1‑75 (MQ=255) cAACGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGAttttt > 2:2242372/1‑73 (MQ=255) aaCGGATGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGtgt > 1:2758429/1‑74 (MQ=255) gtggtgCTGGGGGAGGTGATCCCCCCTTGATGGTGTCCCGCCCCCGACGTCGTGGCGGGGGCGGTATTGTCTATaa < 2:1981008/73‑1 (MQ=255) gATTTGGTGGTGGGGGGGATCCCCTCCGGAGGAGGTGCCGCCCCCGACGTCGTGGCGGGGGGGGTATTGTCTATaa < 2:2192591/76‑1 (MQ=255) aTGTGCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATaaa > 2:2169605/1‑76 (MQ=255) tgtgGTGGGGGGGGTGCTCGCCCCCTGATGATGTGCCGCCCCCGGCGTCGGGTCGGGAGCGTTATTGTCTATaaaa < 2:4702142/76‑1 (MQ=255) tgtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCTATaaaa > 1:3811743/1‑76 (MQ=255) tgtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATaaaa > 2:4135965/1‑76 (MQ=255) ggtgCTGGTGGGGTGGACCCCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATaaaa < 1:2200297/75‑1 (MQ=255) gtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTTTAAAAc > 2:2595760/1‑76 (MQ=255) gtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAc > 2:2986563/1‑76 (MQ=255) gtgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATaaaaa > 2:4511739/1‑75 (MQ=255) tgCTGGTGGAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACt > 1:2530385/1‑76 (MQ=255) gggggAGGTGTCCGCCCCCTGAGGATGTCCCCCCCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTAcc < 2:2492091/73‑1 (MQ=255) ggAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGGATGTCTATTAAAATTTCCCCt > 1:2638054/1‑76 (MQ=255) gAGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCtt > 2:1827736/1‑76 (MQ=255) aGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTa > 1:2300631/1‑76 (MQ=255) aGGTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAAATTAACGCTTa > 2:957634/1‑76 (MQ=255) gTGCTCGCCTCCTGATGATGTGCCGCCCCGGACGTCGTGTCGGGGGCGTTATTGTCTATAAAACTTACCGCTTAtt < 2:4496370/76‑1 (MQ=255) gTGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGGTTGTCTATTAAAATTACCGCTtttt > 1:830798/1‑76 (MQ=255) tGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTTTTTTAAACTTACCGGTTGttt > 1:3859284/1‑76 (MQ=255) tGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCtttttt > 2:1946725/1‑76 (MQ=255) tGATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTAttt > 2:1314659/1‑76 (MQ=255) gctCGCCTCCTGTTGATGGGCCGCCCCCGGGTTGGTGTGGGGGGGGGTATTGTCTATAAAACTTACCGCTTATTTg < 2:2993423/74‑1 (MQ=255) gATCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTTTAAAACTTACCGCTTAtttt > 1:3258527/1‑75 (MQ=255) gATCCCCCCTTGATGATGTGCCGCCCCCGACGCCGTGGGGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTg < 2:871507/76‑1 (MQ=255) tCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTTTAAAACTTACCGCTTATTTgag > 1:3231285/1‑76 (MQ=255) tCGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCCCTTATTTtag > 2:40158/1‑76 (MQ=255) ccccctcctGATGATGTGCCCCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTgag < 1:4586380/73‑1 (MQ=255) ccGCCTCTTGATGATGGCCCGCCCCCGACCTCGGGTAGGAAGGGGTATTGTCTATAAAACTTACCGCTTATTTgag < 2:2436349/75‑1 (MQ=255) cGCCTCCTGATGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTtaga > 2:840114/1‑76 (MQ=255) cctcctGAGGATGTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGata < 2:2269532/76‑1 (MQ=255) ctcctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCTATAAAAATTACCGCCTATTTGAGatat > 2:3221108/1‑76 (MQ=255) ctcctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTTAGatat > 2:2267312/1‑76 (MQ=255) ctcctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCCCTTATTTTAGatat > 1:1613349/1‑76 (MQ=255) ccctGATGATGCCCCCCCCCGGACGTGGGGGCGGGAGGGTTATTGTCTATAAAACTTACCGCTTATTTGAGATAtt < 2:170314/75‑1 (MQ=255) ctGATGATGTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTTTTTTAGATATTCa > 2:1719436/1‑76 (MQ=255) tgatgaTGTCCCGCCCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCAt < 1:712441/76‑1 (MQ=255) ttgatgTGCCGCCCCCGCCGTGGTGTCGGGCGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCg < 1:5121328/75‑1 (MQ=255) atgatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTTAGATATTCATCg > 1:5045877/1‑76 (MQ=255) atgatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAAATTACCGCTTATTTGAGATATTCATCg > 2:2980165/1‑76 (MQ=255) tgttgtgCCCCCCCCGTCGTCGTGGGGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGa < 1:678068/73‑1 (MQ=255) tgatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGa > 2:1147728/1‑76 (MQ=255) tgatgTGCCGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGa > 2:3587781/1‑76 (MQ=255) gatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAAATTACCGCTTATTTGAGATATTCATCGaa > 1:5039092/1‑76 (MQ=255) gatgTGCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaa > 2:3342513/1‑76 (MQ=255) atgTGCCGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaa > 2:4722697/1‑76 (MQ=255) tgtgCCGCTCCCGACGTCGTGTCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGaaaa > 2:3364014/1‑76 (MQ=255) gtcccGCTCCCGACGTGGTGTGGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAAt < 2:1365814/73‑1 (MQ=255) gtcGCTCCCGACGTTCTGCCGGGCGGGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATgg < 2:3739507/74‑1 (MQ=255) gCCGCTCCCGGCGTCGTGGCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATgg < 1:669407/76‑1 (MQ=255) ccGCTCCCGACGTCGTGTCGGGGGGGGGATTGTCTTTAAAACTTTCCGCCTATTTTAGGTATTTCTCGAAAAAggg > 2:3112311/1‑75 (MQ=255) ccGCTCCCGACGTCGTGTCGGGGGCGGGATTGTTTATAAAACTTACCGCTTATTTTAGATATTCATCGAAAATGGc > 2:1308662/1‑76 (MQ=255) ccGCTCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGc > 2:4658273/1‑76 (MQ=255) ccGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCCTATTTGAGATATTCATCGAAAATGGc > 2:2212012/1‑76 (MQ=255) cGCTCCCGACGTCGTGTCGGGGGGGGGGTTGTTTCTAAAAATTAACGCCTATTTTTGATATTCATCGAAAATAGCg > 2:1564034/1‑76 (MQ=255) cGCTCCCGACGTCGTGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCg > 2:3693878/1‑76 (MQ=255) gCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTTAGATATTCATCGAAAATGGCGa > 2:3285865/1‑76 (MQ=255) gCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGa > 2:2117540/1‑76 (MQ=255) gCTCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAAATTACCGCTTATTTGAGATATTCATCGAAAAAGGCGa > 2:2638825/1‑76 (MQ=255) ccTCCCGCCGTGGTGTCGGGGTGTGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGa < 1:579873/75‑1 (MQ=255) tCCGGACGTCGTGTGGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGt < 1:217143/76‑1 (MQ=255) tCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGt > 2:971742/1‑76 (MQ=255) tCCCGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGt > 1:2003526/1‑76 (MQ=255) tCCCGACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAga > 1:3492211/1‑75 (MQ=255) ccgccgTGGTGTGGGGGGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTaa < 1:3264727/76‑1 (MQ=255) ccGACGTCGTGTCGGGGGCGGTATTGTCTATAAAACTTACCGCTTTTTTTAGATATTTATCGAAAATAGCGAGCaa > 2:3521849/1‑76 (MQ=255) gACGTCGTGTCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAAtt < 1:1138496/76‑1 (MQ=255) tgcgtGTGGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCtt < 2:200802/74‑1 (MQ=255) gtcgtGTCGGGAGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATTGAAAATGGGGGGAAATtttt > 1:2429835/1‑76 (MQ=255) gtgtCGGGGGCGGGATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGGAATTCTTGAt > 2:1735135/1‑76 (MQ=255) tgtCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATg > 2:2200281/1‑76 (MQ=255) gtCGGGGGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt < 2:1143238/76‑1 (MQ=255) gtCGGGAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt < 2:1621966/76‑1 (MQ=255) gtCGGGAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt < 2:5025156/76‑1 (MQ=255) gcagggAGCGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGt < 2:1448569/73‑1 (MQ=255) gggAGGGTTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGtata < 2:623113/76‑1 (MQ=255) ggAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATac < 1:1259066/76‑1 (MQ=255) gAGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATaca > 2:1476444/1‑76 (MQ=255) aGCGGTATTGTCTATAAAACTTACCGCTTATTTGAGATATTCATCGAAAATGGCGAGTAATTCTTGATGTATacac > 1:125910/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 < 26 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |