Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,294,730 | 1 | . | G | 33.3% | 8.9 / 13.9 | 12 | coding (1503/2487 nt) | napA | periplasmic nitrate reductase subunit alpha |
Reads supporting (aligned to +/- strand): ref base . (5/3); new base G (2/2); total (7/5) | |||||||||||
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.20e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGT > NZ_CP009273/2294650‑2294817 | cacaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATc < 1:350794/86‑1 (MQ=255) cacaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATc > 2:350794/1‑86 (MQ=255) acaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGGATCgcgc > 1:164648/1‑90 (MQ=255) acaTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGGATCgcgc > 1:363990/1‑90 (MQ=255) cGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATCGCGCCAGCCCg > 1:31915/1‑90 (MQ=255) gATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATCGCGCCAGCCCGGCATACGctc < 1:110707/90‑1 (MQ=255) tgactgTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCt > 1:204298/1‑90 (MQ=255) tACGGATCGGAGACGATGATGAAGTTGCGCGGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATgttgtt < 2:164648/90‑1 (MQ=255) tACGGATCGGAGACGATGATGAAGTTGCGCGGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATgttgtt < 2:363990/90‑1 (MQ=255) ggAGACGATGATGAAGTTGCGC‑GGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACAt < 2:266429/90‑1 (MQ=255) gTTGCGC‑GGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAtt > 2:326696/1‑90 (MQ=255) tGCGC‑GGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCa > 2:51507/1‑90 (MQ=255) cgc‑gGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGt > 1:48331/1‑90 (MQ=255) | CACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGC‑GGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGT > NZ_CP009273/2294650‑2294817 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |