Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,041,298 | 0 | G | . | 26.7% | 33.0 / 10.1 | 15 | coding (1103/2874 nt) | gcvP | aminomethyl‑transferring glycine dehydrogenase |
Reads supporting (aligned to +/- strand): ref base G (5/6); new base . (2/2); total (7/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCA > NZ_CP009273/3041228‑3041382 | caAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGtt > 1:97831/1‑90 (MQ=255) caAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGtt > 2:139685/1‑90 (MQ=255) cGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATAc < 2:192885/90‑1 (MQ=255) tGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAgg < 2:138363/90‑1 (MQ=255) tGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATA‑‑‑‑‑‑‑‑‑‑‑‑TGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCg > 1:33779/1‑90 (MQ=255) tGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATA‑‑‑‑‑‑‑‑‑‑‑‑TGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCg > 1:351290/1‑90 (MQ=255) gcgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATA‑‑‑‑‑‑‑‑‑‑‑‑TGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAAc < 2:33779/90‑1 (MQ=255) gcgcAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATA‑‑‑‑‑‑‑‑‑‑‑‑TGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAAc < 2:351290/90‑1 (MQ=255) aGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCg < 2:23924/90‑1 (MQ=255) aGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATaa > 2:24390/1‑90 (MQ=255) tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATaaa > 1:221810/1‑86 (MQ=255) tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATaaa < 2:221810/86‑1 (MQ=255) tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATa > 1:257292/1‑84 (MQ=255) tttttGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATa < 2:257292/84‑1 (MQ=255) tatCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCcagca < 2:329813/90‑1 (MQ=255) | CAAGGTGTCGAAATAGTGCGCATGGCGCAGTTTCAGACCTTTTTGTTGCAGGCCCGCCGCCAGGATATCGGTCAGACGGTGAATGCGGTTAGCGATACGTTTCAGGCCAACCGGGCCGTGATAAACGGCATACAGGCTGGCGATGTTTGCCAGCA > NZ_CP009273/3041228‑3041382 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |