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% | 12.3 / 19.0 | 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. |
ATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAA > NZ_CP009273/3273741‑3273913 | aTCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAaca > 1:558455/1‑90 (MQ=255) gACGCTAAAAAAATTTATCAAAAACATTCCCCCCGCCGCCGATCGAAAAAAAATCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGAc < 2:524152/90‑1 (MQ=255) cGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGc > 1:412760/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:541186/1‑90 (MQ=255) acaactTTCCTCGCGCCGCCAATCGACAGAAAATCCTGCTGGTTTGTAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTc < 1:593336/90‑1 (MQ=255) acaTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCgg > 1:182779/1‑90 (MQ=255) acaTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCgg > 2:406015/1‑90 (MQ=255) ccccgcgCCGCCGATCGAAAAAAAATCCTGCTGGTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAAt < 1:286688/87‑1 (MQ=255) tCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCg > 2:603708/1‑90 (MQ=255) ccccgccgACCGAAAAAAAATCCTGTTGTTTTTAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAt < 2:179138/88‑1 (MQ=255) gccgccGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAtt > 2:174014/1‑90 (MQ=255) ccGACAAAAAACCCTGCTGTTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg < 1:557881/89‑1 (MQ=255) gACAGAAAACCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCa < 1:227697/90‑1 (MQ=255) gAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat < 1:115587/90‑1 (MQ=255) aaaaaCCCTGCTGTTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat < 1:573997/89‑1 (MQ=255) aTCCTGTTGTTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCAc < 1:453480/90‑1 (MQ=255) ctgctgGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTAc > 2:363181/1‑90 (MQ=255) aaaaaCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCaa < 2:541186/89‑1 (MQ=255) | ATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAA > NZ_CP009273/3273741‑3273913 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |