Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,368,619 | 0 | A | C | 54.3% | 31.9 / 50.8 | 35 | V48G (GTA→GGA) | ampC | BlaEC family class C beta‑lactamase |
Reads supporting (aligned to +/- strand): ref base A (16/0); new base C (0/19); total (16/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.46e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAG > NZ_CP009273/4368533‑4368700 | cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑accg > 1:402751/1‑90 (MQ=255) aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑accgcc > 2:993494/1‑90 (MQ=255) aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑accgcc > 2:1195530/1‑90 (MQ=255) aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑accgcc > 1:179431/1‑90 (MQ=255) tgtgGCGGGCTTTTTTTTGGCGATGCCCCCAACCCCCCAGGAAAAGTAAAAAGTTTCCCCC‑TAAAAAA‑TT‑CCCCCCA‑CCGCCATACCCgg < 1:1017188/90‑1 (MQ=255) tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCgg > 1:288402/1‑90 (MQ=255) ggaggggCTGTTTTTTGGGGATTCCCGCATCGCCCCGGTTAAAGAAATAGGTTTTCCCC‑TAATAAA‑TT‑CCCGCCC‑CCGCCATACCCGGGa < 2:169463/86‑1 (MQ=255) cGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGGAAAAGTAATAAGGTTTCCCC‑TGAAAAA‑TT‑CCCGCCA‑CCGCCATACCCGGGATCt < 1:322858/90‑1 (MQ=255) ctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTttg > 2:138880/1‑90 (MQ=255) ctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTttg > 2:777332/1‑90 (MQ=255) tttttttgGGCGATCCCCACATACCCCAAGAAAAAAAAAAAAGTTTTCCCC‑TAAAAAA‑TC‑CCCCCCC‑CCGCCATACCCGGGATCTTttgt < 2:560475/87‑1 (MQ=255) gttttttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTCCCC‑TGAAAAA‑TT‑CCCGCCA‑CCGCCGTACCCGGGATCTTttgtt < 1:364334/88‑1 (MQ=255) tttttgGGCGATGCCCGAATACCCCCGGGAAAAGTAATAAGTTTTCCCC‑TGAAAAA‑TT‑CCCGCCA‑CCGCCATACCCGGGATCTTttgttg < 1:79953/89‑1 (MQ=255) tttttGGCGATTCCCCCATCCCCCAAGGTAAATAAATATGTTTCCCCT‑TAAAAAA‑TC‑CCCCCCC‑CCGCCATACCCGGGATCTTTTGTTGc < 2:848966/90‑1 (MQ=255) tttttGGCGATGCCCGCATCCCCCGAGAAAAAGAAATATGTTTTCCCC‑TAAAAAA‑TC‑CCCCCCC‑CCGCCATACCCGGGATCTTTTGTTGc < 2:1180504/90‑1 (MQ=255) ttttgGGCGTTGCCCGCATCCCCCAAGAAAAAAAAAAAAGTTTTCCCT‑TAAAAAA‑TC‑CCCCCCA‑CCGCCATACCCGGGATCTTTTGTTGc < 2:309352/90‑1 (MQ=255) ttttgGGCGTTGCCCCCATA‑CCCCCGGGAAAATAAAAAGGTTTTCCC‑TTAAAAATTT‑CCCCCCA‑CCGCCATACCCGGGATCTTTTGTTGc < 1:138880/90‑1 (MQ=255) gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAg > 1:431816/1‑90 (MQ=255) aTGCCCGCATCCCCCAAGGAAAAGAAAAAAGTTTTCCCC‑TAAAAAA‑TT‑CCCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc < 1:492557/90‑1 (MQ=255) ttgccGCAAGGCCCAAAGTAAAAAAAAAATGTTTTCCCTTAAAAAA‑TC‑CCCCCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc < 2:936382/87‑1 (MQ=255) gCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt > 2:153290/1‑90 (MQ=255) ccgaGGAAAAATAAAAATTTTTCCCC‑TAAAAAA‑TT‑CCCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc < 1:1080376/87‑1 (MQ=255) cAGGAAAATTAATAAGTTTTCCCC‑TAAAAAA‑TT‑CCCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa < 1:9477/90‑1 (MQ=255) aaGGAAAAGTAATAGGTTTTCCCC‑TAAAAAA‑TC‑CCCCCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa < 2:538448/89‑1 (MQ=255) cagaaaaaTAAAAAGGTTTTCCC‑TTAAAAA‑TTCCCCCCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa < 2:589607/86‑1 (MQ=255) gTAAAGAAAAATGTTTTCCCC‑TAAAAAT‑TC‑CCCGCCG‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc < 1:1053642/90‑1 (MQ=255) aaaaaTAAAAAGTTTTTCCC‑TTAAAAA‑AT‑CCCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc < 2:306681/89‑1 (MQ=255) aaaGTAATAAGTTTCACCT‑TAAAAAA‑TC‑CCCCCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac < 1:1111526/90‑1 (MQ=255) aaGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca > 1:167001/1‑90 (MQ=255) aaGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca > 1:1150506/1‑90 (MQ=255) aaGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca > 1:115657/1‑90 (MQ=255) aaGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCAGGGTAATTGTGCGATGcaca > 2:1182208/1‑90 (MQ=255) taaGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg > 1:1198046/1‑90 (MQ=255) ttACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt > 1:425469/1‑90 (MQ=255) ttACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt > 1:227868/1‑90 (MQ=255) ttACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt > 1:275882/1‑90 (MQ=255) ttACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt > 1:98846/1‑90 (MQ=255) ttACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt > 1:365281/1‑90 (MQ=255) gATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg > 1:612841/1‑90 (MQ=255) | CAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TGATAAA‑TT‑ACCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAG > NZ_CP009273/4368533‑4368700 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |