Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 536,014 | 0 | T | G | 46.9% | 65.3 / 10.4 | 49 | intergenic (+50/‑7) | allB/ybbY | allantoinase AllB/uracil/xanthine transporter |
Reads supporting (aligned to +/- strand): ref base T (3/23); new base G (23/0); total (26/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.48e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTG > NZ_CP009273/535925‑536099 | tGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgt < 2:1927999/90‑1 (MQ=255) gcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGgtgtg > 1:656006/1‑90 (MQ=255) gcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGgggtg > 1:1634014/1‑90 (MQ=255) gcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGGTTAGgtggg > 2:1308386/1‑88 (MQ=255) gAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGACCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGGGttt > 1:1714298/1‑90 (MQ=255) aaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTAt > 1:164363/1‑90 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTc < 1:740244/90‑1 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTc < 2:805793/90‑1 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTc < 1:618197/90‑1 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTc < 2:520962/90‑1 (MQ=255) tCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGTGTTTTTG‑TTTAAtttttc > 2:657626/1‑88 (MQ=255) agcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGa < 2:2007161/90‑1 (MQ=255) gTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGTGTGTTTTTT‑TTTAATTTTTCAGTCAGCCGCGaaaac > 2:838517/1‑88 (MQ=255) tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGcc < 2:972319/90‑1 (MQ=255) tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGcc < 2:1668073/90‑1 (MQ=255) tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGcc < 1:838517/90‑1 (MQ=255) tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGcc < 1:1483217/90‑1 (MQ=255) aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCGGGTAAGGGGGGTTTTTT‑TTTAATTTTTCAGTCAGCCGCGAAAGCCt > 2:1477679/1‑90 (MQ=255) aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTTTG‑TTCAATTTTTCAGTCAGCCGCGAAAGCCt > 1:1873655/1‑90 (MQ=255) aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTTTG‑TTCAATTTTTGAGTCCGCCGCGAAAGCCt > 2:822198/1‑90 (MQ=255) tCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGGGTTTTTG‑TTTAATTTTTTAGTCAGCCGCGGAAGCCTGt > 2:1602124/1‑90 (MQ=255) ggCCCCTGCAATGCCCGTCCTTGCGGCGGGTTTTCTCCGGTTGAGGGGGGTTTTTGTTTCAGTTTTTTA‑TAAGCCGCGGAAGCCTGGTAt > 2:1289918/1‑90 (MQ=255) ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGGGTTTTTT‑TTTAATTTTTCAGTCAGCCGCGAAAACCTGtttt > 1:1647194/1‑88 (MQ=255) ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGTTTTTG‑TTCAATTTTTCAGTCAGCCGCGAAAGCCTGtttt > 2:938255/1‑88 (MQ=255) ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTTAGGGGGGTTTTTG‑TTCAATTTTTTAGTCAGGCGCGGAAGCCTGTTAt > 2:1590600/1‑90 (MQ=255) cccTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCGGTTAAGGGGGGGTTATG‑TTCCATTTTTCGGTCAGCCGCGCAAGCCCGTTATCAg > 1:1858625/1‑90 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCGGGTTAGGGGGGGTTTTT‑TTTAATTTTTCAGTCAGCCGCGGAAGCCTGTTTTTAggg > 2:728213/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCGGGTAAGGGGTGTTTTTT‑TTCAATTTTTCGGTCAGCCGCGAAAACCTGTTTTTAGGa > 2:1168204/1‑90 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCGGGTAAGGGGGGTTTTTG‑TTTAATTTTTTAGGCAGGCGCGAAAGCCTGTTTTTAggg > 2:1891323/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGGGTTTTTT‑TTTAATTTTTCAGTCAGCCGCGAAAGCCCGTTTTTAggg > 1:1091281/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGcc > 1:505652/1‑79 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGcc < 2:505652/79‑1 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa > 1:103271/1‑90 (MQ=255) tGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGAt < 1:325309/90‑1 (MQ=255) gCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGAtt > 1:170521/1‑90 (MQ=255) aTGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa > 1:179588/1‑90 (MQ=255) tGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTATG‑TTCAATTTTTCAGTCAGCCGCGAAAACCTGTTATCAGGATTTCAg > 2:325309/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGTGTTTTTG‑TTCAATTTTTCCGTCAGCCGCGAAAGCCTGTTATTAGGGTTTTTGTGGt > 2:2070817/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGTGTTTTTT‑TTTAATTTTTCAGGCAGCCGCGGAAACCTGTTTTTAGGGTTTTTGTGGt > 2:1864612/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATG‑TTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTAGTggg > 1:878801/1‑89 (MQ=255) tCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGttt < 2:1873655/90‑1 (MQ=255) gcggcgGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCtt < 1:1864612/90‑1 (MQ=255) ttCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGc > 2:1719725/1‑90 (MQ=255) tcCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAac < 2:109754/90‑1 (MQ=255) tcCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAac < 1:864596/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 1:393797/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:388195/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:2044706/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 1:728213/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:1102179/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:972851/90‑1 (MQ=255) ggTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACg < 1:1187785/90‑1 (MQ=255) aaGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTtg < 1:1296103/90‑1 (MQ=255) | TGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATG‑TTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTG > NZ_CP009273/535925‑536099 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |