| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,461,183 | 0 | T | . | 42.1% | 66.5 / 61.2 | 38 | pseudogene (772/4605 nt) pseudogene (772/4104 nt) | yfaS yfaS | ECK2220:JW2222+JW2221:b4500; hypothetical protein ECK2220:JW2222:b2228; hypothetical protein, N‑ter fragment |
| Reads supporting (aligned to +/- strand): ref base T (2/20); new base . (16/0); total (18/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.88e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGA > minE/1461118‑1461239 | aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:1199624/71‑1 (MQ=255)aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:365774/71‑1 (MQ=255)aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:360694/71‑1 (MQ=255)aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:3435910/71‑1 (MQ=255)aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:307622/71‑1 (MQ=255)aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:2998061/71‑1 (MQ=255)aCCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCgg < 1:1103300/71‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca < 1:2003418/70‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca < 1:901211/70‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca < 1:2575216/70‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca < 1:3864214/70‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca < 1:3119949/70‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGaaaca < 1:1446756/70‑1 (MQ=255) ttttttACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATAGGaaaca < 1:530734/70‑1 (MQ=255) cgGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgtc < 1:1697796/71‑1 (MQ=255) cgGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgtc < 1:366034/71‑1 (MQ=255) cgGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgtc < 1:783732/71‑1 (MQ=255) cgGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgtc < 1:3127748/71‑1 (MQ=255) ggTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACcgtc < 1:3647783/70‑1 (MQ=255) ccACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGtcgtcg > 1:2278914/1‑71 (MQ=255) ccACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGtcgtcg > 1:1894817/1‑71 (MQ=255) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:3871039/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:4154276/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:43120/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:688712/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:1826398/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:1330697/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:3588025/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:3552736/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:3550410/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:1617794/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:2861316/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:2827943/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:2806028/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:2680697/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:2314927/1‑37 (MQ=25) ctTTGCCTGACACTTTGCTAAGCGCCACAG‑ATCGGaa > 1:2172530/1‑37 (MQ=25) cTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTcgacga < 1:1409311/70‑1 (MQ=255) | ACCCTGTTTTTTACCCGCGGTCCACACCAGAAGCTCTTTGCCTGACACTTTGCTAAGCGCCACGGTATCGGAAACAAACACCACCGTCGTCGCCCGATACCCACCAACCATCGCCTCGACGA > minE/1461118‑1461239 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |