Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,123,160 | 0 | A | G | 52.4% | 22.8 / 27.2 | 21 | E206G (GAA→GGA) | metF | methylenetetrahydrofolate reductase |
Reads supporting (aligned to +/- strand): ref base A (2/8); new base G (11/0); total (13/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGAT > NZ_CP009273/4123075‑4123247 | aaCCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGgaaat > 2:154990/1‑88 (MQ=255) cgATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCGTTGATTTGGGAATTA‑TTcccg > 2:230026/1‑88 (MQ=255) aCTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTGATGTGGGAATTT‑TTTCCGGAAtt > 1:187303/1‑90 (MQ=255) gTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTGATGTGGGAATTA‑TTTCGGGGATTTTTcc > 2:700514/1‑90 (MQ=255) gTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGGAAATA‑TTTCGGGGATTTTTcc > 2:376815/1‑90 (MQ=255) gTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGcc < 1:700514/90‑1 (MQ=255) tcttcttcGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCgg < 1:934920/90‑1 (MQ=255) gATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTGATGTGGGAATTA‑TTTCGGGGATTTTTTCGGGATC‑TAAc > 1:800086/1‑90 (MQ=255) gAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGGTTGATGTGGGAATTA‑TTTCGGGGATTTTTCCGGGATT‑TTACTTTtaa > 2:466276/1‑90 (MQ=255) cTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAAATA‑TTTCCGGGATTTTTCCGGGATCTTAATTTTAAAAAgg > 2:530219/1‑90 (MQ=255) aCCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGGATTGATGTGGGAAATA‑TTTCGGGGATTTTTTCGGGGTC‑TTACTTTTAAAAGGGga > 1:367113/1‑90 (MQ=255) tCGTGGCCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGc < 1:877892/90‑1 (MQ=255) tCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGc < 1:376815/90‑1 (MQ=255) gTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCg < 2:367113/90‑1 (MQ=255) gACCGCTGTGTATCGGCGGGCCTTGGTGTGGGAATTATTTCCGGGATTTTTG‑CGGGATC‑TAACTTTAAAAAGGCGAAGAAAATTTGCGat > 2:569422/1‑90 (MQ=255) gACCGCTGTGTATCGGCGGGCATTGATTTGGGAATTA‑TTTCCGGGATTTTTTCGGGGTC‑TTACTTTTAAAAAGCGGAGAAAATTGCCggt > 1:473677/1‑88 (MQ=255) aCCGCTGTGTATCGGCGGGCATTGATGTGGGAAATA‑TTTCGGGGATTTTTCCGGGATC‑TAAATTTAAAAAAGGGGAGAAAATTTTCGata > 2:678013/1‑90 (MQ=255) gcggGCATTGATGTGGGAATTA‑TTTCGGGGATTTTTCCGGGATC‑TAACTTTTAACAAGGGGAGAAATTTTCCGATATGACCAACGTGCGt > 1:641890/1‑90 (MQ=255) gCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTc < 1:69670/90‑1 (MQ=255) gCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTc < 1:662495/90‑1 (MQ=255) tgGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGa > 2:858236/1‑90 (MQ=255) ggAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGAt < 1:226493/90‑1 (MQ=255) | AACCGCGCGATTACTCAGTTCTTCTTCGATGTCGAAAGCTACCTGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTA‑TTCCGGGAATTTTGCCGGTATC‑TAACTTTAAACAGGCGAAGAAATTTGCCGATATGACCAACGTGCGTATTCCGGCGTGGAT > NZ_CP009273/4123075‑4123247 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |