| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,678,659 | 0 | A | C | 57.1% | 12.5 / 15.8 | 14 | Y70D (TAT→GAT) | glyA | serine hydroxymethyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/8); total (6/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.33e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGC > NZ_CP009273/2678595‑2678748 | gCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATCCTCGCAACCGCCGTAGTAGCGTTTg < 1:5643/90‑1 (MQ=255) ggTAGTCAGCGCCAAACATTTTTTCCGCACGATCGATCGCAAGTTGTCCAACGATATCAACATCCCCGCAACCGCCGTAGTAGCGTTTGc < 1:269317/89‑1 (MQ=255) tAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGccc > 1:287873/1‑90 (MQ=255) gcgcCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATaa > 2:269317/1‑90 (MQ=255) gcgcCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATaa > 2:375267/1‑90 (MQ=255) aaCATTTTTTCCGCCCGATGGACCCCAATTTGTTAAACTATATCAACATCCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCa < 2:231856/90‑1 (MQ=255) cgatcgacCGCCTGTTTTCCAACAATACCAACATCCCCCCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc < 2:6163/90‑1 (MQ=255) cgaccgacCGCCAGTTGTCCAACAATATCAACATCCTCCCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc < 1:375267/90‑1 (MQ=255) gCCATTTTTACAACGATATCAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTgagac < 1:285283/90‑2 (MQ=255) tCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGcc > 1:100238/1‑90 (MQ=255) ccaacaaTACAAAAATCCTCCCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGcc < 2:5284/89‑1 (MQ=255) aCATCCCCGCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAcgcgc < 2:405107/90‑1 (MQ=255) cATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAcgcgcg > 2:367913/1‑90 (MQ=255) tACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCggg > 1:185038/1‑90 (MQ=255) aCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTACGCGCGGGc > 2:61328/1‑90 (MQ=255) | GCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGC > NZ_CP009273/2678595‑2678748 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |