Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 753,125 | 0 | T | A | 40.7% | 30.3 / 17.4 | 27 | intergenic (‑330/+60) | ybgD/gltA | putative fimbrial‑like adhesin protein/citrate synthase |
Reads supporting (aligned to +/- strand): ref base T (14/2); new base A (0/11); total (14/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.98e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAG > NC_000913/752956‑753339 | aGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGGTAGAAATATTTTC‑ACATAGCAATCAAACATCAAAGCATg > 2:19881/1‑239 (MQ=255) cGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTTTAAATTATAAATCAACCCTTAAAGCTTGCAATCGCCTTTTaaatc > 2:190018/1‑237 (MQ=255) tCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGGTA‑AAATATTTTAAACTTAGTAATAAACCATTAACGCTTGATATAG‑CTTTTTAAATcca > 2:315730/1‑237 (MQ=255) ttttaAATTGAAAGTTTTGGGTTGCTGATAATTTGAGCTTTTCTATTCTTTTTAAATTTCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGc < 1:315730/239‑1 (MQ=255) tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATTAACCATTAACGCTTGATATCG‑CTTTTTAAGTCGCGTTTTTCTTATCCCTGT‑TAACGCTGAAGCGGACGGGCaa > 2:708276/1‑239 (MQ=255) aTTGAAATTATTGGTTTGCTGATAATTGGAGTTTTTTTATTTTTTTAAATTTTCTAAATAGGTTTGTTAATGGATTTTTTTTTTAAAT‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCtt < 1:190018/239‑1 (MQ=255) ttatataGTTTTGTAAATGTATTTTA‑TTTTACAATTTTTTTGTTTTT‑AGGCAAAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGt < 2:248120/144‑1 (MQ=255) tatataGCTATGTAAATGATTTTTATTTTTAAAT‑TTTTTTGTTTTT‑AGGCAAAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGc < 2:435047/179‑1 (MQ=255) tatgGGTTTGTAAATGTATTTTTTTTTAAAAT‑TTTTTTTTGTTT‑AGGAAAAAAAAAACCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTATTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCa < 2:663371/185‑1 (MQ=255) tataGGTTTGTTAATGAATTTTATTTTTAAAT‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCa < 1:558902/222‑1 (MQ=255) tataGGTCTTTTAATGGTTTTTTTTTTTACAT‑TTTTTTGTTTTT‑AGGAATAAAAAAACCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGc < 2:199198/177‑1 (MQ=255) gttttttAATTGATTTTTATTTTAAAA‑TTTTTTTGGTTTAAGAAATAAAAAAATCA‑‑CCCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCa < 2:319214/213‑1 (MQ=255) gtttGTAAAGGGTTTTTATTTTAAAAT‑TTTTTTTTTTTT‑AGGAATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGcc < 1:618352/217‑1 (MQ=255) gTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCa > 1:416479/1‑188 (MQ=255) cTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑c > 1:248120/1‑145 (MQ=255) tGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATa > 1:435047/1‑179 (MQ=255) ttAATGGTTTTTATTTTTACAT‑TTTTTTGTTTTT‑AGGAATAAAAAAACCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCat < 2:168811/139‑1 (MQ=255) ttAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACt > 1:663371/1‑185 (MQ=255) ttAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGGTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATa > 1:199198/1‑177 (MQ=255) ttAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAATCGCGGTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAAACTTCATACCGTCACTTGGCATTTTGCTCCAATGGGCGATCCACCCAACGGTACGTGCCATTTCGGAAAt > 2:306089/1‑239 (MQ=255) ggATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTT‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTTAAGTCGCGGTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCACTCTTCATAACGTCACTGTGCATTTTGCTCCAGTGGGCGGTCCAGCCAACCGTACGt > 2:618352/1‑220 (MQ=255) ggATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGGGATCCAGCCAACggga > 1:319214/1‑215 (MQ=255) ttattttTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGtata > 1:168811/1‑139 (MQ=255) ttattttTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAg > 2:193542/1‑142 (MQ=255) ttattttTACAA‑TTTTTTGTGTTT‑AGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAg < 1:193542/142‑1 (MQ=255) tttAAAT‑TTTTTTGTTTTT‑AGGCATAAAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGa < 1:306089/239‑1 (MQ=255) aCAA‑TTTTTTGTGTTA‑AAG‑ATATAAAAAATATCACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCCCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAg < 2:558843/239‑1 (MQ=255) | AGGGGGGTATAGATAGACGTCATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAA‑TTTTTTGTGTTT‑AGGCATATAAAAATCA‑‑ACCCGCCATATGAACGGCGGGTTA‑AAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCG‑CTTTTAAAGTCGCGTTTTTCATAT‑CCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCATTTCGCTCCAGTGGGCGATCCAGCCAACGGTACGTGCCATTGCGAAAATGACGGTGAACATGGAAG > NC_000913/752956‑753339 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |