| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,856,024 | 0 | C | G | 51.3% | ‑5.1 / 23.5 | 44 | pseudogene (940/1617 nt) | BW25113_RS19145 | PTS sugar transporter subunit IIB |
| Reads supporting (aligned to +/- strand): ref base C (14/5); new base G (4/16); total (22/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.23e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAC > NZ_CP009273/3855935‑3856104 | ggCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCTAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc < 2:527810/90‑1 (MQ=255) gtgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGgggg > 2:54447/1‑90 (MQ=255) cGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGGTCGGTAATTCCCAGCAGCATGGCGGTTCAGGTTGTGGGGAGCAg > 2:1780094/1‑90 (MQ=255) gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGTGGGGATCAGc > 1:932258/1‑90 (MQ=255) gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGagcagc > 1:1801500/1‑90 (MQ=255) cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCa < 1:1775706/90‑1 (MQ=255) cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCa < 1:731308/90‑1 (MQ=255) aaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGTGGGAATCAGCAAg > 2:554737/1‑90 (MQ=255) cAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGGAAGcc > 1:1458446/1‑90 (MQ=255) gCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTGCGGGGAACAGCAAAcccc > 1:662169/1‑89 (MQ=255) cAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTGGGGGATCAGCAAAccccc > 2:258431/1‑88 (MQ=255) gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAACCCCCCCCtt > 2:540288/1‑90 (MQ=255) taaaaTGAAAAGGAAGGAAAATCCCGGGGGTTCGGTAATCCCCACCAGCAGGGGGGTTAAGGTGGGGGGAATCAGCAAGCCCGCCACttt < 2:246661/88‑1 (MQ=255) gAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGGGGGAATCAGCAAACCCGCCACtttt > 2:1031948/1‑90 (MQ=255) gAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCCCtttt > 1:246661/1‑90 (MQ=255) aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCggg < 1:1653934/66‑1 (MQ=255) aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCggg > 2:1653935/1‑66 (MQ=255) aaTAAAAAGAAAGGTAAATCCCGGCGTTTGGTTAATCCCCCCCAGCAGGGCGTTTAAGGTGGGGGGAACCACCAACCCCGCCACTTTTAc < 2:794739/90‑1 (MQ=255) taaaCAGGAAGGAAATTTCCGGGGTTTGGTTACTCCCCCCCAGCGGGGGGGTTAGGGTTGGGGGAACCAGCAAGCCCGCCACTTTTAccc < 1:258431/88‑1 (MQ=255) aaCAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGGGGGGAACAGCAAACCCGCCACTTTTTACCGa > 1:720599/1‑90 (MQ=255) ggAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑tttt > 2:1176216/1‑90 (MQ=255) ggAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCCCTTTTACCCGactgt > 2:1679949/1‑86 (MQ=255) aaGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAGCAGGAAACCCGCCACCTTTTCCCCA‑tttttc > 1:1437230/1‑88 (MQ=255) aaGGTAAATCCCGGGGGTTCGGTATTTCCCACCAGAAGGGCGGTTAGGGTTGGGGAAATCAGAAAGCCCGCCACTTTTACCCGA‑TTTTcc < 2:1530540/90‑1 (MQ=255) aaaTCCCGGGGGTTGGGTAATCCCCCCCAGCAGGGGGGTTAAGGTGGGGGGAATCAGCAACCCCGCCACTTTTACCCGA‑TTTTCCGGTGc < 2:662169/90‑1 (MQ=255) aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGTGGGGATCAGGAAACCCCCCCCCTTTACCCCC‑TTTTTCGGGGCg > 1:77305/1‑90 (MQ=255) ttCCAGCGGTGGGGAATTCCCCCCCAGCAGGGGGTTTGGGGTGGGGGAAAACAGCAACCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGc < 1:252658/90‑1 (MQ=255) tcccAGCGTTTCGATAATCCCCCCCAGCGGGGGGTTTAGGGTGGGGGAAAAAAACACCCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGc < 1:1649634/88‑1 (MQ=255) ccAGGGGTTCGGTATTCCCCCCCAGCAGGGGGGTTAAGGTGGGGGAAACCAGAAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAg < 2:720599/90‑1 (MQ=255) gccGGTTCGATAATCCCCACCAGCATGGGGGTTAGGTTCGGGGGAACCACAACGCCCGCCACTTTTACCCGA‑TTTTCCGATGCGGCAGta < 2:453643/88‑2 (MQ=255) gCGGGTCGTTAATCCCCCCCAGCGGGGGGGTTAGGGTGGCGGAAACCAACAACCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGa < 1:495707/90‑1 (MQ=255) ggTTCGGTAATTCCCACCAGCATGGTGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAg > 2:1194950/1‑90 (MQ=255) tgggTATTCCCCCCCAGCAGGGGGTTTAGGGTGGGGGAAATCAGAAACCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTAc < 1:1739682/88‑1 (MQ=255) tCGGTAATCCCCCCCAGGATGGGGTTTAAGGTGGGGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGTAGTGAAGTAc < 2:1089010/90‑1 (MQ=255) tCGGAAATCCCCCCCAGAATGGGGTTTAGGGTGGGGGGAACCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTAc < 2:1437231/90‑1 (MQ=255) ttAATCCCCCCCAGCAGGGGGGTTAAGGTGGGGGAAATCAGCAACCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATc < 2:1024107/89‑1 (MQ=255) attccccACCAGCGGGGGGGTAAGGGTGGGGGGAACCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGc < 2:932259/86‑1 (MQ=255) actccccACCAGCGGGGGGGTAAGGGTGGGGGAAACCAGAAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGc < 2:77305/86‑1 (MQ=255) aaTCCCCCCCAGCAGGGGGTTTAGGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGc < 2:1801501/90‑1 (MQ=255) aTTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGGAGTGGAGGACATCGCt > 2:66173/1‑90 (MQ=255) tCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCCCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGGGAAGGACATCGCTaa > 2:1709222/1‑90 (MQ=255) caccaGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGa > 1:630545/1‑90 (MQ=255) caccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCCCTTTTACCCGA‑TTTTCCGGTGGGGCAGTGAAGTACATCGCGaaaaa > 2:1586095/1‑88 (MQ=255) ccaGCATGGCGGTTAAGGTTTCGGGGATCAACAAACCCCCCCCCTTTACCCCACTTTTCCGG‑GCGGGAGGGGAGGACACAGCTAaaaaaa > 1:331021/1‑90 (MQ=255) ggtttAAGGTCAAGGGAGCCAGCAACCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACgg < 2:377621/87‑1 (MQ=255) ttAGGGTGGGGGAAAACAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAc < 2:465300/90‑1 (MQ=255) | GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGA‑TTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAC > NZ_CP009273/3855935‑3856104 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |