Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,435,727 | 0 | G | A | 34.3% | 46.8 / 19.3 | 35 | intergenic (‑108/+33) | uspF/ompN | universal stress protein F/outer membrane porin N |
Reads supporting (aligned to +/- strand): ref base G (21/2); new base A (0/12); total (21/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.09e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAAC > NC_000913/1435482‑1435939 | gtagGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGaaaaa < 1:6937/251‑1 (MQ=255) gggCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGcc > 1:154106/1‑251 (MQ=255) aGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGcc > 2:19567/1‑251 (MQ=255) gAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCAt > 2:40528/1‑250 (MQ=255) gTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTACTCAACGTGGCTAATCACGCGTTGAGTTACTTCTGAATCGGCAATATCGATAGGGACAAGAATCGTTCTGTTCATAACAACTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGAagg > 1:385528/1‑251 (MQ=255) cccctgaaccatcaATTTTTCCCTCTCCCTCAACGTGGCAAACCACGCTTTGAGTTAATTCGGAATCAAAAATATCGATAGGGAAAAAAATCGTTCTGTTCAAAAACCCTCTT‑TTTTTAGTATCCGAATAAAGTGTACCGCCATATAACATTTTTTTT‑GTAATGACGGAGTGT‑ACATTTTTATTTTTTAAAAAAGGA‑G‑AAGAAATAAGAAATAAGACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACt < 1:381883/242‑1 (MQ=255) ctGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACa > 1:259563/1‑226 (MQ=255) ttGCCTCTCCCTAAACGTGGCGAATCGCGCGTTGATTTATTTCTAAACCGAAAATATCGATAGGGAAAAAAATCGTTCTTTTAATAAAACCTCCT‑GTTTTAGTACCCGCAAAAAGTGTAACCCCAGATGACTCTTTTTGT‑GAAAGGAGGGAGTTC‑ACATTTTAATTTTAGAACAAAGGA‑GAAAAAA‑TAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt < 1:57025/247‑1 (MQ=255) ctTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACa > 1:507531/1‑250 (MQ=255) ccccaaaggTGGCTAATCGCGCGTTGATTAAATTCTGAATCGAAAATATCGATAGGGACAAGATTCTTTCTGTT‑AAAAAACCCCCTTGTTTTAGTATCCGAATAAAGTGTACCGCAAGATGACTTTTTTTGT‑GTAAGGAGGGAGTTC‑ACATTTTTAATTTAGAAAAAAGGA‑GAAAAAA‑TAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat < 1:222285/242‑1 (MQ=255) ccTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGACATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAatat > 2:101080/1‑251 (MQ=255) ccccgctttGAGTTAATTTAAAATCGAAAATATCGAAGGGGAAAAGACTCTTTCTTTTAATAAAACCTCCT‑TTTTTAGTACCCAAAAAAAGGGCACCGCCAGATGACTTTTTTTGT‑AAAATGGCGGAGTCC‑ACTTTTTTAATTAAAATAAAAGGA‑GGAAAAA‑AAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGc < 1:19567/244‑1 (MQ=255) aTTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt > 2:439982/1‑248 (MQ=255) aTTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt > 1:285139/1‑248 (MQ=255) cTGAATCGGAAATATCGATAGGGGCAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGc > 1:88250/1‑251 (MQ=255) atgtgtatgggacaaatattCGTTGTGTTCAAAAAACCCCCT‑GTTTTAGTCTCCACATAAAGTGTAACGCAAGATGACTTTTTTTGT‑GTAAGGCGGGATTAC‑ATATTTTTAATTTAGA‑AAAAGGAAGAAAAAA‑AAAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCa < 1:408563/232‑1 (MQ=255) cGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATcc > 1:134363/1‑250 (MQ=255) gATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGGA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa > 2:165022/1‑250 (MQ=255) aaGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAt > 1:416698/1‑251 (MQ=255) gAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAttt > 2:303111/1‑251 (MQ=255) aatctttctttaaaTAAAACCTCCT‑TTTTTAGTATCCACATAAAGTGAACCGCCATATGAC‑TTTTTTGTTGAAATGACGGATTTC‑ACATTTTTAATTTAGAAAAAAGGA‑GAAAAAA‑AAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt < 1:159592/238‑1 (MQ=255) tttAGTATCCGAATAAAGTGTACCCGCATATGACTCTTTTTTT‑GAAGTGACGGAGTCC‑ACTTTTTTAATTAAGATAAAAGGA‑GAAAAAA‑TAAAAAGAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCATCGTAGGTGGACATGtttt < 1:512103/250‑1 (MQ=255) tATCCACATAAAGGGTAACGCCAGATGACTCTTTTTTG‑GTAAGGAGGGAGTTC‑ACATTTTTAATTAAGATCAAAGGA‑GAAAAAA‑TAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTg < 1:375600/250‑1 (MQ=255) aaGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCACAGGC‑GGAAGAA‑TAAGCAGAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGt > 1:80786/1‑251 (MQ=255) agggaaaCGCCAGATGATTTTTTTTGT‑GTAGTGAGGGTGTAC‑AATTTTTTAATTTAAAACAAAGGA‑GGAAGAA‑TAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTc < 1:40528/246‑1 (MQ=255) gtAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAAt > 2:319522/1‑240 (MQ=255) gtAACGCCAGATAACCTTTTTTGG‑GAAAAGACGGAGTTC‑ACTTTTTAAATTAAAAAAAAAGGA‑GAAAAAA‑AAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAAt < 1:319522/240‑1 (MQ=255) aaccccaaaTGACACTTTTTGT‑GAAAGGAGGGAGTTC‑ACATTTTTAATTAAAAAAAAAGGG‑GAAAAAA‑AAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGcc < 1:60558/243‑1 (MQ=255) aCGCCAGGTGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAACAAGCCAGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCCTCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGcc > 2:446019/1‑250 (MQ=255) gATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCg > 1:359765/1‑251 (MQ=255) tGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCCGAACATGCCCGCCCTAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGCCAATATCATCGGTCGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGATGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGc > 1:377962/1‑251 (MQ=255) tttggaaTGACGGA‑TTCTACATTTTAAATTAAGAAAAAAGGA‑GGAAAAA‑TAAAAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTATAAGCTTTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGTTCATGTGGTTATGGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAAc < 1:568141/245‑1 (MQ=255) aCGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCttt > 1:225235/1‑250 (MQ=255) gggcgTTC‑ACATTTTTAATTTAAAGAAAAGGA‑GAAAAAA‑AAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTat < 1:104252/247‑1 (MQ=255) acaTTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAac > 1:243192/1‑250 (MQ=255) | GTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGCTAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCCT‑GTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGT‑GTAATGACGGAGTTC‑ACATTTTTAATTTAGATCAAAGGA‑GGAAGAA‑TAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAAC > NC_000913/1435482‑1435939 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |