Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,631,772 | 0 | A | C | 24.1% | 130.0 / 10.3 | 54 | intergenic (‑267/‑141) | rpoE/nadB | RNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase |
Reads supporting (aligned to +/- strand): ref base A (37/4); new base C (0/13); total (37/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.15e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGCA > REL606/2631697‑2631847 | tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca > 1:355507/1‑76 (MQ=255) aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCGTTTTTATTTCACCACCTAATCCTccac > 2:1044395/1‑76 (MQ=255) tCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacc < 2:1679633/76‑1 (MQ=255) gtgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagc > 1:581485/1‑76 (MQ=255) gtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcca > 1:680702/1‑76 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag > 1:962990/1‑76 (MQ=255) tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag > 2:1758599/1‑76 (MQ=255) gAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTa > 2:898070/1‑76 (MQ=255) gAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTa > 1:1344084/1‑76 (MQ=255) aaaTGCCCATTTAAATGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATGCTCCACCAGCCAGTaa > 2:992277/1‑76 (MQ=255) aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAAc > 1:1363599/1‑76 (MQ=255) aTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACt > 2:834328/1‑76 (MQ=255) aTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACt > 1:1187927/1‑76 (MQ=255) tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTCCCCCCCCTAATCCCCCCCCAGCCAGTAACtt < 2:1271872/76‑1 (MQ=255) ttACCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCtttt < 2:624405/76‑1 (MQ=255) ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAACCCTCCCCCAGCCAGTAACTTCTCtttt < 2:1725322/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCCCCACCTAATCCTCCCCCAGCCAGTAACTTCTCttttt < 2:157334/76‑1 (MQ=255) tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAACCCTCCACCAGCCAGTAACTTCTCttttt < 1:131380/76‑1 (MQ=255) tGCCGGAAGAGTACCCCACCGGCCTTTTTATTTCACCACCTAACCCTCCACCAGCCAGTAACTTCTCTTTTTCTcg < 2:731865/76‑1 (MQ=255) gCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc > 2:235431/1‑76 (MQ=255) ctaaaGAGACACCAAACGCCCTTTTTATTTCCCCCCCTAACCCCCCACCAGCCAGTAACTTCTCTTTTTCTcgccg < 2:1186384/73‑1 (MQ=255) cTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccg < 2:1647510/76‑1 (MQ=255) tGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc > 1:758701/1‑76 (MQ=255) gAAGAGTACCCCAACGCCCTTTTTATTCCCCCCCCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc < 1:126826/76‑1 (MQ=255) aaGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGccc > 1:261025/1‑76 (MQ=255) agagTAACCCACCGGCCTTTTTATTTCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 1:472803/76‑1 (MQ=255) agagTAACCCAACGGCCTTTTTATTTCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt < 2:1564990/76‑1 (MQ=255) gagTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg > 1:1595944/1‑76 (MQ=255) tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt > 1:1547962/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCTCCGCCCTGCGTc > 2:1622369/1‑76 (MQ=255) aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc > 2:415903/1‑76 (MQ=255) ccAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGc > 1:140633/1‑76 (MQ=255) cAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCg > 2:891183/1‑76 (MQ=255) aaCGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt > 2:1444159/1‑76 (MQ=255) aCGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtg < 1:629362/76‑1 (MQ=255) cgcccTTTTTATTCCCCCACCGAAGCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtgt < 1:89059/73‑1 (MQ=255) gCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt > 1:1389533/1‑76 (MQ=255) cTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg > 1:885374/1‑76 (MQ=255) tttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa > 1:695892/1‑76 (MQ=255) ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 1:1034385/1‑76 (MQ=255) ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc > 1:543266/1‑76 (MQ=255) ttATTCCCCCACCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc < 1:1210815/76‑1 (MQ=255) tATTCCCCCCCCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt < 1:1337947/76‑1 (MQ=255) tATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACt < 2:580754/76‑1 (MQ=255) tttCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt > 1:1670148/1‑76 (MQ=255) tccaccccCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa < 1:1560797/74‑1 (MQ=255) tCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa > 1:1056980/1‑76 (MQ=255) caccacCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAAc > 1:494601/1‑76 (MQ=255) caccacCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAAc > 1:1168079/1‑76 (MQ=255) ccaccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaa > 2:675073/1‑76 (MQ=255) ccccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaaa < 2:1691467/74‑1 (MQ=255) caccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaaa > 1:1264501/1‑76 (MQ=255) cccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAat < 1:392944/75‑1 (MQ=255) ccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata > 1:819040/1‑76 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat > 2:718175/1‑76 (MQ=255) cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat > 1:1199601/1‑76 (MQ=255) tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt > 1:624590/1‑76 (MQ=255) aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGCGTTTAGCCACTGTAACAAATATTa > 2:1648071/1‑76 (MQ=255) tcctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaa > 2:1503653/1‑76 (MQ=255) cccaccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATa < 1:1737109/75‑1 (MQ=255) accaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGCa > 2:1170725/1‑76 (MQ=255) | TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAGCA > REL606/2631697‑2631847 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |