Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,812,800 | 1 | . | C | 68.0% | 48.3 / 27.4 | 25 | intergenic (‑43/+62) | celF/chbR | 6‑phospho‑beta‑glucosidase/transcriptional regulator ChbR |
Reads supporting (aligned to +/- strand): ref base . (4/4); new base C (7/10); total (11/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCA > NZ_CP009273/1812718‑1812877 | tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑ccccccc > 2:691038/1‑90 (MQ=255) tCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑ccccccc > 1:373612/1‑90 (MQ=255) cccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑cccccccc < 1:691038/90‑2 (MQ=255) ccGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCgg < 1:568994/90‑1 (MQ=255) acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATGAATAATCCa < 1:211785/90‑1 (MQ=255) acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa < 2:635024/90‑1 (MQ=255) acgacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCa < 1:682272/90‑1 (MQ=255) gacTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCa > 2:69602/1‑90 (MQ=255) cTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCaa > 2:300929/1‑90 (MQ=255) cTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAAc < 1:28868/90‑1 (MQ=255) tAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAgg < 2:299418/90‑1 (MQ=255) aTTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGAc < 1:832222/90‑1 (MQ=255) cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg > 1:358018/1‑83 (MQ=255) cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAg < 2:358018/83‑1 (MQ=255) cTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGat < 2:194562/90‑1 (MQ=255) aTAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt < 1:613965/90‑1 (MQ=255) tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt > 1:701207/1‑90 (MQ=255) tAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAAt < 1:267562/90‑1 (MQ=255) tttCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGt < 1:838784/90‑1 (MQ=255) cTTCAGTACAGAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtat > 2:77349/1‑90 (MQ=255) tCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAAc > 2:224184/1‑90 (MQ=255) gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtata > 1:208022/1‑81 (MQ=255) gAATACTGATATCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGtata < 2:208022/81‑1 (MQ=255) cTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGAc < 2:222235/52‑1 (MQ=255) cTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGAc > 1:222235/1‑52 (MQ=255) atCTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAAtt > 1:329703/1‑90 (MQ=255) atCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTc > 2:679051/1‑90 (MQ=255) cTGGCATATCTGCCCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCa > 1:732713/1‑90 (MQ=255) | TCCCGCCACCAATAGTGACGACTTTTAATTTCTGGCTCATAATTTCTCCCTTCAGTACAGAATACTGATATCTGGCATATCTG‑CCCCCCCGGACATAAATAATCCAGCAACAGGACAGATATGTGAATTGTCAGGTATAACGACTTACTGATTAAATTCA > NZ_CP009273/1812718‑1812877 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |