Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,273,830 | 0 | A | C | 54.5% | 38.8 / 49.4 | 33 | T85T (ACA→ACC) | agaV | PTS N‑acetylgalactosamine transporter subunit IIB |
Reads supporting (aligned to +/- strand): ref base A (15/0); new base C (0/18); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.64e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAA > NZ_CP009273/3273741‑3273917 | aTCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAaca > 1:646698/1‑90 (MQ=255) tCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAacac > 1:1438995/1‑90 (MQ=255) cGTCCGTTTGTGGCCGCTGAAAAAATTTATCAACAACATTCACCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGAAAACCCCccgc < 1:1492282/90‑1 (MQ=255) tGGACGCTGAAAAAAGTTATCGACAACTTTCACCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTg < 2:326655/90‑1 (MQ=255) gggccGCTGAAAAAAGTTATCAACAACATTCACCGCCCCCCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTg < 2:226598/86‑1 (MQ=255) gACGCTAAAAAAATTTATCAAAAACATTCCCCCCGCCGCCGATCGAAAAAAAATCCTGCTGTTTTGAAAAACCCCCGCCGATTTCCTGAc < 1:434834/90‑1 (MQ=255) cGCTGAAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTGACGc < 2:642053/90‑1 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 2:26868/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:144056/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 2:1454389/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:799969/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:803243/1‑90 (MQ=255) cTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCAGGTTTGTAAAACACCCGCCGATTTCCTGACGCTg > 1:1122708/1‑90 (MQ=255) tttATCGACAACACTCACCGCCCCGCCGATCGAAAAAAAATCCTGTTGTTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGt < 1:1454389/89‑1 (MQ=255) tttATCAAAAACACCCACCGCCCCGCCGACCGACAAAAAATCCTGCTGTTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGt < 1:1310997/89‑1 (MQ=255) ttATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTg > 2:1497415/1‑90 (MQ=255) tATCGACAACATTCACCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTgg < 2:618081/90‑1 (MQ=255) aTCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGc > 1:1133843/1‑90 (MQ=255) aTCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGc > 1:639863/1‑90 (MQ=255) aTCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGc > 1:46977/1‑90 (MQ=255) caacaTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTcc > 2:1565518/1‑90 (MQ=255) ccccgccCCCCCGATCGAAAAAAAATCCTGCTGGTTTTAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAAt < 2:1240832/87‑1 (MQ=255) cccgccCCGCCGACCGAAAGAAAATCCTGTTGGTTTTAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATc < 1:1483692/88‑1 (MQ=255) gccgccGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCAtt > 1:907332/1‑90 (MQ=255) ccgccgATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTa > 1:647578/1‑90 (MQ=255) gATCGAAAAAAAACCCTGGGGTTTTTAAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGtt < 1:493477/90‑1 (MQ=255) tagACAAAAAATCTTGTTGTTTTGTAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg < 2:1384072/88‑1 (MQ=255) tCGACAAAAAATCCTGTTGTTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg < 2:379908/90‑1 (MQ=255) tCGAAAAAAAATCCTGCTGTTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg < 2:647578/90‑1 (MQ=255) ccGACAAAAAATCCTGCTGTTTTGTAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg < 2:335440/89‑1 (MQ=255) gAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat > 1:1459449/1‑90 (MQ=255) aaaaaTCCTGTTGTTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat < 2:947497/89‑1 (MQ=255) aaaaaTCCTGCGGGTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAatat < 2:1122708/89‑1 (MQ=255) gTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAAt < 2:144056/90‑1 (MQ=255) tttGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATg > 2:1325656/1‑90 (MQ=255) acacCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAaacaa > 2:1223272/1‑90 (MQ=255) | ATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAA > NZ_CP009273/3273741‑3273917 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |