Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 97,986 | 0 | T | C | 72.5% | 70.4 / 43.0 | 51 | R245R (CGT→CGC) | murC | UDP‑N‑acetylmuramate‑‑L‑alanine ligase |
Reads supporting (aligned to +/- strand): ref base T (9/5); new base C (18/19); total (27/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.63e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGATGCGCGTCACCCTGAATGCGCCAGG > NZ_CP009273/97857‑98120 | ctaccAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCa > 2:776621/1‑137 (MQ=255) ctaccAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCa < 1:776621/137‑1 (MQ=255) ccAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACt > 2:627287/1‑139 (MQ=255) gggCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCa < 2:839297/131‑1 (MQ=255) gggCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCa > 1:839297/1‑131 (MQ=255) gggCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCa > 1:264415/1‑131 (MQ=255) gggCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCa < 2:264415/131‑1 (MQ=255) gACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGAc > 1:298203/1‑125 (MQ=255) gACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGAc < 2:298203/125‑1 (MQ=255) aCTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCtt > 2:742958/1‑139 (MQ=255) cTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCtt > 2:92983/1‑138 (MQ=255) cTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCtt > 2:842674/1‑138 (MQ=255) tttGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCtt > 1:117920/1‑137 (MQ=255) agaATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGa < 1:70380/139‑1 (MQ=255) aaTTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAg < 1:543297/139‑1 (MQ=255) aaTTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAg < 2:282042/139‑1 (MQ=255) aCAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAc < 2:692377/139‑1 (MQ=255) aCAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAc < 1:128605/139‑1 (MQ=255) aCAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAc < 1:863992/139‑1 (MQ=255) cAGTCTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAcg < 1:122918/139‑1 (MQ=255) aGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAcgt < 2:374460/139‑1 (MQ=255) aCTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAcgtgc > 1:352627/1‑139 (MQ=255) ttttATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAcgtgcgt > 2:381271/1‑139 (MQ=255) tCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTAt > 1:668344/1‑139 (MQ=255) ccGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGcc < 1:627287/139‑1 (MQ=255) tACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAg < 1:92983/138‑1 (MQ=255) cGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAggg < 2:120856/138‑1 (MQ=255) tGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGCGTGGTGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACtt > 2:632021/1‑138 (MQ=255) tGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACtt > 2:451248/1‑138 (MQ=255) aTGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACgctgc > 1:620118/1‑139 (MQ=255) aTGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACgctgc > 1:114542/1‑139 (MQ=255) aTGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACgctgc < 2:196881/139‑1 (MQ=255) gtgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGctgctg > 2:818256/1‑139 (MQ=255) tgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATcag > 2:196840/1‑105 (MQ=255) tgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATcag < 1:196840/105‑1 (MQ=255) tgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCtgctgc > 1:465366/1‑139 (MQ=255) tgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAcgt > 2:467094/1‑86 (MQ=255) tgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGAcgt < 1:467094/86‑1 (MQ=255) tgtgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCtgctgc > 2:420068/1‑139 (MQ=255) tgtTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTgcgc > 2:566661/1‑139 (MQ=255) atCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCa < 2:272945/45‑1 (MQ=255) atCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCa > 1:272945/1‑45 (MQ=255) ccGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACa > 1:590421/1‑136 (MQ=255) ggTGATCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAagag > 1:662535/1‑139 (MQ=255) tCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGAt < 2:620118/139‑1 (MQ=255) tCCGCGAATTGTTACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGAt < 1:451248/139‑1 (MQ=255) ttACCGCGAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGATGCGCGTCAcc < 1:818256/138‑1 (MQ=255) tagcgcgAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGATGCGCGTCACCCt > 1:446399/4‑139 (MQ=255) gcgAGTGGGGCGCCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGATGCGCGTCACCCTGAAt < 2:114542/139‑1 (MQ=255) ggCGCCAGACCACGGCTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGATGCGCGTCACCCTGAATGCGCCAgg < 1:720578/139‑1 (MQ=255) ggCGCCAGACCACGACTTCCGGCTTCAGCGAAGATGCCGACGGGCGTGTAGAAGATTATCAGCGGATTGGCCCGCAGGGGCCCTTTACGCTGCTGGGCCAGGACAAAGAGCGGATGGGCGCCACCCTGAATGCCCCAgg > 2:720578/1‑139 (MQ=255) | CTACCAGGGCGACTTTGAGAATTTAAAACAGACTTTTATTAATTTTCTGCACAACCTGCCGTTTTACGGTCGTGCGGTGATGTGTGTTGATGATCCGGTGATCCGCGAATTGTTACCGCGAGTGGGGCGTCAGACCACGACTTACGGCTTCAGCGAAGATGCCGACGTGCGTGTAGAAGATTATCAGCAGATTGGCCCGCAGGGGCACTTTACGCTGCTGCGCCAGGACAAAGAGCCGATGCGCGTCACCCTGAATGCGCCAGG > NZ_CP009273/97857‑98120 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |