Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 536,014 | 0 | T | G | 59.1% | 10.0 / 23.5 | 22 | intergenic (+50/‑7) | allB/ybbY | allantoinase AllB/uracil/xanthine transporter |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (13/0); total (14/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.81e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT > NZ_CP009273/535932‑536103 | aaaGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGGGtttt > 1:17485/1‑89 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTc < 2:706580/90‑1 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGAATTCTCCGGGTAAGGGGTGTTTTTTTTc > 2:156672/1‑90 (MQ=255) ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGTTTTAAGGGTGTTTATGTTCCATTTTTCa > 2:476913/1‑90 (MQ=255) tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTGTTTAATTTTTCAGTCAGCcgc > 1:299338/1‑90 (MQ=255) tAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGcc < 2:126298/90‑1 (MQ=255) ggCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGTTTTTGTTTAATTTTTCAGTCAGCCGCGAAAGCCTGTTAt > 2:317426/1‑90 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCCGTTAAGGGGGGGTTATGTTTAATTTTTCAGTCAGCCGCGAAAGCCCGTTTTCAggg > 1:499592/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGGGTTTTTGTTTAATTTTTCAGTCAGCCGCGGAAGGCTGTTATCAggg > 2:317356/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGTGTTTTTTTTTAATTTTTCAGTCAGCCGCGGAAGCCCGTTTTTAggg > 1:35439/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGGTTTTGTTCAATTTTTCAGTCAGCCGCGGAAGACCGTTTTCAggg > 2:437810/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGa > 2:297904/1‑90 (MQ=255) gCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGt < 2:309884/90‑1 (MQ=255) cGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGTGTTTTTTTTCAATTTTTCAGTCAGCCCCGAAAGCCTGTTGTTAGGGTTTTAGTGGt > 1:376145/1‑90 (MQ=255) cGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTTTGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTTAGTGGt > 2:375408/1‑90 (MQ=255) tCCTTGCGGCGGGGATTCTCCGGGTAAGGGGGGTTTTTGTTCAATTTTTTAGTCAGCCGCGAAAGCCCGTTATCAGGATTTTTGTGGttt > 1:462155/1‑90 (MQ=255) gggCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttt > 2:184473/1‑90 (MQ=255) ttCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGc < 2:218614/90‑1 (MQ=255) ttCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGc < 1:634978/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:414016/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:40519/90‑1 (MQ=255) tAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGtt < 2:499592/90‑1 (MQ=255) tgtgtTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt > 1:177666/1‑90 (MQ=255) | AAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGT > NZ_CP009273/535932‑536103 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |