| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,273,830 | 0 | A | C | 58.8% | 11.8 / 16.4 | 17 | T85T (ACA→ACC) | agaV | PTS N‑acetylgalactosamine transporter subunit IIB |
| Reads supporting (aligned to +/- strand): ref base A (7/0); new base C (0/10); total (7/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.14e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCC > NZ_CP009273/3273747‑3273905 | gTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACccgcc > 2:251053/1‑90 (MQ=255) cGCTGAAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGc < 1:223031/90‑1 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:726483/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 2:672765/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:212918/1‑90 (MQ=255) acaTTCCCCGCGCCGCCGATCGAAAAAAAATCCTGTTGGTTTAAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCgg < 2:737025/90‑1 (MQ=255) acaTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCgg > 1:737025/1‑90 (MQ=255) cccgcgCCCCCGATCGAAAAAAAATCCTGTTGTTTTTAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATc < 2:290586/88‑1 (MQ=255) cccgccCCGCCGACCGAAAAAAAATCCTGTTGTTTTAAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATc < 2:726483/88‑1 (MQ=255) cgccgcCGATCGACAAAAAATCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAt < 1:231121/90‑1 (MQ=255) gccgccGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAtt > 1:226514/1‑90 (MQ=255) gaccGAAAGAAAATCCTGCTGTTTTGTAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGtt < 2:446031/87‑1 (MQ=255) tCGAAAAAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg < 1:436852/90‑1 (MQ=255) aaaaaTCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat < 1:42834/89‑1 (MQ=255) ttCCTGTTGTTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCAc < 1:430074/89‑1 (MQ=255) ctgctgGTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTAc < 1:251053/90‑1 (MQ=255) ctgGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGcc > 1:279251/1‑90 (MQ=255) acacCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGtt > 1:89354/1‑60 (MQ=255) | GTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCC > NZ_CP009273/3273747‑3273905 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |