Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 380,012 | 1 | . | G | 20.0% | 46.5 / 13.6 | 20 | intergenic (‑131/+57) | frmR/yaiO | regulator protein that represses frmRAB operon/outer membrane protein |
Reads supporting (aligned to +/- strand): ref base . (12/4); new base G (3/1); total (15/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGG > NC_000913/379915‑380126 | aGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc > 2:261531/1‑139 (MQ=255) gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAAGAAAGAGGAAAGGGGGGGGGGGATAATCGTAAGGGAAGGCCGCCCCGGATAgcg > 2:112887/1‑139 (MQ=255) gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc > 1:42245/1‑139 (MQ=255) tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGgccgcc > 2:74966/1‑123 (MQ=255) tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGgccgcc < 1:74966/123‑1 (MQ=255) tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc > 1:367843/1‑139 (MQ=255) gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc > 1:177153/1‑139 (MQ=255) aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa > 1:337890/1‑139 (MQ=255) aaTGCATACCCCCATATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa < 2:367843/139‑1 (MQ=255) aTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGaa > 1:178573/1‑139 (MQ=255) ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCa > 1:76726/1‑139 (MQ=255) ccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCa > 2:184238/1‑139 (MQ=255) tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGa < 1:261531/139‑1 (MQ=255) acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGTTAGCGGGCGGCAGAAGGAACCAGAATTTCCAGGTCAGACGGGCTGCAAGt > 2:484634/1‑139 (MQ=255) tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg > 2:484080/1‑139 (MQ=255) gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg > 2:292954/1‑139 (MQ=255) gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg > 2:320587/1‑139 (MQ=255) agaATAAAAATTATCCGGTGAATGGGGGCGGAAAATAAAGAGGAAAGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg < 1:484080/139‑1 (MQ=255) aaTAAAATTTATCCGATGAATGTGGTCGGAAAACAAACAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCACCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGAcc > 1:335610/1‑139 (MQ=255) aTGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg < 1:246060/139‑1 (MQ=255) | AGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGG > NC_000913/379915‑380126 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |