Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 419,492 | 1 | . | T | 50.0% | 33.4 / 46.5 | 24 | coding (1091/1176 nt) | ubiF | 2‑octaprenyl‑3‑methyl‑6‑methoxy‑1,4‑benzoquinol oxygenase |
Reads supporting (aligned to +/- strand): ref base . (9/3); new base T (12/0); total (21/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.17e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATA > minE/419423‑419561 | aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:302409/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:989122/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:4209152/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:403597/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:3249841/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:2644253/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:2104427/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:1886727/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:1740246/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:2092522/1‑71 (MQ=255) aTAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCAGGATTCAGCAATAATCTGCCACCACTGCG‑t > 1:464081/1‑71 (MQ=255) cAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAAtct > 1:2632579/1‑71 (MQ=255) cAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAAtct > 1:312269/1‑71 (MQ=255) cAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAAtct > 1:3200062/1‑71 (MQ=255) aaGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCg > 1:537050/1‑71 (MQ=255) gATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCGGGTTAATgg > 1:2544356/1‑70 (MQ=255) ttttATGCTGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:4251146/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATggcg > 1:3715461/1‑67 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:3366916/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:915348/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:4273469/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:1466003/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:3966077/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:155622/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:1162266/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:2932114/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:2561538/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCGTTTTTATGCGTAATCTCGGGTTAATGgcggcg > 1:2551501/1‑70 (MQ=255) ttttATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑ttttta > 1:499108/1‑42 (MQ=255) ttATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGc > 1:227462/1‑71 (MQ=255) ttATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGc > 1:2183966/1‑71 (MQ=255) ttATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGc > 1:1196352/1‑71 (MQ=255) ttATGCCGGATTCAGCAATAATCTGCCAACACTGCG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGc > 1:66583/1‑71 (MQ=255) cG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAata < 1:3701217/71‑1 (MQ=255) cG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAata < 1:1398236/71‑1 (MQ=255) cG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAata < 1:280919/71‑1 (MQ=255) g‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGGTGAAACGTCAGGCGCTGAAata < 1:399721/70‑1 (MQ=255) | ATAACTTCATTATGCAAAGCGGTATGGATCTGTTTTATGCCGGATTCAGCAATAATCTGCCACCACTGCG‑TTTTATGCGTAATCTCGGGTTAATGGCGGCGGAGCGTGCTGGCGTGTTGAAACGTCAGGCGCTGAAATA > minE/419423‑419561 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |