| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,020,300 | 0 | A | C | 53.3% | 14.2 / 18.4 | 15 | E391D (GAA→GAC) | guaD | guanine deaminase |
| Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/8); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTC > NZ_CP009273/3020215‑3020366 | cGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCa > 1:33211/1‑90 (MQ=255) cTCGGCGGAGCAAAATCTCGGGCCTTTGACGATTTGTTTGGCATCTTTTCACCTGGAAAAGAGGCTGTTTTCGTGGTGATGAACcccccc < 2:58375/90‑1 (MQ=255) aGCAAAATCTCGGGGCCTTGACGTTTTGATTGGCAACTTTTTCCCGGGCAAAGAGGCTGATTTCGTGGTGATGGACCCCACCGCCACTcc < 2:33211/90‑1 (MQ=255) ctcgGGGCTTTGACGATTTGATTGGCAATTTTTTCCCTGGAAAAGAGGCTGTTTTCGTGGTGATGGACCCCCCCGCCACTCCGCTAcagc < 1:380705/90‑1 (MQ=255) tGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGTTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACagcagc < 1:63782/90‑1 (MQ=255) ttGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt > 1:278630/1‑90 (MQ=255) ttGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt > 1:372551/1‑90 (MQ=255) tGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTa > 1:111342/1‑90 (MQ=255) tGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTa > 1:360438/1‑90 (MQ=255) ttttGATTGGCACCTTTTTACCGGGCAAAGAGGCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACa < 2:492843/89‑1 (MQ=255) aaGAGGTTTTTTTCGTGGTGATGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGt < 2:210934/90‑1 (MQ=255) gagGCTGTTTTGGTGGTGATGAACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc < 2:194264/90‑1 (MQ=255) gagGCTGTTTTCGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTc < 1:512497/90‑1 (MQ=255) ggtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt < 1:446233/41‑1 (MQ=255) ggtgATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCt > 2:446233/1‑41 (MQ=255) | CGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTC > NZ_CP009273/3020215‑3020366 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |