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% | 11.5 / 15.3 | 13 | Y70D (TAT→GAT) | glyA | serine hydroxymethyltransferase |
Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/7); total (5/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCG > NZ_CP009273/2678579‑2678745 | gTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAAcc < 2:481523/90‑1 (MQ=255) gTGCGGCTGGACGTTAGCGTAGCCAGCGCCGAACATTTCTTCCGCACGATCGACCGCCAGTTGTTCAACGATATCAACATCCCCGCAAcc < 1:188416/90‑1 (MQ=255) ggACGTTGGGGTAGTCGCCGCCAAACATTTTTTCCCCACGATCGACCGCCAGTTTTTCAACAATATAAACATCCCCCCAACCGCCgtagt < 1:64553/90‑1 (MQ=255) ggACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCgtagt > 1:227846/1‑90 (MQ=255) tAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGtt > 2:453514/1‑90 (MQ=255) caccgcCGAACTGTTTTTCCCCCCGATCGACCGCCTGTTTTCAAACAATATCAACATCCCCGCACCCGCCGTAGTAGCGTTTGCCCGGAt < 1:376928/87‑1 (MQ=255) cAGCGCCAAACTTTTTTTCCCCACGACGGATCGCCATTTGTACAACAATACCAACATCCCCACAACCGCCGTAGTAGCGTTTGCCCGGAt < 1:490155/90‑1 (MQ=255) gAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTc > 1:534664/1‑90 (MQ=255) aaCATTTTTTCCGCACGATCGCCCGCCTTTTTTTAAACAATATCAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCa < 1:190090/90‑1 (MQ=255) gCAGGATCGATCGCCATTTTTTCAACAATATCAACATCCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATttgttg < 2:365600/90‑1 (MQ=255) cACGATCGACCGCAAGTTGTTCAACAATATCAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTgg < 2:255303/90‑1 (MQ=255) atCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAc > 1:38005/1‑90 (MQ=255) cATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGACCCCTGCGCCTGCATTAcgcgcg > 1:587992/1‑90 (MQ=255) | GTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTCAACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCG > NZ_CP009273/2678579‑2678745 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |