Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,360,813 | 0 | C | T | 55.0% | 9.1 / 20.9 | 20 | F88F (TTC→TTT) | arnC | undecaprenyl‑phosphate 4‑deoxy‑4‑formamido‑L‑arabinose transferase |
Reads supporting (aligned to +/- strand): ref base C (2/7); new base T (11/0); total (13/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.64e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAATTATTA‑CCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCG > NZ_CP009273/2360728‑2360884 | aaGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTc < 1:177717/90‑1 (MQ=255) tCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCAC‑ggt > 1:314634/1‑88 (MQ=255) tCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGGCAC‑Ctt > 1:745857/1‑90 (MQ=255) aaCAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTa < 2:681886/90‑1 (MQ=255) aGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAatt < 1:135909/90‑1 (MQ=255) tgtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTTACGGTTCT‑TGGGGCCTTATTTTTTA‑CCttt > 1:299700/1‑90 (MQ=255) gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTTGGCAC‑CTTAC‑TGGGGGCCTAATTTTTA‑TCCTTgt > 1:904520/1‑89 (MQ=255) gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCAC‑GTTTC‑TGGCGGCTTAATTATTA‑CCCttta > 2:148476/1‑88 (MQ=255) gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCAC‑CTTAT‑TGGCGGCTTAATTTTTA‑TCCttta > 1:105135/1‑88 (MQ=255) gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTTAGTCCC‑AGTAC‑TGGCGGCTTAATTTTTATCCCtta > 2:643974/1‑89 (MQ=255) gtCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTTAGTCAC‑GTTAC‑TGGCGGCTTAATTATTT‑TCCttta > 2:236690/1‑88 (MQ=255) tttGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAATTATTA‑CCCTTGATGCCGa < 1:643974/90‑1 (MQ=255) tttGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAATTATTA‑CCCTTGATGCCGa < 2:122816/90‑1 (MQ=255) tttGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAATTATTA‑CCCTTGATGCCGa < 2:629362/90‑1 (MQ=255) tGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCAC‑GTTACTTGGGGCTTTATTTTTTA‑CCCTT‑TTGCCGAtc > 2:921088/1‑90 (MQ=255) aCCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCCGTAAC‑AATAC‑TGGCGGCTTAATTATTA‑ACCTTTTTTTCGCTCTCCAGa > 2:92927/1‑90 (MQ=255) cGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAATTATTA‑CCCTTGATGCCGATCTCCAGAAt < 1:830038/90‑1 (MQ=255) gggCAACATTCAGCGATTATGGCGGGGTTTAGTCAC‑GTTAC‑TGGCGGCCTAATTTTTT‑CCCTTTTTTCCGCTCTCCCGAAACCCCCcgaa > 1:55879/1‑90 (MQ=255) cAGCGATTATGGCGGGGTTTTGTGAC‑CTTAC‑TGGGGGCGTAAATTTTA‑ACCTCGCTGGCGCGCCCCACAAACCCCCCGAAGAAAAccccc > 1:593829/1‑90 (MQ=255) aGCGATTATGGCGGGAGTCCGTCCC‑GTTAC‑TGGCGACGTAATTATTT‑CCCCTTGTGCCGATTCCCCGAAACCGCCCGAAGAAAAcccccc > 2:639535/1‑89 (MQ=255) | AAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCAC‑GTTAC‑TGGCGACTTAATTATTA‑CCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCG > NZ_CP009273/2360728‑2360884 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |