Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 169,172 | 0 | C | G | 53.3% | 11.4 / 19.7 | 15 | G408G (GGC→GGG) | fhuB | Fe(3+)‑hydroxamate ABC transporter permease FhuB |
Reads supporting (aligned to +/- strand): ref base C (0/7); new base G (8/0); total (8/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACC‑GGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGT > NZ_CP009273/169083‑169252 | aGCGGGGCGTTGCTCGAGGATTTAATGCGCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGc < 2:15396/90‑1 (MQ=255) aGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGc < 1:107674/90‑1 (MQ=255) aGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGc < 2:245476/90‑1 (MQ=255) ggggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGt < 2:7347/90‑1 (MQ=255) gggCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGta > 2:45537/1‑90 (MQ=255) tCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGc < 2:398050/90‑1 (MQ=255) aTTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGTATTTTTTAGCGACTGa > 2:246531/1‑90 (MQ=255) tGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACC‑GGaa < 2:139111/90‑1 (MQ=255) gCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGGGTCATGCTGGCGGGGGGGGGGTGTTTTTTTTAGCGGCTGACC‑GGgaa > 1:175598/1‑90 (MQ=255) tGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGCGGGGTGTATTATTCAGCGACTGACC‑GGAAAACCg > 2:199167/1‑90 (MQ=255) gcgcTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGTATTATTCAGCGGCCGACC‑GGAAAccccgt > 2:249042/1‑87 (MQ=255) tGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACC‑GGAAACCCGATGGCa < 1:87325/90‑1 (MQ=255) aggCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACC‑GGAAACCCGATGGCa < 2:50873/89‑1 (MQ=255) gcgAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGCGGGGTGTATTATTCAGCGACTGACC‑GGAAACCCCATGGGAAAccc > 2:319486/1‑90 (MQ=255) tATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACC‑GGAAACCCGATGGCAAGCCCGGAAGt < 1:83397/90‑1 (MQ=255) aTGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGTATTTTTTAGCGACTGACCGGGAAACCCCTGGGCAGCCCCGGAGGt > 1:50873/1‑90 (MQ=255) gcgcTGTTTGCGGGCGTCATGCTGGCGGGGGGGGGGTGGATTTTTTAGCGGCTGACC‑GGAAACCCGATGGCCAAGCCGGGAGGGGTgggg > 1:92580/1‑90 (MQ=255) ggtggcggGCTGTATTATTCAGCGACTGACC‑GGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCgt < 2:319748/90‑1 (MQ=255) | AGCGGGGCGTTGCTCGAGGATTTAATGCCCTGGCGCTGGCCGCGAATTATGGCGGCGCTGTTTGCGGGCGTCATGCTGGCGGTGGCGGGCTGTATTATTCAGCGACTGACC‑GGAAACCCGATGGCAAGCCCGGAAGTGCTGGGGATTAGCTCCGGCGCGGCGTTTGGCGT > NZ_CP009273/169083‑169252 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |