Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,911 | 0 | G | A | 35.3% | 42.9 / 14.0 | 34 | intergenic (‑122/+1) | fliA/fliX | RNA polymerase sigma factor FliA/small RNA FliX |
Reads supporting (aligned to +/- strand): ref base G (16/6); new base A (0/12); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.72e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT > NC_000913/2001664‑2002155 | gTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGaa > 2:427273/1‑250 (MQ=255) gTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTTTCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAAc < 1:272367/251‑1 (MQ=255) cACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCcc‑g > 1:121485/1‑251 (MQ=255) cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCcc‑gcc < 1:362828/251‑1 (MQ=255) aggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATATCGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTgg > 2:47252/2‑249 (MQ=255) gcttcgggcccccccaggggaaCATAACGCTGCACCAGCG‑GTTTTTTTCCCATTACCCCTTCGGGGGTAAAGAGTGAATTCACAATAAACACCCCGGCTTTATATGAGTTATCGGCATTATTTCCCGTTT‑CGGCAGGGTTTTT‑AA‑TCGGACTATTGGTGGGGAAAATGAGGGGTTATTGGGGGG‑TTAAAGGTAAATTCCAGGAAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc < 1:245788/230‑1 (MQ=255) gcgggACGCCCCAGCGGGACATAACGTTGCCACAGGGATTTTTTACCC‑TTTCCCCTTTCAGGGGTATAGAGTGAATACACGAAAAACAGCCCTGTTTTATATGATTTACCGGCATGATTATCCTTTT‑CTGC‑GGGTTTTT‑AAAGCGGACGATTAGGGGGAGAAAGGGGGGTTTATTGGGGGG‑TTAGGGGTAAATTCAGGGAAGAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc < 1:321653/246‑1 (MQ=255) cGCACCAGGGGGACATACCGCTGCCCCAGCGATTTTTTACCC‑ATTCCCCCTTCGGGGTTATAGAGTAAATTCACGAAAAACAGCCCTGCTTTATATGAGTTATCGGAATGTTTTTCCTTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTGGGGGGTGAAATGGGGGTTTATTGGGGGG‑TTACAGGTAAATCCCAGGCAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt < 1:7369/251‑1 (MQ=255) gCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt > 2:341042/1‑251 (MQ=255) gaacataccggtgcaaagaggtttttttACCC‑ATTACCCCTCCAGCGGTATAGGGTAAATTCACGATAAACAGCCCGGCTTTATATGTGTTATGGGAATGTTTATCCGTTT‑CGGCAGGGTTTTA‑AA‑TCGGACTATTAGTGGGTGAAATGGGGGTTTTTTTGGGGG‑TTACAGGAAAATTCCGGGAAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc < 2:97911/226‑1 (MQ=255) ccACAGCGAGGGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGGTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt > 2:58856/1‑251 (MQ=255) cacaGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt > 1:361517/1‑251 (MQ=255) gaggttttttttcccttcctccTTCAGCGTTAGAGAGTGATTTCACAATAAACAGCCCTGCGTTATATTATTTATCGGCATGATTATC‑GTTTGTTGCGGGGTTTTT‑AA‑TGGGAGGATGGGGGGGGGAAA‑GGGGGGTTTTTTGGGGGGTAACTGTTAAATCCCAGGCAAAAAAACACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc < 2:473693/231‑1 (MQ=255) tttttccc‑ATAACACCTTCGGCGGTGAAGAGTCAATTCAGGATAAACAGCCTGGGGTTATTTGGGTTAGCGGTATGATTCTCCGTTC‑CTGAAGGGTTTTA‑AT‑TCGAAGGATTGGTGGG‑GAAATTGGGGGTTTTTTGGGGGGTTACAGTAAAATCCCAGGCAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt < 1:369329/244‑1 (MQ=255) gcgggggtttagagtaaagacaggATAAACAGCCTTGTGTTTTATGATTTAGGGGCATGTTTATCGGTTT‑CTG‑AGGGTTTTTTAA‑TGGGACGGTTAGTGGGTAAAAGGGGGGTTTTTTTGGGGT‑TTAGAAAAAAATTCCCGCCAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc < 2:93291/228‑1 (MQ=255) gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGGAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATTTGCAAATTCAATAAATTGc > 2:402959/1‑250 (MQ=255) gatttCACGAAAAACGGCCCTGGTTTATATGATTTATCGGCATGTTTATCCTTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGTTTAGTGGGTGAAATGAGGGGTTTTTGGGGGG‑TTAAAGTAAAATTCCGGGCAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt < 1:261562/247‑1 (MQ=255) aaaTTCACAAAAAACAGCCCGGGGTTAAATGAGTTATCGGCATGTTTGTCCGTTT‑C‑GGAGGGTTTTT‑AAATGGGAGGTTTAGGGGGTGAAATGGGGGGTTATTGGGGGT‑TTACAGGAAAATCCCAGAC‑GAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc < 1:141459/250‑1 (MQ=255) cACGATAAACAGCCCTGCTTTATAGGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGAGAAATGAGGGTTTTTTGGGGGG‑TTACAGGTAAATCCCAGGCAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt < 1:625093/250‑1 (MQ=255) cAGCCCTGCGTTATATGCGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg > 2:172066/1‑201 (MQ=255) cTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa > 2:464890/1‑249 (MQ=255) tatGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc > 1:617917/1‑251 (MQ=255) gAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc > 1:623709/1‑251 (MQ=255) gAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGGTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc > 2:198527/1‑251 (MQ=255) tATCGGCATGATTATCCTTTT‑CTGCGGGTTTTTT‑AA‑TCGGAGAATTGGTGGGTAAAATGAGGGTTTTTTTGGGGT‑TTACAGAAAAATTCCAGACAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgn < 2:285839/250‑2 (MQ=37) tCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa > 2:297372/1‑250 (MQ=255) cATGATTATCCGATT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCCGCGCATCAGGCAATTTg > 2:554571/1‑251 (MQ=255) tttcctttt‑CTG‑AGGGTTTTTTAA‑TCGGACGTTTAGTGGGTGAAA‑GGGGGGTTTTTTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGtt < 1:364587/244‑1 (MQ=25) gTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGGGGGTGAACTGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAACACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTAGCCGGATGCGGCGTAAACTCCTTATCCGGCCTACAAAAATGCGCAAATTCAATAAATTGCAATTCCACTTGTCGGCCTGATAAGCGCAGTGCATCAGGCAATTAGGCGTTGCCGTc > 2:452216/1‑251 (MQ=255) tttggCAGGTTTTTT‑AA‑GCGGACGATGAGTGGGTAAAAGGAGGGTTTATTTGGGGA‑TGAAAGGTAAATCCCAGGAAAAAAAACACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCa < 2:95553/246‑1 (MQ=37) ggaggtttttt‑AA‑GCGGATGATTGGTGGGTAAAAGGGGGTGTTTTTTGGGTG‑TTACGGGAAATTTCGAGGAAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGtctc < 2:190789/245‑1 (MQ=25) gAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGGTGCTGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCACTTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAAccccgc > 1:159732/1‑248 (MQ=255) tttttgggggg‑TTACAGGAAAATTCCGGGAAGAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACa < 1:427273/244‑1 (MQ=255) cAGGCAAAAAAAACCCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt < 1:341042/251‑1 (MQ=255) | GTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCC‑ATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT‑CTGCAGGGTTTTT‑AA‑TCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCC‑GCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCT > NC_000913/2001664‑2002155 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |