Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,092,675 | 0 | A | G | 54.5% | 2.2 / 12.3 | 11 | noncoding (121/137 nt) | BW25113_RS23065 | RtT sRNA |
Reads supporting (aligned to +/- strand): ref base A (0/5); new base G (6/0); total (6/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.16e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCA > NZ_CP009273/1092591‑1092703 | tcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAc < 2:588462/90‑1 (MQ=33) ccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTGACtcgcct > 2:1774677/1‑88 (MQ=11) cccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTGACtcgcttc > 1:848502/1‑90 (MQ=11) tgtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctcg < 1:1269082/90‑1 (MQ=33) gtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGTTGAATCCCTtccctcgc > 2:2950207/1‑90 (MQ=11) tgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGcc < 2:2151355/90‑1 (MQ=33) gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCCCTTCGCTCcccc > 2:1275383/1‑90 (MQ=11) gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAt < 2:1945815/71‑1 (MQ=34) aTACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTGAATCCCTTCGCTCGCCCCGCgg > 2:1540069/1‑90 (MQ=11) aTACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTTGCTCGCCCCGCgg > 2:946546/1‑90 (MQ=11) cacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa < 1:265894/90‑1 (MQ=31) | TCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCA > NZ_CP009273/1092591‑1092703 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |