| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 377,474 | 0 | C | T | 73.9% | 19.1 / 10.7 | 23 | intergenic (+163/+61) | yaiL/frmB | DUF2058 domain‑containing protein YaiL/S‑formylglutathione hydrolase FrmB |
| Reads supporting (aligned to +/- strand): ref base C (1/5); new base T (12/5); total (13/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.15e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACT > NC_000913/377254‑377711 | gCAGGATGAAGACGACCCGTAGGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTTCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCg < 1:137955/251‑1 (MQ=255) aTGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACggg > 1:167035/1‑105 (MQ=255) aTGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACggg < 2:167035/105‑1 (MQ=255) gggtacaCCGGATAAGGCGTCCCCGCCCGC‑TCCGGCAGGCGGGCCC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTaa < 2:147800/245‑1 (MQ=255) atttgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTATTGATAtt > 1:534026/24‑251 (MQ=18) tgatgtggtaaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCAt > 1:513249/21‑251 (MQ=18) cccccccacccGGCCGTCGTGCGC‑CATGCCTGATGCGACGTTTTTCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGTGTTCACGCCGCATCCGGCATTCGTGCTATTATCAACGCATATTCAGTTTATTGTCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTTATAACGGATTAACCTCTTGATATTCATCTCCTgg < 2:246382/242‑1 (MQ=255) taaacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc > 1:144981/13‑251 (MQ=21) aacatagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGACGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGACAATATTCTCGCCAATAAAACTAGAGACAAAATAATAGCTGAGATCATAACCCTCGTAATAACGGCTTAACGACTTGGTATTCCTCTCCttggcg > 1:81996/11‑246 (MQ=17) catagataaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGc > 1:36749/9‑247 (MQ=25) taaGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCCCGGGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc > 1:352085/3‑251 (MQ=25) aaaggggtgcgcagaTGCCTTATGCGACGTTATGCGCGTTTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGGGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTc < 2:473379/237‑1 (MQ=25) aGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGGAGTGTAGACCGGCTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGAAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCCTAACCCTCGTGATAAAGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTcc > 2:187822/1‑249 (MQ=25) gTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTTTGATCCTAACCCTCGTGATAACGGATTAACGTCTTGATATTAATCTCCTGGCAGATCTTTTCa > 2:359808/1‑249 (MQ=25) cGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGAtt < 2:23646/251‑1 (MQ=255) gcgcagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTg < 1:323473/245‑1 (MQ=255) gcgcagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGACAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCACTAAAACTGGAAACAATATAATAACTGTGATCAATACCCTCGTGATAAAGGTTTAACGTCTTAATATTCATCTCCTGGCATATCCTTTCAAGAtttt > 2:212093/7‑250 (MQ=21) cagaTGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGACTTGATATTAATCTCCTGGCAGATCTTTTCAAGATTTgg > 1:32312/4‑249 (MQ=37) gACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGGGGTAGGCAAAAAGCTCGCCACTAAAACAGGAGACAAAATAATAGCAGTGATCATAACCCTCGTGATAACGGATTAACGGCTTGATATTCATCTCCTGGCAGGTCTTTTAAAAATTTTGAGTCCGCAGCGGtttt > 2:393610/1‑249 (MQ=25) ccgctattggcgttttttCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTTGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATACCCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTg < 2:470939/234‑1 (MQ=38) gTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAatca < 1:250256/251‑1 (MQ=255) ttagccggCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAatcatc < 2:257692/244‑1 (MQ=255) aTCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACt > 1:117054/1‑251 (MQ=255) | GCAGGATGAAGACGACCCGTATGCCGACTTCAAAGTGCCTGATGATTTGATGTGGTAAACATAGATAAGTCGTGCGCAGATGCCTGATGCGACGCTATGCGCGTCTTATCAGGCCTACGGTTTATGGGCGAAGTGTAGACCGGATAAGGCGTTCACGCC‑GCATCCGGCAGTCGTGCGC‑CATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCTACGGCTTACGGGCGAAGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCTATTATCAACGCATATTCAGTTTATTGGCGTGGTAGGCAATATGCTCGCCAATAAAACTGGAGACAAAATAATAGCTGTGATCATAACCCTCGTGATAACGGATTAACGTCTTGATATTCATCTCCTGGCAGATCTTTTCAAGATTTGGAGTCCGCAGCTGTTCTGCGTAAAAATCATCACT > NC_000913/377254‑377711 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |