| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 2,626,314 | 0 | A | C | 56.2% | 14.9 / 19.7 | 16 | G372G (GGT→GGG) | guaB | IMP dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base A (7/0); new base C (0/9); total (7/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.74e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAG > NZ_CP009273/2626225‑2626401 | tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCa > 2:642185/1‑90 (MQ=255) aTCGCGCCCAGGGACCCCATCCCCCGGTAAGATTTGTAAGACCGGCCCTGGTAGAGTTCGTTTTCCCCCGGAGATTCTTCAGTACCCGcc < 1:55125/90‑1 (MQ=255) aTCGCGCCCAGGGACCCCATCCCACGGTAAGATTTGTAAGACCGGCCCTGGTAGAGTTCGATTTCCCCCGGAGATTCTTCAGTACCCGcc < 1:286080/90‑1 (MQ=255) gcCCAGGGAACAAATCCCACGGTAAGATTTGTAAGACCGCCCCTGGTAGAGTTCGATTTCCCCCGGAGATTCTTCAGTCCCCGCCAGCAt < 1:657460/90‑1 (MQ=255) gggACCCCATCCCACGGTAAGATTTGTAAGACCGGCCCTGGTAGAGTTCGTTTTCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAAc < 1:49018/90‑1 (MQ=255) gggAACCCATCCCACGGTAAGTTTTGTAAGACCGGCCCTGGTAGAGTTCGATTTCCCCCGGAGATTCTTCAGTCCCCGCCAGCATGGAAc < 1:382186/90‑1 (MQ=255) gggAACCCATCCCACGGTAAGTTTTGTAAGAACGGCCCTGGTAGAGTTCGTTTTCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAAc < 1:642185/90‑1 (MQ=255) aCCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAc > 2:206967/1‑90 (MQ=255) ggTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCgcgc > 1:406565/1‑90 (MQ=255) tgaaaGACCGGCCCTGGTAGAGTTCTTTTTCCCCCGGAGTTTCTTCAGTCCCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAg < 1:486694/87‑1 (MQ=255) tGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAg > 1:453211/1‑90 (MQ=255) aaCGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGa > 2:16800/1‑90 (MQ=255) aaCGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGa > 2:399193/1‑90 (MQ=255) ccGGCCCTGGTAGAGTTCGATTCCCCCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGAt < 1:546789/89‑1 (MQ=255) gCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGc > 1:276902/1‑90 (MQ=255) agTTCGTTTCCCCCGGGAGATTCTCCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATg < 1:167253/90‑1 (MQ=255) tCACCCGGAGATTCTTCAGTACCCGCCACCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGgag > 1:482325/1‑90 (MQ=255) | TAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCCAGCGGCGATAGCTTTGGCGATGTCGCCGGAG > NZ_CP009273/2626225‑2626401 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |