Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,351,348 | 0 | T | G | 54.2% | 10.2 / 26.2 | 24 | H135P (CAC→CCC) | cadC | lysine decarboxylation/transport transcriptional activator CadC |
Reads supporting (aligned to +/- strand): ref base T (10/1); new base G (0/13); total (10/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.61e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAAAACCATACCCAAAAGGTAGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGCTGTACCA > NZ_CP009273/4351260‑4351436 | aaaaaaCCATACCCAAAAGGTAGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACtgtg < 2:84344/90‑1 (MQ=255) aaaaaaCCATACCCAAAAGGTAGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACtgtg < 1:250438/90‑1 (MQ=255) aGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAgg > 2:276086/1‑90 (MQ=255) aGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAgg > 1:488365/1‑90 (MQ=255) aGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAgg > 1:48725/1‑90 (MQ=255) aGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAgg > 1:301766/1‑90 (MQ=255) aGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGAAGGGAGAATCTGTGGCAgg > 2:358896/1‑90 (MQ=255) gAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAGGAAc > 1:57758/1‑90 (MQ=255) gTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAGGAACCGCct > 1:47596/1‑90 (MQ=255) ttGCTTTTAACGGGGGATTGTTCAGGTGGCGTT‑GCGGTTTTTTGAATGTTAAGACTGGGGGAGGGAAAATCTGTGGCAGGAACCGCctct < 1:276086/90‑1 (MQ=255) tgtttttAACGGGGGTTTGTTCGGGGGGCGTT‑GGGGTTTTTTGAATGTTAAGACTGGGGGAGGGGGAATCTGTGGCAGGAACCGCCTCTg < 2:202434/87‑1 (MQ=255) tttAACGGGGGATT‑TTTAGGGGGGGTTGGGGGTTTTTTGAATTTTAAGACTGTGGGGGGGAAAATCTGGGGCAGAACCCGCCTCTGGtat < 1:491515/90‑1 (MQ=255) gggggATTGTTCAGGTGGCGTT‑GGGGTGTTTTGAATGTTAAGACGGGGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAggggg < 2:148965/89‑1 (MQ=255) ggggATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAgggggg > 2:174392/1‑89 (MQ=255) tCAGGTGGCGTT‑GCGGTTTTTTGAATGTTAAGTCTGGGGGGGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGAt < 1:478924/90‑1 (MQ=255) gTGGGGTG‑GCGGTTTTTTGAATGTTAAGCCGGGGGGAGGGAGAATCTGTGGCAGCCCCCGCCTCTGGTATAGGGGGAGGCGAAGATAGCa < 2:414187/90‑1 (MQ=255) gggTTTTTTGATTGTTAAGCCTGGGGAGGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCct < 1:358896/89‑1 (MQ=255) ttttCAATTTAAAGACGGGGGGGGGAAAACCCTGGGGGAGGTCCCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTccc < 2:298086/90‑1 (MQ=255) aaTTTTAAGACGGGGGGGGGAAGAATTTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCtcttc < 2:142207/90‑1 (MQ=255) ttAAGACTGTGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATa > 1:77748/1‑64 (MQ=255) ttAAGACTGTGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCg > 1:210081/1‑90 (MQ=255) ttAAGACTGGGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATa < 2:77748/64‑1 (MQ=255) tAAGCCTGGGGGAGGGAAAATCTGTGGCAGGACCCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGc < 2:181569/90‑1 (MQ=255) tAAGAGTGGGGGGGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGc < 2:256532/90‑1 (MQ=255) tAAGACTGGGGGAGGGAGAATTTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGc < 2:476977/90‑1 (MQ=255) cTGGGGGGGGGAAAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGCTGTAc < 2:215507/90‑1 (MQ=255) gtgGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGCTGTACCa > 1:1483/1‑90 (MQ=255) | AAAAAACCATACCCAAAAGGTAGTGAATCGTTTGCTTTTAACTGGGGATTGTTCAGGTGGCGTT‑GCGGTGTTTTGAATGTTAAGACTGTGGGAGGGAGAATCTGTGGCAGGAACCGCCTCTGGTATAGGGGGAGGCGAAGATAGCATTATTTCCTCTCCCTCTTCTTCGCTGTACCA > NZ_CP009273/4351260‑4351436 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |