Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,622,815 | 0 | C | A | 46.0% | 39.3 / 21.2 | 38 | intergenic (+74/+145) | yneM/dgcZ | inner membrane‑associated protein/diguanylate cyclase, zinc‑sensing |
Reads supporting (aligned to +/- strand): ref base C (19/1); new base A (1/16); total (21/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.14e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG > NC_000913/1622675‑1622956 | cGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAAc > 1:1079555/1‑151 (MQ=255) gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAAcc > 2:687624/1‑151 (MQ=255) cTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGcc > 2:1098071/1‑151 (MQ=255) gggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCa < 1:258755/150‑1 (MQ=255) ttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGa > 2:299659/1‑150 (MQ=255) ttctggCCGTTTTTTCCACCCAAAAAGGGAATAATAAAAGACGGGAAATAACCCCCCCCCAACCCGCCCCTTTT‑TT‑ACTT‑TTTTGTACAGGGGG‑CCTTT‑TTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc < 2:1137844/150‑1 (MQ=255) tctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt > 2:917781/1‑151 (MQ=255) tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg > 2:759798/1‑151 (MQ=255) ggCCGCTTTTTCCAGCCAAAAGGGGGAGAACAAATGAGGGAAAATACCCCCCCCCCAACCCGCCCCTTTG‑TT‑ACAG‑TTGTGTAAAGGGGG‑CCTTT‑TTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt < 2:946541/148‑1 (MQ=255) ggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg > 2:412953/1‑150 (MQ=255) cgcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa > 2:750155/1‑148 (MQ=255) gcgTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt > 2:703908/1‑148 (MQ=255) ttttcaacccAAAAAGGGGATAAAAAATAAACGGAAATACCCCCCCCCCAACCCGCCCCTTTT‑TT‑ACGG‑TTGTAAACAGGGGG‑CCTGT‑TTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt < 2:569536/142‑1 (MQ=255) tATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑CTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt > 1:1023607/1‑151 (MQ=255) aTTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGCTTGCATGGATAGAttt > 2:699437/1‑150 (MQ=255) gggATGAAAAATGAACGGAGATACCCCCTCCCCTACCCGCCCCTTTT‑TT‑ACAT‑TTGTGTAAAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt < 2:772306/148‑1 (MQ=255) gggATAATTAATGAAGGGAAAT‑ACCCCCCCCCCAACCCGCCCCTTT‑TT‑TCAGTTTGTGAAAAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt < 2:785273/148‑1 (MQ=255) ggtggaaaagtGAAGGGATATACCCCCCCCCCTACCCGCCCCCTTT‑TT‑ACTG‑TTG‑GCAAAAGGGGTTTTGT‑TTTTTATGGGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGc < 1:703494/141‑1 (MQ=255) aCTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc > 1:968943/1‑150 (MQ=255) cGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt > 2:989082/1‑150 (MQ=255) gaTAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc > 2:704536/1‑151 (MQ=255) ccccccccTACCCGCCCCTTTG‑TT‑ACAG‑TTGTTTACAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa < 2:436960/144‑1 (MQ=255) cccccaccaaccctccccTTGTTT‑AAAG‑TTGGTTACAGGGGT‑CTTTA‑TTTTTAGGGGGAAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTg < 1:418999/135‑1 (MQ=255) ccTAACCGCCCCTTTG‑TT‑ACAT‑TTGTGAACAGGGGG‑CTTGT‑TTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa < 2:817180/151‑1 (MQ=255) cTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGaa > 2:550006/1‑150 (MQ=255) gCCCCTTT‑TTAACGG‑TTGTGAAAAGGGGG‑CTTTT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg < 2:894433/150‑1 (MQ=255) gCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGac > 2:838453/1‑151 (MQ=255) ccTTG‑TT‑ACAG‑TTTTGTACAAGGGG‑CCTGA‑TTTTTATGACGGGGCTAAAAAAACGCCAGTCAACCGGCGGTGAATGCTAGTACGGAAAGATTTGTGTTTTGCTTTTACGCCAACAGGGATTTTCCTCCACTGATAACGAATCGTTTactc > 1:626010/1‑148 (MQ=255) ttttaaaatttttgtAAAGGGGG‑CTTGT‑TTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca < 2:542073/137‑1 (MQ=255) tt‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGc > 1:111429/1‑151 (MQ=255) gttttgaagaaGGGG‑CTTGT‑TTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt < 1:1007457/139‑1 (MQ=255) ttgtaaaaaaGGGC‑CTTGT‑TTTTTAAGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt < 2:26329/139‑1 (MQ=255) ttgggaaaaGGGG‑CCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc < 2:563835/141‑1 (MQ=255) aggggg‑CCTTT‑TTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCa < 2:3807/148‑1 (MQ=255) gggg‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGGATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCCGTTTTCTCa > 1:858683/1‑148 (MQ=255) ggg‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatg > 2:687958/1‑150 (MQ=255) ttttAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTaaa < 2:750953/150‑1 (MQ=255) ttATGACGGAGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGAGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg > 1:409342/1‑151 (MQ=255) | CGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTG‑TT‑ACAG‑TTGTGTACAAGGGG‑CCTGA‑TTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG > NC_000913/1622675‑1622956 |
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: ‑ |