Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,678,659 | 0 | A | C | 61.5% | 8.6 / 10.6 | 13 | Y70D (TAT→GAT) | glyA | serine hydroxymethyltransferase |
Reads supporting (aligned to +/- strand): ref base A (4/1); new base C (0/8); total (4/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.99e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTC‑AACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGC > NZ_CP009273/2678575‑2678748 | cGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTC‑AACGATATCAACATACTCGc > 1:967841/1‑90 (MQ=255) gTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTC‑AACGATATCAACATACTCGCAAcc > 1:733547/1‑90 (MQ=255) gTGCGGCTGGACGTTAGCGTAGTCAGCGCCAAACATTTCTTCCGCACGATCGACCGCCAGTTGTTC‑AACGATATCAACATCCCCGCAAcc < 1:408406/90‑1 (MQ=255) gTAGTCAGCGCCGAACATTTCTTCCGCACGATCGATCCCCAGTTTTTC‑AACGATACCAACATCCTCGCAACCGCCGTAGTAGCGTTTGcc < 2:608846/90‑1 (MQ=255) gcCGAACCTTTCTTCCGCACGACCGCTCG‑CATTTTTTCAAACAATAACAACATCCCCCCAACCGCCGTAGTAGCGTCTGCCCGGATAAcc < 1:726056/90‑1 (MQ=255) gAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTC‑AACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTc > 1:553513/1‑90 (MQ=255) aaaCATTTCTTCCCCCCGTCCGATCGCAGTTTGTCC‑AACAATAACAACATACCCGCACCCGCCGTAGTAGCGTTTGCCCGGATAACCTTc < 1:518091/89‑1 (MQ=255) cgatcgatCGCCATTTGTTC‑AACGATATCAACATCCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc < 2:463658/90‑1 (MQ=255) cgatcgatCGCCAGTTGTCC‑AACGATAAAAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTc < 2:733547/90‑1 (MQ=255) accgatcgCCATTTTTCC‑AACGATAACAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAg < 2:553513/88‑1 (MQ=255) gttgttC‑AACAATATCAACATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCt < 2:862493/90‑1 (MQ=255) gttC‑AACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTgcg > 1:243260/1‑90 (MQ=255) aCATCCCCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTAcgcgc < 2:967841/90‑1 (MQ=255) aCTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGc > 2:655938/1‑90 (MQ=255) | CGGAGTGCGGCTGGACGTTAGCGTAGTCAGCGCCGAACAGTTCTTTCGCACGATCGATCGCCAGTTGTTC‑AACGATATCAACATACTCGCAACCGCCGTAGTAGCGTTTGCCCGGATAACCTTCAGCATATTTGTTGGTCAGCTGAGAACCCTGCGCCTGCATTACGCGCGGGC > NZ_CP009273/2678575‑2678748 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |