Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 546,689 | 0 | T | G | 55.2% | 30.3 / 35.9 | 30 | S265R (AGT→AGG) | ybcF | carbamate kinase |
Reads supporting (aligned to +/- strand): ref base T (2/11); new base G (16/0); total (19/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.46e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACC > NZ_CP009273/546605‑546773 | tGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTa < 2:284664/90‑1 (MQ=255) ccATGCCACACCGGGTGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGAGTGGGGATGTGAG‑aaa > 1:561256/1‑89 (MQ=255) aTGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑Gcc < 1:536525/90‑1 (MQ=255) aTGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑Gcc < 1:553117/90‑1 (MQ=255) tGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGGTATGTCAG‑AA‑GCCg > 1:375512/1‑90 (MQ=255) ggATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTGTGGCAGAAA‑GCCGGGGAAAAc > 2:536525/1‑90 (MQ=255) ggATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAAcc < 1:337131/90‑1 (MQ=255) ggATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGGTATGTTAG‑AAGGCCGGGGAAAAc > 1:677271/1‑90 (MQ=255) aTGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGTGGGGGTGTCAG‑AA‑GGCGTGGGAAAcccc > 1:580271/1‑89 (MQ=255) gTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGTTATGTTAG‑AA‑GCCGTGGGAAAACcccggg > 2:878290/1‑88 (MQ=255) tAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTTAA‑AA‑ACCGGGGGGAAACCGCGCggg > 1:451510/1‑89 (MQ=255) tAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGGGGTTTTGTCAG‑AA‑GCCCTGGGAAAACCCCggggg > 2:234159/1‑89 (MQ=255) gcgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTCAG‑AA‑ACCGGGGGAAAACCGCGTGGGtt > 2:3477/1‑90 (MQ=255) gcgcCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGAtt < 1:672020/90‑1 (MQ=255) aTTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGGGGGTATGTCAG‑AA‑GCCGTGGTAAACCCCCGGGGGTTggggg > 2:273529/1‑89 (MQ=255) aTTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGGTTTGTTAG‑AA‑GCCGGGGGAAAACCCCGtggggtggggg > 1:927898/1‑89 (MQ=255) aTTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGGGGGGGGGTATGTCAG‑AA‑ACCGGGGGAAAACCGCGGGGGTTggggg > 1:190851/1‑89 (MQ=255) ttGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGt < 2:78883/90‑1 (MQ=255) aaaGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATcg < 1:775554/90‑1 (MQ=255) gCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGGGGGTTTGTCAG‑AA‑GCCGTGGGAAAACCGCGGGGGTTGGGGGGTTTTTGCGa > 2:861809/1‑90 (MQ=255) ttCGATGGGGCCGAATGTAACGGCGGGGGGGGGGTATTTCAA‑AA‑GCCGGGGGAAACCCCCGGGGGTTGGGGGGTTATTGGGAATTGAaga > 1:506117/1‑90 (MQ=255) tCGATGGGGCCGAATATAACGGCGGTGGGGGGGTATGTCAG‑AA‑GCCGGGGGAAAACCCCGGGGGTTGGGGGGGTATTGCGAAATGAagag > 2:310526/1‑90 (MQ=255) cGATGGGGCCGAATGTAACGGCGGTGGGGGGGTATGTCAG‑AA‑GCCGGGGTAAAACCGCGGGGGTTGGGGGGTTATCGCGAATTGAAgaga > 2:373726/1‑90 (MQ=255) ggCCGAATGTAACGGCGGGGGGGGGGTATGTCAG‑AA‑ACCGGGGGGAAACCGCGGGGGTTGGGGGGGTTTTGCGAATTGAAAAGACGCCgg > 1:708602/1‑90 (MQ=255) ggCCGAATGTAACGGCGGAGGGAGGTTATTTGAG‑AA‑GCCGTGGCAAACCCGCGGGGATTGGGGCGGTTTGGCGAATTTAAGGGGGGCTgg > 2:580230/1‑90 (MQ=255) gCCGAATGTAACGGCGGTGGGGGGGTATGTTAG‑AA‑GGCGTGGGAAAACCGCGGGGGTTGGGGGGGTTTTGCGAATTGAAAAGACGCCggg > 2:286698/1‑89 (MQ=255) cggcggTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCgg < 2:555414/90‑1 (MQ=255) gcggTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCgggg < 2:646795/90‑1 (MQ=255) gTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGAcc < 2:34537/90‑1 (MQ=255) gTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGAcc < 2:2647/90‑1 (MQ=255) | TGCCATTCGCCATGCCACACCGGATGAGTTAGCGCCATTTGCCAAAGCCGATGGTTCGATGGGGCCGAATGTAACGGCGGTGAGTGGTTATGTCAG‑AA‑GCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACC > NZ_CP009273/546605‑546773 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |