| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,678,659 | 0 | A | C | 53.8% | 10.4 / 13.9 | 13 | Y70D (TAT→GAT) | glyA | serine hydroxymethyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/7); total (6/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGC > NZ_CP009273/2678574‑2678731 | ccGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGGTCTTTCGCACGATCGATCGCCAGTTGTTCAGCGATATCAACATACTCg > 2:350765/1‑90 (MQ=255) gTGCGGCTGGACGTTAGCGTAGTCAGCGCCAAACATTTTTTTCCCACGATCGATCGCCAGTTGTTCAACGATACCAACATCCTCGCAAcc < 1:534843/90‑1 (MQ=255) aCGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGtagtag > 1:469755/1‑90 (MQ=255) aCGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGtagtag > 2:489762/1‑90 (MQ=255) gTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCgg > 1:204184/1‑90 (MQ=255) gcgcCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATaa > 1:365139/1‑90 (MQ=255) gAACATTTCTTTCCCCCGATCGATCGCCAGTTGTTCAACGATATCAAAATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTc < 1:454946/90‑1 (MQ=255) cACGTTCGATCGCTATTTTTCAAATAATATAAACATCCCCGCAGCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTgg < 2:44380/90‑1 (MQ=255) cgatcgatCGCCATTTTTTCAACGATATAAACATCCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc < 1:468201/90‑1 (MQ=255) cgatcgatCGCCAGTTGTTCAACGATATCAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc < 1:167938/90‑1 (MQ=255) gttCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTgcg > 1:551013/1‑90 (MQ=255) tccAACGATAACACCATCCCCCCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTgcgc < 2:469755/88‑1 (MQ=255) tacaacaaTAACAACACCCCCGCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTgcgc < 2:204184/88‑1 (MQ=255) | CCGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGC > NZ_CP009273/2678574‑2678731 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |