| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,911,898 | 1 | . | T | 37.5% | 62.2 / 82.7 | 48 | intergenic (+193/+13) | mgtA/yjgF | magnesium transporter/ketoacid‑binding protein |
| Reads supporting (aligned to +/- strand): ref base . (0/30); new base T (18/0); total (18/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.37e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCT > minE/2911851‑2911945 | gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:112495/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:3296037/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:3075965/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:3008333/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:330230/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:2857285/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:3334353/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:3387628/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:2639172/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:3468706/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:512081/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:2005144/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:1932774/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:587219/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:598076/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:1349278/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:851674/1‑60 (MQ=38)gCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGATTTTCCATCAAg > 1:1162353/1‑60 (MQ=38) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:736451/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:756786/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:3535352/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:760680/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:804204/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:3347990/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:806947/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:886500/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1079095/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:963371/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:3223294/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1922761/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1160717/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1238175/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1270784/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1302813/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1682105/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1789888/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1813783/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1830669/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:1832255/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:3159394/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2233163/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2343232/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2357714/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2652043/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2659400/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2839060/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:2988526/62‑1 (MQ=255) gCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCt < 1:3140944/62‑1 (MQ=255) | GCGATTTTGCGTTTGTCATCAGACTTAATCCGGGCATGATAGCCCGGA‑TTTCCATCAAGATTAGCGACGAACAGCGATCGCTTCGATCTCAATCT > minE/2911851‑2911945 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |