Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 380,012 | 1 | . | G | 40.0% | 1.6 / 44.7 | 25 | intergenic (‑131/+57) | frmR/yaiO | regulator protein that represses frmRAB operon/outer membrane protein |
Reads supporting (aligned to +/- strand): ref base . (4/11); new base G (9/1); total (13/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.61e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACC > NC_000913/379874‑380136 | cTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaaa > 1:642997/1‑139 (MQ=255) tataCTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGccc > 1:62343/1‑139 (MQ=255) taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggtggggttaatcggaaggaaaggccgccccggac > 2:1036057‑M1/1‑99 (MQ=255) taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATa > 1:967210/1‑139 (MQ=255) gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggggggctaatcggcagggaaggccgccccggatag > 1:645315‑M1/1‑96 (MQ=255) gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGAGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAgcgg > 1:809618/1‑139 (MQ=255) gggTATTCTCTATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggcaaaccggcaggaaagcccgcccgggaagcggg > 2:933445‑M1/1‑94 (MQ=255) gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccgcggatagcgg > 1:607862‑M1/1‑94 (MQ=255) gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgg > 1:178462‑M1/1‑94 (MQ=255) gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAgcggg > 1:437705/1‑139 (MQ=255) gTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgggc > 1:266048‑M1/1‑92 (MQ=255) tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg > 1:400248/1‑139 (MQ=255) atttgtcAAGGCATACCCCCCTATAGTATATAACACTGGAGAAAAAAATTTATCCGGTGAATGTGGGCGGAAAAAAAAAAGGAAAgggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcaga < 1:732183‑M1/136‑55 (MQ=255) gtaaatgcATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGCTGAATGTGGTCGGAAAACAAAGAGGAAAcggggggggggggctaatgggcggggaaggccgccccggaaagggggcggcagaagga > 2:700416‑M1/4‑81 (MQ=255) gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt > 2:920519/1‑139 (MQ=255) aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGTGGCAGAAGGAATCAg > 1:192108/1‑139 (MQ=255) tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAt > 1:882985/1‑139 (MQ=255) ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg < 2:62343/139‑1 (MQ=255) ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggt > 1:226949‑M1/1‑66 (MQ=255) tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgcccgggatagcgggcggcagaaggaatcagaatttccaggtca > 2:453755‑M1/1‑64 (MQ=255) tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg < 2:401510/138‑1 (MQ=255) tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg > 1:401510/1‑138 (MQ=255) gacagataaaaaactGGAGAATAAAATTTATCCGGGGAATTGGGGCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcag < 2:767075‑M1/134‑77 (MQ=255) taggttaaaaaactggggaaaaaaATTTATCCGGGGAATGGGGTGGGAAAACCAAG‑GGGGAgggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcaga < 2:226949‑M1/128‑79 (MQ=255) aaataacACGGGAAAATAAAATTTATCCGGGGAATGGGGTGGGAAAAAAAAGAGGGAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacggg < 2:48757‑M1/137‑82 (MQ=255) tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggc > 1:554752‑M1/1‑57 (MQ=255) ataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc < 2:809618/139‑1 (MQ=255) acacTGGAGAATAAAATTTATCCGGGGAATGTGGGCGGAAAACAAAGAGGAAAgggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggctgc < 1:700416‑M1/139‑87 (MQ=255) cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggggggcaaatcggaaggaaaggcgcccccggatagcgggcggcagaagaaataagaatttcaagtccggacgggctgaa > 2:387684‑M1/1‑52 (MQ=255) cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggctgcaa > 1:526447‑M1/1‑52 (MQ=255) tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggggggctaatgggcggggaaggccgccccggatagggggcggaggaaggaatcagaatttccaggtcagacgggccgcaagt > 1:363769‑M1/1‑49 (MQ=255) tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggcccccccggatagcggggggaagagggaataagaattcccggggcagccgggcggcaagtt > 2:921014‑M1/1‑49 (MQ=255) tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGTAAATGGGAAGGGAAGCCCCCCCGGGATAGGGGGGGCAAGAGGGAAACAAAATTCCCGGGTCAGGGGGGCTGCAAttt > 1:79285/1‑139 (MQ=255) tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg > 2:204895/1‑139 (MQ=255) ggAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa < 2:882985/139‑1 (MQ=255) gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggcccccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggctgcaactgtc > 1:767075‑M1/1‑47 (MQ=255) tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCg > 1:636051/1‑139 (MQ=255) aTCCGGTGAATGTGGTGGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatc < 1:920519/139‑1 (MQ=255) aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatc > 2:1202381/1‑139 (MQ=255) gAATGTGGTCGGAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggaaagcgggcggcagaaggaataagaatttccagggcagacgggctgaaagttgaagaccgttaaaaacatcggt > 2:961529‑M1/1‑27 (MQ=255) gAAAACAAAGAGGAAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggctgcaagttgcagaccgttaaaatcatcggttggggtgtcgt > 1:449571‑M1/1‑16 (MQ=255) aaaaaaAGAGGGAAggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggctgcaagttgcagaccgttaaaatcatcggttggggtgtcgtac < 1:11586‑M1/139‑126 (MQ=255) aaCAAAGAGGGAAgggggggggggggctaatcggcagggaaggccgccccggatagcgggcggcagaaggaatcagaatttccaggtcagacgggctgcaagttgcagaccgttaaaatcatcggttggggtgtcgtac < 1:387684‑M1/139‑127 (MQ=255) | CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCT‑AATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACC > NC_000913/379874‑380136 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |