Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,622,816 | 0 | G | A | 41.0% | 31.8 / 35.3 | 40 | intergenic (+75/+144) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
Reads supporting (aligned to +/- strand): ref base G (22/1); new base A (0/16); total (22/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.51e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA > NC_000913/1622668‑1622964 | ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGCTAATCCCTCACCTAACCGGCCCCTTGTTAAAGTTGGGTACAAGGGGC‑C‑TGATTTTTATAACGGCGa > 1:915119/1‑150 (MQ=255) ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGa > 2:573020/1‑150 (MQ=255) aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGaaa > 2:158264/1‑150 (MQ=255) ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGaaaaa > 2:426570/1‑150 (MQ=255) gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGaaaaaa > 1:459843/1‑150 (MQ=255) aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCg < 1:136018/150‑1 (MQ=255) ttttttttgtttttctggCCGCGTTTTCCACCCAAAAAGGGAATGACAAAGAAGCGGATAAAACCCCCCACCAACCCGGCCTTTTTTCATTGTAGTGTA‑AAGGGGCTG‑TATTTTTTGTGCGGGCAAAAAAAAACCGCCAGTAAACcggcg < 2:607327/148‑1 (MQ=255) tattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg > 1:774454/1‑151 (MQ=255) attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg > 1:202241/1‑150 (MQ=255) tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAAt > 1:132208/1‑151 (MQ=255) tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg > 1:607291/1‑150 (MQ=255) ttttcctggCCGCGTATTTCACCCAAAAGTGGAATGATTAATGAACGGATATACTCCCTCGCCTACCCGCCACCTTGAAACATTTATTCAAAAGCGCC‑C‑TGTTTTTTATGCCGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAAt < 1:403872/146‑1 (MQ=255) ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc > 1:879419/1‑150 (MQ=255) tctggCCGCTTTTTTCACCCAAAAAGGGGATGACTAAGGAAGGGAAAAACCCCCTCCCCTACCCGCCCCTTTTTTACATTTGT‑AAAAAGGGGC‑CTTTTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa < 2:283967/150‑1 (MQ=255) gCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt > 2:964597/1‑151 (MQ=255) tttttacggCCAAAAAGGGGATAAATAATGAACGAAAAAAACCCCTCCCCTACCCGCCCCTTTTTTACTGTTGTGAAAAGGGGGC‑T‑TGTTTTTTACGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa < 1:12653/139‑1 (MQ=255) ttccggCCAAAAAGGGGAGAACAAAGGACGGGAAATACCCCCCCACCAAACCGCCCCCTTGTAAGAGTGGTGAAAGAGGGTC‑T‑TGTTTTTGGTGCAGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt < 1:935595/144‑1 (MQ=255) ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgtg > 2:590588/1‑151 (MQ=255) aaaaggggtaaaaaaaTGAAGGGAGAAACCCCCTCACCTACCCGCCCCTTTTTTAAGGTTGTGTAGAGGGGGT‑T‑TTTTTTTTGGGGGGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt < 1:762800/137‑1 (MQ=255) aaaGGGGATGAAAAGTAAACGGATATAACCCCTCCCAAACCCGCCCCCTTTTCACAGTTGTG‑AAAAGGGGC‑CTTGTTTTTTAAAAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt < 1:203612/149‑1 (MQ=255) tGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAc > 1:116466/1‑151 (MQ=255) tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCG‑AAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc > 2:277484/1‑151 (MQ=255) tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCACCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg > 1:463996/1‑147 (MQ=255) tcccccctcccctcCCCGGCCTTTTGTTACGTTTGTGAAAAAGGGCC‑T‑TTTTTTTTATGAGGAAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc < 1:639408/145‑1 (MQ=255) cccTCCCCTACCCGCCCCTTTTTTACATTGGTGAACAGGGGGC‑T‑TTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt < 1:136961/148‑1 (MQ=255) ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa > 1:698773/1‑150 (MQ=255) ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAAGCTAACAGGCATTTTCCTGCACTGATAACGa > 1:456335/1‑150 (MQ=255) cGGCCCCTTTTTTCGGTTGTGAACAAGGGCT‑C‑TTTTTTTTAGGAGGGC‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt < 1:783831/149‑1 (MQ=255) cGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGGTAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg > 2:284288/1‑151 (MQ=255) ggCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGATAACAGGCATTTTCCTGCACTGATAACGAATCg > 1:134191/1‑147 (MQ=255) gcccctTTTTTAAATTTGTGAACAGGGGGC‑T‑TGATTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt < 1:716480/147‑1 (MQ=255) tttttACATTGGTTTACAAGGCGG‑C‑TTTTTTTTATGACGACAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca < 1:421100/146‑1 (MQ=255) gTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc > 1:769005/1‑149 (MQ=255) ttGTGTAAAGGGGGC‑T‑TGTTTTTAAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt < 1:944114/147‑1 (MQ=255) tgtACAGGGGCC‑C‑TGTTTTTTATGAGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt < 2:710878/150‑1 (MQ=255) tgtACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt > 1:484911/1‑151 (MQ=255) ggggctttttttttAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat < 2:132478/141‑1 (MQ=255) tgttttttATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa < 1:966390/144‑1 (MQ=255) tttttAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt < 1:665221/147‑1 (MQ=255) cGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTaaa > 1:204985/1‑150 (MQ=255) | TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA > NC_000913/1622668‑1622964 |
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: ‑ |