Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,622,802 | 0 | A | T | 42.9% | 61.5 / 34.6 | 49 | intergenic (+61/+158) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
Reads supporting (aligned to +/- strand): ref base A (21/7); new base T (0/21); total (21/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.24e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATG > NC_000913/1622652‑1622948 | ggTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGa < 2:214745/151‑1 (MQ=255) aaTATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑ttt < 2:935988/151‑1 (MQ=255) aaTATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑tt < 2:639882/150‑1 (MQ=255) atatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑tttt > 2:550454/1‑151 (MQ=255) atatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑tttt < 2:560536/151‑1 (MQ=255) tatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑ttttt > 1:667346/1‑151 (MQ=255) atGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑ttttt > 2:268439/1‑150 (MQ=255) gTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACg < 1:104302/150‑1 (MQ=255) ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGa > 1:594321/1‑150 (MQ=255) tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGaa > 1:242965/1‑150 (MQ=255) ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGaaaaa > 1:679536/1‑150 (MQ=255) gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑cc > 2:357729/1‑151 (MQ=255) gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑cc < 2:871019/151‑1 (MQ=255) tGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGcc > 2:365767/1‑150 (MQ=255) ggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGt > 1:140035/1‑151 (MQ=255) ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACcgg > 2:914901/1‑151 (MQ=255) tttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACcgg > 1:339189/1‑150 (MQ=255) tttttttctgtttttctggCCGCTTTTTCCACCAAAAAAGGGGATGACTAATAAAGGAAAATACCCCCCCCCCTACCCGCCCCTTT‑TTTACA‑GTTGTGTACAAGGGGCT‑TGT‑TTTTTATGACGG‑CAAAAAAAAA‑CCGCCAGTAAACcggc < 1:168591/147‑1 (MQ=255) ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGt > 2:305902/1‑150 (MQ=255) tgttgtgtttttcgggCCGATTATTACACACAAAAGGGGGATGAAAAATGAGGGGAAAAAACCCCTCA‑CAAACCGGCCCCTTGGTTAGA‑TTTG‑GGAAAGGGGGCCTTGT‑TTTTTAGGACGA‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGt < 2:113674/142‑1 (MQ=255) tttttgtttttcgggCCGCTTATTCCACCCAAAAAGGGGATGAATAATGAAGGGA‑AAAAACCCCCCCCCAACCCGCCCCTT‑TTTAAATGTTGTGAACAGGGGCCC‑TGT‑TTTTTATGCCGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTg < 1:909710/150‑1 (MQ=255) tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGaa > 1:546514/1‑150 (MQ=255) gtttttcgggCCGGTTTTTTCACCAAAAAAGGGGATAAATAATGAACGGAAATAACCCCTCCCCTACCCGCCCCCTT‑GTACCG‑GTGGTGCACAGGGGGCT‑TGT‑TTTTTGGGACGG‑CAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGc < 1:554708/148‑1 (MQ=255) gttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCtt > 1:732063/1‑151 (MQ=255) ttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCtt > 2:389026/1‑150 (MQ=255) tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATgg > 1:643774/1‑151 (MQ=255) gcccGCTTTTTTCACCAAAAAAGGGGATGACAAATGAAGGGAAATACCCCCTCCCCTACCCGCCCCTTT‑TTTACG‑TTTGTGTAAAGGGGGCC‑TGT‑TTTTTAGGACGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCAt < 1:99563/146‑1 (MQ=255) gcccGCTTATTCCACCCAAAAAGGGGATGTCAAATGAAGGGAGATACCCCCCCCCCAACCCGCCCCTTT‑TTTAAT‑GTTGTGAAAAAGGGCCT‑TGT‑TTTTTATGAGGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa < 2:89593/149‑1 (MQ=255) ttATTCCAGCAAAAAAGGGGAAGACAAAGGAGAGGAGAAATTCCCCCACCTACCCGCCCTTTT‑TTTAGG‑GTTGTAAAAAAGGGGCT‑TGT‑TTTTTAGGGGGGAAAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg < 1:22816/148‑1 (MQ=255) acccAAAAAGGGGATAACTAAGGAACGGAAATACCCCCCCCCCAACCCGCCCCTTT‑TTTAAA‑GTTGTGTACAGGGGGCT‑TTT‑TTTTTATGACGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtt < 1:119749/148‑1 (MQ=255) aaaaGGGGATGACTAATAAACGAAAATACCCCCTCCCCAACCCGCCCCTTT‑TTTACG‑TTTGTGTAAAGGGGGCT‑TTT‑TTTTTAGGACGG‑CAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg < 1:747251/147‑1 (MQ=255) aagtGGGATAACAAATGAACGAGAAAAACCCCTCCCCTACCCCGCCCCTT‑TTTAAA‑TTTGTTTACAAGGGCCT‑TTT‑TTTTTAGCACGG‑AAAAAAAAAT‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt < 2:758904/148‑1 (MQ=255) aaaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt > 1:157548/1‑151 (MQ=255) aaggggATGAAAAATGACGGGAAAAAACCCCCCCCCTACCCGCCCCTTT‑TTTACG‑GTTGTGTACAGGGGGCT‑TTT‑TTTTTATGACGG‑GAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt < 1:304861/148‑1 (MQ=255) gggAAGAATAATAAACGGATAAACCCCCCCACCTACCCCCCCCTTT‑GTTATG‑GTGGTTAACAAGGGGCT‑TGT‑TTTTGATGAGGA‑CAAAAAAAAT‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt < 2:506328/148‑1 (MQ=255) aTGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt > 1:756864/1‑150 (MQ=255) aaTAAAGGGAAATAACCCCTCCCCAACCCGCCCCTTT‑TTTACA‑TTTTTG‑AAAAGGGGCCTTTT‑TTTTTAGGCGGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg < 1:692412/149‑1 (MQ=255) aaCGGAAATACCCCCCCCCCTACCCGCCCCTTT‑GTTAGG‑TTTGTGTACAGGGGGGG‑TGT‑TTTTTAGGAGGG‑AAAAAAAAAT‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa < 2:37782/144‑1 (MQ=255) aaaaaccccccCACCTACCGGCCCCCTT‑TTTACA‑G‑TGTTTAACAGGGGCCTTGT‑TTTTTATGGCGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt < 1:937981/140‑1 (MQ=255) caccgccccGTCCCTTT‑GTTGCA‑GTTGTGTACAAGGGGCT‑TGT‑TTTTTTTGACGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg < 1:598467/144‑1 (MQ=255) ccaacccGGCCCTTT‑TTTACT‑GTTGTGTAAAAGGGTCT‑TTT‑TTTTGATGAGGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg < 2:772904/144‑1 (MQ=255) ccGCCCCTTT‑TTTACA‑TTTGTGTAAAGGGGCCC‑TTT‑TTTTTAGGGGGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt < 1:680077/147‑1 (MQ=255) ttGTGTAGAAGGGGCT‑TGT‑TTTTTAGGACGA‑CAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc < 2:620522/148‑1 (MQ=255) gtgtACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt > 1:256410/1‑151 (MQ=255) aaaaGGGGCC‑TTTTTTTTTAGGACGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct < 1:684662/149‑1 (MQ=255) aCAGGGGGCT‑TGT‑TTTTTAGGGCGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc < 2:643766/149‑1 (MQ=255) cAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc > 2:587303/1‑150 (MQ=255) acaGGGGGC‑CGA‑TTTTTTAGACGG‑AGAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc < 2:944703/147‑1 (MQ=255) aGGGGCCTTTA‑TTTTTAGGAGGG‑AAAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt < 1:608465/146‑1 (MQ=255) cc‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg > 1:825493/1‑151 (MQ=255) | GGTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTT‑GTTACA‑GTTGTGTACAAGGGGCC‑TGA‑TTTTTATGACGG‑CGAAAAAAAA‑CCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATG > NC_000913/1622652‑1622948 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |