| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 1,292,259 | 0 | T | C | 61.5% | 17.0 / 30.2 | 27 | intergenic (+100/+44) | galU/hns | glucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS |
| Reads supporting (aligned to +/- strand): ref base T (9/1); new base C (0/16); total (9/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.20e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAT > REL606/1292184‑1292333 | aaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑t > 2:52235/1‑76 (MQ=255)aaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑t > 2:1794827/1‑76 (MQ=255) cGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑Tcc < 1:1323531/76‑1 (MQ=255) gACGCTTTTTTTTTATAGTTTTTTTTTA‑TAAA‑TTGTTCTAAAC‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGAt < 2:627904/76‑1 (MQ=255) gcgcGTTTTTTTTTAATTTTTTTTTTTTAAATTTG‑TCTAAAC‑CCGGAAAAAAAAAAA‑CCCCGCCGATGGCGGGAtt < 2:1888897/73‑1 (MQ=255) gctgTTTTTTTATAGCTTATTTTTATAAAA‑TTG‑TCTAAAA‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGAtttt < 2:369791/73‑1 (MQ=255) ctttttttttATATTTTTTTTTTA‑TAAA‑TTT‑TTTTAAACCCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTa < 1:1721795/73‑1 (MQ=255) tttttttttATAGTTTTTTCTTATTAAA‑TTG‑TCTTAAC‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTaa < 2:1578484/74‑1 (MQ=255) tttttttttATAGTTTATTTTTATAAAA‑TTG‑TCTAAAA‑CCGGAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTaa < 2:1541475/74‑1 (MQ=255) tttttATATTTTTTTTTTATAAAT‑TTT‑TCTAAAC‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCaa < 2:1395509/76‑1 (MQ=255) tttttATATTTTTTTTTTATAAAA‑TTG‑T‑TTAAA‑CCCGAAAAAAAAAAACCCCCGCCGATGGCGGGATTTTAAGCaa < 1:88146/76‑1 (MQ=255) ttATAGCTTATTCTTATTAAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑TCCCGCCGATGGCGGGATTTTAAGCAAGTg > 1:32234/1‑76 (MQ=255) tatttttattAAA‑TTG‑TCTAAAA‑CCGGAAAATAAAAAC‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTAc < 2:1763457/76‑1 (MQ=255) tttattaAA‑TTG‑TCTAAAA‑CCGGAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACaaaa < 2:215162/75‑1 (MQ=255) attaAA‑TTG‑TCTAAAA‑CCGGAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGat < 1:36493/76‑1 (MQ=255) ttaAT‑TTT‑TCTTAAA‑CCGGAAAATAAAAAA‑GCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt < 2:45455/76‑1 (MQ=255) ttaAT‑TTG‑TCAAAAC‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt < 1:643868/76‑1 (MQ=255) taaaa‑TTG‑TTTAAAC‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt < 2:999743/74‑1 (MQ=255) taAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑TCACCCCGGTGGCGGGTTTTTTATCAAATGCAATCTTCCAAAGatta > 2:155189/1‑76 (MQ=255) aaa‑TTT‑TCTAAAC‑CCGGACAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattat < 2:1643717/76‑1 (MQ=255) aaa‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑TCCCGCCGATGGCGGGATTTTAAGCAATTGCAATCTACAAAAGattat > 2:138146/1‑76 (MQ=255) tCTAAAC‑CCGAAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg < 1:346347/76‑1 (MQ=255) tCTAAAA‑CCGGACAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg < 2:398822/76‑1 (MQ=255) tAAA‑CCGGACAATAAAAAA‑TCCCCCCCGTGGCGGGGGTTTTAACAAAGGCAATCTACAAAACAATATTTCTTTAtt > 1:599248/1‑75 (MQ=255) aaaa‑CCCGAAAAAAAAAAACCCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt < 1:1084831/75‑1 (MQ=255) aaa‑CCGGACAATAAAAAA‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCa > 2:266991/1‑76 (MQ=255) aaa‑CCGGAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCa < 1:1520075/76‑1 (MQ=255) aa‑CCGGAAAAAAAAAAA‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAg < 1:1659373/76‑1 (MQ=255) a‑CCGAAAAAAAAAAAC‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg < 1:1095979/76‑1 (MQ=255) ggAAAAAAAAAAA‑CCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaaa < 2:1205033/76‑1 (MQ=255) aaaaa‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg > 2:1920915/1‑76 (MQ=255) aaa‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAggg > 2:194963/1‑76 (MQ=255) aaa‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAggg > 1:31465/1‑76 (MQ=255) actCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa < 2:1732828/74‑1 (MQ=255) a‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAATATTATTGCTTGATCAGGAAATCGTCGAGGGAt > 2:1026514/1‑76 (MQ=255) a‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAt > 2:1026538/1‑76 (MQ=255) | AACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAA‑TTG‑TCTTAAA‑CCGGACAATAAAAAA‑TCCCGCCGATGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAT > REL606/1292184‑1292333 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |