Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,571,647 | 0 | T | C | 55.0% | 8.2 / 24.1 | 20 | E147G (GAA→GGA) | gntR | gluconate operon transcriptional repressor GntR |
Reads supporting (aligned to +/- strand): ref base T (7/2); new base C (0/11); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.64e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCG > NZ_CP009273/3571514‑3571784 | gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAATTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTccac < 1:151052/138‑1 (MQ=255) gtcgAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTccac > 2:151052/1‑138 (MQ=255) ggaaaTGGGGCGATGCCCGCGACAAATAATGGCGGGGGTCATCTGGCGTGCTGCTTAAAAGTTATCAAACCCGACGGCTATACCAGGGCAT‑GGCGCTTTGCTGCCCATCATTCCCCCCACGGGAATCCCCGCCACTTCa < 1:47406/136‑1 (MQ=255) cattgtgGCGAGGCCCGCGAGAAATAAGGGCGGTGGTCATCGGGCGTGCTGCTTAAAAGTTAAAAAACCCGACGGCGATATAAGGGCAG‑GGCGATTGGCTGTCCATCAGTCCCACCAGGGGAATACCCGCCACTTCAat < 1:240579/136‑1 (MQ=255) tgtgGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAatcat > 2:240579/1‑139 (MQ=255) gATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCttt > 1:380716/1‑137 (MQ=255) cgcgAAAAATAAGGGAAGTGGTCATCTGGCGTGCTCCTTAAAAGTTACCAAACCCGACGGCAATATCAAGGCAG‑GGCGATTTGCTGTCCATCAGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTgcg < 1:240611/139‑1 (MQ=255) ggagaAATAGGGGCATTGGTCATCTGGCGTGCGTCTAAAAATTAAAAAAACCCGACGGAGATA‑CAAGGCCTGGGCGATTTGCTGTCCATCAGTCCCCCCGGGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg < 2:399322/138‑1 (MQ=255) ggagaAATAATGGCGTTGGTTATCGGGCGTGCTGCTTAAAATTTACAAAAACCGACGGCGATATCAGGGCAT‑GGCGTCTTGCTGTCCATCGTTCCCACCCGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg < 1:396868/138‑1 (MQ=255) ggagaAATAATGGCAGTGGTCCTCGGGCGTGCTGCTCAAAAGTTATAAAAACCGACGGCGATATCAGGGCAG‑GGCGCCTTGCTGCCCAGCGGTCCCACCGGGGGAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg < 2:55848/138‑1 (MQ=255) caagaaataagGCCAGTGGTCATCTGGCTTGCTTCTACAAATTTAACAAACCCGACGGCGAAAACAGGCCTG‑GGCGACTTGCTGTCCATCATTCCCACCACGGAAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg < 1:34404/134‑1 (MQ=255) caaGCAATAGTGCCGGTGGTCTTGTGGCGGTATGCTTAAATTTTAGCAAACCCGCCGGCAATATCAGGTCGT‑GGCGCCTTGCTGTCCAGCGGTCCCACCCCGGGAATACCCGCCACTTCAATCATCTTTAAGGTGcgcg < 1:63929/137‑1 (MQ=255) ataaataatGGCGTTGGTCGTCTGGCGTGCTGTTAAAAATTTAAAAAACCCGCCGGCGATAACAAGCCAT‑GGCGATTTGCGTTCCATCTGTCCCCCCCCGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGcgcg < 2:413316/134‑1 (MQ=255) gTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCgg > 2:142641/1‑139 (MQ=255) ttATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTc > 2:305527/1‑138 (MQ=255) atCAAGGCGT‑GGCGTCTTGCTGCCCATCAGTTCCCCCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAgg < 1:44737/139‑1 (MQ=255) t‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGc > 2:323710/1‑139 (MQ=255) cTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCAtt > 2:32749/1‑139 (MQ=255) tCCACCAGTCCCCCCAGGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt < 2:275384/139‑1 (MQ=255) cATCGGTCCCCCCACGGGAATCCCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGtt < 2:100281/137‑1 (MQ=255) ttCCACCACGGGAATACCCGCCACTTCAATCATCTTTACGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGCGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCg > 2:93796/1‑139 (MQ=255) | GTCGAGACGTGCGCCGAGATAGGCAATGTGGCGATGCCCGCGAGCAATAATGGCAGTGGTCATCTGGCGTGCTGCTTCAAAGTTATCAAAACCGACGGCGATATCAAGGCAT‑GGCGACTTGCTGTCCATCAGTTCCACCACGGGAATACCCGCCACTTCAATCATCTTTAAGGTGCGCGGCGTGTGGGTACGTTCGGTGAGGATCAGGCCGTCGATATTCCAGGAGAGCATGGATTCGAGGCGTTCTTGCTCCATTTCCGGTTTATAACCG > NZ_CP009273/3571514‑3571784 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |